%《數字信號處理教程——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')