example8_8.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

%例8-8程序:example8_8
Wp=0.2*pi;Ws=0.3*pi;Rp=1;Ar=15; T=1;    %数字滤波器设计指标
OmegaP=Wp*T;                            %模拟原型滤波器通带截止频率
OmegaS=Ws*T;                            %模拟原型滤波器阻带截止频率
%
[cs1,ds1]=afd_butt(OmegaP,OmegaS,Rp,Ar);    %巴特沃思模拟滤波器设计
[b1,a1]=imp_invr(cs1,ds1,T);                  %用脉冲响应不变法实现模拟滤波器到数字滤波器的变换
[db1,mag1,pha1,grd1,w1]=freqz_m(b1,a1);        %计算频率响应和相位响应
[C1,B1,A1]=dir2par(b1,a1)                    %直接型到并联型转换
%画图
figure(1)
subplot(3,4,1)
plot(w1/pi,mag1);
axis([0 1 0 1.1]);
ylabel('(a)')
title('幅度响应');
subplot(3,4,2);
plot(w1/pi,db1);
axis([0 1 -40 5]);
title('幅度响应（dB）');
subplot(3,4,3);
plot(w1/pi,pha1/pi);
axis([0 1 -1 1]);
title('相位响应');
subplot(3,4,4);
plot(w1/pi,grd1);
axis([0 1 0 10]);
title('群延迟');
%
[cs2,ds2]=afd_chb1(OmegaP,OmegaS,Rp,Ar);            %未归一化切比雪夫Ι型模拟低通滤波器设计
[b2,a2]=imp_invr(cs2,ds2,T);                       %用脉冲响应不变法实现模拟滤波器到数字滤波器的变换
[db2,mag2,pha2,grd2,w2]=freqz_m(b2,a2);             %计算频率响应和相位响应
[C2,B2,A2]=dir2par(b2,a2)                         %直接型到并联型转换
%画图
subplot(3,4,5)
plot(w2/pi,mag2);
axis([0 1 0 1.1]);
ylabel('(b)')
subplot(3,4,6);
plot(w2/pi,db2);
axis([0 1 -40 5]);
subplot(3,4,7);
plot(w2/pi,pha2/pi);
axis([0 1 -1 1]);
subplot(3,4,8);
plot(w2/pi,grd2);
axis([0 1 0 15]);
%
[cs3,ds3]=afd_elip(OmegaP,OmegaS,Rp,Ar);            %未归一化切比雪夫Ⅱ型模拟低通滤波器设计
[b3,a3]=imp_invr(cs3,ds3,T);                       %用脉冲响应不变法实现模拟滤波器到数字滤波器的变换
[db3,mag3,pha3,grd3,w3]=freqz_m(b3,a3);             %计算频率响应和相位响应
[C3,B3,A3]=dir2par(b3,a3)                         %直接型到并联型转换
%画图
subplot(3,4,9)
plot(w3/pi,mag3);
axis([0 1 0 1.2]);
ylabel('(c)')
subplot(3,4,10);
plot(w3/pi,db3);
axis([0 1 -40 5]);
subplot(3,4,11);
plot(w3/pi,pha3/pi);
axis([0 1 -1 1]);
subplot(3,4,12);
plot(w3/pi,grd3);
axis([0 1 0 45]);