数字信号处理第2版张小虹8数字信号处理2课件教学.pptVIP

Fnm Wm Wn Wa Wb X Y 系数Fnm的敏感度为： 如图8.3-5所示任意数字系统，设a为输入节点， Wa为 输入节点变量的z变换； b为输出节点， Wb为输出节 由伯德敏感度定义，定义数字网络H(z)对某一支路传递 点变量的z变换；系统函数H(z)为Tab=Wb / Wa。 (8.3-12) SFnm = Tab ?Tab ?Fnm =TanTmb 式中：Tan——a点输入、 n点输出的传递函数； Tmb——m点输入、 b点输出的传递函数； Tab——a点输入、 b点输出的传递函数； Fnm——任一支路n、m的传递系数。 与极点位置敏感度类似，系统的敏感度越低，当系数有 量化误差时，系统函数的变化越小。同样的，不同结构 的系统敏感度是不同的，高阶系统的级联、并联形式 比直接形式的敏感度低。 3 2 1 4 例：如图所示二阶数字网络，传递系数已在图中标明。 ?a1 b1 x(n) y(n) b0 z?1 z?1 ?a2 求支路41传递系数?a1变化时网络的敏感度 SF41 T12 解：由敏感度的一般公式 所以 SF41 = T12 ?T12 ?F41 =T14T12 SFnm = Tab ?Tab ?Fnm =TanTmb T14 = 1+a1z?1+a2z?2 z?1 T12 = 1+a1z?1+a2z?2 b0+b1z?1 SF41 = T12 ?T12 ?F41 =T14T12 (1+a1z?1+a2z?2)2 (b0+b1z?1)z?1 = 因为T14 = ?a1，且T12=H(z)，而 与前面结果相同。 所以 直接得到系统的敏感度， 此例也可以利用 SF41 = T12 ?T12 ?F41 H(z) = T12 = 1+a1z?1+a2z?2 b0+b1z?1 (1+a1z?1+a2z?2)2 (b0+b1z?1)z?1 = ? H(z) ?(?a1 ) = 1+a1z?1+a2z?2 b0+b1z?1 ? ?(?a1 ) 8.3.4系数量化对IIR系统频响的影响 系统的系数量化对零、极点位置的改变，将导致系 统频响特性的改变。 在以下讨论中为分析简便，只考虑系数量化误差，不 涉及运算误差。 (8.3-13) △Tab的推导留作习题。 △Tab TanT mb △ Fnm 1? Tnm △ Fnm = 当△Fnm较大时，引起系统函数Tab的变化△Tab为 无限精度的理想系统函数一般形式为 有效字长b位、定点舍入量化，△bk、△ak的误差范围 设量化后的系数为 ?q/2 △bk，△ak ≤q/2 (8.3-15) 为 ak=ak+ ?ak (8.3-14) H(z)= ? k=0 M bkz?k ? k=1 N akz?k 1? = B(z) A(z) △bk、△ak是系数量化误差，如果系数bk、 ak采用 ^ ^ bk=bk+ ?bk ^ 量化后的实际系统函数为 式中 (8.3-16) ? k=0 M bkz?k ? k=1 N akz?k 1? H(z)= ^ ^ ^ ? k=0 M bkz?k ? ? k=1 N akz?k ? 1? ? k=0 M ?bkz?k ? k=1 N ?akz?k = ? k=0 M ?bkz?k ， ?B(z)= ? k=0 N ?akz?k ， ?A(z)= B(z) +?B(z) = A(z) ??A(z) 将H(z)代入上式并整理 由上式得到实际系统与理想系统的偏差 HE (z) B(z) +?B(z) A(z) ??A(z) H(z)= ^ B(z) A(z) ?B(z) + ?A(z)B(z)/A(z) A(z) ??A(z) = + =H(z)+ ?B(z) +?A(z)H(z) A(z) ??A(z) (8.3-17) HE (z)= H(z)? H(z) ^ (8.3-18) ?B(z) +?A(z)H(z) A(z) ??A(z) = 实际系统，可以由如图8.3-7所示的理想系统与偏差系统 系统的频响误差为 HE(z)的并联组合而成。 x(n) y(n) H(z) ? (n)= y(n) +e(n) HE(ej?) = ? (ej?)?H (ej?) e(n) HE(z) 但实际上计算时，很难知道△bk、△ak的准确值。为了 频响的均方误差为 将HE(z)代入，可以计算频响的均方误差。 理论上可以频响的均方误差即误差能量度量频响偏差。 为了估计? 2的大小，假设△bk、△ak是独立均匀等概分 布的随机变量。 (8.3-19) 2? 1 ? ? ?? ?HE(ej?)?2d? ? 2 = j2? 1 c HE (z)HE (1/z) z?1 dz = 若△bk系数共有M1项，△ak系数共有N1项，则 考虑△bk、△ak是统计独立的