%《數字信號處理教程——MATLAB釋義與實現》第四章例4.5.3程序hc453
% 連續信號頻譜計算中對稱雙邊序列的算例
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
T=1;M=5/T;N=2*M+1;n=1:M;        % 原始數據
DOmega=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*DOmega,H)     % 繪圖
set(gcf,'color','w');              % 置圖形背景色為白
grid on,
title('T=1,N=11')