ms18_1_1.m

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

%巴特沃斯模拟滤波器原型

for n=1:8;              %由使用者输入滤波器阶数N
%计算n阶模拟低通原型，得到左半平面零极点
[z0,p0,k0]=buttap(n);   %由滤波器阶数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,1.1]),
title('巴特沃斯模拟滤波器原型幅频特性(N=1:8)');
end