复信号的傅里叶变换是什么?频谱是什么样子的?3D频谱图长啥样子?

1.复信号的数学表达式 

 

大家都知道,复数是由实数与虚数构成。同理,复信号也可以有一个实信号和一个虚信号构成。

数学表达式可以表示为:

S(t)=Sr(t)+jSi(t) (1)

这里我们还可以回想起经典的欧拉公式:

e^{jwt}=coswt+jsinwt (2)

这个公式将复变函数,三角函数以及指数函数巧妙的结合在了一起。如果定义一个复平面,其横坐标就是实数,纵坐标就是虚数,诸如此类的函数我们叫它复变函数,并且它实际上是绕原点旋转的圆,如下图:

 其中θ=wt=2*pi*t/T,该复变函数可以看做是一个角速度为 w 周期为T在复平面上绕原点旋转的半径为R的圆,如果R=1,则该圆为单位圆。

2.复信号的傅里叶变换

如果我们直接给公示(1)套用CFT公式,我们会发现摸不着头脑,不知道该如何求复信号的傅里叶变换。这时我们其实可以类比于公示(2)发现所有的复信号都可以写成实部是余弦函数,虚部是正弦函数的组合形式,那么我们有:

S(t)=Acosw_{_{s}}t+Ajsinw_{s}t(3)

我们这里就架设A=1,则:

S(t)=cosw_{_{s}}t+jsinw_{s}t(4)

给S(t)求CFT,得到:

F[S(t)]=\int_{-\infty }^{\infty }(cosw_{s}t+jsinw_{s}t)e^{-jwt}dt=\int_{-\infty }^{\infty }(cosw_{s}t)e^{-jwt}dt+j\int_{-\infty }^{\infty }(sinw_{s}t)e^{-jwt}dt(5)

我们可以清楚的发现复信号的傅里叶变换可以变为正弦函数的傅式变换和余弦函数的傅式变换之和,我们知道正弦函数和余弦函数的傅式变换分别为:

-j\pi \delta (w-w_{s})+j\pi \delta (w+w_{s})(6)

\pi\delta (w-w_{s})+\pi \delta (w+w_{s})(7)

这里就不帮大家回顾该怎么推导了,那我们就可以很容易得到复信号S(t)的傅式变换为(6)式与(7)式之和,即:

F[S(t)]=\pi\delta (w-w_{s})+\pi \delta (w+w_{s})+j[-j\pi \delta (w-w_{s})+j\pi \delta (w+w_{s})]=2\pi \delta(w-w_{s})(8)

这样我们就得到了复信号S(t)的傅式变换。

3.复信号的频谱图

从公式(8)我们可以看到复信号S(t)的傅式变换就Delta冲激函数的向右频移的2π倍,也就是频点Ws对应的冲激的2π倍,那么它的频谱其实就是一根简单的谱线,如下图:

4.复信号3D频谱图 

如果我们要把频谱画到三维图里面,在Matlab是可以用mesh函数来实现的。如下图:

 但是我们需要注意mesh的X,Y,Z参数必须设置正确,否则会造成矩阵维度不一致的情况出现。

5.Matlab仿真验证

%%=======================================================================================
% 程序说明:本程序验证复信号的傅式变换并且仿真复信号的频谱
% 程序名称:复信号的频谱分析
% 作者: hill5678
% 当前版本:1.0
%%=======================================================================================
clear all;
close all;
clc;
N=128; % 采样点数
n=0:N-1;
fm=1000; % 载波频率
fs=1e4; % 采样频率 
t=1/fs*n; % 时间向量
s3=cos(2*pi*fm*t)+1j*sin(2*pi*fm*t); % 构造复信号
s3_1=repmat(s3,N,1);
y1=zeros(length(s3_1)/2-1,length(s3_1));
y2=zeros(length(s3_1)/2,length(s3_1));
s3_1(1:end/2-1,:)=y1;
s3_1(end/2+1:end,:)=y2;
f=[0:N-1]*fs/N; % 真实频率
mag1=abs((fft(s3_1,N,2))); % FFT点数与采样点数相同
mag2=abs((fft(s3,N,2)));
figure(1)
plot(f,mag2);
xlabel('频率f/Hz');
ylabel('幅度');
title('复信号s3的频谱');
grid on;
figure(2)
plot(t,s3);
xlabel('时间');
ylabel('幅度');
title('复信号的时域波形');
grid on;
figure(3)
mesh(f,n,mag1)
set(gca, 'color', [202 / 255, 234 / 255, 206 / 255]);
xlabel('频率/Hz','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'r');
ylabel('FFT点数','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'r');
zlabel('幅值','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'r');
title('复信号s3的3D频谱图','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'k')

*如果上述内容有误,请大家留言指正,我会及时修改*

  • 21
    点赞
  • 84
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

一枚雷达信号处理领域的小白

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值