以求解下面的傅里叶变换为例

核心代码：

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)
>>

实验截图