有限脉冲响应数字滤波器的设计ppt课件

1、第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计第第7章章 有限脉冲呼应数字滤波器的设计有限脉冲呼应数字滤波器的设计 7.1 线性相位线性相位FIR数字滤波器的条件和特点数字滤波器的条件和特点 7.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器7.3 利用频率采样法设计利用频率采样法设计FIR滤波器滤波器7.4 利用切比雪夫逼近法设计利用切比雪夫逼近法设计FIR滤波器滤波器7.5 IIR和和FIR数字滤波器的比较数字滤波器的比较第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计7.1 线性相位线性相位FIR数字滤波器的条件和特点数字滤波器的条件和特

2、点 FIR滤波器具有线性相位的特性，它的幅度有特殊的对称性，零点和网络构造也很特别。 1. 线性相位的条件 对于长度为N的h(n)，传输函数为)(10)()()(jgNnnjjeHenheH第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 式中，式中，Hg()称为幅度特性，称为幅度特性，()称为相位特性。称为相位特性。留意，这里留意，这里Hg()不同于不同于|H(ej)|, Hg()为为的实函数，的实函数，能够取负值，而能够取负值，而|H(ej)|总是正值。总是正值。H(ej)线性相位是线性相位是指指()是是的线性函数，即的线性函数，即 ()=-, 为常数为常数 (7.1.

3、3) 假设假设()满足下式也称为线性相位，满足下式也称为线性相位， ()=0-, 0是起始相位是起始相位 (7.1.4) 严厉地说，此时严厉地说，此时()不具有线性相位。但以上两种情不具有线性相位。但以上两种情况都满足群时延是一个常数，即况都满足群时延是一个常数，即( )dd 第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 所以两种情况都称为线性相位。普通称满足(7.1.3)式是第一类线性相位；满足(7.1.4)式为第二类线性相位。 FIR滤波器具有第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即 h(n)=h(N-n-1) FIR滤波器具有第二类线性

4、相位的条件是：h(n)是实序列且对(N-1)/2奇对称，即 h(n)=-h(N-n-1) 2. 线性相位FIR滤波器幅度特性Hg()的特点第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 3. 线性相位FIR滤波器零点分布特点 第一类和第二类线性相位的系统函数分别满足(7.1.7)式和(7.1.10)式，综合起来用下式表示：(1)1( )()NH zzH z 图7.1.1 线性相位FIR滤波器零点分布 第第7章章 有限脉冲响应数字滤波器的设

5、计有限脉冲响应数字滤波器的设计 4. 线性相位线性相位FIR滤波器网络构造滤波器网络构造 设设N为偶数，那么有为偶数，那么有11120021122(1)00( )( )( )( )( )( )(1)( )(1)NNNnnnNnmnNNnN mnmH zh n zh n zh n zH zh n zh Nmzh nh Nn 令令m=N-n-1,那么有那么有第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计12(1)0(1) 121(1)20( )( )( )( )(1)2NnN nnNNnN nnH zh n zzNH zh n zzhz 假设假设N为奇数，那么将中间项为奇数，

6、那么将中间项h(N-1)/2单独列出，单独列出， 从上面的公式可以看出，线性相位从上面的公式可以看出，线性相位FIR滤波器比滤波器比FIR滤波器的直接型构造节省乘法器近一半。线性相滤波器的直接型构造节省乘法器近一半。线性相位滤波器的网络构造图如下：位滤波器的网络构造图如下： 第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计图7.1.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 奇数第第7章章 有限脉冲响应数字滤波器的设计有

7、限脉冲响应数字滤波器的设计图7.1.3 第二类线性相位网络构造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 奇数111111111第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计7.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器 设希望设计的滤波器传输函数为设希望设计的滤波器传输函数为Hd(ej)，hd(n)是是它的单位脉冲呼应，因此它的单位脉冲呼应，因此 对于理想的低通滤波器，它的传送函数为对于理想的低通滤波器，它的传送函数为()( )1

8、( )()2jjddnjj nddHeh n eh nHeed第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 相应的单位取样呼应相应的单位取样呼应hd(n)为为,()0,j acjdceHe1sin()( )2()ccj aj ncdnah needna 为了构造一个长度为为了构造一个长度为N的线性相位滤波器，只需将的线性相位滤波器，只需将hd(n)截取一段，并保证截取的一段对截取一段，并保证截取的一段对(N-1)/2对称。设对称。设截取的一段用截取的一段用h(n)表示，即表示，即 h(n)=hd(n)RN(n)第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波

9、器的设计 我们实践实现的滤波器的单位取样呼应为我们实践实现的滤波器的单位取样呼应为h(n),长长度为度为N，其系统函数为，其系统函数为H(z)，10( )( )NnnH zh n z图图7.2.1 理想低通的单位脉冲呼应及矩形窗理想低通的单位脉冲呼应及矩形窗第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 以上就是用窗函数法设计以上就是用窗函数法设计FIR滤波器的思绪。另滤波器的思绪。另外，我们知道外，我们知道Hd(e j)是一个以是一个以2为周期的函数，可为周期的函数，可以展为傅氏级数，即以展为傅氏级数，即()( )jj nddnHeh n e根据复卷积定理，得到：根据复

10、卷积定理，得到：(1()()()2jjjdNH eHeRed 式中，式中，Hd(e j)和和RN(e j)分别是分别是hd(n)和和RN(n)的傅里叶变换，即的傅里叶变换，即111(1)200sin(/2()( )( )sin(/2)NNjNjj nj njaNNNnnNReRn eeeRe第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计sin(/21( ),sin(/2)2NNNRRN()称为矩形窗的幅度函数；将称为矩形窗的幅度函数；将Hd(ej)写成下式：写成下式：()( )jj addHeHe按照按照(7.2.1)式，理想低通滤波器的幅度特性式，理想低通滤波器的幅度特

11、性Hd()为为1,( )0,cdcH将将Hd(e j)和和RN(e j)代入代入(7.2.4)式，得到：式，得到：()1()( )()21( )()2jj ajadNj adNH eHeRedeHRd 第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 将将H(ej)写成下式：写成下式： 其卷积结果用图表示比较直观。其卷积结果用图表示比较直观。 ()( )1( )( )()2jj adNH eHeHHRd第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 图图7.2.2 矩形窗对理想低通矩形窗对理想低通 幅度特性的影响幅度特性的影响 第第7章章 有限脉冲响应

12、数字滤波器的设计有限脉冲响应数字滤波器的设计 经过以上分析可知，对经过以上分析可知，对hd(n)加矩形窗处置后，加矩形窗处置后，H()和原理想低通和原理想低通Hd()差别有以下两点：差别有以下两点： (1)在理想特性不延续点在理想特性不延续点=c附近附近H()有过渡带。有过渡带。过渡带的宽度，近似等于过渡带的宽度，近似等于RN()主瓣宽度，即主瓣宽度，即4/N。 (2)通带内通带内H()有动摇，最大的峰值在有动摇，最大的峰值在c-2/N处。处。阻带内阻带内H()有余振，最大的负峰在有余振，最大的负峰在c+2/N处。处。 第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 下面

13、引见几种常用的窗函数。设下面引见几种常用的窗函数。设 h(n)=hd(n)w(n) 式中式中w(n)表示窗函数。表示窗函数。 1. 矩形窗矩形窗(Rectangle Window) 2. 三角形窗三角形窗(Bartlett Window) 3. 汉宁汉宁(Hanning)窗窗升余弦窗升余弦窗 4. 哈明哈明(Hamming)窗窗改良的升余弦窗改良的升余弦窗 5. 布莱克曼布莱克曼(Blackman)窗窗第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计图图7.2.4 常用的窗函数常用的窗函数第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 图图7.2.5

14、常用窗函数的幅度特性常用窗函数的幅度特性 (a)矩形窗；矩形窗；(b)巴特利巴特利特窗特窗(三角形窗三角形窗)；(c)汉宁窗；汉宁窗；(d)哈明窗；哈明窗；(e)布莱克曼窗布莱克曼窗 第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 图图7.2.6 理想低通加窗后的幅度特性理想低通加窗后的幅度特性(N=51,c=0.5) (a)矩形窗；矩形窗；(b)巴特利特窗巴特利特窗(三角形窗三角形窗)；(c)汉宁窗；汉宁窗； (d)哈明窗；哈明窗；(e)布莱克曼窗布莱克曼窗第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 下面引见用窗函数设计下面引见用窗函数设计FI

15、R滤波器的步骤。滤波器的步骤。 (1)根据技术要求确定待求滤波器的单位取样呼应根据技术要求确定待求滤波器的单位取样呼应hd(n)。假设给出待求滤波器的频响为假设给出待求滤波器的频响为Hd(ej)，那么单位取样，那么单位取样响运用下式求出：响运用下式求出：1( )()2jjddhnHee d22101( )()MjkjknMMMdkhnHeeM第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 (2)根据对过渡带及阻带衰减的要求，选择窗函数的方根据对过渡带及阻带衰减的要求，选择窗函数的方式，并估计窗口长度式，并估计窗口长度N。设待求滤波器的过渡带用。设待求滤波器的过渡带用表示，

16、它近似等于窗函数主瓣宽度。表示，它近似等于窗函数主瓣宽度。 (3) 计算滤波器的单位取样呼应计算滤波器的单位取样呼应h(n)， h(n)=hd(n)w(n) (4)验算技术目的能否满足要求。设计出的滤波器频率验算技术目的能否满足要求。设计出的滤波器频率响运用下式计算：响运用下式计算：10()( )Njj nnH eh n e第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 例例7.2.1 用矩形窗、汉宁窗和布莱克曼窗设计用矩形窗、汉宁窗和布莱克曼窗设计FIR低通滤低通滤波器，设波器，设N=11,c=0.2rad。 解：解： 用理想低通作为逼近滤波器，按照用理想低通作为逼近滤

17、波器，按照(7.2.2)式，有式，有sin()( ),010()1(1)52sin(0.2 (5)( ),010(5)cddnh nnnNnh nnn第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 用汉宁窗设计：用汉宁窗设计：( )( )( ),0102( )0.5(1cos)10dHnHnh nh nnnnn用布莱克曼窗设计：用布莱克曼窗设计：它们的幅频特性如图它们的幅频特性如图 所示。所示。)()104cos08. 0102cos5 . 042. 0()()()()(11nRnnnnnhnhBlBld第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计

18、图图7.2.7 例例7.2.1的低通幅度特性的低通幅度特性第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计7.5 IIR和和FIR数字滤波器的比较数字滤波器的比较 首先，从性能上来说，首先，从性能上来说，IIR滤波器传输函数的极点滤波器传输函数的极点可位于单位圆内的任何地方，因此可用较低的阶数获可位于单位圆内的任何地方，因此可用较低的阶数获得高的选择性，所用的存贮单元少，所以经济而效率得高的选择性，所用的存贮单元少，所以经济而效率高。但是这个高效率是以相位的非线性为代价的。高。但是这个高效率是以相位的非线性为代价的。FIR滤波器可得到线性相位，阶数较高。滤波器可得到线性相位，阶数较高。 第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 从构造上看，从构造