以求解下面的傅里叶变换为例![](https://img-blog.csdnimg.cn/c15a8d386a694b9ba68ef1f82daab871.png?x-oss-process=image/watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBA6K-35oiR5ZCD54K46bih,size_20,color_FFFFFF,t_70,g_se,x_16)
核心代码:
Fw=fourier(ft)
代码:
clear all;
close all;
ft=sym('exp(-1i*t)*(heaviside(t))')
%ft=sym('sin(pi*t)*(heaviside(t)-heaviside(t-2))')
%ft=sym('sin(2*pi/100*t)*(heaviside(t+50))')
%ft=sym('exp(-i*t)*(heaviside(t))')
Fw=fourier(ft)
subplot(211),ezplot(abs(Fw)),grid on,title('幅度谱')
phase=atan(imag(Fw)/real(Fw));
subplot(212),ezplot(phase);grid on,title('相位谱')
命令行窗口显示:
ft =
exp(-t*1i)*heaviside(t)
Fw =
pi*dirac(w + 1) - 1i/(w + 1)
>>
实验截图