5线性相位FIR数字滤波器

1、数字信号处理数字信号处理by Zaiyue Y ZJU, 2016第第4章章 数字滤波器的结构数字滤波器的结构4.1 引言引言 4.2 用信号流图表示网络结构用信号流图表示网络结构4.3 IIR系统的基本网络结构系统的基本网络结构4.4 FIR系统的基本网络结构系统的基本网络结构4.5 FIR系统的线性相位结构系统的线性相位结构4.6 FIR系统的频率采样结构系统的频率采样结构4.7 数字信号处理中的量化效应数字信号处理中的量化效应2数字滤波器的设计与实现 （1）确定性能指标 （2）求系统函数H(z) （3）确定运算结构 （4）确定实现方法 已知寻求本章内容4.1 引言引言 关键点：同一个H(

2、z)可以写成不同形式，因此可以由不同结构来实现。3 一般时域离散系统或网络可以用差分方程、单位脉冲响应以及系统函数进行描述。如果系统输入输出服从N阶差分方程 0101( )()()( )( )( )1MNiiiiMiiiNiiiy nb x nia y nibzY zH zX za z其系统函数H(z)为 给定一个差分方程，不同的算法有很多种，例如： 1122113111( )10.80.151.52.5( )10.310.511( )10.310.5H zzzHzzzHzzz不同算法直接影响系统运算误差、运算速度以及系统复杂程度和成本5加法器乘法器单位延时基本运算单元方框图流图基本运算单元的

3、方框图及流图表示4.2 用信号流图表示网络结构用信号流图表示网络结构6流图结构 节点 源节点 支路 输出节点 网络节点 分支节点 输入支路 相加器节点的值=所有输入支路的值之和 输出支路支路的值=支路起点处的节点值传输系数7流图的化简 （1）并联支路 abxyyaxbxbaxy()ab xabxyybubaxyuabxyuaxcybcab1xycuxbcaby1（2）串联支路 ()b ax()ab xybu()b axcy（3）反馈支路 8例： 122221221211202( )(1)( )(1)( )( )( )( )( )( )( )nnnnnx nanany nbnbnbn(4.2.1

4、) 图4.2.2 信号流图(a)基本信号流图；(b)非基本信号流图可得120121212( )( )( )1bb zb zY zH zX za za z9基本信号流图 (1) 信号流图中所有支路都是基本的，即支路增益是常数或者是z-1； (2) 流图环路中必须存在延时支路； (3) 节点和支路的数目是有限的。 10 FIR：无反馈支路 差分方程，0( )()Miiy nb x ni单位脉冲响应h(n)有限长，,0( )0, other nbnMh nnFIR网络 v.s. IIR网络IIR：有反馈支路 差分方程，例如： 单位脉冲响应h(n)无限长，例如：( )(1)( )y nay nx n(

5、 )( )nh na u n114.3 IIR系统的基本网络结构系统的基本网络结构 1.直接型 N阶差分方程： 系统函数：01( )()()MNiiiiy nb x nia y niIIR的三种结构：直接型、级联型、并联型01201200( )( )( )11( ),( )1MiiiNiiiMiiNiiiib zH zH z Hza zH zb zHza z12 图4.3.1 IIR网络直接型结构 b0b1b2z1z1z1z1a1a2x(n)x(n 1)x(n 2)y(n)y(n 1)y(n 2)x(n)y(n)b0b1b2z1z1z1z1a1a2w2w1H1(z)H2(z)H2(z)H1(z

6、)x(n)y(n)a1a2b0b1b2z1z1（ a ）（ b ）（ c ）1!1!01( )()()MiiNiiy nb x nia y ni13 例4.3.1 IIR数字滤波器的系统函数H(z)为12312384112( )5311448zzzH zzzz画出该滤波器的直接型结构。 解：由H(z)写出差分方程531( )(1)(2)(3)8 ( )4 (1)44811 (2)2 (3)y ny ny ny nx nx nx nx n14图4.3.2 例4.3.1图x(n)y(n)z1z1z1 4811 245438115直接型特点 （1）简单直观, 运算速度快, 要求的内存少； （2）不能

7、直接调整滤波器系统函数的零、极点；（3）系数的有限字长效应对零、极点位置的影响很大，甚至可能使原设计稳定的滤波器变为不稳定的。 直接型结构多用于低阶（23阶）滤波器。 16 2. 级联型 将H(z)的分子、分母多项式分别因式分解1111(1)( )(1)MrrNrrC zH zAd z(4.3.1) Cr、dr为零、极点。由于它们是实数或共轭成对复数，因此上式可写作：1201212112( )( ),( )1kjjjjjjjjzzH zHzHza za z(4.3.2) 其中，0j、1j、2j、1j和2j均为实数。 17 图4.3.3 一阶和二阶直接型网络结构(a)直接型一阶网络结构；(b)直

8、接型二阶网络结构 x(n)y(n)z1x(n)y(n)z1z1( a )( b )j0j1j2j0j1j2j1j1j0 Hj(z)表示一个一阶或二阶的数字网络的系统函数，可由直接型网络结构表示： 结论：Hj(z)网络级联构成H(z) 网络。18 例4.3.2 设系统函数H(z)如下式： 12312384112( )1 1.250.750.125zzzH zzzz试画出其级联型网络结构。 解：将H(z)分子分母进行因式分解，得到 112112(20.379) (4 1.245.264)( )(1 0.25)(10.5)zzzH zzzzx(n)z12y(n)z14z10.3790.251.245

9、.2640.5 图4.3.4 例4.3.2图 19级联型特点 （1）每个一阶网络决定一个零点、一个极点，每个二阶网络决定一对零点、一对极点；（2）能直接调整滤波器系统函数的零、极点；（3）信号不会回流，运算误差的积累比直接型小；20 Hi(z) 为一阶或二阶网络，1( )( )kiiH zH z(4.3.4) 1011212( )1iiiiizH za za z 0i、1i、1i和2i为实数。 3.并联型将H(z)展成部分分式形式结论：Hi(z)网络并联构成H(z) 网络。21例4.3.3 画出例题4.3.2中的H(z)的并联型结构。 解：将H(z)展成部分分式形式：111281620( )1

10、610.510.5zH zzzz 将每部分用直接型结构实现，然后并联。 x(n )y(n )z 1z 11680.520 16 0.520z 1 图4.3.5 例4.3.3图 22并联型特点：（1）可以直接控制极点； （2）各二阶节的误差互不影响，故误差一般比级联型稍小；（3）有限字长效应的影响小； （4）零点不能独立地调节（二阶节的零点并不一定是系统的零点）； （5）系数较多 乘法次数多。234.4 FIR系统的基本网络结构系统的基本网络结构 FIR网络结构特点是没有反馈支路，即没有环路，其单位脉冲响应是有限长的。设单位脉冲响应h(n)长度为N，其系统函数H(z)和差分方程为1010( )(

11、 )( )( ) ()NnnNmH zh n zy nh m x nm24 1.直接型 按照H(z)或者差分方程直接画出结构图如图4.4.1所示。这种结构称为直接型网络结构或者称为卷积型结构。 图4.4.1 FIR直接型网络结构 x(n)y(n)z1z1z1h(0)h(1)h(2)h(N2)h(N1)25 2. 级联型 将H(z)进行因式分解，并将共轭成对的零点放在一起，形成一个系数为实数的二阶形式，这样级联型网络结构就是由一阶或二阶因子构成的级联结构，其中每一个因式都用直接型实现。 例4.4.1 设FIR网络系统函数H(z)如下式： H(z)=0.96+2.0z-1+2.8z-2+1.5z-

12、3 画出H(z)的直接型结构和级联型结构。 26 解：将H(z)进行因式分解，得到： H(z)=(0.6+0.5z-1)(1.6+2z-1+3z-2) 其直接型结构和级联型结构如图所示。 图4.4.2 例4.4.1图z1z1z1x(n)0.60.51.623y(n)y(n)x(n)z1z1z10.9622.81.5( a )( b )特点比较：（1）级联型的每个基本节控制一对零点，便于控制滤波器的传输零点（2）系数比直接型多，所需的乘法运算多274.5 线性相位线性相位FIR数字滤波器数字滤波器 考虑长度为N的h(n)，系统函数为：1()0()( )( )Njj njgnH eh n eHe

13、什么是线性相位FIR？10( )( )NnnH zh n z 频率响应函数为：jzeHg()称为幅度特性，()称为相位特性。注意，Hg()不同于|H(ej)|，Hg()为的实函数，可能取负值，而|H(ej)|总是正值。(4.5.1)(4.5.2)28 H(ej)线性相位是指： ()是的线性函数，即 为常数 (4.5.3) 或()满足下式： ,0是起始相位 (4.5.4) 严格地说， (4.5.4) 中()不具有线性相位，但以上两种情况都满足群时延是一个常数，即( )dd 0 第一类线性相位第二类线性相位2930第一类线性相位： h(n)是实序列且对(N-1)/2偶对称，即线性相位条件：( )(

14、1)h nh Nn第二类线性相位： h(n)是实序列且对(N-1)/2奇对称，即( )(1)h nh Nn 注意：充分条件(4.5.5)(4.5.6)31第一类线性相位条件证明：1010( )( )( )(1)NnnNnnH zh n zH zh Nnz令m=N-n-111(1)(1)00(1)1( )( )( )( )()NNN mNmmmNH zh m zzh m zH zzH z( )(1)h nh Nn321(1)1(1)01111()222011( )( )()( )221( ) 2NNnNnnNNNNnnnH zH zzH zh n zzzzh nzz z=ej11()20101(

15、)( )cos() 211( )( )cos() ,( )(1)22NNjjnNgnNH eeh nnNHh nnN 010,(1)2N33110011(1)(1)00(1)1( )( )(1)( )( )( )( )()NNnnnnNNN mNmnnNH zh n zh NnzH zh m zzh m zH zzH z 令m=N-n-11(1)1(1)01111222011( )( )()( )221( )2NNnNnnNNNNnnnH zH zzH zh n zzzzh nzz 第二类线性相位条件证明：34112011220()( )1( )sin ()21( )sin ()2jjz eN

16、NjnNNjjnH eH zNjeh nnNeh nn 1011( )( )sin (),( )()222NgnNNHh nn 01,(1)22N 35幅度特性Hg()的特点Case 1：第一类线性相位、N为奇数101( )( )cos() 2NgnNHh nn(3)/2011( )()2 ( )cos() 22NgnNNHhh nnh(n)对(N-1)/2偶对称，余弦项也对(N-1)/2偶对称，可以以(N-1)/2为中心，把两两相等的项进行合并，由于N是奇数，故余下中间项n=(N-1)/236cos(n-)对=0,2皆为偶对称因此Hg() 也对=0,2是偶对称的。可以实现各种（低通、高通、带

17、通、带阻）滤波器10( )( )2 ( )cos() MgnHhh nn 11,22NNM 37幅度特性Hg()的特点Case 2：第一类线性相位、N为偶数101( )( )cos() 2NgnNHh nn1201( )2 ( )cos ()2NgnNHh nn与N=奇数相似，不同点是由于N=偶数，Hg()中没有单独项，相等的项合并成N/2项。0( )2 ( )cos() MgnHh nn 11,22NNM 381cos ()cos ()2cos ()sin ()0222NnnNNnn cos(n-)对=为奇对称对=0, 2皆为偶对称因此Hg()=0， Hg()关于=是奇对称，关于=0, 2偶

18、对称可以实现低通和带通不能实现高通和带阻滤波器因为N是偶数，所以当=时有：3940幅度特性Hg()的特点Case 3：第二类线性相位、N为奇数101( )( )sin() 2NgnNHh nn3201( )2 ( )sin()2NgnNHh nnh(n)奇对称，因此10( )2 ( )sin() MgnHh nn 11,22NNM 1()02Nh41当=0,2 时，sin(n-)=0且sin(n-)对过零点奇对称因此，Hg()关于=0,2是奇对称只能实现带通滤波器不能实现低通、高通和带阻因为N是奇数，所以=(N-1)/2是整数。42幅度特性Hg()的特点Case 4：第二类线性相位、N为偶数1

19、001( )( )sin() 2 ( )sin() 2NMgnnNHh nnh nn 当=0, 2 时，sin(n-)=0；当= 时，sin(n-)=1，为峰值点sin(n-)对过零点奇对称，对峰值点偶对称因此，Hg()关于=0, 2是奇对称，关于=偶对称可以实现高通和带通不能实现低通和带阻因为N是偶数，所以=(N-1)/2=N/2-1/2。4344在第一类和第二类线性相位系统的证明中用到：(1)1( )()NH zzH z 线性相位FIR滤波器零点分布特点第一类取+第二类取如果zi为H(z)的零点： (zi)-1也是零点由于h(n)为实序列，零点共轭成对： zi*和(zi*)-1也是零点45

20、 图4.5.1 线性相位FIR滤波器零点分布 线性相位FIR滤波器零点分布特点46回顾线性相位FIR滤波器网络结构：N为偶数：12(1)0(1) 121(1)20( )( )( )( )(1)2NnN nnNNnN nnH zh n zzNH zh n zzhz N为奇数，则将中间项h(N-1)/2单独列出: 第一类取+第二类取 47图4.5.2 第一类线性相位网络结构x(n)y(n)z1z1z1z1z1z1z1h(0)h(1)h(2)h(N/21)x(n)y(n)z1z1z1z1z1z1h(0)h(1)h(2)h(N1)/2)N 偶数N 奇数48图4.5.3 第二类线性相位网络结构x(n)y

21、(n)z1z1z1z1z1z1z1h(0)h(1)h(2)h(N/21)x(n)y(n)z1z1z1z1z1z1h(0)h(1)h(2)h(N1)/2)N 偶数N 奇数111111111 由DFT可知，H(z)与频域采样值H(k)满足 1101( )( )(1)1NNkkNH kH zzNWz2( )( ),0,1,2,1jkNz eH kH zkN4.6 FIR系统的频率采样结构系统的频率采样结构条件：满足频率域采样定理，即频率域采样点数N大于等于原序列的长度M推论：M有限，因此频率采样结构只使用于FIR，不适用于IIR（4.6.1）2jNNWe49 将(4.6.1)式写成下式： 1011(

22、 )( )( )( )1( )( )1NckkNckkNH zHzHzNHzzH kHzWz (4.6.2) 式中 Hc(z)是一个梳状滤波网络，其零点为2,0,1,2,1jkkNkNzeWkNHc(z) 零点与Hk(z)极点对消50图4.6.1 FIR滤波器频率采样结构 x(n)y(n)z1z1 z NH(0)H(1)H(N 1)0NW1NW1NNWz1N151优点：(1)便于调整：在频率采样点k， ，只要调整H(k)(即一阶网络Hk(z)中乘法器的系数H(k)，就可以有效地调整频响特性(2)便于标准化、模块化：只要h(n)长度N相同，其梳状滤波器部分和N个一阶网络部分结构完全相同，只是各支

23、路增益H(k)不同缺点：(1)系统稳定性脆弱：位于单位圆上的N个零极点对消 (2)硬件实现不方便：H(k)和W-kN一般为复数，要求乘法器完成复数乘法运算 kjH eH k52修正：(1) 将单位圆上的零极点向单位圆内收缩到半径为r的圆上，取r1且r1。此时H(z)为1101( )( )(1)1NNNrkkNH kH zr zNrWz(4.6.3) 53(2) 将Hk(z)和HN-k(z)合并为一个二阶网络，记为Hk(z)1()111101122( )()( )11( )( )1121 2 cos()kkN kNNkkNNkkH kH NkHzrWzrWzH kHkrWzrW zaa zrk

24、zr zN由DFT的共轭对称性知道，如果h(n)是实数序列，则H(k)关于N/2点共轭对称，即H(k)=H*(N-k)；且WN-(N-k)=WNk=(WN-k)*012Re( )1,2,3,12 Re( )2kkkNaH kNkarH k W 其中，实系数为：1k0kz1z1 r 2)2cos(2kNrx(n)y(n)z1H(0)z N r r1/NH1(z)H2(z)z1 rH(N/2)( b )( a )(12zHN54当N为偶数时，H(z)可表示为11201111221()1(0)2( )(1)2111 2cos()NNNkkkNHaa zHH zr zNrzrzk zr zN 其中，H

25、(0)和H(N/2)为实数。频率采样修正结构由(N/2)-1个二阶网络和两个一阶网络并联构成1k0kz1z1 r 2)2cos(2kNrx(n)y(n)z1H(0)z N r r1/NH1(z)H2(z)z1 rH(N/2)( b )( a )(12zHN(4.6.4) 551(1)/201112211(0)( )(1)2112cos()NNNkkkHaa zH zr zNrzk zr zN当N为奇数时，H(0)为实数，H(z)可表示为(4.6.5) 564.7 数字信号处理中的量化效应数字信号处理中的量化效应 信号x(n)值量化后用Qx(n)表示， 量化误差用e(n)表示， e(n)=Qx(

26、n)-x(n)图 4.7.1 量化噪声e(n)的概率密度曲线 (a) 截尾法； (b) 舍入法 q0q1pe(n)e(n)(a)0q1e(n)(b)pe(n)12q12q57 1. A/D变换器中的量化效应 A/D变换器的功能原理图如图 4.7.2(a)所示， 图中 是量化编码后的输出， 如果未量化的二进制编码用x(n)表示， 那么量化噪声为e(n) = -x(n)， 因此A/D变换器的输出 为 x n ( )( )x nx ne n(4.7.1) 那么考虑A/D变换器的量化效应， 其方框图如图 4.7.2(b)所示。 这样， 由于e(n)的存在而降低了输出端的信噪比。 x n x n58 图

27、 4.7.2 A/DC功能原理图(a) A/DC变换器功能原理图； (b) 考虑量化效应的方框图采样量化编码x(n)xa(t)xa(nT)(a)理想A/DCxa(t)x(n)(b)e(n)x(n)59 假设A/D变换器输入信号xa(t)不含噪声， 输出 中仅考虑量化噪声e(n)，信号xa(t)平均功率用 表示， e(n)的平均功率用 表示， 输出信噪比用S/N表示， x n2x2e22xeSN或者用dB数表示 2210lgxeSdBN(4.7.2) A/D变换器采用定点舍入法， e(n)的统计平均值me=0， 方差2221121212beq60 将 代入(4.7.2)式， 得到:2e26.02

28、10.7910lgxSbN(4.7.3) 为充分利用其动态范围， 取 ，代入(4.7.3)式， 得13x6.021.29SbN61 2. 数字网络中系数的量化效应 数字网络或者数字滤波器的系统函数用下式表示：01( )1MrrrNrrrb zH za z 式中的系数br和ar必须用有限位二进制数进行量化， 存贮在有限长的寄存器中,经过量化后的系数用 和 表示，量化误差用br和ar表示， rrba,rrrrrraaa bbb62 对于N阶系统函数的N个系数ar，都会产生量化误差ar，每一个系数的量化误差都会影响第i个极点Pi的偏移。可以推导出第i个极点的偏移Pi服从下面公式：iiiPPP(4.7

29、.4) 11N rNiirNrilll iPPaPP(4.7.5) 63推导过程64 111111111111,1,iiiNNrrriirrNNNllllil il iriiirrrizpzpzpN rN rNiiiirNNrrilillll il iA zA za zPzzaA zzPzzzPPA zA zA zaPPPaaaA zPPPPPaaPPPP 上式表明极点偏移的大小与以下因素有关： (1) 极点偏移和系数量化误差大小有关。 (2) 极点偏移与系统极点的密集程度有关。 (3) 极点的偏移与滤波器的阶数N有关，阶数愈高， 系数量化效应的影响愈大， 因而极点偏移愈大。 系统的结构最好不

30、要用高阶的直接型结构，而将其分解成一阶或者二阶系统，再将它们进行并联或者串联，以减小极点偏移量。65例：设计一带通滤波器，并对其系数用16位字长量化，其中尾数15位。663. 数字网络中的运算量化效应 1) 运算量化效应 考虑定点乘法运算：21211222( )( ),( )( )( )( )( ),( )( )( )v nav nv nQ av ne nQ av nv nv nv ne n67图 4.7.3 考虑运算量化效应的一阶网络结构e1(n)e2(n)x(n)abz1y(n)ef(n)在图 4.7.3 中，有两个乘法支路，采用定点制时共引入两个噪声源，即e1(n)和e2(n)，噪声e2

31、(n)直接输出， 噪声e1(n)经过网络h(n)输出，输出噪声ef(n)为68 ef(n)=e1(n)*h(n)+e2(n) 如果尾数处理采用定点舍入法， 则输出端噪声平均值为121212 ( )( )( )( ) ()( )( )fmmmE e nh nE e nEh m e nmE e nmh mm 上式中E 表示求统计平均值， m1和m2分别表示两个噪声源的统计平均值， 这里m1=m2=0， 因此，0fm 69 由于e1(n)和e2(n)互不相关，求输出端噪声方差时， 可分别求其在输出端的方差，再相加。这里，每个噪声源的方差均为2222222121,212( ) ( )( )( )bef

32、fffffq qE enmE enE enE en输出端的噪声ef(n)的方差为70 式中，ef1(n)和ef2(n)分别表示e1(n)和e2(n)在输出端的输出； 222211100110020022022220( )( )( ) ()( ) ()( ) ( ) () ()( ) ( )()( )( )effmlmlemlemfeemE enE enEh m e nmh l e nlh m h l E e nm e nlh m h lmlhmhm 71 根据帕斯维尔定理， 也可以用下式计算：2212111( )()21( )1feedzH z H zjzbzH zaz72 2) 网络结构对输出噪声的影响 例 4.7.1 已知网络系统函数为1120.40.2( ),0.91 1.70.7zH zzzz 网络采用定点补码制， 尾数处理采用舍入法。 试分别计算直接型、 级联型和并联型结构输出噪声功率。 解 112111110