% 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