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

1/2 -1/2 1/3 -1/3 1/4 -1/4 x(n) y(n) z?1 z?1 z?1 将流图最左边连接上、下部增益为1的直通支路移至最右边； 通路的支路增益乘以?1； z?1放置在各基本单元的右下侧； 除了无延时的直通通路外，其余支路反向；指向无延时的直通 系数按 kM、kM-1、…、 k2、k1从左到右排列，特别的与kM有关 的支路(虚线)可以不要。 1/3 1/3 1/4 1/4 1/2 1/2 x(n) y(n) z-1 z-1 z-1 5.4.3 具有零、极点(IIR)的格型滤波器 具有零、极点的IIR系统的系统函数为 通常 M≤N 。 (5.4-16) H(z)= B(z) A(z) ? i=0 M bM z?i ? i=1 N aN z?i 1+ = (i) (i) M=N IIR系统的格型梯形结构如图5.4-10所示。 图中c1、c2、…、cM-1、cM为确定系统函数零点的梯形 系数。 y(n) c1 cM cM-1 c2 c0 x(n) z -1 z -1 z -1 z -1 p0 pM pM-1 p2 p1 qM qM-1 q2 q1 q0 k1 -k1 -kM kM -kM-1 kM-1 -k2 k2 1、若c0=1，而c1=c2=…=cM-1=cM=0，则图5.4-10是一个 由图5.4-10可见 2、若k1=k2=…=kM?1=kM=0，即所有反射支路开路，则图 k1、k2、…、kM?1、kM、仍按全极点系统的方法得出。 全极点的IIR格型结构。 5.4-10是一个全零点的FIR直接型结构。 3、由上述两点，图5.4-10的上半部分格型实现全极点 系统1/A(z)；下半部分梯形实现全零点系统B(z)。因下 半部分零点系统B(z)无反馈，对上半部分无影响，所以 …、cM-1、cM 与全零点系统时的方法有所不同。 现在的任务是求出系数 ci，i =0,1,2,3,…,M。 直接给出两种递推公式。 i =0,1,2,3,…,M (5.4-17) m=0,1,2,3,…,M?1 (5.4-18b) 而由于上半部分对下半部分有影响，因此求c1、c2、 (2) (5.4-18a) cM = bM cm = bm (m) (1) ci = bi ? cm a m m?i ? m=i+1 M 其中bM是已知的。 式系数{ am }、{ bm }，求出格型梯形结构系数{ km }、 利用MATLAB函数dir2ladr，可以由已知的零、极点形 { cm }。 利用MATLAB函数ladr2dir，还可以由已知的格型梯形 结构系数{ km }、{ cm }求出零、极点形式系数{ am }、 { bm }。 例5.4-5 已知某IIR的系统函数为 求其格型结构系数并作图。 解 求解例5.4-5的MATLAB的程序与结果如下。 b=[1 2 2 1]; a=[1 13/24 5/8 1/3]; [K,C]= dir2ladr(b,a) C =-0.2695 0.8281 1.4583 1.0000 答案 K =0.2500 0.5000 0.3333 1+2z?1 +2z?2 +z?3 H(z)= 1?(13/24)z?1 +(5/8)z?2 +(1/3)z?3 例5.4-5的格型结构如图所示。 1/4 -1/4 1/3 -1/3 1.4583 0.8281 -0.2695 1/2 -1/2 x(n) y(n) z -1 z -1 z -1 即 k1=1/4，k2=1/2，k3=1/3； c0=-0.2695，c1=0.8281，c2=1.4583，c3=1。 a =1.0000 0.5417 0.6250 0.3333 例5.4-6 用MATLAB函数ladr2dir验算例5.4-5的结果。 解 例5.4-6 的MATLAB程序及结果如下 K=[1/4 1/2 1/3]; C=[-0.2695,0.8281,1.4583,1]; [b,a]= ladr2dir(K,C) 答案 b = 1.0000 2.0000 2.0000 1.0000 * 在线教务辅导网： 更多课程配套课件资源请访问在线教务辅导网 敏感，以及一个M阶的格型滤波器可以产生从1阶到M 结构便于高速并行处理外，还具有对有限字长效应不 §5.4格型滤波器结构 广泛应用在功率谱估计、语音处理、自适应滤波、线 性预测、逆滤波等方面的格型滤波器，除了其模块化 阶的M个横向滤波器的输出性能。本节分别讨论全零 点、全极点、以及一般IIR系统的格型滤波器。 一个M 阶直接形式FIR滤波器的系统函数为 5.4.1 全零点(FIR)的格型滤波器 如用FIR直接结