samp5_7.m

MATLAB7.x数字信号处理 光盘内容

%Samp5_7
Rp=0.5;   %设置滤波器的通带波纹为0.5dB
[z,p,k]=cheb1ap(6,0.5);   %设计Chebyshev I型模拟原型滤波器
[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 I |H(jw)|^2');
title('Chebyshev I 低通原型滤波器(wc=1)')
[bt,at]=lp2hp(b,a,0.8);   %由低通原型滤波器转换为截止频率为0.8的高通滤波器
[Ht,wt]=freqs(bt,at,0:0.01:2);   %计算滤波器的复数频率响应
subplot(2,1,2),plot(wt,abs(Ht).^2);  %绘制传递函数的平方幅频响应
xlabel('w/wc');ylabel('Chebyshev I |H(jw)|^2');
title('Chebyshev I 高通滤波器(wc=0.8)')