第九章_数字滤波器的分类及结构.ppt

数字信号处理,数字滤波器的实质是用有限精度算法实现的离散时间LSI系统，以完成对信号进行滤波处理的功能。其输入是一组由模拟信号经过抽样和量化的数字信号，输出是经过处理的另一组数字信号。数字滤波器既可以是一台由数字硬件装配成的用于完成滤波计算功能的专用机，也可以是由通用计算机完成的一组运算程序。本章主要介绍数字滤波器的分类及结构。,引言,第九章数字滤波器的分类及结构,数字滤波器的分类数字滤波器结构的表示方法IIR滤波器的结构FIR滤波器的结构离散时间系统的Lattice结构,根据单位冲激响应h(n)的时间特性分类：无限冲激响应数字滤波器（IIR）有限冲激响应数字滤波器（FIR）根据实现方法和形式分类:递归型数字滤波器非递归型数字滤波器根据频率特性分类：低通数字滤波器高通数字滤波器带通数字滤波器带阻数字滤波器,1.滤波器的分类,数字滤波器可以用一个差分方程来描述：由上式可以看出，实现一个数字滤波器需要几种基本的运算单元：加法器、单位延迟和常数乘法器。这些基本的单元可以有两种表示方法：方框图法和信号流图法，因此一个数字滤波器也可也有两种表示方法：方框图法和信号流图法。,2.数字滤波器结构的表示方法,考虑如下二阶数字滤波器的信号流图：x(n)处称为输入节点或源节点，y(n)处称为输出节点或阱节点，其余节点称为网络节点。节点之间用有向支路连接，每个节点可以有几条输入支路和几条输出支路，节点值等于它所有输入支路的信号之和，而输入支路的信号值等于这一支路起点处的节点信号值乘以之路上的传输系数。延迟算子z-1表示单位延迟。,2.数字滤波器结构的表示方法,源节点没有输入支路，阱节点没有输出支路。如果某节点有一个输入，一个或多个输出，该节点称为分支节点。如果某节点有两个或两个以上的输入，该节点称为相加器。各节点值为：对分支节点2有，故,2.数字滤波器结构的表示方法,IIR滤波器的特点：单位冲激响应h(n)是无限长的；系统函数在有限Z平面上（0|z|）有极点存在；结构上存在着输出到输入的反馈。,3.IIR滤波器的结构,数字滤波器可用差分方程来描述：也可以用系统函数来表示：,3.IIR滤波器的结构,表示为两个系统级联的形式：称为直接型结构。,3.IIR滤波器的结构,H1(z),H2(z),x(n),y(n),y(n),直接型的变型：,3.IIR滤波器的结构,H2(z),H1(z),x(n),y(n),y(n),z-1,-aN-1,b1,x(n),y(n),b0,z-1,z-1,b2,bM-1,bM,z-1,z-1,z-1,-a2,-a1,-aN,y(n),直接型结构（典范型）:,3.IIR滤波器的结构,H2(z),H1(z),x(n),y(n),y(n),习题：已知数字滤波器的系统函数画出该滤波器的直接型结构。解答：如右图所示。直接型结构的特点:所需要的延迟单元最少；系统调整不方便；受有限字长影响较大。,3.IIR滤波器的结构,对系统函数H(z)进行因式分解：式中M=M1+2M2,N=N1+2N2。级联型结构图：,3.IIR滤波器的结构,习题：已知数字滤波器的系统函数画出该滤波器的级联型结构。解答：级联型结构的特点:所用存储单元较少；系数调整方便，便于准确实现系统的零极点；受有限字长影响较大。,3.IIR滤波器的结构,对系统函数H(z)进行因式分解：上式中N=N1+2N2。当MN时，公式中不包含最后一项。当M=N时，最后一项变成G0。一般IIR系统皆满足M0处收敛，极点全部在z=0处；结构上不存在输出到输入的反馈。,4.FIR滤波器的结构,数字滤波器可用一个差分方程来描述：对于FIR滤波器则有：直接型结构：,4.FIR滤波器的结构,把FIR滤波器的系统函数用二阶因子乘积表示：级联型结构：,4.FIR滤波器的结构,前已证明，当FIR系统的单位冲激响应满足时，该系统具有线性相位。类型滤波器：，且N为奇数。,4.FIR滤波器的结构,令m=N-1-n,类型滤波器的结构：,z-1,z-1,4.FIR滤波器的结构,类型滤波器：，且N为偶数。,z-1,z-1,z-1,h(N/2-1),4.FIR滤波器的结构,类型滤波器：，且N为奇数。,z-1,z-1,4.FIR滤波器的结构,类型滤波器：，且N为偶数。,z-1,z-1,-z-1,h(N/2-1),4.FIR滤波器的结构,给定一个FIR系统的单位冲激响应为h(n),n=0,1,N-1,其系统函数和h(n)的DFT分别为：显然，H(k)实际上是H(z)在单位圆上的N个值，即H(k)是H(j)在频域的抽样。因此，我们可以用H(k)来表示H(z)，即:,4.FIR滤波器的结构,H(z)可看做是两个子系统级联。一个是FIR子系统H1(z)，一个是IIR子系统H2(z)。FIR子系统由N个延时单元组成，系统函数为H1(z)，该系统在单位圆上有N个等分的零点：频率响应：,4.FIR滤波器的结构,2/N,0,“梳状滤波器”,IIR子系统由N个一阶系统并联组成，系统函数为该系统有N个极点：IIR系统与FIR系统级联后，N个IIR系统在单位圆上的极点正好和FIR系统在单位圆上的零点相互抵消，所以整个系统是FIR系统，称为FIR系统的频率抽样型结构。,4.FIR滤波器的结构,4.FIR滤波器的结构,频率抽样型结构的优点：系统在频率采样点=2k/N上的响应等于H(k)，因此改变H(k)就等于改变了系统的频率响应；只要h(n)的长度N相同，不论频率响应如何，梳状滤波器以及个N一阶网络的结构相同，便于标准化和模块化。缺点：稳定性差；有限字长效应导致零极点不能完全相消，从而造成系统不稳定。结构中系数为复数，因此运算量大。利用H(k)的对称性可以一定程度上降低计算量。,4.FIR滤波器的结构,快速卷积结构：NNh+Nx-1,N=2m,4.FIR滤波器的结构,h(n)Nh点序列,x(n)Nx点序列,X(k),H(k),y(n),一个M阶的FIR系统的系统函数可以写作：上式中假定H(z)=B(z)的首项系数等于1，表示M阶FIR系统的第i个系数，该系统的Lattice结构如下图：,5.离散时间系统的Lattice结构,Lattice结构（格型结构）的特点：Lattice结构有M个参数，共需2M次乘法，M次延迟。直接型结构有M个参数，共需M+1次乘法，M次延迟。,5.离散时间系统的Lattice结构,Lattice结构,直接型结构,信号的传递是从左至右，中间没有反馈回路，所以这是一个FIR系统。,5.离散时间系统的Lattice结构,信号流图中的基本单元具有如下关系：其中：,5.离散时间系统的Lattice结构,与FFT算法中的蝶形单元做比较：,定义：为输入端x(n)至第m个基本单元后所对应系统的系统函数，对应上端，对应下端。,5.离散时间系统的Lattice结构,如何由给定的系数b(1),b(2),b(M)求出Lattice结构的参数k1,k2,kM呢？两端分别除以P0(z)和Q0(z),得：,5.离散时间系统的Lattice结构,由前面的定义可知：，则有归纳得出：,5.离散时间系统的Lattice结构,将的定义带入上式类似的有,5.离散时间系统的Lattice结构,可按照如下步骤求出参数k1,k2,kM：步骤1由上述关系式，首先得到步骤2由kM及系数求出的系数，或者由下式直接求出，那么。步骤3重复步骤2，即可全部求出。,5.离散时间系统的Lattice结构,习题：一个FIR系统的零点分别在和0.8处，求其Lattice结构。解答：首先写出系统函数可知：,5.离散时间系统的Lattice结构,因此有：可知：进一步：最后：Lattice结构：,5.离散时间系统的Lattice结构,5.离散时间系统的Lattice结构,-km,km,IIR系统（全极点）的Lattice结构图：,5.离散时间系统的Lattice结构,令M=1，则：,5.离散时间系统的Lattice结构,令M=2，则：,5.离散时间系统的Lattice结构,5.离散时间系统的Lattice结构,由此类推，若定义：则有：,对比FIR系统与IIR系统：结论：参数k1,k2,kM的求解方法同FIR系统完全一样。,5.离散时间系统的Lattice结构,km,-km,习题：一个IIR系统的系统函数为，求其Lattice结构。解答：已求出则：,5.离散时间系统的Lattice结构,一个既具有零点又具有极点的IIR系统的Lattice结构：,5.离散时间系统的Lattice结构,若，即为全极点系统。,5.离散时间系统的Lattice结构,y(n),若，则为FIR系统的直接实现形式。,5.离散时间系统的Lattice结构,y(n),根据前面的定义：其中是由至之间的系统函数，令为由至之间的系统函数，那么,5.离散时间系统的Lattice结构,整个系统的的转移函数应该是加权后的并联，即：,5.离散时间系统的Lattice结构,如何求解？令M=3，,5.离散时间系统的Lat