%《數字信號處理教程——MATLAB釋義與實現》第四章例5.5.2程序hc552
% 截取長度及對稱性對頻率特性的影響
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
N0=[6,15];
for k=[1:2] 
    n0=N0(k);N=2*n0+1;n=[0:N-1]+1e-5;   % 設置關于n0對稱的位置向量
    h=sin(n-n0)./(n-n0)/pi;             % 對稱截斷理想濾波器的脈沖響應序列
    [H,w]=freqz(h,1,4096);              % 對稱截斷理想濾波器的頻率響應
    subplot(3,2,2*k-1),plot(w,abs(H))   % 畫幅頻特性
    title('n0=6,N=13')
    line([0,pi],[0,0])
    subplot(3,2,2*k),plot(w,angle(H))   % 畫相頻特性
    title('n0=15,N=31')
    axis([0,pi/2,-4,4])
    line([0,pi],[0,0])
end
n0=6;N=2*n0+2;n=[0:N-1]+1e-5;   % 設置關于n0不對稱的位置向量
h=sin(n-n0)./(n-n0)/pi;         % 不對稱截斷理想濾波器的脈沖響應序列
[H,w]=freqz(h,1);               % 不對稱截斷理想濾波器的頻率響應
subplot(3,2,5),plot(w,abs(H))   % 畫幅頻特性
title('n0=6,N=14')
line([0,pi],[0,0])
subplot(3,2,6),plot(w,angle(H))  % 畫相頻特性
title('n0=6,N=14')
axis([0,pi/2,-4,4])
line([0,pi],[0,0])
set(gcf,'color','w');                % 置圖形背景色為白