FIR-滤波器和-IIR-滤波器的格型-梯形结构.docx

本文讨论全零点格型结构、全极点格型结构以及零极点格型结构4.3.5全零点格型结构1973年，Gray和Markel提出一种新的系统结构形式，即格型结构(latticestructure)<>这是一种很有用的结构，在功率谱估计、语音处理、自适应滤波等方面以得到了广泛的应用。这种结构的优点是，对有限字长效应的敏感度低，且适合递推算法。这种结构有三种形式，即适用于FIR系统的全极点格型结构和适用于IIR系统的全极点和零极点格型结构。下面先介绍图7.10所示的全零点格型结构。其他两种个性结构将留到第4.3节讨论。格型结构是由多个基本单元级联起来的一种极为规范化的结构。图7.11示出其中的第阳极。与FIR滤波器的宜接型结构一样，全零点格型结构也是没有反馈支路的，图7.10全零点格型结构图7.11全零点格型结构的基本单元让我们从一组FIR滤波器的系统函数开始研究全零点格型结构。图7.10中，以x()为输入序列，后接M个格型级，这样就形成M个滤波器：第加(m=1,2,M)个滤波器有两个输出，即上输出力5)和下输出g”,5)。以，5)为输出的滤波器称为前向滤波器；以为输出的滤波器称为后向滤波器。对于M个前向FlR滤波器，它们的系统函数为：Hlfl(z)=Am(zm=1,2,.,M(18)式中，4,(z)是多项式:4(Z)=I+Z4"(%)z”,mM(19)k=这里，为了数学推导的方便，令式子右边第1项为1；下标用代表滤波器序号，也代表滤波器的阶数，例如，给定4(0)=1以及。(I)M(2),(M),则第4个滤波器的系统函数为W4(z)=1+4(l)z,+4(2)z2÷4(3)z-3+d4(4)z4设第团个滤波器的输入、输出序列分别是x()和y()，则y(n)=X(")十Zq,r(%)x5-Z)(21)=其直接型实现如图12所示。x(n)图7.12FlR滤波器的一种直接实现形式m=1阶滤波器的输出可表示为y()=%()+4-1)(22)该输出也可以从图12所示的第一级格型滤波器得到。图中，两个输入端联在一起，激励信号为X5)。从两个输出端得到的信号分别为工()和g5):15)=M)+-0、(23)g1()=x(M)-X(M-I)其次我们考虑二阶FIR滤波器，它的直接型结构输出为y(n)=x(ri)+a2()x(n-l)+a2(2)x(n2)=x(n)x(n-l)x(-2)la1()2(2)t(24)上式将输出y()表示为两个向量的内积，T表示向量转置。相应地，这个二阶滤波器可以用两个级联的格型单元(图10前面的两级)来实现。，图中，第一级的输出为f5)=+kix(n-l),(25)gi(n)=kix(n)-x(n-)25)=E5)+%2g5-Dg2()=&2/l()+g|(T)将式(25)中的工5)代入式(26)中，得f2(n)=x(n)+kxn-1)+k2klx(n-1)+xn-2)=%()+K(1+k2)x(n-1)+k2x(n-2)现在令式(24)和式(27)的系数相等，即O2=k2,%=K(I+&)于是，得二阶格型结构的参数其中，&二%(2)这个结果是很容易理解的。从图7.12看，如果滤波器阶数机=2,则时延为2的输入输出传输值为2(2),而从图7.10看，从输入到上端输出有三条可能的支路,而其中时延为2的支路传输值为占。如果这两个流图等效，则应有G=(2)。因此可以推论，若有加个格型级，则其最右边的支路线与直接型结构的参数册(M)相等：kfn=a<m)为了得到其它支路传输值幻I，Z吁2”“，K与直接型结构的参数之间的关系，我们需要从图7.10所示的M阶格型结构的最右边做起：根据M阶混波器的直接型参数，依次求M-LM-2,M-3,.J阶滤波器的直接型参数。这是降阶递推。只要求出机阶滤波器的系数组4优)«=l,2,.,m,则格型结构的支路传输=4(加)。式(29)表明，二阶格型结构的两个参数仁和期可以根据直接型结构的参数求出。继续这个过程，可以得到一个m阶直接型FlR滤波器和一个机阶或机级格型滤波器之间的等效性。按照图7.10,格型滤波器可用递归方程描述为>5)=go5)=5)(31)fmW=九.|()+先达吁(-1),rn=l,2,.,f-1(32)gtnW=kmfn-5)+gg5-1),m=1,2,.,M-1(33)因此，第M-I级滤波器的输出相当于M-I阶FlR滤波器的输出，即(34)yW=fM_/n)因为FlR滤波器和格型滤波器的输出力.()可以表示为fm(n)=Yjam(k)x(n-k)an(0)=1(35)Ar=O而这个式子是两个序列的卷积和，所以它遵从Z变换关系En(Z)=At(Z)X(Z)现在我们来看二级格型漉波器的另一个输出g2()。由图710得g2(ri)=k2fx(n)+gi(n-1)=2x(h)+kix(n-1)+klx(n-1)+x(n-2)=k2x(n)+kl(l-k2)x(n-1)+x(n2)=a2(2)x(n)+a2(l)x(w-1)+x(n2)=x(n)x(n-l)x(n-2)a2(2)O2(l)1(37)可见，对于g2()为输出的后向滤波器，滤波系数组为%(2)02(l)1，而对于以力()为输出的滤波器，漉波系数组按相反次序排列，为1a2(l)a2(l)o根据以上分析。可见加级格型滤波器的输出gm(n)可以用卷积和形式表示为(38)gfn5)=E<k)x(j-k)K=O式中，滤波系数以(Q与产生输出，的另滤波器有关，只不过操作次序相反。例如，如果7=6,(O)=1,«6(1)=2,4(2)=4,4=7,4(4)=5,4=3,a6(6)=6,则&=1,&=2,4=4,其=7,风=5,风=3,&(O)=60m(k)= %<tn-k),Am5) = 1在Z域中，式(38)变为k= O,.,m(39)GJz)= Bm(z)X(z)(40)BQ)=Gfn(z)X(Z)(41)这里，BjZ)是下输出端相对于输入端的系统函数:Btt,(z)=Yma)z-k(42)A=O因为Sm(k)=afn(m-k),故/nmBm(Z)=Eam(m-k)z"=EaJJ)zi=ZfEam(j)z，k=Qj=OJ=O=Zf4(ZT)<43)这个式子描述前、后向滤波器系统函数之间的关系。现在我们回到式(31)(33)的递推方程组，并把它们变换到Z域，得玲(Z)=GO(Z)=X(Z)(44)工<z)="I(Z)+晨ZTGAl(z),m=1,2,.,M-1(45)G<Z)=kmF,ll(z)+zlGtl(zm=1,2,(46)各式除以X(Z)并利用前面的关系式，可得A)=为=1(47)AZ)=AjnJz)+kmz'iBltlJz),m=1,2,.,M-1(48)BJZ)=KAI+ZTB,I,m=1,2,.,M-I(49)因此，在Z域，一个格型级可用矩阵方程描述为4(Z)Il,T4-(2)&IAlLT4.=(50)利用式(47)(49)可以根据格型滤波器系数，从机=1开始按升阶递推法求出直接型滤波器系数。例给定三级格型滤波器如图13所示。确定与之等效的直接型结构的FIR滤波器系数。图13给定三级格型滤波器解根据式(48),得1(z)= 4(2) + z,(z)=l + z, =l + -z,14因此，对应于单级格型的FlR滤波器系数为4(0)= 1。4(1)= K=1,因&(Z)是4,(Z)的反转多项式，故其次,耳="对于m=2得格型滤波器，根据式(48)得A2(z)= , (z) + z,Bl(z)13_|1 -2=l+-z +-Z82因此,31对应于二级格型的FlR滤波器系数为a2(O)= l, 2(l) = ,2(2) = -o此外 82最后,B2(z)= - + -z-2 + z32 8添上第三个格型级，得出多项式4(Z)= 2(z) + k3zlB2(z)113-15-21-3=1 + z +-Z +-Z2483因此，与给定三级格型滤波器等效的直接型FIR滤波器系数为1351%(0)= L %= 9=;24 o 31351G(O)=1,生=有，。3=-,«3(3)=-Z4oJ假定已知M阶直接型FIR滤波器的系数或者多项式4式z),我们希望确定相应的格型滤波器的系数组化,i=1,2,.,M。对于第M个格型级，可直接得出女M=AM(A/)，所以,只需从M-I开始降阶递推过程。为了得到ZM-,只需求出多项式Aw-.(2)=1+Aw-O)z,+.÷Ai(M-I)Z-(MT)就可以得到ZM-=4WT(M-1)。根据式(48)和式(49),可以得到降阶递推关系：An(Z)=A”(Z)+3Z-RT(Z)=An-I(z)+k”IBKZ)-kmAm,l(Z)J于是，Al(Z)=A"(z)g凡(z),m=M-lM-2,.4(51)1-km例设FIR滤波器的系统函数为/(z)=A,f(z)=l+-ZT+-Z-2+-ZT“2483确定对应于该FIR滤波器的格型系数。解首先，直接得出左3=。3(3)=g,而且r*/、15T13-2-3R(z)=-+-z+-z+z33824在机=4的情况下，利用式(51)降阶递推，得(¾-y=1+3z-+lz-221-片82113因此，女2=生(2)=2和与(2)=5+WZ'÷Zo最后，在相=2的情况下，再降阶递推，得A2一NS(Z)1-1A(Z)=2:2=l+-z个1-好4因此，K=A(I)=I图13示出所得三级格型滤波器。4.3.5IlR系统的全极点格型结构IIR滤波器的全极点系统函数H(Z)为“(Z)=M(12)i-z-A=I与M阶FIR系统函数相比较，可见这两种系统互为逆系统。我们在第节以研究了FlR系统的(全零点)格型结构。现在我们要基于式(12)找出IIR系统的全极点格型结构。最简单的途径就是研究逆系统的信号流图，从中找出规律。给定一阶FIR系统函数为H(z)=-=c+k,z-i(13)X(Z),则差分方程为y(n)=cx(n)+kix(n-1)(14)图19是相应的信号流图工(月)O->9+oJ包拳自OZT图19一阶FIR系统逆系统的系统函数为"'(Z)=Y(Z)_1X(z) c + kiz-l(15)其差分方程为(16)y()=-x()-k、y(n-1)图20示出相应的信号流图。1cX(H)OHpk-OJ色)图20一阶FIR系统的逆系统所以，可以按照图21所示的中间步骤从原系统得到逆系统：