1第四章数字滤波器的原理和设计方法PrinciplesandDesignofDigitalFilters2020/1/202本章作业习题4.24.34.4(1)4.84.104.114.134.183第一节概述Introduction4概述Introduction什么是数字滤波器?数字滤器相对于模拟滤波器的优点?数字滤器的如何分类?5概述Introduction数字滤波器的功能就是把输入序列通过一定的运算,变成符合要求的输出序列。数字滤波器可以用差分方程来表示,即对应的系统函数为knxbknyanyMkkNkk01NkkkMkkkzazbzXzYzH1016概述Introduction与模拟滤波器比较,数字滤波器有以下优点:精度和稳定性高改变系统函数比较容易,因而比较灵活不存在阻抗匹配问题便于大规模集成可以实现多维滤波7概述Introduction数字滤波器一般可以用两种方法实现:一种是根据描述数字滤波器的数学模型或信号流图,用数字硬件构成专用的数字信号处理机;另一种是编写滤波运算程序,在计算机上运行。8概述Introduction由前面的滤波器差分方程看出,实现数字滤波器需要3种基本运算单元,即加法器、单位延迟器和常数乘法器。这些单元有方框图和信号流程图两种表示法。方框图和流程图两种表示法中的3种基本运算单元如下图所示。9概述Introduction方框图和信号流程图nx1nx2nxnx21nx1z1nxnxanxanx1nx2nxnx211znx1nxnxanxa加法器单位延迟器常数乘法器10概述Introduction数字滤波器的结构一般也有方框图和信号流程图两种表示方法,我们主要用信号流程图表示法。设一个二阶的数字滤波器的差分方程为其方框图和信号流程图如下图所示。nxbnyanyany0212111概述Introductionnxbnyanyany02121nx0b1z1z1a2any1z1znx0b1a2any信号流程图结构方框图结构12概述Introduction数字滤波器根据它的单位取样响应是一个有限长序列还是一个无限长序列可以分成有限冲激响应(FIR)和无限冲激响应(IIR)滤波器两种有限冲激响应(FIR)的方框图:DDD00ba30ba10ba20ba()yn()xn(1)xn(3)xn(2)xn13第二节无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters14无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters无限冲激响应(IIR)数字滤波器有哪些基本网络结构?不同的网络结构有什么特点?15无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters用什么结构实现数字滤波器是一个重要的研究课题。运算结构的不同将会影响系统运算的精度、误差、速度和经济性等性能指标。在一般情况下,都要求使用尽可能少的常数乘法器和延迟器来实现系统,并要求运算误差尽可能小。然而,这些要求有时是互相矛盾的,例如,为了获得具有较小运算误差的结构,使用的乘法器和延迟器的数目往往不是最少的。对于同一个系统函数,可以有多种不同的结构。下面讨论IIR系统的几种结构。16无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFiltersIIR数字滤波器是一种递归系统,描述它的差分方程如下式所示下面两图分别是IIR数字滤波器的方框图和流程图,称为IIR数字滤波器的直接I型结构。(假设M=N)knxbknyanyMkkNkk011.IIR系统的直接I型结构17无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters1z1znx0b1a2any1z1NaNa1ny2ny1NnyNny1z1z1z1nx2nx1NnxNnx1b2b1NbNb型方框图系统的直接阶IIIRNnQ18无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFiltersnx0b1a2any1NaNa1z1b2b1NbNb1z1z1z1z1znQ型信号流程图系统的直接阶IIIRN19无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters可以看出,这种结构可以看成是两部分级联组成左边部分对应的是差分方程等号右边的第二项,或系统函数的分子,实现的是系统函数的各零点。右边部分对应的是差分方程等号右边的第一项,或系统函数的分母,实现的是系统函数的各极点。直接I型需要2N个延迟器和2N个乘法器。20无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters2.IIR系统的直接II型结构10()NkkkHzbz我们所讨论的IIR滤波器是线性非移变系统,所以2112HzHzHzHzHz211()1NkkkHzaz21无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFiltersnx0b1a2any1NaNa1z1b2b1NbNb1z1z1z1z1znwnwII调换直接型级联位置后的流程图22无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters上图所示的中间两条延迟线都起着对中间变量进行延迟的作用,因此可以进行合并,于是得到下图所示的直接II型结构。nwnx0b1a2any1NaNa1z1b2b1NbNb1z1z型信号流程图系统的直接IIIIR23无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters比较直接I型结构和直接II型结构可以看出,直接II型结构首先实现系统的各极点,然后实现各零点;对于一个N阶数字滤波器,直接II型结构仅需要N个延迟单元。因此,用硬件实现数字滤波器时,直接II型结构比直接I型结构少用寄存器;用软件实现时,直接II型少用存储单元。24无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters3.IIR系统的级联结构如果将N阶IIR系统函数分解成二阶因式连乘积,则可得到级联结构,即这样,整个系统将由M个二阶系统级联构成,如下图所示zHzHzHzHM21zH1zH2zHMnxny个二阶系统级联系统分解成阶将MIIRN25无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters将系统函数的分子分母多项式都进行因式分解,得到由于的系数都为实数,所以的零点和极点或者是实数,或者是共轭复数。22111111111111101111111NkkkMkkkNkkMkkNkkkMkkkzdzdzhzhzczgAzazbzHzHzH26无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters若将上式每对共轭复数因子合并起来构成一个实系数二阶因子,则得下式2211122111221111111111NkkkMkkkNkkMkkzzzzzczgAzH27无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters如果把实数一阶因子看作二阶实数因子的特列,即看作二次项系数等于零的二阶因子,则可以看成全部由实系数二阶因子连乘来构成即:其中称为滤波器的二阶基本节。zHMkkMkkkkkzHAzzzzAzH1122112211112211221111zzzzzHkkkkk28无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters如果每个二阶基本节都采用直接II型结构来实现,那么整个滤波器的级联结构如下图所示。nxA1121ny1z111z1z1b2b1z1z1z2112221222M1M2M1M2型结构的级联形式每个二阶节都用直接IISeconorderSection29无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters级联结构的特点二阶基本节搭配灵活,可以按实际需要调换二阶节的次序。直接控制系统的零点和极点。因为每个二阶基本节有一对零点和一对极点,所以通过调整系数就可以独立调整第k个二阶节的零点而不影响其它零点,调整系数就可以独立调整第k个二阶节的极点。kk21,kk21,30无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters4.IIR系统的并联结构如果将IIR系统的系统函数化成部分分式之和,则可得到IIR数字滤波器的并联结构。QkkkkkPkkkQkkkkkPkkkNkkkMkkkzzzzcAczdzdzeBzcAczazbzH122111101101111110101111111zH31无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters所以,滤波器可由P个一阶网络、Q个二阶网络和一个常支路并联构成。Q1Q2Q011210c1cPc1z1z1z1z1z1zQ111011APAnxny系统的并联结构IIR32无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFilters网络转置定理如果将原信号流图中所有支路的方向加以反转,并将输入和输出相互交换,则新的信号流程图表示的系统的系统函数不会改变。利用转置定理,可以将前面讨论的各种IIR结构进行转置处理,从而得到各种新的实现结构。下图表示的是其中一种网络的转置处理情况,即将直接II型结构转置得到的结果。33无限冲激响应数字滤波器的基本网络结构BasicStructureofIIRDigitalFiltersnx0b1a2any1NaNa1z1b2b1NbNb1z1znx0b1a2any1NaNa1z1b2b1NbNb1z1z型结构转置处理系统直接IIIIR右图画成习惯的左输入右输出形式