数字信号处理-实验-无限冲激响应(IIR)数字滤波器的三种结构.docx

《数字信号处理-实验-无限冲激响应(IIR)数字滤波器的三种结构.docx》由会员分享，可在线阅读，更多相关《数字信号处理-实验-无限冲激响应(IIR)数字滤波器的三种结构.docx（7页珍藏版）》请在优知文库上搜索。

1、数字信号处理实验报告实验名称：无限冲激响应(IIR)数字滤波器的三种结构学号：姓名:评语：成绩：一、实验目的1、掌握HR滤波器的直接II型、级联型和并联型三种结构的根本原理和特点。2、掌握利用MATLAB实现IIR滤波器的三种结构的程序设计方法，并能够进行三者之间的相互转换。3、掌握滤波器频响特性的绘制方法。二、实验原理与计算方法按照结构划分数字滤波器，有递归式和非递归式两种。递归式数字滤波器的差分方程为jVjWy()+Z)(一口=Zbix(i)(1)k=lr=l其中至少有一个40.非递归式数字滤波器的差分方程为M)=Z3(-i)(2)i=l可以看出递归式数字滤波器的响应)，()不仅与鼓励x(

2、)有关，而且与以前的输出信号M-灯有关；而非递归式数字滤波器的响应y()仅只与鼓励x()有关。按照单位样值响应划分数字滤波器，那么有无限冲激响应(IIR)和有限冲激响应(FlR)之分。IIR滤波器是递归式的，差分方程如I(I)式所示，FIR滤波器一般是非递归式的，差分方程如1(2)式所示。HR滤波器常用的典型结构有直接II型、级联型和并联型，分别介绍如下：1、直接Il型(也称为正准型结构)根据(D式，HR滤波器的传输函数为NSbkZ-kH(Z)=T其中已假设(1)式中的N=M,对于其它情况，那么可令相应的某些系数为零。N1令HKZ)=Ebicz-k,H2(Z)=-N“力1工ek=那么有”(Z)

3、=Hl(Z)”2(Z)(4)由此可以得到相应的时域中鼓励/()与响应M)之间的关系为N为S)=Zak为(-&)+M)其中为()是与式中的“2(Z)相应的中间函数序歹U。x(n)y2(n) bo y(tn)图6-1直接II型结构信号流由式确定的直接II型的信号流图如图6-1所示，其中将中间的两条延时链合并为一条，实际的信号流将按(5)式分成两个延时链独立运行。编程时，用三个数组分别存放系数ak,bk和y2(n-k,k=l,2,N0由图可以看出，乃()沿中间的延时链自上向下传播过程中将逐级向右移位，每一级2(-幻向左边与4.的乘积按A=1,2,N累加，再和()相加，得到的力()沿中间的延时链又自上

4、向下传播逐级向右移位，每一级以(-A)向右边与电的乘积按A=O,12,N累加，其结果就是对应于鼓励x()的响应y(n)o然后y2(n)向右移位一个单位时间，输入鼓励x()计算下一个时刻的响应y(n)o直接II型结构具有简单直观的典型网络结构形式，在计算机上很容易实现。但是它对系数幺,“的量值变化比拟敏感，4，%直接确定了系统零、极点的位置，从而影响到系统的性能。尤其当阶数N较高时，系统对系数的字长效应很敏感，产生的误差也较大。2、级联型结构由于当直接IT型结构传输函数阶数增加时，系数量化引起的误差影响到滤波器的性能，因此要采用其它形式的结构。将(3)式的传输函数分子和分母进行因式分解，即用它的

5、零、极点表示为%M2n(izT)n(厢+%z-z-2)H(Z)=2jj(l-pz-,)(l+z-,+z-2)k=*=1其中N+2N2=N,+22=N,分子和分母中的实系数二阶因子分别对应于共聊零、极点。可以将M 9 I 112 M“FA产二n /A= l + %z +a2kz k= /八 +z,+z2H A(Z) = q-l + htz +a2kz上式分子和5母中单根一阶因子作为二阶因子的一种特例，那么(6)式就可以表示为必个实系数二阶根本节级联的形式(7)(8),、LA卜J.flMr1,IMr,FSMr,历MMn) () 即im(h) Zom n)图6-2级联型结构信号流其中(Z)为滤波器的

6、二阶根本节，M=为取整结果。小2&，力0人-，仇A，么A那么分别是第个根本节的分母、分子系数，它们只关系到滤波器的某一对零、极点。级联型结构如图6-2所示，它是一系列二阶根本节“(z)的级联，每一个Hjt(Z)可以用直接II型结构实现。级联型结构的特点是对滤波器性能的调整比拟方便，调整系数比t2A,只单独涉及到第4级零、极点，而不会影响到其它任i级的零、极点，因而可以独立地控制滤波器的各零、极点的分布。3、并联型结构这种结构将传输函数H(Z)展开为局部分式，即表示为假设干一阶和二阶根本节网络与一个常数B0之和+War+Z”砥Lk=I-PkZA=I1-MZ-Oikz其中必+2=N,同样也可以统一

7、表示为二阶根本节的形式,.+z,、W(z)=B0y上产_-=B.+YHk(z)(10)A=Iakza2kzA=IBo图6-3并联型结构信号流并联型结构信号流如图6-3所示，其中二阶根本节网络可以用直接II型结构实现，程序设计也可参考直接型II结构的方法。并联型结构也可以单独调整极点位置，但却不能象级联型结构那样直接控制零点的分布。因为并联型结构各二阶根本节网络的零点并不是整个系统函数的零点。因此，当要准确传输零点时，以采用级联型结构为宜。不过，由于并联型根本节之间互不影响，所以运算误差比级联型的要小一些。通常，IIR滤波器用系统函数的有理式(直接结构形式)描述，MataIb中提供相应的函数可以

8、把直接型结构转换成级联和并联型结构。-MATLAB,直接型结构由两个行向量描述，b包含系数仇,a包含系数公,可由filter函数实现。要由的直接型结构实现级联型结构，可利用函数dir2cas根据直接型的系数6和公得到系数bo,&,和A。,该函数把矢量b和a转换成KX3维矩阵B和首先计算加，它等于ZWao,(ol),然后通过给较短的矢量添零，使矢量b和a一样长。这就保证每一个双二阶环节的分子和分母均不为零。接着计算多项式8(z)和4z)的根，用Cplxpair函数把这些根以共规复根对的次序排列，最后用poly函数把每一对根再转换成二阶分子或分母多项式。级联型由CaSfntr函数实现，它把每个二阶

9、环节的系数存放在矩阵8和A中，放在一个循环里，而在此循环中采用filter函数。输入被乘以加，每个滤波器的输出作为下一级滤波器的输入，最后一个滤波器的输出即为总的输出。函数cas2dir可以把级联形式转换成直接形式。这是一个包含几次多项式乘法的简单运算，在一个K次的循环中实用ConV函数。类似的，利用dir2par函数可以把直接型系数6和诙转换成并联型系数&,和A6,并在parfiltr函数中使用这些系数以实现并联形式。而par2dir函数那么可用来实现并联形式转换成直接形式。三种结构的IIR滤波器均可通过以下方式计算单位抽样响应。delta=impseq(0,0,7);hdir=f11ter

10、(b,a,delta)%直接型结构hcas=casfiltr(bO,B,A,delta)%直接型结构实现级联型结构hpar=parfiltr(C,B,A,delta)%直接型结构实现并联型结构传输函数为H(Z)的滤波器的频率特性为”(Cm)=(z)I：=/。由于我们所讨论的滤波器都是稳定系统，其所有极点都在Z平面的单位圆内，单位样值响应力()是一个无限长的衰减序列，所以要截取有限长的一段来作频率特性分析。截取的长度越长，逼近的程度越高。如果截取长度为M滤波器的频率特性由Fourier变换得N-IN-N-(11)(12)(13)H(ej)h(n)ejn=h(n)cosn-jh(n)snnW=O1

11、1=071=0那么滤波器的幅频特性和相频特性可表示为COS。r0NT/Zh(n)sinnM=ON-I，()Sinn-arctanh(n)cosn=0三、实验内容设三阶滤波器的传输函数为W1+0.5814zl+0.2114z2H(z)=：；r1-03984z-,+0.2475z2-0.04322z-3鼓励信号为X1(7?)=eancosn0Rfj(n)a=0.5,0=,N=40x2W=3()-20W20根据所给定的滤波器系统函数和参数值，用三种结构实现该IIR数字滤波器结构，分别使鼓励阳()和必()通过该滤波器，求出相应于不同鼓励的响应，作出响应信号的图形和幅频、相频特性曲线，并判断该滤波器为何

12、种滤波器(低通、高通、带通、带阻)。1、1()；实验代码：b=l0.58140.21140a=l-0.39840.2475-0.04322delta=impseq(0,-10,10)hdir=filter(b,a,delta)b,B,A=dir2cas(b,a)hcas=casfiltr(b,B,A,delta)C,B,A=dir2par(b,a)hpar=arfiltr(C,B,A,delta)n=0:40x1=exp(-0.5.*n).*cos(n.*(pi./2)yl=filter(b,a,xl)y2=casfiltr(b,B,A,xl)y3=parfi1tr(C,B,A,xl)subp

13、lot(4t3,1);stem(xl,.,)；title(,x)subplot(4,3,4);plot(yl);titleC直接型)subplot(4,3,5);plot(angle(yl);titleC直接型相频特性)subplot(4,3,6);plot(abs(yl);titleC直接型幅频特性)subplot(4,3,7);plot(y2);title(级联型)subplot(4,3,8);plot(angle(y2);title(级联型相频特性)subplot(4,3,9);plot(abs(y2);title(级联型幅频特性)subplot(4,3,10);plot(y3);tit

14、le(,并联型)subplot(4,3,11);plot(angle(y3);title(,并联型相频特性)subplot(4,3,12);plot(abs(y3);titleC并联型幅频特性)实验截图：2、x2(n)实验公码：b=l0.58140.21140a=l-0.39840.2475-0.04322delta=impseq(0,-10,10)hdir=filter(b,a,delta)b,B,A=dir2cas(bta)hcas=casfiltr(b,B,A,delta)C,B,A=dir2par(b,a)hpar=arfiltr(C,B,A,delta)n=-20:20n,x2=impseq(0,-20,20)yl=filter(b,a,x2)y2=casfiltr(b,B,A,x2)y3=parfiItr(C,B,A,x2)subplot(4,3,1);stem(x2,.,);title(,x2,)subplot(4,3,4);plot(yl);titleC直接型)subplot(4,3,5);plot(angle(yl);titleC直接型相频特性)subplot(4,3,6)：plot(abs(yl);title(直接型幅频