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