% dm10201
% 門信號的傅里葉變換
syms t w ;                                      %定義符號變量
ut = sym('Heaviside(t+0.5)-Heaviside(t-0.5)');  %借助階躍函數定義門函數
subplot(211);  
ezplot(ut);                                     %繪制門信號
hold on;
axis([-1 1 0 1.1]);                             %給定坐標系的范圍
plot([-0.5 -0.5], [0 1]);                       %繪制左跳變沿
plot([0.5 0.5], [0 1]);                         %繪制右跳變沿
title('時域波形')
Fw = fourier(ut,t,w);                           %用fourier函數對門信號求傅里葉變換
FFw = maple('convert',Fw,'piecewise');          %將傅里葉變換的符號表達式轉換為樣條曲線
FFP = abs(FFw);                                 %求幅度頻譜 
subplot(212);
ezplot(FFP,[-10*pi 10*pi]);                     %繪制幅度頻譜
axis([-10*pi 10*pi 0 1.1]);                     %給定坐標系的范圍
title('幅度頻譜')
%End