q712.m

MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全合著

%《MATLAB及在电子信息课程中的应用》第七章例7.12程序q712
% 梳状滤波器幅频特性与极点位置
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear;close all
b=[1,0,0,0,0,0,0,0,-1];
a0=1;
a1=[1,0,0,0,0,0,0,0,-(0.8)^8];          % 给出四种梳状滤波器分母系数
a2=[1,0,0,0,0,0,0,0,-(0.9)^8];
a3=[1,0,0,0,0,0,0,0,-(0.98)^8];
[H,w]=freqz(b,a0);                      % 求它们的频率特性                     
[H1,w1]=freqz(b,a1);
[H2,w2]=freqz(b,a2);
[H3,w3]=freqz(b,a3);
figure(1);                 
set(gcf,'color','w')                    % 设置图形背景色为白色
subplot(2,2,1);zplane(b,a0);            % 画出FIR梳状滤波器的零极点图
title('FIR梳状滤波器零点图');           
subplot(2,2,2);zplane(b,a1);            % 画出第一个IIR梳状滤波器的零极点图
title('IIR梳状滤波器零、极点图,a=0.8');
subplot(2,2,3);plot(w/pi,abs(H));
title('FIR梳状滤波器幅频响应曲线');     % 画出FIR梳状滤波器的幅频特性
ylabel('幅度');xlabel('ω/π');
subplot(2,2,4);plot(w1/pi,abs(H1));     % 画出第一个IIR梳状滤波器的幅频特性
title('IIR梳状滤波器幅频响应曲线,a=0.8');
ylabel('幅度');xlabel('ω/π');
figure(2);
set(gcf,'color','w')                    % 设置图形背景色为白色
subplot(2,2,1);zplane(b,a2);title       % 画出第二个IIR梳状滤波器的零极点图
('IIR梳状滤波器零、极点图,a=0.9');
subplot(2,2,2);zplane(b,a3);            % 画出第三个IIR梳状滤波器的零极点图
title('IIR梳状滤波器零、极点图,a=0.98)');
subplot(2,2,3);plot(w2/pi,abs(H2));     % 画出第二个IIR梳状滤波器的幅频特性
title('IIR梳状滤波器幅频响应曲线,a=0.9');
ylabel('幅度');xlabel('ω/π');
subplot(2,2,4);plot(w3/pi,abs(H3));     % 画出第三个IIR梳状滤波器的幅频特性
title('IIR梳状滤波器幅频响应曲线,a=0.98');
ylabel('幅度');xlabel('ω/π')