hc552.m

来自「matlab数字信号处理代码」· M 代码 · 共 29 行

%《数字信号处理教程——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');                % 置图形背景色为白