example8_5.m

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

%例8-5程序:example8_5
Wp=0.2*pi;Ws=0.3*pi;Rp=1;Ar=16;     %设计指标
wmax=0.5*pi;
%
[b1,a1]=afd_butt(Wp,Ws,Rp,Ar);        %巴特沃思模拟滤波器设计
[db1,mag1,pha1,w1]=freqs_m(b1,a1,wmax);   %计算频率响应和相位响应
[ha1,x1,t1]=impulse(b1,a1);
[C1,B1,A1]=sdir2cas(b1,a1)               %直接型到级联型转换
%画图
figure(1)
subplot(4,3,1)
plot(mag1)
axis([0 500 0 1.05]);
ylabel('(a)')
title('频率响应')
subplot(4,3,2)
plot(pha1)
axis([0 500 -3.2 3.2]);
title('相位响应')
subplot(4,3,3)
plot(ha1)
axis([0 100 -0.08 0.25]);
title('脉冲响应')
%
[b2,a2]=afd_chb1(Wp,Ws,Rp,Ar)        %未归一化切比雪夫Ι型模拟低通滤波器设计
[db2,mag2,pha2,w2]=freqs_m(b2,a2,wmax);   %计算频率响应和相位响应
[ha2,x2,t2]=impulse(b2,a2);
[C2,B2,A2]=sdir2cas(b2,a2)               %直接型到级联型转换
%画图
subplot(4,3,4)
plot(mag2)
axis([0 500 0 1.1]);
ylabel('(b)')
subplot(4,3,5)
plot(pha2)
axis([0 500 -3.2 3]);
subplot(4,3,6)
plot(ha2)
axis([0 110 -0.08 0.22]);
%
[b3,a3]=afd_chb2(Wp,Ws,Rp,Ar)        %未归一化切比雪夫Ⅱ型模拟低通滤波器设计
[db3,mag3,pha3,w3]=freqs_m(b3,a3,wmax);   %计算频率响应和相位响应
[ha3,x3,t3]=impulse(b3,a3);
[C3,B3,A3]=sdir2cas(b3,a3)               %直接型到级联型转换
%画图
subplot(4,3,7)
plot(mag3)
axis([0 500 0 1.05]);
ylabel('(c)')
subplot(4,3,8)
plot(pha3)
axis([0 500 -3 3.2]);
subplot(4,3,9)
plot(ha3)
axis([0 110 -0.4 0.3]);
%
[b4,a4]=afd_elip(Wp,Ws,Rp,Ar)        %模拟椭圆滤波器设计
[db4,mag4,pha4,w4]=freqs_m(b4,a4,wmax);   %计算频率响应和相位响应
[ha4,x4,t4]=impulse(b4,a4);
[C4,B4,A4]=sdir2cas(b4,a4)               %直接型到级联型转换
%画图
subplot(4,3,10)
plot(mag4)
axis([0 500 0 1.05]);
ylabel('(d)')
subplot(4,3,11)
plot(pha4)
axis([0 500 -3 3.2]);
subplot(4,3,12)
plot(ha4)
axis([0 110 -0.08 0.3]);