%《MATLAB在電子信息課程中的應(yīng)用》第七章例7.12程序q712
% 梳狀濾波器幅頻特性與極點(diǎn)位置
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 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梳狀濾波器零點(diǎn)圖');
subplot(2,2,2);zplane(b,a1);title('IIR梳狀濾波器零、極點(diǎn)圖,a=0.8');
subplot(2,2,3);plot(w/pi,abs(H));title('FIR梳狀濾波器幅頻響應(yīng)曲線');
ylabel('幅度');xlabel('ω/π');
subplot(2,2,4);plot(w1/pi,abs(H1));title('IIR梳狀濾波器幅頻響應(yīng)曲線,a=0.8');
ylabel('幅度');xlabel('ω/π');
figure(2);
subplot(2,2,1);zplane(b,a2);title('IIR梳狀濾波器零、極點(diǎn)圖,a=0.9');
subplot(2,2,2);zplane(b,a3);title('IIR梳狀濾波器零、極點(diǎn)圖,a=0.98)');
subplot(2,2,3);plot(w2/pi,abs(H2));title('IIR梳狀濾波器幅頻響應(yīng)曲線,a=0.9');
ylabel('幅度');xlabel('ω/π');
subplot(2,2,4);plot(w3/pi,abs(H3));title('IIR梳狀濾波器幅頻響應(yīng)曲線,a=0.98');
ylabel('幅度');xlabel('ω/π')