数字信号处理

1、第第8章章 数字滤波器的结构数字滤波器的结构 在得到在得到H(z)之后，还应该将其实现。之后，还应该将其实现。H(z)实际上实际上是一个数学表达式，它的不同的数学表示形式是一个数学表达式，它的不同的数学表示形式对应着不同的算法结构。这就是说，同一个系对应着不同的算法结构。这就是说，同一个系统函数统函数H(z)，可以用不同的算法结构来实现。，可以用不同的算法结构来实现。不同的算法结构，使得滤波器具有不同的性能。不同的算法结构，使得滤波器具有不同的性能。与算法结构有关的滤波器性能主要指以下几个与算法结构有关的滤波器性能主要指以下几个方面。方面。 1相乘、相加、延迟等运算的运算量，这关系到运算的复杂

2、程度和运算速度。 2所用的加法器、乘法器、延迟器的数目，这关系到系统的复杂度和系统成本，而延迟器的数目直接影响系统所需存储器的多少。 3系统频率特性对于乘法器系数变化的灵敏度，这是数字滤波器的非常重要的性能。在设计H(z)时，并没有考虑系数的精度问题，即是在无限精度的条件下进行设计的。但是，在实现H(z)时，这些系数都要化为二进制数以有限的字长即有限的精度来进行运算。于是，实现时H(z)产生的系数误差使得所设计的零极点位置发生变化，系统的频率特性也就发生变化。 n在相同的二进制字长的条件下，或者说在相同的系数误差范围的情况下，不同的滤波器结构所引起的频率特性的变化程度不同，频率特性变化大的结构

3、对于乘法器系数变化的灵敏度就高，或者说这种结构的灵敏度特性差，反之则是灵敏度特性好。 n对于IIR滤波器，结构的灵敏度特性不仅影响系统的频率特性，还影响系统的稳定性，这是因为滤波器实现时系数的误差若使得极点的位置变化较大，就有可能从单位园内移出，这将使系统不稳定。 n综上所述，算法结构的选择对于一个系统的实现非常重要。一个确定的数字系统有其确定的差分方程、确定的单位抽样响应h(n) 以及确定的系统函数H(z)，而同一个系统函数H(z)，却可以画出不同的算法结构。用什么样的结构来实现H(z)，也是数字滤波器的一个很重要的问题，而信号流图是数字网络结构的一种非常方便、直观、有效的表示。 8.1 8

4、.1 数字网络的信号流图数字网络的信号流图 8.1.1 8.1.1 信号流图及其有关概念信号流图及其有关概念信号流图是由连接节点的有向线段构成的网络，它是表示信信号流图是由连接节点的有向线段构成的网络，它是表示信号流通的几何图形。信号流图清楚地表示了系统的算法结号流通的几何图形。信号流图清楚地表示了系统的算法结构，通过信号流图，可以对系统进行有效的分析，还可以构，通过信号流图，可以对系统进行有效的分析，还可以方便地求出系统函数。方便地求出系统函数。下面结合图下面结合图8.18.1来介绍信号流图的有关概念。来介绍信号流图的有关概念。 图图8.1 8.1 一个信号流图一个信号流图 1 1节点：节点

5、： 信号流图中每一节点都对应一个变量，或者说信号流图中每一节点都对应一个变量，或者说代表一个信号。节点又叫做节点变量。代表一个信号。节点又叫做节点变量。 2 2支路与支路传输：支路与支路传输： 支路是连接两个节点的有向线段。支路是连接两个节点的有向线段。支路旁边标注的系数或称加权叫做支路传输，它们起支路旁边标注的系数或称加权叫做支路传输，它们起着相乘的作用。一种支路传输就是乘法器系数；另一种实着相乘的作用。一种支路传输就是乘法器系数；另一种实际上表示延迟，如图中支路际上表示延迟，如图中支路X0X1X0X1的传输的传输z-1z-1，有，有X1 = z-1 X1 = z-1 X0 , X0 , 即

6、即x1(n) = x0(n-1)x1(n) = x0(n-1)。 3. 3. 源源( (节节) )点点: : 对于一个节点，流入该节点的信号叫输入，对于一个节点，流入该节点的信号叫输入，流出该节点的信号叫输出。若一个节点只有输出支路与之流出该节点的信号叫输出。若一个节点只有输出支路与之相连接，则称之为源节点，或输入节点。相连接，则称之为源节点，或输入节点。 4. 4. 汇点汇点: : 若一个节点只有输入支路与之相连接，则称之为若一个节点只有输入支路与之相连接，则称之为汇点，或输出节点。汇点，或输出节点。 5. 5. 混合节点混合节点: : 若一个节点既有输入支路与之相连接，又有若一个节点既有输

7、入支路与之相连接，又有输出支路与之相连接，则称之为混合节点。输出支路与之相连接，则称之为混合节点。 6. 6. 开路径开路径: : 也叫通路，即是从某一节点出发，沿支路方向也叫通路，即是从某一节点出发，沿支路方向连续经过一些支路而中止到另一节点上的路径。连续经过一些支路而中止到另一节点上的路径。 留意，留意，每一节点只通过一次。每一节点只通过一次。 7 7闭路径闭路径: : 也叫自环或者环路，即是从某一节点出发，沿也叫自环或者环路，即是从某一节点出发，沿着支路方向，连续经过一些支路又中止在出发节点的路径。着支路方向，连续经过一些支路又中止在出发节点的路径。留意，途中各节点只通过一次。留意，途中

8、各节点只通过一次。 8 8节点变量的值节点变量的值: : 设连接节点设连接节点Xi Xi 与与XjXj的支路的支路XjXiXjXi的支路的支路传输为传输为tjitji，则节点变量，则节点变量XiXi的值为：的值为： (8.1) (8.1) 即节点变量的值等于流入该节点的全部信号的叠加。要特别即节点变量的值等于流入该节点的全部信号的叠加。要特别注意，计算节点变量值时不要考虑输出支路。注意，计算节点变量值时不要考虑输出支路。jjijitXXn图8.1中，节点变量X1、X2、X3、X4之值分别为：n X1 =z-1 (8.2) n X2=aX0+bX1+eX3 (8.3)n X3=cX1+dX2+f

9、X3 (8.4)n X4=X2 (8.5) n而源节点X0 无值可言。 8.1.2 8.1.2 解代数方程组求节点变量之值解代数方程组求节点变量之值将流图中各节点变量的值都表示出来，就组成了一代数方程将流图中各节点变量的值都表示出来，就组成了一代数方程组。如果流图中除了源点之外有组。如果流图中除了源点之外有M M个节点，那末可以得到个节点，那末可以得到含有含有M M个未知数由个未知数由M M个方程组成的线性方程组。但是为了简个方程组成的线性方程组。但是为了简化方程组的求解过程，应当尽量减少需要求解的节点个数。化方程组的求解过程，应当尽量减少需要求解的节点个数。 n一般情况下，这样的线性方程组有

10、且只有一组解，另外，往往也不需要解出每个节点变量，因此，采用克莱姆法则来求解比较方便。n 例8.1 用代数方程组求解法求图8.1的系统函数 H=X4/X0 。n 解：由于X4=X2，因此只需要解出节点变量X2 。 由上面已有的节点之值表示式(8.2)、(8.3)和(8.4)，可以得到方程组的矩阵表示形式： 0 1 10 0 1001321aXXzXXX f-c d -e-b 根据克莱姆法则，可求得： X2 =2/这里为系数矩阵行列式： def11-f de- 1 1 -f d ce- 1 b-0 0 1而 ec1)-b(fXz-1)-(faX 1-f ce- b-1-f 0 e- 1-f 0

11、c e- aX b-0 Xz 1 01 -001000-12XzaX于是可求得系统函数： defeczzfbfaXXXXXH1) 1() 1(11020204n 看得出，从X2 到X4这条传输为1的支路是从X2延伸出来的，实际上，可以直接从X2 获取输出，并且直接求出节点X2与源点X0之间的系统函数。以此类推，可以根据信号流图求出从源点到其它任何一个节点的传输函数也即系统函数。 8.1.3 8.1.3 化简信号流图求系统函数化简信号流图求系统函数如果一个信号流图最终被化简为图如果一个信号流图最终被化简为图8.28.2所示的最简形式，那所示的最简形式，那末支路传输末支路传输H H显然就是系统函数

12、，因为显然就是系统函数，因为Y=HX0Y=HX0，H=Y/X0 H=Y/X0 。信号流图化简的依据是节点变量的表示式和代数方程的恒信号流图化简的依据是节点变量的表示式和代数方程的恒等关系，主要有下面三种情况。等关系，主要有下面三种情况。 图8.2 信号流图的最简形式 YX0 H 8.1.3.1 8.1.3.1 合并支路合并支路串联的支路通过支路传输相乘来合并，并联的支路通过支路串联的支路通过支路传输相乘来合并，并联的支路通过支路传输相加来合并。例如图传输相加来合并。例如图8.38.3中：中： H=a(bdeg+cfh)ij H=a(bdeg+cfh)ij。 图图 8.3 8.3 支路的合并支路

13、的合并 8.1.3.2 8.1.3.2 消除节点消除节点 图图8.4 8.4 节点的消除节点的消除 图图8.48.4中，由于中，由于 X2 =bX1 =b(aX0)=abX0 X2 =bX1 =b(aX0)=abX0 以及以及 X3 =cX1 X3 =cX1 =c(aX0)=acX0 =c(aX0)=acX0 ，因此消去了节点，因此消去了节点X1X1。 8.1.3.3 8.1.3.3 消除自环消除自环 图图8.5 8.5 自环的消除自环的消除n由图8.5左边的流图有：X2 =bX1 =b(aX0+cX2)，于是X2 =abX0+bcX2 ，即为中间的流图；再将含X2的项移到等式左边，得到： ，

14、即为右边的流图，此时消除了自环，流图化为最简。021XbcabXn 例例8.2 如图如图8.6所示所示, 利用自环消除的规则首先消除了左边利用自环消除的规则首先消除了左边流图的自环和节点流图的自环和节点X2，然后经支路合并后得到了最简形，然后经支路合并后得到了最简形式。式。n 此例的流图还可以首先消除节点此例的流图还可以首先消除节点X3来进行简化，如图来进行简化，如图8.7所示。所示。 图图 8.6 8.6 例例8.28.2流图的化简一）流图的化简一） 图图 8.7 8.7 例例8.28.2流图的化简二）流图的化简二）n 例例8.3 8.3 信号流图如图信号流图如图8.8(a) 8.8(a)

15、所示，中间的三个节点所示，中间的三个节点X2X2、X3X3、X4X4以及三个自环按图中所示步骤逐一消除。应注意，以及三个自环按图中所示步骤逐一消除。应注意，在流图在流图(c)(c)中，节点中，节点X3X3后面有一条输入支路的权值要受后面有一条输入支路的权值要受X3X3上自环消除的影响，事实上，在流图上自环消除的影响，事实上，在流图(b)(b)中，有：中，有： X3 = abX1 + bcX3 + eX4X3 = abX1 + bcX3 + eX4n由此式得：由此式得： n这即为流图这即为流图(c)(c)中的情形。后面的化简也类似处理。中的情形。后面的化简也类似处理。41311XbceXbcab

16、X 图图 8.8 8.8 例例8.38.3流图的化简流图的化简n最后得： n显然，D=X5/X1 就是这个系统的系统函数。115)1 (1XbcgfedbcabdfDXX8.1.4 Mason8.1.4 Mason公式公式 在给出在给出MasonMason公式之前，首先解释以下名词。公式之前，首先解释以下名词。1 1通路传输：通路边界间即开始节点与终止节点之间通路传输：通路边界间即开始节点与终止节点之间各支路传输之乘积。各支路传输之乘积。2 2环路传输：绕环路一周各支路传输之乘积。环路传输：绕环路一周各支路传输之乘积。 3不接触：两条通路间或两个环路间或一条通路与一个环路之间若无公共节点则称它

17、们互不接触。Mason公式给出了一个信号流图中，从源点到其余任一节点的传输函数即系统函数的表达式： (8.6) iigHn其中， 为流图的行列式：n = 1 -（所有环路传输之和）+（每两个互不接触的环路传输乘积之和）-（每三个互不接触的环路传输乘积之和）+ n而 gi 是从源点到这一节点的第 i 条通路的通路传输，i 则是此通路流图的余子式： n i = 1-（与此通路不接触的各环路传输之和）+（与此通路不接触的每两个互不接触的环路传输乘积之和）-（与此通路不接触的每三个互不接触的环路传输乘积之和）+ n 例8.4 由一电路网络所得到的信号流图如图8.9所示，试用Mason公式求其系统函数：

18、H=V4/Vg 。 图图 8.9 8.9 一个电路网络的信号流图一个电路网络的信号流图n 解： 这实际上是一个模拟网络的信号流图，但是，从流图的观点，并不需要考虑节点和支路传输的物理意义，也就是说，模拟网络的信号流图可以与数字网络的信号流图同样地处理。 n 图中，每条支路用其编号来代表，如、等等；另外，对于同一条支路，有： ZjYj=1, RjGj=1。 n为求得Mason公式中的分母，应找出流图中所有环路：na.2,13； b.10,11； c.3,12； d.8,9； e.6,7,8,13; f.8,3,4,11。n于 是 = 1 - ( - Z 1 Y 2 - G 3 R 4 + - Z

19、 1 G 3 + Z 1 G 3 -G3R4)+(G3R4Z1Y2-G3R4)n为求得Mason公式中的分子，需要找出由Vg到V4的所有通路。 各条通路的具体情况如下： 1). ,； g1=G3R4； 1=1-(-Z1Y2+Y1Z1)=1+Z1Y2-.2). ,； g2=Y2Z1Y1R4=Y2R4； 2=1.3). ,； g3=-Y2Z2G3R4=-G3R4； 3=1-Y1Z1=1-.4). ,； g4=G3Z1Y1R4=G3R4； 4=1.5). ,； g5=Y2Z1(-G3)R4=-Z1Y2G3R4； 5=1.6). ,； g6=Y2Z2(-G3)Z1Y1R4=-G3R4； 6=1.n故可

20、得到: n n 采用Mason公式可以快捷地从一个信号流图得到系统函数，也可以快捷地由已知的系统函数画出信号流图。3143214342434)1 ()1 (1)1 (GZRGYZRGRYRGgVVHiig8.1.5 8.1.5 信号流图的转置信号流图的转置将信号流图转置是对其形式的一种变换，这种变换包括三项将信号流图转置是对其形式的一种变换，这种变换包括三项操作：操作：1.1.将输入变量源点和输出变量汇点交换位将输入变量源点和输出变量汇点交换位置；置；2.2.将各条支路反向；将各条支路反向；3.3.保持每条支路的支路传输不变。保持每条支路的支路传输不变。 n当信号流图中只有一个源点和一个汇点时

21、，转置后的流图与原流图有相同的系统函数。n 例8.5 将图8.10(a)的信号流图进行转置，并根据Mason公式说明转置前后的流图具有相同的系统函数。 图图 8.10 8.10 信号流图的转置信号流图的转置n解：在图8.10中，对流图(a)进行转置的三项操作就得到了流图(b)，再按习惯将输入画在左边，输出画在右边，就得到流图(c)，图(b)和图(c)完全等价。n根据Mason公式，很容易得到这三个流图的系统函数。在这三个流图中，都只有一个环路，而且环路传输都是acz-1 ；又，三个流图中从源点到汇点都只有一条通路，通路传输为c, 并且这条通路与环路接触。于是，三个信号流图具有相同的系统函数，即

22、为：11)()()(aczczXzYzH8.2 IIR数字滤波器的结构数字滤波器的结构 8.2.1 直接型直接型IIR数字滤波器是递归型的线性时不变因果系统，其差分方数字滤波器是递归型的线性时不变因果系统，其差分方程可以表示为：程可以表示为： （8.7）MiNiiiinybinxany01)()()(n对上式两边进行z变换，可得：n (8.8)n图8.11所示的信号流图是对差分方程(8.7)式或者其z变换式(8.8)式的直接实现，因此这种形式的IIR滤波器结构叫做直接型。MiNiiiiizYzbzXzazY01)()()( 图图 8.11 8.11 直接直接 I I 型型n由(8.8)式可以得

23、到IIR滤波器的系统函数： n n (8.9) n显然，分子的系数ai 确定了H(z)的零点，而分母的系数bi 确定了H(z)的极点。NiiiMiiizbzazXzYzH101)()()(n图8.11所示的直接型结构可以看成是两个独立网络的级联，也即第一个网络的输出y1(n)正是第二个网络的输入。第一个网络是：n Miiinxany01)()(n故可得到第一个网络的系统函数：n (8.10)MiiizazXzYzH011)()()(nH1(z)与H(z)的分子多项式对应，也就是说，第一个网络实现了滤波器的零点。n第二个网络是：n n故其系统函数是：n (8.11)Niiinybnyny11)(

24、)()(NiiizbzYzYzH11211)()()(nH2(z)的分母正是H(z)的分母，因此第二个网络实现了滤波器的极点。n整个滤波器是这两个网络的级联：n (8.12)n 如今，将图8.11中的两个网络的级联次序交换，如图8.12所示。 )()()()()()()()()(2111zHzHzYzYzXzYzXzYzH 图图8.12 8.12 直接直接型型。n这时，前面的网络即反馈网络为：n n对这个式子进行z变换后可以得到其系统函数：n (8.13) Niiinybnxny122)()()()(11)()(212zHzbzXzYNiiin 图8.12后面的网络为： n经z变换后得到系统函

25、数： (8.14)Miiinyany02)()()()()(102zHzazYzYMiiin图8.12的系统函数为： n （8.15） n(8.15)式与(8.9)式最后结果完全相同。NiiiMiiizbzazHzHzYzYzXzYzXzYzH1012221)()()()()()()()()(n 图8.12也是IIR滤波器的直接型结构，叫直接型。上面的讨论同时也说明了，级联网络总的输入输出关系即总的系统函数与各子网络级联的次序无关；级联网络总的系统函数等于各级联子网络系统函数的乘积。n 8.2.2 正准型n将图8.12中两列传输比为z-1的支路合并为一列假设M=N），就得到正准型结构 图图8.

26、13 IIR8.13 IIR滤波器的正准型结构滤波器的正准型结构n观察图8.11、图8.12以及图8.13(a)、(b)，容易看出，这四个流图中环路的情况以及由输入到输出的通路情况都完全相同(设M=N），根据Mason公式，很快就可以知道，这四个流图的系统函数都如(8.9)式所示；并且，四个流图中的前向通路与系统函数H(z)的分子很好地对应，而流图中的反馈回路与H(z)的分母很好地对应。因此，在这样的情况下，我们很容易根据信号流图写出其系统函数，也很容易根据系统函数画出相应结构的信号流图。n 例8.6 已知系统函数： ，用正准型结构实现。n 解：将H(z)化成IIR滤波器系统函数的标准形式：n

27、 n流图如图8.14所示。这里应该特别注意，H(z)的分母的第一项应该化为1，否则容易出错。4 . 06 . 028 . 02 . 43)(22zzzzzH21212 . 03 . 014 . 01 . 25 . 1)(zzzzzH 图图8.14 8.14 例例8.68.6的信号流图的信号流图n IIR滤波器的直接型结构和正准型结构都是对系统差分方程的直接实现，它们的共同缺点是灵敏度特性差，系数误差会引起零极点位置的较大的变化，而正准型比直接型优越的地方是节省一半的存储单元，因此，直接型结构一般不采用，而正准型结构当阶次高N大于3或者4时也不采用。 8.2.3 8.2.3 级联型级联型 将将(

28、8.9)(8.9)式的分子、分母多项式进行因式分解：式的分子、分母多项式进行因式分解： （A=a0A=a0）NiiMiiNiiiMiiizdzcAzbzazH111110)1 ()1 (1)(n因系数ai、bi都为实数，故零点ci与极点di或是实根，或是共轭复根，而每一对共轭因式又可以合并为一个实系数的二次三项式，于是有：n (8.16)122111221112211111)1 ()1 ()1 ()1 ()(NiNiiiiiiMiMiizzzqzzzpAzHn单实根因式可以看作二阶因式的特例，设NM，则上式又可写成：n （8.17）LiiLiiiiizHAzzzzAzH1122112211)(

29、11)(n这里L是由N/2到N范围内的一个整数。 n 子网络Hi(z)的结构如图8.15所示，而整个系统就是由L个这样的二阶子网络级联而成。当然，可能有的子网络中有些系数等于0，此时图8.15所示的结构会得到简化。 图图 8.15 8.15 一个子网络的结构一个子网络的结构n 由(8.17) 式可以清楚地看到，1i和2i确定第i对零点，1i和2i确定第i对极点，也就是说，级联型结构是将零点和极点分散到各个子网络中来实现的，这就使得这种结构的灵敏度特性优于直接型和正准型结构。并且，调整任何一个子网络的零点和极点都不影响其它的零点和极点，即零极点具有独立性。n 应该特别指出，零点对与极点对之间共有

30、L! 种不同的搭配方式，而对于每一种搭配所得到的L个子网络，又有L! 种不同的排列次序。这些不同的方案总的系统函数都相同，但系统特性或者说零极点位置对于系数变化的灵敏度不同。因此，可以选择最优的零极点搭配和级联次序，以降低系统频率特性对于系数变化的灵敏度，提高系统的稳定性。 8.2.4 8.2.4 并联型并联型还可以将还可以将(8.9)(8.9)式的式的H(z)H(z)表示成如下的部分分式展开式：表示成如下的部分分式展开式： （8.198.19）NMiiiNiiiiiNiiizCzdzdzeBzpAzH011*111121)1)(1 ()1 (1)(n将式中每一对共轭的一次因式合并为实系数的二

31、次三项式，并且假设M=N，就得到：n (8.20) 012211110112111)(CzzzzpAzHNiiiiiNiii 图图8.16 IIR8.16 IIR滤波器的并联型结构滤波器的并联型结构n由于并联支路的极点也是整个网络的极点，而并联支路的零点却不一定是整个网络的零点，因此并联型的网络结构，可以独立地调整极点位置，但不能独立地控制零点。当然，与级联型一样，并联型结构的灵敏度特性优于直接型和正准型。另外，并联型结构的运算误差比级联型结构小。 8.3 FIR 8.3 FIR 数字滤波器的结构数字滤波器的结构 FIR FIR数字滤波器的差分方程为：数字滤波器的差分方程为： (8.21) (

32、8.21)并且系统函数为：并且系统函数为： (8.22) (8.22)10)()()(Niinxihny10)()(NnnznhzH 8.3.1 8.3.1 横截型横截型 由于横截型结构实际上是对差分方程由于横截型结构实际上是对差分方程(8.21)(8.21)式的直接实现，式的直接实现，所以横截型又叫做直接型。另外，由于差分方程：所以横截型又叫做直接型。另外，由于差分方程： 所以横截型还可以叫做卷积型。所以横截型还可以叫做卷积型。)()()()()( 10nxnhinxihnyNi图图 8.17 8.17 横截型结构一）横截型结构一）图图 8.18 8.18 横截型结构二）横截型结构二）n运用

33、Mason公式很容易将系统函数H(z)与图8.17或者图8.18对应起来。显然，这两个流图中都没有环路，因此H(z)的分母为1；而分子是从输入到输出的N条通路传输之和。n 如果FIR滤波器是线性相位的，那末其冲激响应h(n)满足对称条件，这时，图8.17和图8.18都可以简化。1. 偶对称的情形, 有: h(n)=h(N-1)-n 0nN-1 当N为偶数时，有： (8.23)信号流图如图8.19(a) 所示，此结构与FIR滤波器一般的直接型结构相比，乘法器减少了一半。12/0)1(12/0)1(12/010)()1()()()(NnnNnNnnNNnnNnnzznhznNhznhznhzH 图

34、图8.19 8.19 线性相位线性相位FIRFIR滤波器的直接型结构滤波器的直接型结构 n当N为奇数时，有：n n (8.24)n信号流图如图8.19(b) 所示，乘法器数目也减少了近一半。211210)1(211210)1(12101021)(21)1()()()(NNnnNnNNnnNNnnNnnzNhzznhzNhznNhznhznhzH2. 奇对称的情形, 有: h(n)=-h(N-1)-n 0nN-1当N为偶数时，可得： （8.25）只需在图8.19(a)中，在由z-(N-1-n) 来的信号旁写一个减号即可得到这种情况下的网络结构。 12/0)1()()(NnnNnzznhzHn当N

35、为奇数时，可得： (8.26)n此时的网络结构，除了在图8.19(b)中，在由z-(N-1-n) 来的信号旁写一减号以外，还应去掉乘法器系数为h(N-1)/2)的支路。1210)1()()(NnnNnzznhzH8.3.2 8.3.2 级联型级联型一般情况下，冲激响应一般情况下，冲激响应h(n)h(n)为实数，所以将为实数，所以将(8.22)(8.22)式中的多式中的多项式分解因式后，项式分解因式后，H(z)H(z)可以写成若干实系数的二次三项式可以写成若干实系数的二次三项式的乘积，即：的乘积，即： (8.27) (8.27)这里这里K K是是(N-1)/2(N-1)/2到到N-1N-1范围内

36、的某一整数，于是可以得到范围内的某一整数，于是可以得到FIRFIR滤波器的级联型结构。滤波器的级联型结构。 KiiiiNnnzzznhzH12211010)()()( 图图8.20 FIR8.20 FIR滤波器的级联型结构滤波器的级联型结构n这种形式的结构中，每一个子网络控制一对零点，即零点可以独立调整，而且零点位置变化的灵敏度优于横截型结构。但是这种结构所需的乘法运算次数比横截型多。 8.3.3 8.3.3 频率抽样型频率抽样型 非递归型的非递归型的FIRFIR滤波器也可以用递归算法来实现，这就滤波器也可以用递归算法来实现，这就是频率抽样型结构。事实上，如果按照用频率抽样法设计是频率抽样型结

37、构。事实上，如果按照用频率抽样法设计的系统函数的系统函数(7.56)(7.56)式来实现式来实现FIRFIR滤波器，那末得到的就是滤波器，那末得到的就是频率抽样型结构。频率抽样型结构。 n显然，这样的结构中出现了反馈回路，当然同时也出现了极点 。这又如何解释FIR滤波器的极点都集中在z=0呢？n (8.28)n显然，这个结构主要由两个子网络级联而成，它们的流图分别如图8.21和图8.22所示。kNjkNeWz2)()(11)(1)(101zHzHNzWkHNzzHbaNkkNN图图 8.21 Ha(z) 8.21 Ha(z)的信号流图的信号流图 图图 8.22 Hb(z) 8.22 Hb(z)

38、的信号流图的信号流图n 可以看到，Ha(z)在单位圆上均匀分布的N个零点正是Hb(z)的N个极点，因此，当这两个网络前后级联时，它们正好相互抵消；另外，级联结构还使得Ha(z)在z=0处的极点抵消了Hb(z)在z=0的一阶零点。因此，最终的结果是保留了FIR滤波器原有的零极点，即在z = 0的N-1阶极点和有限z平面上的N-1个零点。n但是实际上，Ha(z)在单位园上的N个零点是靠延时来实现的，因此能够准确实现；而Hb(z)在单位圆上的极点却是靠复数乘法来实现的，故不能准确实现。既然单位圆上的极点不能被零点完全抵消，滤波器就会出现不稳定现象，因此，应当对上面所述的网络结构进行修正。n令0r1并

39、且r1, 用rz-1 来代替(8.28)式的H(z)中的z-1，即有：n （8.31）n 于是，的零点移到了半径为r的圆上： 1011011)()(11)()1 ()(NkkNeNkkNNNzrWkHzHNzrWkHNzrzH)1 ()(NNezrzH1, 1 , 0,2NkrezkNjkn而并联网络的极点也移到了半径为r的圆上：n n这样，即使极点不能完全被零点抵消，由于是在单位圆内，也不会导致整个系统不稳定。n 下面，还要将并联网络部分所含的运算都化为实数运算，为此要利用以下周期性：n , 1, 1 , 0,2NkrerWzkNjkNkkNkNNWW)()()(kNHkHn还可以得到： n

40、设0kN-1，那么也有0N-kN-1，此时有：H*(k)=H(N-k)； 又令H(k)=|H(k)|ej(k)。n 现在将并联网络中的第k及第(N-k)子网络合并为一个二阶网络，即有： )()()()()(*10*10kNHkHWnhWnhkHNnnkNNnnkN (8.33)这里 (8.34) )()(2)/2cos(21/2)(cos2)(cos2)(1)(*1)(1)(1)( 2211111)(1zHkHzrzNkrzNkkrkkHzrWkHzrWkHzrWkNHzrWkHkkNkNkNNkN2211)/2cos(21/2)(cos)(cos)(zrzNkrzNkkrkzHkn 若N为偶

41、数，那么(8.33)和(8.34)式中k=1,2,N/2-1，故N-k=N-1,N-2,N/2+1。n当k = 0： (8.35)n当k = N/2： (8.36)11001)0(1)0()(rzHzrWHzHN112/2/1)2/(1)2/()(rzNHzrWNHzHNNNn 若N为奇数，那么(8.33)和(8.34)式中k=1,2,(N-1)/2，故N-k=N-1,N-2,(N+1)/2，而H0(z)如(8.35)式所示。n 因此，由(8.31)式，有：n当N为偶数， n (8.37)n当N为奇数， n (8.38) )()()()(2)(1)(12120NkNkezHzHzHkHzHNz

42、H )()()(2)(1)(2110NkkezHzHkHzHNzHn 当N为偶数时，FIR滤波器的整个频率抽样型的网络结构就如图8.27所示。在这样的结构中，He(z)的N个零点也即并联网络的N个极点都移到了半径为rr h=1,-1,0,0,1; h=1,-1,0,0,1;n C,B,A=fir_freq_structure(h); C,B,A=fir_freq_structure(h); n程序运行结果为：n C= 4.2979 3.0867 1.0000n B= 0.4653 -0.9856 0.6479 0.0765n A= 1.0000 -0.6180 1.0000 1.0000 1.

43、6180 1.0000 1.0000 -1.0000 0n于是有： 11618. 110765. 06479. 00867. 3618. 019856. 04653. 02979. 451)(12112115zzzzzzzzzH 第第9章章 数字信号处理中的有限字长效应数字信号处理中的有限字长效应9.1 概述概述 在这部分之前，我们讨论的实际上是抽样数据信号的处理，在这部分之前，我们讨论的实际上是抽样数据信号的处理，因为前面的内容，都只是涉及信号在时间上是离散的这一因为前面的内容，都只是涉及信号在时间上是离散的这一特征，并没有涉及数值上离散的特征。而对于真正的数字特征，并没有涉及数值上离散的特

44、征。而对于真正的数字信号的处理，只需要在前面所讨论的离散时间信号处理的信号的处理，只需要在前面所讨论的离散时间信号处理的原理和方法的基础上，加入字长效应的影响。原理和方法的基础上，加入字长效应的影响。 9.1.1 数字系统与有限字长效应数字系统与有限字长效应 一个线性、时不变、因果系统的差分方程为：一个线性、时不变、因果系统的差分方程为： NiiMiiinybinxany10)()()(n到目前为止，在讨论问题时，对于系统的输入序列x(n)、输出序列y(n)，对于方程中的系数ai、bi等等，都认为它们的数值有无限精度。但是当具体实现一个离散系统时，无论用软件方式还是硬件方式，实际处理的数都是有

45、限字长的二进制数。当认为一个离散系统的数据是无限精度的，这样的系统就叫做抽样数据系统；若离散系统的数据是有限字长的，则此系统就是数字系统。n 对于一个数字系统，由于本应为无限精度的数据变为有限字长来进行处理，因此肯定会对系统的特性产生一定影响，这就是有限字长效应问题。这个问题本来是数字信号处理中的一个重要问题，但是，随着计算机和微处理器技术的飞速发展，运算速度和运算精度都在不断提高，使得有限字长效应的重要性已逐渐降低。不过，在数字信号处理的一些实际应用中，这个问题还是存在的，因此，有必要了解它的影响以及降低影响的一些方法。 9.1.2 关于数的表示关于数的表示 进行数字信号处理时，数的表示有定

46、点制和浮点制两种。进行数字信号处理时，数的表示有定点制和浮点制两种。浮点制运算比定点制运算的动态范围大，处理精度高，但浮点制运算比定点制运算的动态范围大，处理精度高，但实现较复杂而且运算速度较慢，因而常用于计算机上的软实现较复杂而且运算速度较慢，因而常用于计算机上的软件实现，进行非实时处理。件实现，进行非实时处理。 n在实时处理中定点制运算得到广泛应用，因为它运算速度较快而且硬件实现较经济，但是由于定点运算的动态范围和处理精度受限制较大，因而有限字长效应问题比较突出。本章主要讨论定点制算法的有限字长效应。 9.1.3 量化误差量化误差 数的定点表示：设寄存器长数的定点表示：设寄存器长L+1位，则除了一位符号位外，位，则除了一位符号位外，可表示的最小数为可表示的最小数为q=2-L，这个值称为量化间距。若要处，这个值称为量化间