数字信号处理 第5章.ppt

《数字信号处理 第5章.ppt》由会员分享，可在线阅读，更多相关《数字信号处理 第5章.ppt（109页珍藏版）》请在装配图网上搜索。

1、5.1 引言 5.2 用信号流图表示网络结构 5.3 无限长脉冲响应基本网络结构 5.4 有限长脉冲响应基本网络结构 5.5 线性相位结构 5.6 频率采样结构 5.7 格型网络结构 习题与上机题,第5章 时域离散系统的网络结构,5.1 引 言 一般时域离散系统或网络可以用差分方程、单位脉冲响应以及系统函数进行描述。如果系统输入、输出服从N阶差分方程： （5.1.1） 则其系统函数H(z)为 (5.1.2),为了用计算机或专用硬件完成对输入信号的处理（运算），必须把（5.1.1）式或者（5.1.2）式变换成一种算法，按照这种算法对输入信号进行运算。其实(5.1.1)式就是对输入信号的一种直接算

2、法，如果已知输入信号x(n)以及ai、bi和n时刻以前的y(ni)，则可以递推出y(n)值。但给定一个差分方程，不同的算法有多种，例如：,可以证明以上H1(z)=H2(z)=H3(z)，但它们具有不同的算法。不同的算法直接影响系统运算误差、运算速度以及系统的复杂程度和成本等，因此研究实现信号处理的算法是一个很重要的问题。我们用网络结构表示具体的算法，因此网络结构实际表示的是一种运算结构。这一章是第9章数字信号处理实现的必要基础。 在介绍数字系统的基本网络结构之前，先介绍网络结构的表示方法。,5.2 用信号流图表示网络结构 观察（5.1.1）式可知，数字信号处理中有三种基本算法，即乘法、加法和单

3、位延迟。三种基本运算框图及其流图如图5.2.1所示。,图5.2.1 三种基本运算的流图表示,z1与系数a作为支路增益写在支路箭头旁边，箭头表示信号流动方向。如果箭头旁边没有标明增益，则认为支路增益是。两个变量相加，用一个圆点表示（称为网络节点），这样整个运算结构完全可用这样一些基本运算支路组成，图5.2.2所示的就是这样的流图，该图中圆点称为节点，输入x(n)的节点称源节点或输入节点，输出y(n)称为吸收节点或输出节点。每个节点处的信号称节点变量，这样信号流图实际上是由连接节点的一些有方向性的支路构成的。和每个节点连接的有输入支路和输出支路，节点变量等于所有输入支路的输出之和。在图5.2.2中

4、，,图5.2.2 信号流图,不同的信号流图代表不同的运算方法，而对于同一个系统函数可以有多种信号流图与之相对应。从基本运算考虑，满足以下条件，称为基本信号流图。 （1） 信号流图中所有支路都是基本支路，即支路增益是常数或者是z1； （2） 流图环路中必须存在延迟支路； （3） 节点和支路的数目是有限的。,图5.2.2（a）是基本信号流图，图中有两个环路，环路增益分别为a1z1和a2z2，且环路中都有延时支路，而图5.2.2（b）不是基本信号流图，它不能决定一种具体的算法，不满足基本信号流图的条件。 根据信号流图可以求出网络的系统函数，方法是列出各个节点变量方程，形成联立方程组，并进行求解，求出

5、输出与输入之间的z域关系。 【例5.2.1】求图5.2.2（a）信号流图决定的系统函数H(z)。 解 图5.2.2（a）信号流图的节点变量方程为(5.2.1)式， 对（5.2.1）式进行Z变换，得到：, 经过联立求解得到： 当结构复杂时，上面利用节点变量方程联立求解的方法较麻烦，不如用梅逊（Masson）公式直接写H(z)表示式方便。关于梅逊公式请参考本书附录A。,一般将网络结构分成两类，一类称为有限长单位脉冲响应网络，简称FIR（Finite Impulse Response）网络，另一类称为无限长单位脉冲响应网络，简称IIR（Infinite Impulse Response）网络。FIR

6、网络中一般不存在输出对输入的反馈支路，因此差分方程用下式描述： （5.2.2） 其单位脉冲响应h(n)是有限长的，按照（5.2.2）式，h(n)表示为,另一类IIR网络结构存在输出对输入的反馈支路，也就是说，信号流图中存在反馈环路。这类网络的单位脉冲响应是无限长的。例如，一个简单的一阶IIR网络的差分方程为 其单位脉冲响应h(n)=anu(n)。这两类不同的网络结构各有不同的特点，下面分类叙述其网络结构。,5.3 无限长脉冲响应基本网络结构 IIR网络的基本网络结构有三种，即直接型、级联型和并联型。 1 直接型 将N阶差分方程重写如下： 对应的系统函数为 ,设M=N=2，按照差分方程可以直接画

7、出网络结构如图5.3.1(a)所示。图中第一部分系统函数用H1(z)表示，第二部分用H2(z)表示，那么H(z)=H1(z)Hz(z)，当然也可以写成H(z)=H2(z)H1(z), 按照该式，相当于将图5.3.1(a)中两部分流图交换位置，如图5.3.1(b)所示。该图中节点变量w1=w2，因此前后两部分的延时支路可以合并，形成如图5.3.1(c)所示的网络结构流图，我们将图5.3.1(c)所示的这类流图称为IIR直接型网络结构。,M=N=2时的系统函数为 对照图5.3.1(c)的各支路的增益系数与H(z)分母分子多项式的系数可见，可以直接按照H(z)画出直接型结构流图。,图5.3.1 II

8、R网络直接型结构,【例5.3.1】 设IIR数字滤波器的系统函数H(z)为 画出该滤波器的直接型结构。 解 由H(z)写出差分方程如下： 按照差分方程画出如图5.3.2所示的直接型网络结构。,图5.3.2 例5.3.1图,上面我们按照差分方程画出了网络结构，也可以按照H(z)表达式，直接画出直接型网络结构，这里需要用到Masson公式。 下面讲述直接型的MATLAB的表示与实现。 在MATLAB中，直接型结构由2个行向量B和A表示，B和A与数字滤波器系统函数的关系如下： A=a0, a1, a2, , aN, B=b0, b1, b2, , bM,则直接型系统函数为 调用1.4.2节介绍的MA

9、TLAB 信号处理工具箱函数filter就是按照直接型结构实现滤波器。 如果滤波器输入信号向量为xn，输出信号向量为yn，则yn=filter(B, A.xn) 按照直接型结构实现对xn的滤波，计算系统对输入信号向量xn的零状态响应输出信号向量yn，yn与xn长度相等。,2 级联型 在（5.1.2）式表示的系统函数H(z)中，分子、分母均为多项式，且多项式的系数一般为实数。现将分子、分母多项式分别进行因式分解，得到： （5.3.1） 式中, A是常数; Cr和dr分别表示H(z)的零点和极点。由于多项式的系数是实数，Cr和dr是实数或者是共轭成对的复数，将共轭成对的零点（极点）放在一起，形成一

10、个二阶多项式，其系数仍为实数；再将分子、分母均为实系数的二阶多项式放在一起，形成一个二阶网络Hj(z)。,图5.3.3 一阶和二阶直接型网络结构,【例5.3.2】 设系统函数H(z)如下式： 试画出其级联型网络结构。 解 将H(z)的分子、分母进行因式分解，得到： 为减少单位延迟的数目，将一阶的分子、分母多项式组成一个一阶网络，二阶的分子、分母多项式组成一个二阶网络，画出级联结构图如图5.3.4所示。,图5.3.4 例5.3.2图,级联型结构中每一个一阶网络决定一个零点、一个极点，每一个二阶网络决定一对零点、一对极点。在（5.3.2）式中，调整0j、1j和2j三个系数可以改变一对零点的位置，调

11、整1j和2j可以改变一对极点的位置。因此，相对直接型结构，其优点是调整方便。此外，级联结构中后面的网络输出不会再流到前面，运算误差的积累相对直接型也小。,3 并联型 如果将级联形式的H(z)展成部分分式形式，则得到： （5.3.4） 对应的网络结构为这k个子系统并联。上式中，Hi(z)通常为一阶网络或二阶网络，网络系统均为实数。二阶网络的系统函数一般为 ,式中，0i、1i、1i和2i都是实数。如果1i=2i=0，则构成一阶网络。由（5.3.4）式，其输出Y(z)表示为 上式表明将x(n)送入每个二阶（包括一阶）网络后，将所有输出加起来得到输出y(n)。 【例5.3.3】 画出例题5.3.2中的

12、H(z)的并联型结构。 解 将例5.3.2中H(z)展成部分分式形式： 将每一部分用直接型结构实现，其并联型网络结构如图5.3.5所示。,图5.3.5 例5.3.3图,在这种并联型结构中，每一个一阶网络决定一个实数极点，每一个二阶网络决定一对共轭极点，因此调整极点位置方便，但调整零点位置不如级联型方便。另外，各个基本网络是并联的，产生的运算误差互不影响，不像直接型和级联型那样有误差积累，因此，并联形式运算误差最小。由于基本网络并联，可同时对输入信号进行运算，因此并联型结构与直接型和级联型比较，其运算速度最高。,MATLAB信号处理工具箱提供了14种线性系统网络结构变换函数，实现各种结构之间的变

13、换。可惜缺少并联结构于其他结构之间的变换函数，参考文献10，18中开发了直接型与并联型的相互变换函数tf2par和par2tf。本书涉及的3种常用结构（直接型、级联型、格型）之间的变换函数有如下4种： (1) tf2sos 直接型到级联型结构变换。 (2) sos2tf 级联型到直接型网络结构的变换。 (3) tf2latc 直接型到格型结构变换。 (4) latc2tf 格型到直接型结构变换。,下面先简要介绍变换函数tf2sos和sos2tf及其调用格式，tf2latc 和 latc2tf在5.7节介绍。 (1) S, G = tf2sos(B, A): 实现直接型到级联型的变换。B和A分别

14、为直接型系统函数的分子和分母多项式系数向量，当A=1时，表示FIR系统函数。返回L级二阶级联型结构的系数矩阵S和增益常数G。,S为L6矩阵，每一行表示一个二阶子系统函数的系数向量，第k行对应的2阶系统函数为 级联结构的系统函数为 H(z)=H1(z)H2(z)HL(z),例5.3.2的求解程序如下： B=8，4，11，2； A=1，1.25，0.75，0.125； S，G=tf2sos(B, A) 运行结果： S = 1.0000 0.1900 0 1.0000 0.2500 0 1.0000 0.3100 1.3161 1.0000 1.0000 0.5000 G = 8,该结果与例5.3.

15、2所得结果等价，但本程序结果更标准。 (2) B, A =sos2tf(S, G)：实现级联型到直接型网络结构的变换。B、A、S和G的含义与S, G = tf2sos(B, A)中相同。,5.4 有限长脉冲响应基本网络结构 FIR网络结构特点是没有反馈支路，即没有环路，其单位脉冲响应是有限长的。设单位脉冲响应h(n)长度为N，其系统函数H(z)和差分方程分别为 1 直接型 按照H(z)或者卷积公式直接画出结构图如图5.4.1所示。这种结构称为直接型网络结构或者称为卷积型结构。,图5.4.1 FIR直接型网络结构,2 级联型 将H(z)进行因式分解，并将共轭成对的零点放在一起，形成一个系数为实数

16、的二阶形式，这样级联型网络结构就是由一阶或二阶因子构成的级联结构，其中每一个因式都用直接型实现。 【例5.4.1】 设FIR网络系统函数H(z)如下式： 画出H(z)的直接型结构和级联型结构。 解 将H(z)进行因式分解，得到： 其级联型结构和直接型结构如图5.4.2所示。,图5.4.2 例5.4.1图,例5.4.1的求解程序如下： B=0.96, 2, 2.8, 1.5;A=1; S, G=tf2sos(B, A) 运行结果： S =1.0000 0.8333 0 1.0000 0 0 1.0000 1.2500 1.8750 1.0000 0 0 G = 0.9600,级联结构的系统函数为

17、 H(z)=0.96(1+0.833z1)(1+1.25z1+1.875z2) 级联型结构每一个一阶因子控制一个零点，每一个二阶因子控制一对共轭零点，因此调整零点位置比直接型方便，但H(z)中的系数比直接型多，因而需要的乘法器多。在例5.4.1中直接型需要四个乘法器，而级联型则需要五个乘法器。分解的因子愈多，需要的乘法器也愈多。另外，当H(z)的阶次高时，也不易分解。因此，普遍应用的是直接型。,5.5 线性相位结构 线性相位结构是FIR系统的直接型结构的简化网络结构，特点是网络具有线性相位特性，比直接型结构节约了近一半的乘法器。第7章将证明， 如果系统具有线性相位，它的单位脉冲响应满足下面公式

18、:,（5.5.1）,观察（5.5.2）式，运算时先进行方括号中的加法（减法）运算，再进行乘法运算，这样就节约了乘法运算。按照这两个公式，第一类线性相位网络结构的流图如图5.5.1所示，第二类线性相位网络结构的流图如图5.5.2所示。和直接型结构比较，如果N取偶数，直接型需要N个乘法器，而线性相位结构减少到N/2个乘法器，节约了一半的乘法器。如果N取奇数，则乘法器减少到(N1)/2个，也近似节约了近一半的乘法器。,图5.5.1 第一类线性相位网络结构流图,图5.5.2 第二类线性相位网络结构流图,5.6 频率采样结构 我们已经知道，频率域等间隔采样，相应的时域信号会以采样点数为周期进行周期性延拓

19、。如果在频率域采样点数N大于等于原序列的长度M，则不会引起信号失真，此时原序列的Z变换H(z)与频域采样值H(k)满足下面关系式： （5.6.1）,设FIR滤波器单位脉冲响应h(n)长度为M，系统函数H(z)=ZTh(n)，则(5.6.1)式中H(k)用下式计算: 要求频率域采样点数NM。（5.6.1）式提供了一种称为频率采样的网络结构。由于这种结构是通过频域采样得来的，存在时域混叠的问题，因此不适合IIR系统，只适合FIR系统。但这种网络结构中又存在反馈网络，不同于前面介绍的FIR网络结构，下面进行分析。,将（5.6.1）式写成下式： （5.6.2） 式中Hc(z)是前面学习过的梳状滤波器，

20、Hk(z)是IIR的一阶网络。这样，H(z)是由梳状滤波器Hc(z)和N个一阶网络Hk(z)的并联结构进行级联而成的，其网络结构如图5.6.1所示。我们看到该网络结构中有反馈支路，它是由Hk(z)产生的，其极点为 ,即它们是单位圆上有等间隔分布的N个极点，第2章已学过Hc(z)是一个梳状滤波网络，其零点为 刚好和极点相同，也是等间隔地分布在单位圆上。理论上，极点和零点相互抵消，保证了网络的稳定性， 使频率域采样结构仍属FIR网络结构。,图5.6.1 FIR滤波器频率采样结构,频率域采样结构有两个突出优点： （1） 在频率采样点k处， , 只要调整H(k)（即一阶网络Hk(z)中乘法器的系数H(

21、k)），就可以有效地调整频响特性，使实践中的调整方便，可以实现任意形状的频响曲线。 （2） 只要h(n)长度N相同，对于任何频响形状，其梳状滤波器部分和N个一阶网络部分结构完全相同，只是各支路增益H(k)不同。这样，相同部分便可以标准化、模块化。各支路增益可做成可编程单元，生产可编程FIR滤波器。,然而，上述频率采样结构亦有两个缺点： （1） 系统稳定是靠位于单位圆上的N个零极点相互对消保证的。实际上，因为寄存器字长都是有限的，对网络中支路增益 量化时产生量化误差，可能使零极点不能完全对消，从而影响系统稳定性。 （2） 结构中，H(k)和 一般为复数，要求乘法器完成复数乘法运算，这对硬件实现是

22、不方便的。 为了克服上述缺点，对频率采样结构作以下修正。,首先将单位圆上的零极点向单位圆内收缩一点，收缩到半径为r的圆上，取r1且r1。 此时H(z)为 （5.6.3） 式中，Hr(k)是在r圆上对H(z)的N点等间隔采样之值。由于r1，因此可近似取Hr(k)H(k)。这样，零极点均为 。如果由于实际量化误差，零极点不能抵消时，极点位置仍处在单位圆内，保持系统稳定。,另外，由DFT的共轭对称性知道，如果h(n)是实数序列，则其离散傅里叶变换H(k)关于N/2点共轭对称，即H(k)=H*(Nk)。而且，我们将Hk(z)和HNk(z)合并为一个二阶网络，并记为Hk(z)，则,图5.6.2 频率采样

23、修正结构,N等于奇数的修正结构由一个一阶网络和(N1)/2个二阶网络结构构成。 由图5.6.2可见，当采样点数N很大时，其结构显然很复杂，需要的乘法器和延时单元很多。但对于窄带滤波器，大部分频率采样值H(k)为零，从而使二阶网络个数大大减少。所以频率采样结构适用于窄带滤波器。,5.7 格型网络结构 5.7.1 全零点格型网络结构 1. 全零点格型网络的系统函数 全零点格型网络结构的流图如图5.7.1所示。该流图只有直通通路，没有反馈回路，因此可称为FIR格型网络结构。观察该图，它可以看成是由图5.7.2的基本单元级联而成。,图5.7.1 全零点格型网络结构,图5.7.2 基本单元,2. 由FI

24、R直接型网络结构转换成全零点格型网络结构 假设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。,下面仅给出转换公式，推导过程请参考文献19： (5.7.11) (5.7.12) (5.7.13) 式中, l=N, N1, , 1。,解释 公式中的下标k（或l）表示第k(或l)个系数，这里FIR结构和格型结构均各有N个系数; (5.7.13)式是一个递推公式，上标（带圆括弧）表示递推序号，从（N）开始，然后是N

25、1, N2, , 2；注意(5.7.12)式，当递推到上标圆括弧中的数字与下标相同时，格型结构的系数kl刚好与FIR的系数相等。下面举例说明。 【例 5.7.1】 将下面三阶FIR系统函数3(z)转换成格型网络，要求画出该FIR直接型结构和相应的格型网络结构流图。 ,解 例题中N=3, 按照(5.7.11)式，有 由(5.7.12) 式，得到: 按照(5.7.13) 式，递推得到:,l=3, k=1时, l=3, k=2时, ,l=2, k=1时, 最后按照算出的格型结构的系数，画出三阶FIR直接型结构和三级格型网络结构流图如图 5.7.3所示。 ,图5.7.3 例5.7.1图,略去由全零点格

26、型网络结构转换到FIR直接型网络结构的公式，如需要了解该内容，请参考文献19。 实际上，调用MATLAB函数实现直接型网络结构与格型网络结构之间的相互转换非常容易。tf2latc实现直接型到格型结构变换，latc2tf 实现格型到直接结型结构变换。 K=tf2latc(hn): 求FIR格型结构的系数向量K=k1, k2, , kN, hn为FIR滤波器的单位脉冲响应向量，并关于hn(1)=h(0)归一化。应当注意，当FIR系统函数在单位圆上有零极点时，可能发生转换错误。,hn=latc2tf(K) 将FIR格型结构转换为FIR直接型结构。K为FIR格型结构的系数向量K，hn为FIR滤波器的单

27、位脉冲响应向量，即FIR直接型结构系数向量。显然，该函数可以用于求格型结构的系统函数的系数。 例 5.7.1的求解程序如下： hn=1, 0.9, 0.64, 0.576; K=tf2latc(hn) 运行结果： K=0.6728 0.1820 0.5760 与上面的递推结果相同。,5.7.2 全极点格型网络结构 全极点IIR系统的系统函数用下式表示: (5.7.14) (5.7.15) 式中, A(z)是FIR系统，因此全极点IIR系统H(z)是FIR系统A(z)的逆系统。下面先介绍如何将H(z)变成A(z)。 假设系统的输入和输出分别用x(n)、y(n)表示，由 (5.7.17) 式得到全

28、极点IIR滤波器的差分方程为,图5.7.4 全极点IIR系统的直接型结构,由于重新定义了输入输出，将el(n)按降序运算，rl(n)不变，即 ,图5.7.5 全极点IIR格型结构,当x(n)和y(n)分别作为输入和输出时，（5.7.27）式就是一个全极点的差分方程，由（5.7.24）（5.7.27）式描述的结构就是一阶的单极点格型网络，如图5.7.6(a)所示。如果N2，可得到下面方程组:,图5.7.6 单极点和双极点IIR格型网络结构,（5.7.35）,由上面分析知道, 全极点网络可以由全零点格型网络形成，这是一个求逆的问题。对比全零点格型结构和全极点结构，可以归纳出下面的一般求逆准则： （

29、1） 将输入到输出的无延时通路全部反向，并将该通路的常数支路增益变成原常数的倒数（此处为1）； （2） 将指向这条新通路的各节点的其它节点的支路增益乘以1； （3） 将输入输出交换位置。,调用MATLAB 转换函数可以实现全极点系统的直接型和格型结构之间的转换。 K = tf2latc(1, A): 求IIR全极点系统格型结构的系数向量K，A为(5.7.14)式给出的IIR全极点系统函数的分母多项式A（z）的系数向量。 具有零点和极点的IIR格型网络称为格梯型网络结构，这部分内容请参考文献19。 K, V = tf2latc(B, A): 求具有零点和极点的IIR格型网络系数向量K，及其梯型网

30、络系数向量V。应当注意，当IIR系统函数在单位圆上有极点时，可能发生转换错误。,B，A = latc2tf(K, allpole): 将IIR全极点系统格型结构转换为直接型结构。K为IIR全极点系统格型结构的系数向量，A为IIR全极点系统系数函数的分母多项式A（z）的系数向量。显然，该函数可以用于球格型结构的系统函数，这时分子为常数1，所以B=1。 B，A = latc2tf(K, V)：将具有零点和极点的IIR格梯型网络结构转换为直接型结构。 例如： ,则求IIR全极点系统格型结构系数向量K的程序为 A=1, 13/24, 5/8, 1/3; K=tf2latc(1, A) 运行结果： K

31、=0.2500 0.5000 0.3333 对上面所求格型结构的系数向量K，调用latc2tf求其对应的格型结构的系统函数的程序如下： K =0.2500，0.5000，0.3333； B, A=latc2tf(K, allpole) 运行结果： B = 1 0 0 0 A =1.0000 0.5417 0.6250 0.3333,对应的系统函数为 下面再推导全极点网络结构的传输函数，将(5.7.21)和(5.7.22)式进行Z变换，得到:,（5.7.36）,（5.7.37）,习题与上机题 1. 已知系统用下面差分方程描述: 试分别画出系统的直接型、级联型和并联型结构。式中x(n)和y(n)分

32、别表示系统的输入和输出信号。 2 设数字滤波器的差分方程为 试分别画出系统的直接型、级联型和并联型结构。,3. 设系统的差分方程为 式中, |a|1，|b|1, 试画出系统的直接型、级联型和并联型结构。x(n)和y(n)分别表示系统的输入和输出信号。 4. 设系统的系统函数为 试画出各种可能的级联型结构，并指出哪一种最好。,5 题5图中画出了四个系统，试用各子系统的单位脉冲响应分别表示各总系统的单位脉冲响应，并求其总系统函数。 6 题6图中画出了10种不同的流图，试分别写出它们的系统函数及差分方程。,题5图,题6图,7. 假设滤波器的单位脉冲响应为 求出滤波器的系统函数，并画出它的直接型结构。

33、 8. 已知系统的单位脉冲响应为 试写出系统的系统函数，并画出它的直接型结构。 9. 已知FIR滤波器的系统函数为 试画出该滤波器的直接型结构和线性相位结构。,10 已知FIR滤波器的单位脉冲响应为: (1) N=6 h(0)=h(5)=15 h(1)=h(4)=2 h(2)=h(3)=3 (2) N=7 h(0)=h(6)=3 h(1)=h(5)=2 h(2)=h(4)=1 h(3)=0 试画出它们的线性相位型结构图，并分别说明它们的幅度特性、相位特性各有什么特点。,15 写出题15图中系统的系统函数和单位脉冲响应。 16. 画出题15图中系统的转置结构，并验证两者具有相同的系统函数。 17

34、. 用b1和b2确定a1、a2、 c1和c0，使题17图中的两个系统等效。,题15图,题17图,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 画出系统的零极点分布图，并检验系统的稳定性。,题18图,19*. 假设滤波器的系统函数为 在单位圆上采样六点，选择r0.95，试画出它的频率采样结构，并在计算机上，用DFT求出频率采样结构中的有关系数。 20. 已知FIR滤波器的系统函数为： （1） H(z)=1+0.8z1+0.65z2 （2） H(z)=10.6z1+0.825z20.9z3 试分别画出它们的直接型结构和格形结构，并求出格型结构的有关参数。,21. 假设FIR格型网络结构的参数k1=0.08, k2=0.217, k3=1.0, k4=0.5，求系统的系统函数并画出FIR直接型结构。 22. 假设系统的系统函数为 要求: （1） 画出系统的直接型结构以及描述系统的差分方程; （2） 画出相应的格形结构，并求出它的系数; （3） 判断系统是否是最小相位。,