IIR数字滤波器设计原理.docxVIP

IIR 数字滤波器设计原理 利用双线性变换 设计 IIR滤波器（只介绍 巴特沃斯数字低通滤波器 的设计），首先要设计出满足指标要求的模拟滤波器的传递函数 H a (s) ，然后由 H a ( s) 通过双线性变换可得所要设计的 IIR 滤波器的系统函数 H (z) 。 如果给定的指标为数字滤波器的指标， 则首先要转换成模拟滤波器的技术指 标，这里主要是边界频率 w p和 ws 的转换，对 p 和 s 指标不作变化。边界频率的 2 tan( 1 w) 转换关系为 T 2 。接着，按照模拟低通滤波器的技术指标根据相应 设计公式求出滤波器的阶数 N 和 3dB 截止频率 c ；根据阶数 N 查巴特沃斯归一 s 化低通滤波器参数表， 得到归一化传输函数 H a ( p) ；最后，将 p c 代入 H a ( p) 去归一，得到实际的模拟滤波器传输函数 H a (s) 。之后，通过双线性变换法转换 2 1 z 1 s z 1 ，得到所要设计的 IIR 滤波器的系统函数 H ( z) 。 公式T1 步骤及内容 用双线性变换法设计一个巴特沃斯 IIR 低通数字滤波器。设计指标参数为：在通带内频率低于 0.2 时，最大衰减小于 1dB ；在阻带内 0.3 , 频率区间上，最小衰减大于 15dB 。 以 0.02 为采样间隔，绘制出数字滤波器在频率区间 0, / 2 上的幅频响应特性曲 线。 程序及图形 程序及实验结果如下： %%%%%%%%%%%%%%%%%% % %lskyp %%%%%%%%%%%%%%%%%% rp=1;rs=15; wp=.2*pi;ws=.3*pi; wap=tan(wp/2);was=tan(ws/2); [n,wn]=buttord(wap,was,rp,rs,s); [z,p,k]=buttap(n); [bp,ap]=zp2tf(z,p,k); [bs,as]=lp2lp(bp,ap,wap); [bz,az]=bilinear(bs,as,.5); [h,f]=freqz(bz,az,256,1); plot(f,abs(h)); title( 双线性 z 变换法获得数字低通滤波器，归一化频率轴 ); xlabel(\omega/2\pi); ylabel(低通滤波器的幅频相应 );grid; figure; [h,f]=freqz(bz,az,256,100); ff=2*pi*f/100; absh=abs(h); plot(ff(1:128),absh(1:128)); title( 双线性 z 变换法获得数字低通滤波器 ,频率轴取 [0,\pi/2]); xlabel(\omega); ylabel(低通滤波器的幅频相应 );grid on; 运行结果： 双 线 性 z 变 换 法 获 得 数 字 低 通 滤 波 器 ， 归 一 化 频 率 轴 1 0.9 0.8 0.7 应 相 频 0.6 幅 的 0.5 器 波 0.4 滤 通 低 0.3 0.2 0.1 0 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 /2 双 线 性 z 变 换 法 获 得 数 字 低 通 滤 波 器 , 频 率 轴 取 [0, /2] 1 0.9 0.8 应 0.7 相 频 0.6 幅 的 0.5 器 波 0.4 滤 通 低 0.3 0.2 0.1 0 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6