%《數(shù)字信號處理教程——MATLAB釋義與實現(xiàn)》第三章例3.5.2程序hc352
% 序列的離散傅立葉變換
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
xm=[2,-1,1,1,zeros(1,60)];      % 輸入序列，
lxm=length(xm);nxm=0:lxm-1;
Xdm=fftshift(fft(xm));       % 求FFT并把它移位到對稱位置
%Xdm=fft(xm);
kxm=floor(-lxm/2+0.5:lxm/2-0.5); % 把頻率下標也移位到對稱位置
plot(kxm*2*pi/lxm,abs(Xdm),'o')    % 畫DFT幅頻特性
set(gcf,'color','w')              % 置圖形背景色為白