1.用Matlab符号运算求解法求单边指数信号f(t)=e−2tu(t)f(t)=e^{-2t}u(t)f(t)=e−2tu(t)的FT
MATLAB源程序为:
ft=sym('exp(-2*t)*heaviside(t)');
fw=fourier(ft)
运行结果为:
fw =
1/(2 + w*1i)
2.用Matlab符号运算求解法求F(jw)=11+w2F(jw)=\frac{1}{1+w^2}F(jw)=1+w21的IFT。
解:matlab源程序为
syms t
fw=sym('1/(1+w^2)');
ft=ifourier(fw,t)
运行结果为:
ft =
exp(-abs(t))/2
3.用Matlab命令绘出信号f(t)=e−2tu(t)f(t)=e^{-2t}u(t)f(t)=e−2tu(t)的频谱图。
Matlab源程序为:
ft=sym('exp(-2*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('相位谱')
调用heaviside函数需要与我们定义的函数名一致,一开始h大写了,就报错,后来改为小写就正确了。