samp5_8.m

来自「MATLAB7.x数字信号处理 光盘内容」· M 代码 · 共 16 行

%Samp5_8
Rs=20;    %滤波器的阻带衰减为20dB
[z,p,k]=cheb2ap(6,Rs);   %设计6阶阻带衰减为20dB的Chebyshev II型原型滤波器
[b,a]=zp2tf(z,p,k);    %将零点极点增益形式转换为传递函数形式
[H,w]=freqs(b,a,0:0.01:2);   %计算模拟原型滤波器的复数频率响应
subplot(2,1,1),plot(w,abs(H).^2);  %绘制平方幅频响应
xlabel('w/wc');ylabel('Chebyshev II |H(jw)|^2');
title('Chebyshev II 型低通原型滤波器(wc=1)')
w1=0.8;w2=1.4;   %给定欲设计滤波器通带的下限和上限频率
w0=sqrt(w1*w2);bw=w2-w1; %计算中心点频率和频带宽度
[bt,at]=lp2bp(b,a,w0,bw);  %频率转换
[Ht,wt]=freqs(bt,at,0:0.01:2);  %计算滤波器的复数频率响应
subplot(2,1,2),plot(wt,abs(Ht).^2); %绘制平方幅频响应
xlabel('w/wc');ylabel('Chebyshev II |H(jw)|^2');
title('Chebyshev II型带通滤波器(wc=0.8~1.4)')