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

1)、第一步先求出FIR DF的单位脉冲响应h(n) 换为y(n)。可见是相同的分析方法，可以得到相同的设 上式与最小平方滤波分析时的误差公式相比只是将x(n) 计结果。 将v(n)= y(n) ? h(n)代入式(6.7-26) E=? [d(n)?y(n)?h(n)]2 n=0 ? =? [d(n)? ? h(k)y(n?k) ]2 n=0 ? k=0 N i=0,1,2,…,N 可以得到N+1个线性方程，表示为矩阵形式为 或表示为 ryy h =ryd (6.7-28) 即由 ?E ?h(i) =0 h(N) h(0) h(1) h(2) ??? ryy (0) ryy (N) ryy (N?1) ryy (N?2) ryy (0) ryy (1) ryy (1) ryy (2) ryy (N) ryy (N?1) ryy (0) ryy (1) ??? ??? ??? ??? ??? ??? ??? ??? ryd (0 ) ryd (1) ryd (2) ryd (N ) ??? = (6.7-28) ryy h =ryd 解此矩阵方程，求出滤波器的h(n) (n=0,1,2,…,N) 。 h(n)是有限项的，其z变换 (6.7-29) 式中ryy阵是一个(N+1) ? (N+1)的对称、正定(所有系 数大于等于零)方阵； h阵是一个N ? 1的列矩阵； ryd阵是一个N ? 1的列矩阵。 ? n=0 N h(n) z?n H(z)= 2)、第二步根据求出的h(n)再求出IIR DF的系统函数 H(z)。 (6.7-30) 要让(6.7-29)、(6.7-30)相等，从而找出由h(n)求 ak (k=1,2,…,N1)、bk (k=1,2,…,M)的关系式。 为了得到IIR DF的系统函数 H(z)= ? k=0 M bkz?k ? k=1 N akz?k 1+ 上式两边z-i项(0≤i≤M)的系数应相等，于是有 或 (6.7-32) (6.7-33) (0≤i≤M) 令 ? n=0 N h(n)z?n H(z)= (6.7-31) ? k=0 M bkz?k ? k=1 N1 akz?k 1+ = ? n=0 N h(n)z?n ? k=1 N1 akz?k 1+ ? k=0 M bkz?k = ? k=1 N1 ak h(i?k) =h(i)+ h(i)+ ? k=1 N1 akz?k bi= 其中N1 =N-M 由于bk (k=1,2,…,M)在k M后为零，即bi =0在(i M)， 代入式(6.7-32)可得到： 式中M为分子系数的个数； N1 =N? M为分母系数的 (M +1≤i≤N) (6.7-34) ? k=1 N1 ak h(i?k) bi= h(i)+ ? k=1 N1 ak h(i?k) = ?h(i) 个数。 由(6.7-34)式可以求出ak，然后代入(6.7-33)中求出bk。 最终求出H(z)的M+1个bk系数，N-M个ak系数。 时系数个数要少。 一般M的值较小，N的值较大。如果希望设计的系数再 少些，可以进一步压缩ak的个数（降低系统阶数）。 由此设计的系统函数系数总个数为N+1个，比E=0 令ak的个数为K个(K N?M)，即 (6.7-35) 为使式(6.7-35)逼近(6.7-29)式，考虑到(6.7-34)式，以 (l=0,1,2,…,N)， H(z)= ? k=0 M bkz?k ? k=1 K akz?k 1+ 最小，令 ? k=1 K ak h(i?k) h(i)+ i=M+1 及使误差 E= ? N 2 ?E ?al =0 由上式可列K个方程，解出K个ak系数。将K个ak系数 代入(6.7-33)，再求出M+1个bk系数。不过ak系数不能 随意减少，否则有可能导致系统不稳定。 ?E ?al ? k=1 K ak h(i?k) h(i?l) =0 2 h(i)+ i=M+1 ? N = 即 ? k=1 K ak h(i?k)h(i?l) =? h(i)h(i?l) i=M+1 ? N 得到 例6.7-2：已知hd(n)的长度为2，hd(0)=1，hd(1)=1/2，用 此时的 d(n)应为? (n)，与y(n)相比误差较大。 解： hd(n)相当于有一次回声的干扰滤波器，要求IIR 时域最小平方准则设计设计一个有两项系数的IIR DF。 反滤波，当输入为? (n)时，输出是 DF的系统函数H(z)，必须知道h (0)、h(1)。若不经过 y(n)=hd (n)=? (n)+0.5? (n?1) 计算托布里兹(Toeplitz)矩阵的各元素 经过反滤波器的输出为v(0)、v(1),…，v(n)应逼近? (n)。 又由于y