dm10101.m

来自「各种周期信号的傅里叶变换」· M 代码 · 共 17 行

% dm10101
% 单边指数信号的傅里叶变换
syms t v w x phase im re;                 %定义符号变量
f = exp(-2*t)*sym('Heaviside(t)');        %f(t)=exp(-2*t)*u(t)
Fw = fourier(f);                          %求傅里叶变换
subplot(311);  
ezplot(f);                                %绘制f(t)的时域波形
axis([-1 2.5 0 1.1]);
subplot(312);
ezplot(abs(Fw));                          %绘制幅度谱 
im = imag(Fw);                            %计算F(w)的虚部 
re = real(Fw);                            %计算F(w)的实部
phase = atan(im/re);                      %计算相位谱 
subplot(313);
ezplot(phase);                            %绘制相位谱
%End