example8_10.m

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

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