《FIR-滤波器和-IIR-滤波器的格型-梯形结构.docx》由会员分享,可在线阅读,更多相关《FIR-滤波器和-IIR-滤波器的格型-梯形结构.docx(14页珍藏版)》请在优知文库上搜索。
1、本文讨论全零点格型结构、全极点格型结构以及零极点格型结构4.3.5全零点格型结构1973年,Gray和Markel提出一种新的系统结构形式,即格型结构(latticestructure)这是一种很有用的结构,在功率谱估计、语音处理、自适应滤波等方面以得到了广泛的应用。这种结构的优点是,对有限字长效应的敏感度低,且适合递推算法。这种结构有三种形式,即适用于FIR系统的全极点格型结构和适用于IIR系统的全极点和零极点格型结构。下面先介绍图7.10所示的全零点格型结构。其他两种个性结构将留到第4.3节讨论。格型结构是由多个基本单元级联起来的一种极为规范化的结构。图7.11示出其中的第阳极。与FIR滤
2、波器的宜接型结构一样,全零点格型结构也是没有反馈支路的,图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项为
3、1;下标用代表滤波器序号,也代表滤波器的阶数,例如,给定4(0)=1以及。(I)M(2),(M),则第4个滤波器的系统函数为W4(z)=1+4(l)z,+4(2)z24(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)
4、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)
5、=%()+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=a5)=go5)=5)(31)fmW=九.|()+先达吁(-1),rn=l,2,.,f-1(32)gtnW=kmfn-5)
6、+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
7、)+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)=Ek)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
8、-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)G9+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所示的中间步骤从原系统得到逆系统: