%《數字信號處理教程——MATLAB釋義與實現》第四章例4.3.2程序hc432
% 用FFT計算離散時間傅立葉變換DTFT
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
C=[20,1024];T=0.5;      % 設定兩種補零后長度
for i=[1,2]
    N=C(i);D=2*pi/(N*T);                    % 根據N值確定頻率分辨率
    x=[1 1 1 1 1 1 zeros(1,N-11) 1 1 1 1 1];% 根據N值確定序列
    X=fftshift(fft(x,N));                   % 作DFT，移到對稱位置
    k=floor(-N/2+0.5:N/2-0.5);              % 確定對稱的頻率位置向量
    subplot(1,2,i),plot(k*D,X)              % 繪制幅頻特性
    xlabel('頻率'),ylabel('幅度')
    line([-10,10],[0,0])
end
set(gcf,'color','w')              % 置圖形背景色為白