ms18_7_1.m

《数字信号处理实验(MATLAB版)》这

%椭圆模拟滤波器原型

for n=1:5;               %由使用者输入滤波器阶数N
[z0,p0,k0]=ellipap(n,1,20);%计算n阶模拟低通原型，得到左半平面零极点
b0=k0*real(poly(z0))     %求滤波器系数b0
a0=real(poly(p0))        %求滤波器系数a0
[h,w]=freqs(b0,a0);      %显示系统的频率特性
plot(w,abs(h),'K'),hold on
axis([0,5,-0.02,1.1]),
title('椭圆模拟滤波器原型幅频特性(N=1:5)');
end