hc453a.m

某本matlab书籍中的代码

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