hc453.m

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

%《数字信号处理教程——MATLAB释义与实现》第四章例4.5.3程序hc453
% 连续信号频谱计算中对称双边序列的算例
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
T=1;M=5/T;N=2*M+1;n=1:M;        % 原始数据
D=2*pi/(N*T);                  % 频率分辨率
hp=sin(2*n*T)./(pi*n*T);        % 生成正序列
h=[2/pi hp fliplr(hp)];         % 构成循环对称（偶）序列
H=T*fftshift(fft(h));H(1)       % 求x的FFT，移到对称位置
k=floor(-(N-1)/2:(N-1)/2);      % 乃奎斯特频率下标向量
subplot(1,2,1),plot(k*D,H)     % 绘图
set(gcf,'color','w');              % 置图形背景色为白
grid on,
title('T=1,N=11')