q712.m

这是matlab在电子信息课程中的应用

%《MATLAB在电子信息课程中的应用》第七章例7.12程序q712
% 梳状滤波器幅频特性与极点位置
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

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);
subplot(2,2,1);zplane(b,a0);title('FIR梳状滤波器零点图');
subplot(2,2,2);zplane(b,a1);title('IIR梳状滤波器零、极点图,a=0.8');
subplot(2,2,3);plot(w/pi,abs(H));title('FIR梳状滤波器幅频响应曲线');
ylabel('幅度');xlabel('ω/π');
subplot(2,2,4);plot(w1/pi,abs(H1));title('IIR梳状滤波器幅频响应曲线,a=0.8');
ylabel('幅度');xlabel('ω/π');
figure(2);
subplot(2,2,1);zplane(b,a2);title('IIR梳状滤波器零、极点图,a=0.9');
subplot(2,2,2);zplane(b,a3);title('IIR梳状滤波器零、极点图,a=0.98)');
subplot(2,2,3);plot(w2/pi,abs(H2));title('IIR梳状滤波器幅频响应曲线,a=0.9');
ylabel('幅度');xlabel('ω/π');
subplot(2,2,4);plot(w3/pi,abs(H3));title('IIR梳状滤波器幅频响应曲线,a=0.98');
ylabel('幅度');xlabel('ω/π')