hc352.m

来自「matlab数字信号处理代码」· M 代码 · 共 12 行

%《数字信号处理教程——MATLAB释义与实现》第三章例3.5.2程序hc352
% 序列的离散傅立叶变换
% 电子工业出版社出版  陈怀琛编著 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')              % 置图形背景色为白