ms18_7.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 12 行

%椭圆模拟滤波器原型

n=input('N= ');             %由使用者输入滤波器阶数N
Rp=input('Rp= ');           %由使用者输入滤波器通带衰减指标
As=input('As= ');           %由使用者输入滤波器阻带衰减指标
[z0,p0,k0]=ellipap(n,Rp,As);%计算n阶模拟低通原型，得到左半平面零极点
b0=k0*real(poly(z0))        %求滤波器系数b0
a0=real(poly(p0))           %求滤波器系数a0
freqs(b0,a0);               %求系统的频率特性
Pb=poly2str(b0,'s')         %给出b0决定的关于s多项式
Pa=poly2str(a0,'s')         %给出a0决定的关于s多项式