ms18_2.m

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

%巴特沃斯模拟滤波器
fp=2000;Omgp=2*pi*fp;   %输入滤波器的通带截止频率
fs=5000;Omgs=2*pi*fs;   %输入滤波器的阻带截止频率
Rp=1;As=20;             %输入滤波器的通阻带衰减指标
%计算滤波器的阶数和3dB截止频率
[n,Omgc]=buttord(Omgp,Omgs,Rp,As,'s')
%计算n阶模拟低通原型，得到左半平面零极点
[z0,p0,k0]=buttap(n);   %由滤波器阶数N求模拟滤波器原型
b0=k0*real(poly(z0))     %求滤波器系数b0
a0=real(poly(p0))        %求滤波器系数a0
freqs(b0,a0);            %求系统的频率特性