1、第五章学习目标第五章学习目标理解数字滤波器结构的表示方法掌握IIR滤波器的基本结构掌握FIR滤波器的直接型、级联型、线性相位结构，理解频率抽样型结构了解数字滤波器的格型结构本章作业练习本章作业练习 P226: 1234678 (1)第五章第五章 数字滤波器的基本结构数字滤波器的基本结构 一、数字滤波器结构的表示方法 数字滤波器的系统函数：01( )( )( )1MkkkNkkkb zY zH zX za z10( )()()NMkkkky na y nkb x nk常系数线性差分方程:aa1z1z10( )()()NMkkkky na y nkb x nk加法器常数乘法器单位延时基本运算单元方
2、框图流图例：二阶数字滤波器120( )(1)(2)( )y na y na y nb x n方框图结构流图结构流图结构节点 源节点支路 阱节点 网络节点 分支节点 输入支路 相加器节点的值=所有输入支路的值之和 输出支路支路的值=支路起点处的节点值 传输系数二、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差分方程： II
3、R数字滤波器的基本结构： 直接型 直接型（典范型） 级联型 并联型1、直接型差分方程:10( )()()NMkkkky na y nkb x nk需N+M个延时单元2、直接型（典范型）NM只需实现N阶滤波器所需的最少的N个延时单元，故称典范型。（ ） 直接型的共同缺点：kakb系数 ， 对滤波器的性能控制作用不明显极点对系数的变化过于灵敏，易出现不稳定或较大误差运算的累积误差较大3、级联型将系统函数按零极点因式分解:121211*101111*1111(1)(1)(1)( )1(1)(1)(1)MMMkkkkkkkkNNNkkkkkkkkb zp zq zq zH zAa zc zd zd z
4、A为常数*,kkkkq qd d和分别为复共轭零、极点kkpc和 分别为实数零、极点122MMM122NNN将共轭成对的复数组合成二阶多项式，系数即为实数。为采用相同结构的子网络，也将两个实零点/极点组合成二阶多项式121212121( )( )1kkkkkkkzzH zAAHzzz20k当零点为奇数时： 有一个20k当极点为奇数时： 有一个12NMN当时，共有节121212121( )( )1kkkkkkkzzH zAAHzzz1!2N 各二阶基本节的排列次序有 种1!2N 当M=N时，二阶因子配对方式有 种级联型的特点：调整系数 ， 能单独调整滤波器的第k对零点，而不影响其它零极点1k2k
5、运算的累积误差较小具有最少的存储器便于调整滤波器频率响应性能1k2k调整系数 ， 能单独调整滤波器的第k对极点，而不影响其它零极点4、并联型将因式分解的H(z)展成部分分式：111220100121112( )( )1NNkkkkkkkzH zGGHzzz210kk当N为奇数时，有一个()MN1210101121112( )11NNkkkkkkkkAzH zGc zzz122NNN组合成实系数二阶多项式：111220100121112( )( )1NNkkkkkkkzH zGGHzzz并联型的特点：通过调整系数 ， 可单独调整一对极点位置，但不能单独调整零点位置1k2k各并联基本节的误差互相不
6、影响，故运算误差最小可同时对输入信号进行运算，故运算速度最高转置定理：原网络中所有支路方向倒转，并将输入x(n)和输出y(n)相互交换，则其系统函数H(z)不改变。例：设IIR数字滤波器差分方程为：试用四种基本结构实现此差分方程。( )8 ( )4 (1)11 (2)2 (3)y nx nx nx nx n531(1)(2)(3)448y ny ny n 123123841125311448zzzH zzzz解：对差分方程两边取z变换，得系统函数： 123123841125311448zzzH zzzz得直接型结构：典范型结构： 11211220.37941.245.264111142zzzH
7、 zzzz1121128 10.1910.311.32111142zzzzzz将H(z)因式分解：得级联型结构： 11128162016111142zH zzzz将H(z)部分分式分解：得并联型结构：三、FIR数字滤波器的基本结构1）系统的单位抽样响应 h(n)有限长，设N点 FIR数字滤波器的特点：0z 2）系统函数H(z)在 处收敛，有限z平面只有零点，全部极点在 z = 0 处（因果系统）3）无输出到输入的反馈，一般为非递归型结构10( )( )NnnH zh n z系统函数：z=0处 是N-1阶极点有N-1个零点分布于z平面10( )()()NMkkkky na y nkb x nk0
8、1( )( )( )1MkkkNkkkb zY zH zX za z10( )( )NnnH zh n z10( )( ) ()Nmy nh m x nm1、横截型（卷积型、直接型）差分方程:10( )( ) ()Nmy nh m x nm2、级联型/211201201( )( )()NNnkkknkH zh n zzzN为偶数时，其中有一个 （N-1个零点）20k将H(z)分解成实系数二阶因式的乘积形式:级联型的特点系数比直接型多，所需的乘法运算多每个基本节控制一对零点，便于控制滤波器的传输零点3、频率抽样型1101( )( )(1)1NNkkNH kH zzNWz101( )( )Nckk
9、HzHzNN个频率抽样H(k)恢复H(z)的内插公式：2jkNkze0,1,.,1kN()1jj NcH ee 22sin2NjNje( )1NcHzz 子系统： 是N节延时单元的梳状滤波器在单位圆上有N个等间隔角度的零点：频率响应：222NNNjjjeee2jkkNkNzWe单位圆上有一个极点：2kN与第k个零点相抵消，使该频率 处的频率响应等于H(k)1( )( )1kkNH kHzWz谐振器子系统： 频率抽样型结构的优缺点调整H(k)就可以有效地调整频响特性若h(n)长度相同，则网络结构完全相同，除了各支路增益H(k)，便于标准化、模块化有限字长效应可能导致零极点不能完全对消，导致系统不
10、稳定系数多为复数，增加了复数乘法和存储量 修正频率抽样结构1101( )( )(1)1NNNrkkNH kH zr zNrWz2jkNkzre极点： 0,1,.,1kN11rr且将零极点移至半径为r的圆上：为使系数为实数，将共轭根合并*N kkzz()*()N kkkNNWWW由对称性：*( )()( )NNH kHNkRk又h(n)为实数，则1()1( )()( )11kkN kNNH kH NkHzrWzrWz*1*1( )( )11()kkNNH kHkrWzr Wz101122212 cos()kkzzrkr zN11,2,.,21,2,.,12NkNNkN为奇数为偶数02Re( )k
11、H k其中：12 Re( )kkNrH k W 将第k个和第(N-k)个谐振器合并成一个实系数的二阶网络：当N为偶数时，还有一对实数根01(0)( )1HHzrz/21(/2)( )1NH NHzrzzr k=0, N / 2处：/2 10/211( )1( )( )( )NNNNkkH zr zHzHzHzNN为奇数时(1)/2011( )1( )( )NNNkkH zr zHzHzN只有一个实数根在 k = 0处：z = r4、快速卷积结构5、线性相位FIR滤波器的结构01nNFIR滤波器单位抽样响应h(n)为实数，且满足：( )(1)h nh Nn 偶对称：( )(1)h nh Nn 或奇对称：即对称中心在 (N-1) / 2处则这种FIR滤波器具有严格线性相位。N为奇数时10( )( )NnnH zh n z11112210121( )( )2NNNnnNnnNh n zhzh n z1112(1)201( )2NNnNnnNh nzzhz 1nNm 令h(n)偶对称，取“+”102Nhh(n)奇对称，取“ ”，且N为偶数时10( )( )NnnH zh n z11202( )( )NNnnNnnh n zh n z12(1)0( )NnNnnh nzz
