当前位置:   article > 正文

复信号的傅里叶变换是什么?频谱是什么样子的?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仿真验证

  1. %%=======================================================================================
  2. % 程序说明:本程序验证复信号的傅式变换并且仿真复信号的频谱
  3. % 程序名称:复信号的频谱分析
  4. % 作者: hill5678
  5. % 当前版本:1.0
  6. %%=======================================================================================
  7. clear all;
  8. close all;
  9. clc;
  10. N=128; % 采样点数
  11. n=0:N-1;
  12. fm=1000; % 载波频率
  13. fs=1e4; % 采样频率
  14. t=1/fs*n; % 时间向量
  15. s3=cos(2*pi*fm*t)+1j*sin(2*pi*fm*t); % 构造复信号
  16. s3_1=repmat(s3,N,1);
  17. y1=zeros(length(s3_1)/2-1,length(s3_1));
  18. y2=zeros(length(s3_1)/2,length(s3_1));
  19. s3_1(1:end/2-1,:)=y1;
  20. s3_1(end/2+1:end,:)=y2;
  21. f=[0:N-1]*fs/N; % 真实频率
  22. mag1=abs((fft(s3_1,N,2))); % FFT点数与采样点数相同
  23. mag2=abs((fft(s3,N,2)));
  24. figure(1)
  25. plot(f,mag2);
  26. xlabel('频率f/Hz');
  27. ylabel('幅度');
  28. title('复信号s3的频谱');
  29. grid on;
  30. figure(2)
  31. plot(t,s3);
  32. xlabel('时间');
  33. ylabel('幅度');
  34. title('复信号的时域波形');
  35. grid on;
  36. figure(3)
  37. mesh(f,n,mag1)
  38. set(gca, 'color', [202 / 255, 234 / 255, 206 / 255]);
  39. xlabel('频率/Hz','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'r');
  40. ylabel('FFT点数','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'r');
  41. zlabel('幅值','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'r');
  42. title('复信号s3的3D频谱图','FontSize', 15 , 'FontWeight', 'bold', 'Color', 'k')

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

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/Li_阴宅/article/detail/967890
推荐阅读
相关标签
  

闽ICP备14008679号