samp5_6.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 14 行

%Samp5_6
Rp=2;Rs=30;   %模拟原型滤波器的通带波纹为2dB,阻带衰减为30dB。
[z,p,k]=ellipap(4,Rp,Rs);   %设计椭圆滤波器
[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('椭圆 |H(jw)|^2');
title('原型低通椭圆滤波器(wc=1)')
[bt,at]=lp2lp(b,a,0.5);  %将模拟原型低通滤波器的截止频率变换为0.5
[Ht,wt]=freqs(bt,at,0:0.01:2);  %给出复数频率响应
subplot(2,1,2),plot(wt,abs(Ht).^2); %绘出平方幅频函数
xlabel('w/wc');ylabel('椭圆 |H(jw)|^2');
title('低通椭圆滤波器(wc=0.5)')