数字信号处理-格型网络结构

第5章时域离散系统的网络结构5.7格型网络结构5.7.1全零点格型网络结构1.全零点格型网络结构的流图如图5.7.1所示。该流图只有直通通路，没有反馈回路，因此可称为FIR格型网络结构。观察该图，它可以看成是由图5.7.2的基本单元级联而成。第5章时域离散系统的网络结构图5.7.1全零点格型网络结构第5章时域离散系统的网络结构图5.7.2基本单元第5章时域离散系统的网络结构将上式进行Z变换，得到按照图5.7.2写出差分方程如下：llllknrnene)1()()(11(5.7.1)(5.7.2))1()()(11nrknenrllllllllkzRzzEzE)()()(111)()()(111zRzkzEzRllll(5.7.3)(5.7.4)第5章时域离散系统的网络结构再将上式写成矩阵形式(5.7.5)将N个基本单元级联后，得到:(5.7.6))()(1)()(1111zRzEzkkzzRzEllllll)()(111)()(001111111111zRzEzkkzzkkzzkkzzRzENNNNNN第5章时域离散系统的网络结构令Y(z)=EN(z)，X(z)=E0(z)=R0(z),其输出为由上式得到全零点格型网络的系统函数为只要知道格型网络的系数kl，l=1,2,3,…,N,由上式可以直接求出FIR)(11101)()(01)(111zXzkkzzRzEzYNlllNN（5.7.7)（5.7.8)11101)()()(111NlllzkkzzXzYzH第5章时域离散系统的网络结构2.由FIR直接型网络结构转换成全零点格型网络结构假设N阶FIR型网络结构的系统函数为(5.7.9)式中,h(0)=1;h(n)是FIR网络的单位脉冲响应。令ak=h(k)，得到:(5.7.10)式中，a0=h(0)=1;kl为全零点格型网络的系数,l=1,2,…,NNnnznhzH0)()(NkkkzazH0)(第5章时域离散系统的网络结构转换公式(5.7.11)(5.7.12)(5.7.13)式中,l=N,N－1,…,1kNkaa)()(lllak)1(,,3,2,112)()()1(lkkakaallklllklkNkkkzazH0)(第5章时域离散系统的网络结构【例5.7.1】将下面三阶FIR系统函数Ｈ3(z)转换成格型网络，要求画出该FIR直接型结构和相应的格型网络结构流图。3213576.064.09.01)(zzzzH解例题中N=3,按照(5.7.11)式，有由(5.7.12)式，得到:按照(5.7.13)式，递推得到:576.0,64.0,9.0)3(3)3(2)3(1aaa576.0)3(33ak第5章时域离散系统的网络结构l=3,k=1时,l=3,k=2时,45182795.0)576.0(164.0576.09.01223)3(23)3(1)2(1kakaa91974181.0)576.0(19.0576.064.01223)3(13)3(2)2(2kakaa91974181.0)2(22ak第5章时域离散系统的网络结构l=2,k=1时,最后按照算出的格型结构的系数，画出三阶FIR直接型47757672.0)91974181.0(1)91974181.01(45182795.01222)2(12)2(1)1(1kakaa47757672.01k第5章时域离散系统的网络结构第5章时域离散系统的网络结构实际上，调用MATLAB函数实现直接型网络结构与格型网络结构之间的相互转换非常容易。tf2latc实现直接型到格型结构变换，latc2tf实现格型到直接结型结构K=tf2latc(hn):求FIR格型结构的系数向量K=［k1,k2,…,kN］,hn为FIR滤波器的单位脉冲响应向量，并关于hn(1)=h(0)归一化。应当注意，当FIR系统函数在单位第5章时域离散系统的网络结构hn=latc2tf(K)将FIR格型结构转换为FIR直接型结构。K为FIR格型结构的系数向量K，hn为FIR滤波器的单位脉冲响应向量，即FIR直接型结构系数向量。例5.7.1hn=［1,－0.9,0.64,－0.576］;K=tf2latc(hn)K=［－0.67280.1820－0.5760第5章时域离散系统的网络结构5.7.2全极点IIR系统的系统函数用下式表示:(5.7.14)(5.7.15)式中,A(z)是FIR系统，因此全极点IIR系统H(z)是FIR系统A(z)假设系统的输入和输出分别用x(n)、y(n)表示，由(5.7.17)式得到全极点IIR滤波器的差分方程为)(111)(1zAzazHNkkkNkkkzazA11)(第5章时域离散系统的网络结构如果将x(n)、y(n)的作用相互交换，差分方程则变成下式:则(5.7.17))()()(1nxknyanyNkk)()()(1nyknxanxNkkNkkknxanxny1)()()((5.7.16))(111)(1zAzazHNkkkNkkkzazA11)(第5章时域离散系统的网络结构图5.7.4全极点IIR系统的直接型结构第5章时域离散系统的网络结构基于上面的事实，我们将FIR格型结构通过交换公式中的输入输出作用，形成它的逆系统，即全极点格型IIR再将FIR格型结构的基本公式(5.7.1)、(5.7.2)重写如下：（5.7.18）（5.7.19）)()(),()(0nenynenxNllllknrnene)1()()(11)1()()(11nrknenrlllldefdef第5章时域离散系统的网络结构由于重新定义了输入输出，将el(n)按降序运算，rl(n)不变，即（5.7.20）（5.7.21）（5.7.22）（5.7.23）)()(nenxNllllknrnene)1()()(11－＝1,,1,NNl)1()()(11nrknenrllll1,,1,NNl)()()(00nrneny第5章时域离散系统的网络结构图5.7.5全极点IIR格型结构第5章时域离散系统的网络结构令N=1,得到方程为（5.7.24）（5.7.25)（5.7.26）（5.7.27）)()(1nenxlknrnene)1()()(010－＝)1()()(001nrknenrl)1()()()(10nyknxneny第5章时域离散系统的网络结构单极点和双极点IIR格型网络结构第5章时域离散系统的网络结构全极点网络可以由全零点格型网络形成，可以归纳出下（1）将输入到输出的无延时通路全部反向，并将该通路的常数支路增益变成原常数的倒数（此处为1）；（2）将指向这条新通路的各节点的其它节点的支路增益乘以－1（3）第5章时域离散系统的网络结构调用MATLAB转换函数可以实现全极点系统的直接K=tf2latc(1,A):求IIR全极点系统格型结构的系数向量K，A为IIR全极点系统函数的分母多项式A（z）的系数第5章时域离散系统的网络结构［B，A］=latc2tf(K,′allpole′):将IIR全极点系统格型结构转换为直接型结构。K为IIR全极点系统格型结构的系数向量，A为IIR全极点系统系数函数的分母多项式A（z）的系数向量。显然这时分子为常数1，所以B=1。321z31z85z24131A(Z)则求IIR全极点系统格型结构系数向量KA=［1,13/24,5/8,1/3］;K=tf2latc(1,A)第5章时域离散系统的网络结构全极点格型网络同样存在稳定问题，可以证明稳定的充分必要条件是|kl|≤1,l=1,2,…,N第5章时域离散系统的网络结构1.已知系统用下面差分方程描述:试分别画出系统的直接型、级联型和并联型结构。式中x(n)和y(n)2．设数字滤波器的差分方程为)1(31)()2(81)1(43)(nxnxnynyny＋－＝)2(41)1(31)1()()(nynynxnxny第5章时域离散系统的网络结构3.设系统的差分方程为式中,|a|1，|b|1,试画出系统的直接型、级联型和并联型结构。x(n)和y(n)分别表示系统的输入和输出信号。4.设系统的系统函数为试画出各种可能的级联型结构，并指出哪一种最好。abnxbanxnabynybany)1()()2()2()1()()()81.09.01)(5.01()414.11)(1(4)(211211zzzzzzzH第5章时域离散系统的网络结构5．题5图中画出了四个系统，试用各子系统的单位脉冲响应分别表示各总系统的单位脉冲响应，并求6．题6图中画出了10种不同的流图，试分别写出第5章时域离散系统的网络结构题5图第5章时域离散系统的网络结构题6图第5章时域离散系统的网络结构7.假设滤波器的单位脉冲响应为求出滤波器的系统函数，并画出它的直接型结构。8.已知系统的单位脉冲响应为试写出系统的系统函数，并画出它的直接型结构。9.已知FIR试画出该滤波器的直接型结构和线性相位结构。10)()(anuanhn)5(5.0)3(5.2)2(3.0)1(2)()(nnnnnnh)9.01.29.01(101)(4321zzzzzH第5章时域离散系统的网络结构10．已知FIR滤波器的单位脉冲响应为:(1)N=6h(0)=h(5)=15h(1)=h(4)=2h(2)=h(3)=3(2)N=7h(0)=h(6)=3h(1)=－h(5)=－2h(2)=－h(4)=1h(3)=0试画出它们的线性相位型结构图，并分别说明它们的第5章时域离散系统的网络结构11.已知FIR滤波器的16个频率采样值为：H(0)=12，H(3)~H(13)=0H(1)=－3－j，H(14)=1－jH(2)=1＋j，H(15)=－3＋j试画出其频率采样结构，选择r=1，可以用复数乘法12.已知FIR滤波器系统函数在单位圆上16个等间隔采样点为:H(0)=12H(3)~H(13)=0H(1)=－3－jH(14)=1－jH(2)=1+jH(15)=－3＋j3333第5章时域离散系统的网络结构13．已知FIR滤波器的单位脉冲响应为试用频率采样结构实现该滤波器。设采样点数N=5，要求画出频率采样网络结构，写出滤波器参数的计算14.3211864.0414.16.01)(zzzzH)4()1()()(nnnnh3212898.09.098.01)(zzzzH)(/)()(213zHzHzH第5章时域离散系统的网络结构15．写出题15图中系统的系统函数和单位脉冲响应。16.画出题15图中系统的转置结构，并验证两者具有17.用b1和b2确定a1、a2、c1和c0，使题17图中的两个系统等效。第5章时域离散系统的网络结构题15图第5章时域离散系统的网络结构题17图第5章时域离散系统的网络结构18.对于题18（1）（2）（a）b0=b2=1,b1=2,a1=1.5,a2=－0.9（b）b0=b2=1,b1=2,a1=1,a2=－2第5章时域离散系统的网络结构题18图第5章时域离散系统的网络结构19*.假设滤波器的系统函数为在单位圆上采样六点，选择r＝0.95，试画出它的频率采样结构，并在计算机上，用DFT求出频率采样结构20.已知FIR（1）H(z)=1+0.8z－1+0.65z