数字滤波器的原理ppt课件

1、数字滤波器的原理与结构14-1 数字滤波器的原理数字滤波器的原理24-2 数字滤波器的分类数字滤波器的分类34-3 数字滤波器的表示数字滤波器的表示44-4 IIR数字滤波器的基本网络结构数字滤波器的基本网络结构54-5 FIR数字滤波器的基本网络结构数字滤波器的基本网络结构 4-1 数字滤波器的原理 滤波的定义：对输入信号通过一定的处理得输出信号，这个处理通常是滤波的定义：对输入信号通过一定的处理得输出信号，这个处理通常是滤除输入信号的某些频率成分；保留信号中某些频率范围内的有用信号成分。滤除输入信号的某些频率成分；保留信号中某些频率范围内的有用信号成分。所以把这种处理的过程称为滤波。所以把

2、这种处理的过程称为滤波。 滤波器的定义：实现滤波处理的运算电路、或设备称为滤波器。滤波器的定义：实现滤波处理的运算电路、或设备称为滤波器。 数字滤波的定义：对数字信号序列通过一定的算法进行处理，提取信数字滤波的定义：对数字信号序列通过一定的算法进行处理，提取信号中某频率范围内的信号成分的过程称为数字滤波。号中某频率范围内的信号成分的过程称为数字滤波。 数字滤波器是数字信号处理中使用得最广泛的一种线性系统环节，是数字信号处理的重要基础。 数字滤波器的本质是将一组输入的数字序列通过一定的运算后转变为另一组输出的数字序列。 4-2 数字滤波器的分类按系统冲击响应或差分方程分类：按计算方法分类: 递归

3、系统 ，非递归系统按冲击响应长度分类：IIR ，FIR按频带分类： 低通, 高通, 带通, 带阻 可以分成无限冲击响应 IIR 和有限冲击响应FIR滤波器两类。这两种滤波器都可以现实各种基本滤波器频率特性要求，但它们在计算流程、具体特性逼近等方面是有差别的。4-3 数字滤波器结构的表示 数字滤波器的系统函数：01( )( )( )1MkkkNkkkb zY zH zX za z10( )()()NMkkkky na y nkb x nk常系数线性差分方程:aa1z1z加法器常数乘法器单位延时基本运算单元基本运算单元方框图方框图流图流图10( )()()NMkkkky na y nkb x nk

4、例：二阶数字滤波器120( )(1)(2)( )y na y na y nb x n方框图结构方框图结构流图结构流图结构节点 源节点源节点 支路支路 阱节点阱节点 网络节点网络节点 分支节点分支节点 输入支路输入支路 相加器相加器节点的值=所有输入支路的值之和 输出支路输出支路支路的值=支路起点处的节点值 传输系数数字滤波器的设计步骤:1按照实际需要确定滤波器的性能要求； 2用一个因果稳定的系统函数传递函数去逼近这个性能要求，这种传递函数可分为两类：IIR和FIR。 3用一个有限精度的运算去实现这个传递函数。包括选择运算结构：如级联型、并联型、卷积型、频率采样型以及快速卷积FFT型等，及选择合

5、适的字长和有效的数字处理方法等。4-4 IIR数字滤波器的基本网络结构1系统的单位抽样相应系统的单位抽样相应h(n)无限长无限长IIR数字滤波器的特点：数字滤波器的特点：3存在输出到输入的反馈，递归型结构存在输出到输入的反馈，递归型结构2系统函数系统函数H(z)在有限在有限z平面（平面（ ）上有极点存在）上有极点存在0z 01( )( )( )1MkkkNkkkb zY zH zX za z系统函数： 10( )()()NMkkkky na y nkb x nk差分方程： IIR数字滤波器的基本结构： 直接型 直接型典范型） 级联型 并联型差分方程:10( )()()NMkkkky na y

6、nkb x nk需需N+M个个延时单元延时单元上述结构缺点：上述结构缺点： 需要需要N+M个延迟器个延迟器z-1），太多。），太多。 系数系数ak、bk对滤波器性能的控制不直接，对极、对滤波器性能的控制不直接，对极、 零点的控制难，一个零点的控制难，一个ak、bk的改变会影响系统的的改变会影响系统的 零点或极点分布。零点或极点分布。 对字长变化敏感对对字长变化敏感对ak、bk的准确度要求严格）。的准确度要求严格）。 易不稳定，运算的累积误差较大，阶数高时，上述易不稳定，运算的累积误差较大，阶数高时，上述影响更影响更 大。大。NM只需实现只需实现N阶滤波器所需的最少的阶滤波器所需的最少的N个延时

7、单元，故称典范型。个延时单元，故称典范型。（ ） 直接直接II型优缺点：型优缺点： 优点：延迟线减少优点：延迟线减少M个，为个，为N个，可节省寄存个，可节省寄存器或存储单元。器或存储单元。 缺陷：同直接型。缺陷：同直接型。 通常在实际中很少采用上述两种结构实现高通常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合的低阶阶系统，而是把高阶变成一系列不同组合的低阶系统一、二阶来实现。系统一、二阶来实现。例1、用直接I型及典范结构实现以下系统函数： 121234.20.820.60.4zzH zzz=解：根据IIR滤波器的系统函数标准式 12121.52.10.410.30.2

8、zzH zzz=将系统函数整理为： 011MmmmNnnnb zY zH zX za z=12121.52.10.410.30.2zzzz 10.3a 20.2a 得 ，01.5b 12.1b 20.4b 直接I型结构：典范型结构：12121.52.10.4( )10.30.2zzH zzz 3、级联型将系统函数按零极点因式分解:121211*101111*1111(1)(1)(1)( )1(1)(1)(1)MMMkkkkkkkkNNNkkkkkkkkb zp zq zq zH zAa zc zd zd zA为常数*,kkkkq qd d和分别为复共轭零、极点kkpc和 分别为实数零、极点12

9、2MMM122NNN将共轭成对的复数组合成二阶多项式，系数即为实数。为采用相同结构的子网络，也将两个实零点/极点组合成二阶多项式121212121( )( )1kkkkkkkzzH zAAHzzz20k当零点为奇数时：当零点为奇数时： 有一个有一个20k当极点为奇数时：当极点为奇数时： 有一个有一个12NMN当时，共有节121212121( )( )1kkkkkkkzzH zAAHzzz1!2N 各二阶基本节的排列次序有 种1!2N 当M=N时，二阶因子配对方式有 种级联型结构的优缺点：级联型结构的优缺点：优点：优点： 简化实现，用二阶基本节，通过变换系数就可实现整个简化实现，用二阶基本节，通

10、过变换系数就可实现整个 系统；系统； 极、零点可单独控制、调整，调整极、零点可单独控制、调整，调整 、 可单独调整第可单独调整第 k 对零点，调整对零点，调整 、 可单独调整第可单独调整第k对极点，从而便于调对极点，从而便于调 整滤波器频率响应性能整滤波器频率响应性能 各二阶基本节零、极点的搭配可互换位置，优化组合以减小各二阶基本节零、极点的搭配可互换位置，优化组合以减小 运算误差；运算误差； 可流水线操作。可流水线操作。 运算的累积误差较小运算的累积误差较小 具有最少的存储器具有最少的存储器缺陷：缺陷： 二阶节电平难控制，电平大易导致溢出，电平小则使信噪二阶节电平难控制，电平大易导致溢出，电

11、平小则使信噪比减小。比减小。1k2k1k2k解： 1121124 11 1.410.510.90.8zzzzzz 1212121211kkkkkzzH zAzz111210110.5210121.4 221220.8 120.9 4A那么例例2 2、用级联型结构实现以下系统函数：、用级联型结构实现以下系统函数： 22411.410.50.90.8zzzH zzzz=试问一共能构成几种级联型网络。考虑分子分母的组合及级联的次序，共有以下四种级联型网络：111210110.5210121.4 221220.8 120.9 将因式分解的H(z)展成部分分式：111220100121112( )( )

12、1NNkkkkkkkzH zGGHzzz210kk当N为奇数时，有一个()MN1210101121112( )11NNkkkkkkkkAzH zGc zzz122NNN组合成实系数二阶多项式：111220100121112( )( )1NNkkkkkkkzH zGGHzzz27并联型的特点： 通过调整系数通过调整系数 ， 可单独调整一对极点位置。可单独调整一对极点位置。1k2k 各并联基本节的误差互相不影响，故运算误差各并联基本节的误差互相不影响，故运算误差最小最小 可同时对输入信号进行运算并行运算），故可同时对输入信号进行运算并行运算），故运算速度最高运算速度最高 缺陷：缺陷： 不能直接调整

13、零点，因多个二阶节的零不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。传输零点时，级联型最合适。例例3 3、给出以下系统函数的并联型实现：、给出以下系统函数的并联型实现： 1231125.21.581.411.610.510.90.8zzzH zzzz=解：对此函数进行因式分解并展成部分分式，得 1231125.21.581.411.610.510.90.8zzzH zzzz11120.210.341 0.510.90.8zzzz那么04G 110.5210010.2110120.9 220.8 021

14、120.3 1 /21010121121NkkkkkzH zGzz则并联结构：04G 110.5210010.2110120.9 220.8 021120.3转置定理： 原网络中所有支路方向倒转，并将输入x(n)和输出y(n)相互交换，则其系统函数H(z)不改变。例例4 4：设：设IIRIIR数字滤波器差分方程为：数字滤波器差分方程为：试用四种基本结构实现此差分方程。试用四种基本结构实现此差分方程。( )8 ( )4 (1) 11 (2)2 (3)y nx nx nx nx n531(1)(2)(3)448y ny ny n 123123841125311448zzzH zzzz解：对差分方程

15、两边取z变换，得系统函数： 123123841125311448zzzH zzzz得直接型结构：典范型结构： 11211220.37941.245.264111142zzzH zzzz1121128 10.1910.311.32111142zzzzzz将H(z)因式分解：得级联型结构： 11128162016111142zH zzzz将H(z)部分分式分解：得并联型结构：4-5 FIR数字滤波器的基本网络结构1系统的单位抽样响应系统的单位抽样响应 h(n)有限长，设有限长，设N点点 FIR数字滤波器的特点：数字滤波器的特点：0z 2 2系统函数系统函数H(z)H(z)在在 处收敛，有限处收敛，

16、有限z z平面只有零点，平面只有零点，全部极点在全部极点在 z = 0 z = 0 处因果系统）处因果系统）3无输出到输入的反馈，一般为非递归型结构无输出到输入的反馈，一般为非递归型结构10( )( )NnnH zh n z系统函数：系统函数：z=0z=0处处 是是N-1N-1阶极点阶极点有有N-1个零点分布于个零点分布于z平面平面10( )()()NMkkkky na y nkb x nk01( )( )( )1MkkkNkkkb zY zH zX za z10( )( )NnnH zh n z10( )( ) ()Nmy nh m x nm1、横截型卷积型、直接型）差分方程:10( )(

17、)NnnH zh n z10( )( ) ()Nmy nh m x nm例例1 1 用横截型结构实现以下系统函数：用横截型结构实现以下系统函数： 1111111116121126H zzzzzz=解： 1111111116121126H zzzzzz=1121121111216126zzzzzzz=1212153711126zzzzz=12345820520581312123zzzzz 1234582052058( ) 1312123H zzzzzz=则横截型结构：/211201201( )( )()NNnkkknkH zh n zzzN为偶数时，其中有一个为偶数时，其中有一个 （N-1个零点

18、）个零点）20k将将H(z)分解成实系数二阶因式的乘积形式分解成实系数二阶因式的乘积形式:级联型的特点每个基本节控制一对零点，便于控制滤波器的传输零点系数比直接型多，所需的乘法运算多3、频率抽样型1101( )( )(1)1NNkkNH kH zzNWz101( )( )NckkHzHzNN个频率抽样个频率抽样H(k)恢复恢复H(z)的内插公式：的内插公式：2jkNkze0,1,.,1kN()1jj NcH ee 22sin2NjNje( )1NcHzz 子系统：子系统： 是是N节延时单元的梳状滤波器节延时单元的梳状滤波器在单位圆上有在单位圆上有N个等间隔角度的零点：个等间隔角度的零点：频率响

19、应：频率响应：222NNNjjjeee2jkkNkNzWe单位圆上有一个极点：2kN与第k个零点相抵消，使该频率 处的频率响应等于H(k)1( )( )1kkNH kHzWz谐振器子系统： 频率抽样型结构的优缺点调整H(k)就可以有效地调整频响特性 若h(n)长度相同，则网络结构完全相同，除了各支路增益H(k)，便于标准化、模块化 有限字长效应可能导致零极点不能完全对消，导致系统不稳定 系数多为复数，增加了复数乘法和存储量例例2 2 用频率抽样结构实现以下系统函数：用频率抽样结构实现以下系统函数： 3615231zzH zz=抽样点数 6N ，修正半径 0.9r 解：由N = 6，得频率抽样型

20、结构： 266031116kkH zr zHzHzHz=312531zzz112 11zzz 3315311zzH zz=又233531jkjkj keee 024H即 122 3Hj 20H 32H 40H 522 3Hj那么 011024110.9HHzrzz 31132110.9HHzrzz 312 531H zzzz由得 2zk NH kH z然后求 kHz时 1k 101111122212coszHzz rr zN=4 111620.9Re13.6HW 111243.6 10.90.81zHzzz 101122212 coskkkzHzzrkr zN=其中 02RekH k 12 RekkNrH k W 122 3Hj 012Re1H 101122212 cosk