hc351.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第三章例3.5.1程序hc351
% 序列的离散时间傅立叶变换(DTFT)和离散傅立叶变换(DFT)
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
x=[2,-1,1,1];Xd=fft(x);            % 输入序列，求FFT
Xd1=fftshift(Xd);                  % 把FFT移位到对称位置
K=64;dw=2*pi/K;                  % 求频率分辨率
k=floor((-K/2+0.5):(K/2-0.5));      % 确定频率下标向量
X=x*exp(j*dw*nx'*k);                % 求离散傅立叶变换DFT
subplot(2,1,1),plot(k*dw,abs(X)),hold on      % 画DTFT幅频特性
plot([0:3]*2*pi/4,abs(Xd),'o')          % 画DFT幅频特性
subplot(2,1,2),plot(k*dw,abs(X)),hold on   % 画DTFT相频特性
plot([-2:1]*2*pi/4,abs(Xd1),'x')      % 画DFT相频特性
set(gcf,'color','w')              % 置图形背景色为白