实验四 IIR数字滤波器设计.docx

《实验四 IIR数字滤波器设计.docx》由会员分享，可在线阅读，更多相关《实验四 IIR数字滤波器设计.docx（5页珍藏版）》请在第一文库网上搜索。

1、试验四IIR数字滤波器的设计与MATLAB实现一、试验目的：1、要求把握R数字滤波器的设计原理、方法、步骤。2、能够依据滤波器设计指标进行滤波器设计。3、把握数字巴特沃斯滤波器和数字切比雪夫滤波器的设计原理和步骤。二、试验原理：R数字滤波器的设计方法：频率变换法、数字域直接设计以及计算机帮助等。这里只介绍频率变换法。由模拟低通滤波器到数字低通滤波器的转换，基本设计过程：1、将数字滤波器的设计指标转换为模拟滤波器指标2、设计模拟滤波器G（S）3、将G（S）转换为数字滤波器H （Z）在低通滤波器设计基础上，可以得到数字高通、带通、带阻滤波器的设计流程如下：1、给定数字滤波器的设计要求（高通、带通、

2、带阻）2、转换为模拟（高通、带通、带阻）滤波器的技术指标3、转换为模拟低通滤波器的指标4、设计得到满意3步骤中要求的低通滤波器传递函数5、通过频率转换得到模拟（高通、带通、带阻）滤波器6、变换为数字（高通、带通、带阻）滤波器三、标准数字滤波器设计函数MATLAB供应了一组标准的数字滤波器设计函数，大大简化了滤波器设计过程。1 butter例题1设计一个5阶Butterworth数字高通滤波器，阻带截止频率为250Hz ,设采样频率为IKHz.I k hJ-h t er (5. 250/500. high)L z, , kJ but i er(5t 250 500, , i glt)f r eq

3、z (b15 I 2, I 000)50100150200250300350400450500Frequency (Hz)oo o oopo1 ? 3in3 3w=e250100150200250300350400450500Fequetcy (Hz)- A A oo o o o o oo o 力 o o1 -23 utter(5,250/500, high)：frm (b a 512,1000)100150200250300350400450500Frequency (H(S3 苔 llq图27阶chebyshevll型数字低通滤波器四、冲激响应不变法一般来说,在要求时域冲激响应能仿照模拟滤

4、波器的场合,一般使用该方法。冲激响应不变法一个重要的特点是频率坐标的变换时线性的，因此假如模拟滤波器的频响带限于折叠频率的话，则通过变换后滤波器的频率响应可不失真的反映原响应与频率的关系。例题3设计一个中心频率为500Hz,带宽为600 Hz的数字带通滤波器，采样频率为iKHz。解：源代码如下：zlp k=buttap(3);b, azp2tf (zl p, k);bttat=lp2bp(b,av 500*2*pi,600*2*pi);bz, azimpinvar (bt, at, 1000)； %将模拟滤波器变换成数字滤波器frez(bz.az. 512. * whole*. 1000)p)

5、 3pn=uww1002003004005006007008009001000Frequency (Hz)1002003004005006007008009001000Frequency (Hz)o o o o oo o o o2 1 12se60p) =seqd图3数字带通滤波器五、双线性变换法与冲激响应不变法比较，双线性变换法的主要优点是靠频率的非线性关系得到S平面与Z平面的单值一一对应关系，整个值对应于单位圆一周，所以从模拟传递函数可直接通过代数置换得到数字滤波器的传递函数。例题4设计一个截止频率为200 Hz的数字低通滤波器，采样频率为IK Hzo解，源代码如下工zt pf kJ h

6、t t p(3) ；ht n=z2tf(zt , k);ht# at = lp21p(K a, 200*2*pi)；bzv az =bi linear (bti.1000);freqz (bz azs 512, 1000)结果如图64所示:0 07530 22590. 22590. 0753az =O. 515 1-O 0865L OOOO-O 266例题5基于Butterworth模拟滤波器原型，使用双线性转换设计数字滤波器，其中参数指标为：通带截止频率p=O.25通带波动值：Rp=ldB阻带截止频率Qs=0.311；阻带波动值：s=15dB0解：先确定滤波器阶数N,同时依据Qp和Qs确定Q

7、c=0.5U,接着使用bilinear进行双线性变换，最终绘制在频域上的各种图像，其源代码如下：外数字滤波泮指标wp=O.2j*piws=0.3*pi;Rp=l;As= 15;%将数字港波制指标反转变化为模拟滋波器的参数T=l;fs= I/T;omegap=(2T)*tan( wp/T);omegas=(2/T)*tan(ws/T);epsqrc( I OA(Rp/l O)-1);Ripp)e = sqrt( 1/( 1 +ep*ep);Attn= 1/(1 O*(As/2O); butterworth原型模拟滤波器的设计cs,dsJ=afd-butt(onegap.omegas,Rp.As)

8、%双线性变换ba= bilinear(cswds,Twfs);%频域图像的绘制freqx(ba)程序运行后，产生4阶Butterwor由数字滤波器，频率响应如图所示：20OJ 020 3040 506070 8091N(idlzed Fequency ( ample)0204060p- i-s0 10.20.3040.50.60.70809Nu3zed Fequfcy ( . djple)o o Oo o O3 2 1-swslE)oSBLd六、编程练习1、基于chebyshevl型模拟滤波器原型使用冲激不变转换方法设计数字滤波器,要求具有下面参数指标：通带截止频率p=O.25通带波动值：Rp=ldB阻带截止频率Qs=O.3阻带波动值：As=15dBo2、一个椭圆数字滤波器的设计，要求采纳双线性变换方法，指标参数如下：通带截止频率Qp=0.21L通带波动值：Rp=1.5dB阻带截止频率s=O.35阻带波动值：As=20dBo