数字信号处理-3(IIR补充).ppt

IIR数字滤波器设计步骤 巴特沃斯滤波器设计 (1)求滤波器的最小阶数N [N,Wc]=buttord(Wp,Ws,Rp,Rs,’s’) 根据滤波器设计指标计算巴特沃斯模拟滤波器的阶数N及频率参数Wc。 Wp：通带截止频率 Ws：阻带截止频率 Rp：通带衰减（分贝） Rs：阻带衰减（分贝） N：滤波器的阶数 Wc：3dB衰减时的截止频率 ‘s’：表示这是计算模拟滤波器的参数 巴特沃斯滤波器设计 (2)设计模拟低通滤波器的原型 [z,p,k]=buttap(N) 根据最小阶数N计算滤波器原型的零极点及增益。 巴特沃斯滤波器设计 (3)频率转换（变为高通、带通、带阻等） [Bt,At]=lp2lp(B,A,Wc) 把单位截止频率的模拟低通滤波器系数B、A 转换为另一截止频率Wc的低通滤波器系数Bt和At B、A及Bt、At 的具体含义为： 巴特沃斯滤波器设计 (4)将模拟滤波器转换为数字滤波器 [Bd,Ad]=bilinear(B,A,Fs) 把模拟滤波器系数B、A变为近似等价的数字滤波器系数Bd、Ad，它使用的是双线性变换法。 Fs是采样频率，变换公式为： 巴特沃斯滤波器设计 (5)把2、3、4步和为一步 [B,A]=butter(N,Wc) 设计低通数字滤波器系数B、A，它使用的是双线性变换法。 [B,A]=butter(N,Wc,’high’) 设计高通数字滤波器系数B、A，它同样使用的是双线性变换法。 [B,A]=butter(N,Wc,’stop’) 设计低通数字滤波器系数B、A，它同样使用的是双线性变换法。 巴特沃斯滤波器设计－－示例 function [b,a]=IIR(fs,Fp,Rp,Fs,Rs) Rs=40;Fs=1500; Rp=0.5;Fp=1000; fs=8000; Wp=2*pi*Fp;Ws=2*pi*Fs; [N,Wn]=buttord(Wp,Ws,Rp,Rs,s); [b1,a1]=butter(N,Wn,s); [b,a]=bilinear(b1,a1,fs); [H,W]=freqz(b,a); plot(W*fs/(2*pi),abs(H));grid on;title(频率响应); xlabel(频率);ylabel(幅值); 巴特沃斯滤波器设计－－示例结果 巴特沃斯滤波器设计－－示例结果 * * 设计指标 fs：采样频率 Fp,Rp：通带参数 Fs,Rs：阻带参数 模拟低通滤波器原型设计 Buttap,cheb1ap,cheb2ap ellipap 频率变换 lp2lp,lp2hp,lp2bp,lp2bs 模拟数字转换 bilinear impinvar 求最小阶数N buttord,cheb1ord, cheb2ord,ellipord 合为一步的设计函数 Butter,cheb1,cheb2,ellip 滤波器 系数B,A N,Wc