数字信号处理_第七章_第1页
数字信号处理_第七章_第2页
数字信号处理_第七章_第3页
数字信号处理_第七章_第4页
数字信号处理_第七章_第5页
已阅读5页,还剩62页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 7.1 线性相位线性相位FIR数字滤波器的条件和特点数字滤波器的条件和特点 7.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器7.3 利用频率采样法设计利用频率采样法设计FIR滤波器滤波器7.4 利用等波纹最佳逼近法设计利用等波纹最佳逼近法设计FIR滤波器滤波器7.5 IIR和和FIR数字滤波器的比较数字滤波器的比较7.6 几种特殊类型滤波器简介几种特殊类型滤波器简介7.7 滤波器分析设计工具滤波器分析设计工具FDATool7.1 线性相位线性相位FIR数字滤波器的条件和特点数字滤波器的条件和特点 本节主要介绍本节主

2、要介绍FIR滤波器具有线性相位的条件、滤波器具有线性相位的条件、 FIR滤波器的幅度特性、零点分布特点和网络结构的特点。滤波器的幅度特性、零点分布特点和网络结构的特点。 1. 线性相位条件线性相位条件 对于长度为对于长度为N的的h(n),传输函数为:,传输函数为:式中,式中,Hg()称为幅度特性,称为幅度特性,()称为相位特性。这里称为相位特性。这里Hg()不同于不同于|H(ej)|,Hg()为为的实函数,可能取负值,而的实函数,可能取负值,而|H(ej)|为为的正实函数。的正实函数。(7.1.1) (7.1.2) H(ej)线性相位:线性相位: ()是是的线性函数的线性函数,即:,即: ()

3、= -, 为常数为常数 (7.1.3) 若若()满足下式:满足下式: ()=0 -, 0是起始相位是起始相位 (7.1.4) 也称这种情况为线性相位。也称这种情况为线性相位。 以上两种情况都满足以上两种情况都满足群时延是一个常数群时延是一个常数,即,即一般称:满足一般称:满足(7.1.3)式是式是第一类第一类线性相位;线性相位; 满足满足(7.1.4)式为式为第二类第二类线性相位。线性相位。 2. 线性相位线性相位FIR的时域约束条件的时域约束条件线性相位线性相位FIR滤波器的时域约束条件是指满足线性滤波器的时域约束条件是指满足线性相位时,对相位时,对h(n)的约束条件。的约束条件。1) 第一

4、类线性相位对第一类线性相位对h(n)的约束条件的约束条件 相位函数相位函数()=,由式(,由式(7.1.1)和()和(7.1.2)得)得到到:1jjjg0(e)( )e( )e NnnHh nH(7.1.5) 1g0( )(cosjsin)( )(cosjsin) Nnh nnnH(7.1.6)将(将(7.1.6)式中两式相除得到:)式中两式相除得到:1010( )coscossin( )sinNnNnh nnh nn1g01g0( )cos( )cos ( )sin( )sinNnNnHh nnHh nn即即 移项并用三角公式化简得到移项并用三角公式化简得到:(7.1.7)函数函数h(n)s

5、in(n)关于求和区间的中心关于求和区间的中心(N1)/2奇对称,奇对称,是满足(是满足(7.1.7)式的)式的一组解一组解。 因为因为sin(n)关于关于n=奇对称,如果取奇对称,如果取=(N1)/2,则要,则要求求h(n)关于关于(N1)/2偶对称,所以要求偶对称,所以要求和和h(n)满足如下条件满足如下条件: 10( )sin ()0Nnh nn1100( )cossin( )sincosNNnnh nnh nn(7.1.8) 1( ) , 2( )(1), 01Nh nh NnnN 表表7.1.1 线性相位线性相位FIR数字滤波器的时域和频域特性一览数字滤波器的时域和频域特性一览 2)

6、 第二类线性相位对第二类线性相位对h(n)的约束条件的约束条件相位函数相位函数()=/2,由式(,由式(7.1.1)和()和(7.1.2),),可得到可得到:(7.1.9) 函数函数h(n)cos(n)关于关于求和区间的中心求和区间的中心(N1)/2奇奇对称对称,是满足式(,是满足式(7.1.9)的)的一组解一组解, 因为因为cos(n)关于关于n=偶对称,所以要求偶对称,所以要求和和h(n)满满足如下条件:足如下条件:1jjj( /2)g0(e)( )e( )eNnnHh nH10( )cos ()0Nnh nn(7.1.10) 1( ) , 22( )(1), 01Nh nh NnnN 2

7、 线性相位线性相位FIR滤波器幅度特性滤波器幅度特性Hg()的特点的特点实质上实质上,幅度特性的,幅度特性的特点特点就是线性相位就是线性相位FIR滤波滤波器的器的频域约束条件频域约束条件。 引入两个参数符号:引入两个参数符号: 1,2N12NM情况情况1: h(n)=h(Nn1), N为奇数。为奇数。将时域约束条件将时域约束条件h(n)=h(Nn1)和和()=代入代入式(式(7.1.1)和()和(7.1.2),得到),得到:所以所以(7.1.11) 因为因为cos(n-)关于关于=0, , 2三点三点偶对称偶对称,所以由,所以由式(式(7.1.11)可以看出,)可以看出,Hg()关于关于=0,

8、 , 2三点三点偶对称偶对称。因此情况因此情况1可以实现各种(低通、高通、带通、带阻)滤可以实现各种(低通、高通、带通、带阻)滤波器。波器。1g0( )( )2 ( )cos ()MnHhh nn情况情况2: h(n)=h(Nn1), N为偶数。为偶数。仿照情况仿照情况1的推导方法得到的推导方法得到:1jjjjg00(e)( )e =( )ee2 ( )cos( ()NMnnnHHh nh nn(7.1.12)g0( )2 ( )cos ()MnHh nn式中,。因为式中,。因为是偶数,所是偶数,所以当时以当时(1)/2/21/2NN cos ()cossin022NNnnn 而且而且cos(

9、n)关于关于过零点奇对称过零点奇对称,关于,关于=0和和2偶对称。所以偶对称。所以Hg()=0,Hg()关于关于=奇对奇对称,关于称,关于=0和和2偶对称偶对称。因此,情况。因此,情况2不能实现高通不能实现高通和带阻滤波器。和带阻滤波器。情况情况3: h(n)=h(Nn1),N为奇数。为奇数。102Nh 将时域约束条件将时域约束条件h(n)=h(Nn1)和和()=/2代入式(代入式(7.1.1)和()和(7.1.2),), 并考虑并考虑,得到,得到:1g0( )2 ( )sin ()MnHh nn当当=0,, 2时,时,sin(n)=0, 且且sin(n)关于过零点奇对称关于过零点奇对称Hg(

10、)关于关于=0, , 2三三点奇对称点奇对称情况情况3只能实只能实现带通滤波器现带通滤波器情况情况4: h(n)=h(Nn1), N为偶数。为偶数。用情况用情况3的推导过程可以得到的推导过程可以得到:(7.1.13) N是偶数,是偶数,=(N1)/2=N/21/2。所以,当。所以,当=0, 2时,时,sin(n)=0;当;当=时,时,sin(n)=(1)nN/2, 为峰值点。而且为峰值点。而且sin(n)关于过零点关于过零点=0和和2两点奇对称,关于峰值点两点奇对称,关于峰值点=偶对称。因此偶对称。因此Hg()关于关于=0和和2两点奇对称,关于两点奇对称,关于=偶对称。由此可见,情况偶对称。由

11、此可见,情况4不能实现低通和带阻滤波器。不能实现低通和带阻滤波器。g0( )2 ( )sin ()MnHh nn 3. 线性相位线性相位FIR滤波器零点分布特点滤波器零点分布特点线性相位的系统函数满足:线性相位的系统函数满足:(7.1.21) “+”和和“-”分别对应第一类和第二类线性相位。分别对应第一类和第二类线性相位。一般情况:一般情况: 线性相位线性相位FIR零点零点分布的特点:分布的特点:互为互为倒数的共轭对倒数的共轭对图图7.1.1 线性相位线性相位FIR数字滤波器的零点分布数字滤波器的零点分布当然,也有当然,也有一些特殊情一些特殊情况,如图中况,如图中z1、z2和和z4情情况。况。

12、7.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器7.2.1 窗函数法设计原理窗函数法设计原理设计思想设计思想:从时域出发,设计:从时域出发,设计h(n),逼近,逼近hd(n) 。 设希望设计的滤波器传输函数为设希望设计的滤波器传输函数为Hd(ej),hd(n)是与其对应是与其对应的单位脉冲响应,因此的单位脉冲响应,因此 :ccde )e (21)(e )()e (jjddjdjdnnnHnhnhH一旦一旦Hd(ej)给定,就可得到给定,就可得到hd(n) 。一般情况下,一般情况下, Hd(ej)逐段恒定,在边界频率处有不连续点,因而逐段恒定,在边界频率处有不连续点,因而hd(n)是是

13、无限时宽无限时宽的,且的,且是非因果是非因果的,为得到一有限长滤波器的,为得到一有限长滤波器h(n),最直接,最直接的方法是截断的方法是截断 hd(n),即用一个窗口函数,即用一个窗口函数RN(n)对对hd(n)进行加窗处理:进行加窗处理: h(n)= hd(n) RN(n)。并保证截取的一段关于并保证截取的一段关于n=(N1)/2偶对称。偶对称。 相应的单位取样响应相应的单位取样响应hd(n)为:为:(7.2.1)它是一个无限长的非因果序列。波形如下图所示。它是一个无限长的非因果序列。波形如下图所示。设理想低通滤波器的传输函数设理想低通滤波器的传输函数Hd(ej) 为:为:(7.2.2) 实

14、际设计的滤波器的单位脉冲响应为实际设计的滤波器的单位脉冲响应为h(n),长度为,长度为N,其系统函数为其系统函数为H(z),10)()(NnnznhzH用一个有限长的序用一个有限长的序列列h(n)去代替去代替hd(n),肯定会引起误差肯定会引起误差图图7.2.1 理想低通的单位脉冲响应及矩形窗理想低通的单位脉冲响应及矩形窗)()()(dnRnhnhN吉布斯(吉布斯(Gibbs)效应:)效应:用一个有限长的序列用一个有限长的序列h(n)去代替去代替hd(n) 会引起误差,表现在频域就是会引起误差,表现在频域就是吉布斯效应吉布斯效应。该效应引。该效应引起通带内和阻带内的波动性,尤其使阻带的衰减不足

15、,从起通带内和阻带内的波动性,尤其使阻带的衰减不足,从而不能满足技术上的要求。而不能满足技术上的要求。 吉布斯效应是由于将吉布斯效应是由于将hd(n)直接截断引起的,因此,也直接截断引起的,因此,也称为称为截断效应截断效应。图图7.2.2 吉普斯效应吉普斯效应Hd(ej)是一个以是一个以2为周期的函数,可以展为傅里叶级数,为周期的函数,可以展为傅里叶级数,即即傅里叶级数的系数为傅里叶级数的系数为hd(n),当然就是,当然就是Hd(ej)对应的对应的单位脉冲响应。设计单位脉冲响应。设计FIR滤波器就是根据要求找到滤波器就是根据要求找到N个傅个傅里叶级数系数里叶级数系数h(n),n=1, 2, ,

16、 N1,以,以N项傅氏级数去项傅氏级数去近似代替无限项傅氏级数,这样在一些频率不连续点附近似代替无限项傅氏级数,这样在一些频率不连续点附近会引起较大误差,这种误差就是前面说的近会引起较大误差,这种误差就是前面说的截断效应。截断效应。nnnhHjdjde )()e ( 因此,从这一角度来说,窗函数法也称为因此,从这一角度来说,窗函数法也称为傅氏级数法傅氏级数法。显然,选取傅氏级数的项数愈多,引起的误差就愈小,但显然,选取傅氏级数的项数愈多,引起的误差就愈小,但项数增多即项数增多即h(n)长度增加,也使成本和滤波计算量加大,长度增加,也使成本和滤波计算量加大,应在满足技术要求的条件下,尽量减小应在

17、满足技术要求的条件下,尽量减小h(n)的长度。的长度。)( jRjdjd)e ()e (21)e (WHH(7.2.4)根据傅里叶变换的时域卷积定理,得到(根据傅里叶变换的时域卷积定理,得到(7.2.3)式的傅里叶变换:式的傅里叶变换:(7.2.5)jRg) 1(21j10j10jR10jRe )()2/sin()2/sin(e eee )()(WNwnWeWNNnnNnnNnnjR第一过零点第一过零点内为主瓣内为主瓣旁瓣旁瓣将将Hd(ej)写成写成Hd(ej)=Hdg()ej , 则按照(则按照(7.2.1)式,)式,理想低通滤波器的幅度特性函数为理想低通滤波器的幅度特性函数为|0|1)(c

18、cdg,H将将Hd(ej)和和WR(ej)代入(代入(7.2.4)式,得到:)式,得到:Rgdgj)( jRgjdgjd)()(21ede )(e )(21)e (WHWHH将将H(ej)写成写成H(ej)=Hg()ej ,则,则(7.2.6)gdgRg1( )( )()d2HHW加窗后的滤波器的幅度特性等于理想低通滤波器的幅度特性加窗后的滤波器的幅度特性等于理想低通滤波器的幅度特性Hdg()与与矩形窗幅度特性矩形窗幅度特性WRg()的卷积。的卷积。图图7.2.3 矩形窗加窗效应矩形窗加窗效应当当=c时,如时,如(c)所示,当所示,当c 2/N时,积分近似为时,积分近似为WRg()一半波一半波

19、形的积分,对形的积分,对Hg(0)归一化后的值近归一化后的值近似为似为1/2当当=0时,时,Hg(0)等于图等于图(a)与与(b)两波两波形乘积的积分,相当于对形乘积的积分,相当于对WRg()在在c之间一段波形的积分,当之间一段波形的积分,当c2/N时,近似为时,近似为之间波形之间波形的积分的积分当当=c2/N时,如时,如(d)所示,所示,WR()主瓣完全在区间主瓣完全在区间c, c之之内,而最大的一个负旁瓣移到区间内,而最大的一个负旁瓣移到区间c, c之外,因此之外,因此Hg(c2/N)有一个最大的正峰。有一个最大的正峰。当当=c+2/N时,如时,如(e)所示,所示,WRg()主瓣完全移到积

20、分区间外边,主瓣完全移到积分区间外边,由于最大的一个负旁瓣完全在区间由于最大的一个负旁瓣完全在区间c, c内,因此内,因此Hg(c+2/N)形成最大的负峰。形成最大的负峰。Hg()最大的正峰与最大的负峰对最大的正峰与最大的负峰对应的频率相距应的频率相距4/N对对hd(n)加矩形窗处理后,加矩形窗处理后,Hg()与原理想低通与原理想低通Hdg()的差别有以下两点:的差别有以下两点:在理想特性不连续点在理想特性不连续点=c附近形成过渡带。过渡带的宽附近形成过渡带。过渡带的宽度近似等于度近似等于WRg()主瓣宽度主瓣宽度4/N。通带内产生了波纹,最大的峰值在通带内产生了波纹,最大的峰值在c2/N处。

21、阻带内处。阻带内产生了余振,最大的负峰在产生了余振,最大的负峰在c+2/N处。通带与阻带中波处。通带与阻带中波纹的情况与窗函数的幅度谱有关,纹的情况与窗函数的幅度谱有关, WRg()旁瓣幅度的大旁瓣幅度的大小直接影响小直接影响Hg()波纹幅度的大小。波纹幅度的大小。v 这两点就是对这两点就是对hd(n)截断后,在频域的反映,即截断后,在频域的反映,即吉布斯吉布斯效应效应。 如何减少吉布斯效应的影响,设计一个满足要求的如何减少吉布斯效应的影响,设计一个满足要求的FIR滤波器呢?滤波器呢? 直观上,增加矩形窗口的宽度(即加大直观上,增加矩形窗口的宽度(即加大N)可以减少吉布斯效应)可以减少吉布斯效

22、应的影响。的影响。N 时,时, 在主瓣附近,在主瓣附近, WRg()近似为:近似为: 该函数的性质是:该函数的性质是: 随随x (N ),主瓣幅度加高,同时旁瓣也加高,保持主瓣,主瓣幅度加高,同时旁瓣也加高,保持主瓣和旁瓣幅度相对值不变;波动的频率加快,当和旁瓣幅度相对值不变;波动的频率加快,当x 时,时,sinx/x趋近于趋近于 函数,因此,当函数,因此,当N加大时,加大时, H()的波动幅度没有多大的波动幅度没有多大改善。改善。 N加大带来的最大好处是过渡带变窄加大带来的最大好处是过渡带变窄(过渡带过渡带:4/N)。 因此因此加大加大N并不是减少吉布斯效应的有效方法。并不是减少吉布斯效应的

23、有效方法。Rgsin(/2)sin( )/2NxWNx图图7.2.4 矩形窗函数长度的影响矩形窗函数长度的影响所以有如下结论:所以有如下结论:v 调整窗口长度调整窗口长度N可以有效地控制过渡带的宽度。可以有效地控制过渡带的宽度。v 减少带内波动以及加大阻带的衰减只能从窗函数的形状减少带内波动以及加大阻带的衰减只能从窗函数的形状上找解决方法。上找解决方法。为了消除吉布斯效应,取得较好的频率特性,一般采用其为了消除吉布斯效应,取得较好的频率特性,一般采用其它类型的窗口,它类型的窗口,这些窗口函数在靠近两端的这些窗口函数在靠近两端的h(n)逐步衰减到零,逐步衰减到零,使窗口频谱的主瓣包含更多的能量。

24、使窗口频谱的主瓣包含更多的能量。从而降低通带内的纹波并从而降低通带内的纹波并加大阻带的衰减。加大阻带的衰减。下面是几种常用的窗函数下面是几种常用的窗函数1. 矩形窗矩形窗(Rectangle Window) 2. 三角形窗三角形窗(Bartlett Window) 3. 汉宁汉宁(Hanning)窗窗-升余弦窗升余弦窗 4. 哈明哈明(Hamming)窗窗改进的升余弦窗改进的升余弦窗 5. 布莱克曼布莱克曼(Blackman)窗窗 6. 凯塞凯塞贝塞尔窗贝塞尔窗(Kaiser-Basel Window) 7.2.2. 几种常用的窗函数几种常用的窗函数(7.2.7) 为了描述方便,定义窗函数的几

25、个为了描述方便,定义窗函数的几个参数参数: 旁瓣峰值旁瓣峰值 n窗函数的幅频函数窗函数的幅频函数|Wg()|的最大旁瓣的最大旁瓣的最大值相对主瓣最大值的衰减值(的最大值相对主瓣最大值的衰减值(dB););过渡带宽度过渡带宽度Bg用该窗函数设计的用该窗函数设计的FIR数字滤波器数字滤波器(FIRDF)的过渡带宽度;)的过渡带宽度;阻带最小衰减阻带最小衰减 s用该窗函数设计的用该窗函数设计的FIRDF的阻带最的阻带最小衰减。小衰减。Rgsin(/2)( )sin(/2)NW1 矩形窗(矩形窗(Rectangle Window)wR(n)=RN(n)其幅度函数为其幅度函数为图图7.2.4 矩形窗函数

26、长度的影响矩形窗函数长度的影响2 三角形窗(三角形窗(Bartlett Window)(7.2.8) 其频谱函数为其频谱函数为1) 1(21122) 1(21012)(BNnNNnNnNnn,21j2jBe)2/sin(4/sin2)e (NNNW其幅度函数为其幅度函数为 (7.2.10)2Bg2sin(/4)( )sin(/2)NWN图7.2.5 三角窗的四种波形主瓣主瓣能量能量更集更集中,中,但过但过渡带渡带宽了宽了很多很多3 汉宁(汉宁(Hanning)窗)窗升余弦窗升余弦窗(7.2.11) HnN2( )0.5 1 cos( )1nwnRnN1jj2RNRgjHnHn1j2RgRgRg

27、1j2Hng(e )( )( )e(e )( )22 0.5( )0.25e11 ( )eNNNWFT RnWWFT WnWWWNNW当当N1时,时, N1NHngRgRgRg22( )0.5( )0.25WWWWNN汉宁窗的幅度函数汉宁窗的幅度函数WHng()由三部分相加,旁瓣互相对由三部分相加,旁瓣互相对消,使能量更集中在主瓣中。汉宁窗的四种波形如图消,使能量更集中在主瓣中。汉宁窗的四种波形如图7.2.6所示,参数为所示,参数为: n=31 dB; Bg=8/N; s=44 dB。图图7.2.6 汉宁窗的四种波形汉宁窗的四种波形4 哈明(哈明(Hamming)窗)窗改进的升余弦窗改进的升余

28、弦窗 (7.2.12)其频谱函数其频谱函数WHm(ej)为为)(12cos46. 054. 0)(NHmnRNnn)e (23. 0)e (23. 0)e (54. 0)e (12R12jRERjHmNENWWWW其幅度函数其幅度函数WHmg()为为当当N时,其可近似表示为时,其可近似表示为1223. 01223. 0)(54. 0)(RgRgRgHmgNWNWWWNWNWWW223. 0223. 0)(54. 0)(RgRgRgHmg图图7.2.7 哈明窗的四种波形哈明窗的四种波形这种改进这种改进的升余弦的升余弦窗,能量窗,能量更加集中更加集中在主瓣,在主瓣,但其主瓣但其主瓣宽度和汉宽度和汉

29、宁窗的相宁窗的相同同(7.2.13)5 布莱克曼(布莱克曼(Blackman)窗)窗其频谱函数为其频谱函数为)(14cos08. 012cos5 . 042. 0)(BlnRNnNnnN)()()(14jR14jR12jR12jRjRjBlee04. 0 e)e (25. 0)e (42. 0)e (NNNNWWWWWW其幅度函数为其幅度函数为141404. 0 121225. 0)(42. 0)(RgRgRgRgRgBlgNWNWNWNWWW幅度函数由五部分组成,它们都是移幅度函数由五部分组成,它们都是移位不同,且幅度也不同的位不同,且幅度也不同的WRg()函数,函数,使旁瓣再进一步抵消。旁

30、瓣峰值幅度使旁瓣再进一步抵消。旁瓣峰值幅度进一步增加,其幅度谱主瓣宽度是矩进一步增加,其幅度谱主瓣宽度是矩形窗的形窗的3倍倍图7.2.8 布莱克曼窗的四种波形旁瓣进旁瓣进一步减一步减小,但小,但过渡带过渡带变宽变宽6 凯塞凯塞贝塞尔窗(贝塞尔窗(Kaiser-Basel Window)以上五种窗函数都称为参数固定窗函数,每种窗以上五种窗函数都称为参数固定窗函数,每种窗函数的旁瓣幅度都是固定的。凯塞函数的旁瓣幅度都是固定的。凯塞贝塞尔窗是一种贝塞尔窗是一种参数可调的窗函数,是一种最优窗函数。参数可调的窗函数,是一种最优窗函数。 (7.2.15)式中式中10)()()(00kNnIIn,21121

31、NnI0()是零阶第一类修正贝塞尔函数,可用下面级数计算:是零阶第一类修正贝塞尔函数,可用下面级数计算:一般一般I0()取取1525项,便可以满足精度要求。项,便可以满足精度要求。 参数可以参数可以控制窗的形状。一般控制窗的形状。一般 加大,主瓣加宽,旁瓣幅度减小,加大,主瓣加宽,旁瓣幅度减小,典型数据为典型数据为4 9。当。当 =5.44时,窗函数接近哈明窗。时,窗函数接近哈明窗。 =7.865时,窗函数接近布莱克曼窗。在设计指标给定时,时,窗函数接近布莱克曼窗。在设计指标给定时,可以调整可以调整 值,使滤波器阶数最低,所以其性能最优。凯值,使滤波器阶数最低,所以其性能最优。凯塞(塞(Kai

32、ser)给出的估算)给出的估算和滤波器阶数和滤波器阶数N的公式如下的公式如下:2011( )1! 2kkIk (7.2.17) 式中,式中,Bt=|sp|, 是数字滤波器过渡带宽度。应当注是数字滤波器过渡带宽度。应当注意,因为式(意,因为式(7.2.17)为阶数估算,所以必须对设计结)为阶数估算,所以必须对设计结果进行检验。另外,凯塞窗函数没有独立控制通带波纹果进行检验。另外,凯塞窗函数没有独立控制通带波纹幅度,实际中通带波纹幅度近似等于阻带波纹幅度。凯幅度,实际中通带波纹幅度近似等于阻带波纹幅度。凯塞窗的幅度函数为塞窗的幅度函数为(7.2.16) ss0.4ssss0.112(8.7),50

33、 dB0.5842(21)0.07886(21),2150 dB0,21st82.285NB(1)/2kgkk1( )(0)2( )cosNnWww nn(7.2.18)对对 的的8种典型值,将凯塞窗函数的性能列于表种典型值,将凯塞窗函数的性能列于表7.2.1中,中,供设计者参考。由表可见供设计者参考。由表可见, 当当 =5.568时时, 各项指标都各项指标都好于哈明窗。好于哈明窗。表表7.2.1 凯塞窗参数对滤波器的性能影凯塞窗参数对滤波器的性能影 表表7.2.2 6种窗函数的基本参数种窗函数的基本参数 7.2.3 用窗函数法设计用窗函数法设计FIR滤波器的步骤滤波器的步骤(1) 根据对过渡

34、带及阻带衰减的指标要求,选择窗根据对过渡带及阻带衰减的指标要求,选择窗函数的类型,并估计窗口长度函数的类型,并估计窗口长度N。(2) 构造希望逼近的频率响应函数构造希望逼近的频率响应函数Hd(ej),即,即jj(1)/2ddg(e)( )eNHH一般取一般取2spc(3) 计算计算hd(n)。 如果给出待求滤波器的频响函数为如果给出待求滤波器的频响函数为Hd(ej),那,那么单位脉冲响应用下式求出:么单位脉冲响应用下式求出:(7.2.21)jjddde )e (21)(nHnh2jdd( )(e)kMMHkH如果如果Hd(ej)较复杂,或者不能用封闭公式表示,则不较复杂,或者不能用封闭公式表示

35、,则不能用上式求出能用上式求出hd(n)。可以对。可以对Hd(ej)从从=0到到=2采样采样M点,采样值为点,采样值为 ,k=0 ,1, 2, , M1,进行,进行M点点IDFT(IFFT),得到:,得到: (7.2.22)根据频域采样理论,根据频域采样理论,hdM(n)与与hd(n)应满足如下关系:应满足如下关系:dd( )IDFT( )MMMhnHkdd( )()( )MMrhnh nrM Rn因此,如果因此,如果M选得较大,可以保证在窗口内选得较大,可以保证在窗口内hdM(n)有效有效逼近逼近hd(n)。 对(对(7.2.19)式给出的线性相位理想低通滤波器作)式给出的线性相位理想低通滤

36、波器作为为Hd(ej),由(,由(7.2.2)式求出单位脉冲响应)式求出单位脉冲响应hd(n):为保证线性相位特性,为保证线性相位特性, =(N1)/2。(4) 加窗得到设计结果:加窗得到设计结果:h(n)=hd(n)w(n)。cdsin()( )()nh nn例例7.2.1 用窗函数法设计线性相位高通用窗函数法设计线性相位高通FIRDF,要求,要求通带截止频率通带截止频率p=/2 rad,阻带截止频率,阻带截止频率s=/4 rad,通带最大衰减通带最大衰减 p=1 dB,阻带最小衰减,阻带最小衰减 s=40 dB。解解 (1) 选择窗函数选择窗函数w(n),计算窗函数长度,计算窗函数长度N。

37、已知阻带最小衰减已知阻带最小衰减 s=40 dB,由表(,由表(7.2.2)可知汉宁)可知汉宁窗和哈明窗均满足要求,我们选择窗和哈明窗均满足要求,我们选择汉宁窗汉宁窗。 过渡带宽度过渡带宽度Btps=/4, 汉宁窗的精确过渡带汉宁窗的精确过渡带宽度宽度Bt=6.2/N,所以要求,所以要求Bt=6.2/N/4,解之得,解之得N24.8。对高通滤波器。对高通滤波器N必须取奇数,取必须取奇数,取N=25。25( )0.5 1cos( )12nw nRn(2) 构造构造Hd(ej):式中式中jjcdce,(e)0,0Hspc1312, 22N(3) 求出求出hd(n):jjdd1( )(e)ed2nh

38、 nHccjjjj1 eedeed2nncsin()sin () ()()nnnn将将=12代入得代入得 dsin3 (12)/8( )(12)(12)nh nnn(4) 加窗加窗:d( )( ) ( )h nh n w n25sin3 (12)/8(12)0.50.5cos( )(12)12nnnRnn全通滤全通滤波器波器理想低通滤理想低通滤波器的单位波器的单位脉冲响应脉冲响应例例7.2.2 对模拟信号进行低通滤波处理,要求通带对模拟信号进行低通滤波处理,要求通带0f 1.5kHz内衰减小于内衰减小于1 dB, 阻带阻带2.5kHzf 上衰减上衰减大于大于40 dB。希望对模拟信号采样后用线

39、性相位。希望对模拟信号采样后用线性相位FIR数字滤数字滤波器实现上述滤波,采样频率波器实现上述滤波,采样频率Fs=10 kHz。用窗函数法设。用窗函数法设计满足要求的计满足要求的FIR数字低通滤波器,求出数字低通滤波器,求出h(n),并画出损,并画出损耗函数曲线。为了降低运算量,希望滤波器阶数尽量低。耗函数曲线。为了降低运算量,希望滤波器阶数尽量低。解解 (1) 确定相应的数字滤波器指标确定相应的数字滤波器指标: 通带截止频率为通带截止频率为3 . 010000150022sppFf 阻带截止频率为阻带截止频率为 阻带最小衰减为阻带最小衰减为 s =40dB(2) 用窗函数法设计用窗函数法设计

40、FIR数字低通滤波器,为数字低通滤波器,为了降低阶数选择凯塞窗。根据式(了降低阶数选择凯塞窗。根据式(7.2.16)计算凯塞)计算凯塞窗的控制参数为窗的控制参数为5 . 010000250022sssFf0.4ss0.5842(21)0.07886(21)3.3953指标要求过渡带宽度指标要求过渡带宽度Bt=sp=0.2,根据式(,根据式(7.2.17)计算滤波器阶数为计算滤波器阶数为取满足要求的最小整数取满足要求的最小整数M=23。所以。所以h(n)长度为长度为N=M+1=24。但是,如果用汉宁窗,。但是,如果用汉宁窗,h(n)长度为长度为N=40。理想低通滤波器的通带截止频率理想低通滤波器

41、的通带截止频率c=(s+p)/2=0.4,所,所以由式(以由式(7.2.2)和式()和式(7.2.3), 得到得到:式中,式中,w(n)是长度为是长度为24( =3.395)的凯塞窗函数。)的凯塞窗函数。st840822.28872.2852.2850.2MBdsin0.4 ()1( )( ) ( )( ) , 11.5()2nNh nh n w nw nn7.3 利用频率采样法设计利用频率采样法设计FIR滤波器滤波器1 用频率采样法设计滤波器的基本思想用频率采样法设计滤波器的基本思想设希望逼近的滤波器的频响函数设希望逼近的滤波器的频响函数Hd(ej) ,对它,对它=0到到2之间等之间等间隔采

42、样间隔采样N点,得到点,得到Hd(k):1210| )e ()(2jddNkHkHkN,再对再对Hd(k)进行进行N点点IDFT,得到,得到h(n): (7.3.2)h(n)作为所设计的作为所设计的FIR滤波器的单位脉冲响应,其系统函数滤波器的单位脉冲响应,其系统函数H(z)为为 (7.3.3)1210e )(1)(102jdNnkHNnhNkknN,10)()(NnnznhzH根据根据频率域采样频率域采样理论,利用频率域采样值恢复原信号理论,利用频率域采样值恢复原信号Z变换式,得到变换式,得到H(z)的内插表示形式:的内插表示形式:(7.3.4)此式就是直接利用频率采样值此式就是直接利用频率

43、采样值Hd(k)形成滤波器的系统函形成滤波器的系统函数,(数,(7.3.3)式适合)式适合FIR直接型网络结构,(直接型网络结构,(7.3.4)式)式适合频率采样结构。适合频率采样结构。1012jde1)(1)(NkkNNzkHNzzH 下面讨论两个问题:下面讨论两个问题:一个一个是为了设计线性相位是为了设计线性相位FIR滤滤波器,频域采样序列波器,频域采样序列Hd(k)应满足的条件;应满足的条件;另一个另一个是逼近是逼近误差问题及其改进措施。误差问题及其改进措施。2 设计线性相位滤波器时对设计线性相位滤波器时对Hd(k)的约束条件的约束条件FIR滤波器具有线性相位的条件是滤波器具有线性相位的

44、条件是h(n)为实序列,为实序列,且满足且满足h(n)=h(Nn1),其频响函数应满足的条件是:,其频响函数应满足的条件是:jj ()ddg(e )( )eHH (7.3.5)(7.3.6)21)(N(7.3.7)dgdg( )(2)HHN = 奇数 (7.3.8)dgdg( )(2)HH N = 偶数 在在=02区间上区间上N个等间隔的采样频点为个等间隔的采样频点为将将=k代入(代入(7.3.5)(7.3.8)式中,并写成)式中,并写成k的函数:的函数:kNk2k = 0,1,2,N 1 (7.3.9))(jgde )()(kkHkH(7.3.10)kNNkNNk1221)((7.3.11)

45、)()(ggkNHkHN = 奇数 (7.3.12) )()(ggkNHkHN = 偶数 Hg(k)关于关于N/2点偶对称点偶对称Hg(k)关于关于N/2点奇对称,且点奇对称,且Hg(N/2)=0 设用理想低通作为希望逼近的滤波器设用理想低通作为希望逼近的滤波器Hd(ej),截止频,截止频率为率为c,采样点数为,采样点数为N,Hg(k)和和(k)用下列公式计算:用下列公式计算:(7.3.13) 12101)(1210)(2101)()(cccgcggNkkNNkkNkkkkHkkkNHkH,N=偶数时,偶数时,ggcccgc()1 0,1,2,()0 1,2,1()1 1,2,1() 0,1,

46、2,1cHkkkHkkkkNkHNkkkNkkkNN N=奇数时,奇数时,kc是通带内最后一是通带内最后一个采样点的序号,个采样点的序号,所以所以kc的值取不大的值取不大于于cN/(2)的的最大整数。对于最大整数。对于高通和带阻滤波高通和带阻滤波器,器,N只能取奇数。只能取奇数。3 逼近误差及其改进措施逼近误差及其改进措施时域角度分析:时域角度分析: 待逼近的滤波器为待逼近的滤波器为Hd(ej),对应的单位脉冲响应为,对应的单位脉冲响应为hd(n),在频域,在频域02范围等间隔采样范围等间隔采样N点,利用点,利用IDFT得到的得到的h(n)应是应是hd(n)以以N为周期的周期延拓的主值区序列,

47、即为周期的周期延拓的主值区序列,即d( )()( )Nmh nh nmN Rn由于时域混叠由于时域混叠及截断,使及截断,使h(n)与与hd(n)有有偏差偏差频域角度分析:频域角度分析:H(ej)=FTh(n)的内插表示形式:的内插表示形式:10j2)()e (NkkNkHH式中式中21je)2/sin()2/sin(1)(NNN在采样频点,在采样频点,(2k/N)=1,因此采样,因此采样点处与点处与H(k)相等,逼近误差为相等,逼近误差为0。在采样点之间,在采样点之间,由由N项项H(k)(2k/N)之和形成。之和形成。j(e)kHj(e)kH图图7.3.1 频域幅度采样序列频域幅度采样序列Hg

48、(k)及其内插波形及其内插波形Hg()dg()的误差与的误差与Hdg()特性的平滑特性的平滑程度有关,程度有关,Hdg()特性愈平滑的区特性愈平滑的区域,误差愈小;特性曲线间断点处,域,误差愈小;特性曲线间断点处,误差最大。误差最大。间断点变成倾斜下降的过渡带曲线,间断点变成倾斜下降的过渡带曲线,过渡带宽度近似为过渡带宽度近似为2/N。通带和阻。通带和阻带内产生震荡波纹,且间断点附近带内产生震荡波纹,且间断点附近振荡幅度最大,使阻带衰减减小,振荡幅度最大,使阻带衰减减小,往往不能满足技术要求。往往不能满足技术要求。 增加增加N可以使过渡带变窄,但是通可以使过渡带变窄,但是通带最大衰减和阻带最小

49、衰减随带最大衰减和阻带最小衰减随N的的增大并无明显改善。且增大并无明显改善。且N太大,会太大,会增加滤波器的阶数,即增加了运算增加滤波器的阶数,即增加了运算量和成本。量和成本。提高阻带衰减提高阻带衰减的具体方法是在的具体方法是在频响间断点附近区间内插一个频响间断点附近区间内插一个或几个过渡采样点,使不连续或几个过渡采样点,使不连续点变成缓慢过渡带,这样,虽点变成缓慢过渡带,这样,虽然加大了过渡带,但阻带中相然加大了过渡带,但阻带中相邻内插函数的旁瓣正负对消,邻内插函数的旁瓣正负对消,明显增大了阻带衰减。明显增大了阻带衰减。表表7.3.1 过渡带采样点的个数过渡带采样点的个数m与滤波器阻带最小衰减与滤波器阻带最小衰减 s的经验数据的经验数据 随着随着m, s明显明显s4 频率采样法设计步骤频率采样法设计步骤(1) 根据阻带最小衰减根据阻带最小衰减 s选择过渡带采样点的个数选择过渡带采样点的个数m。(2) 确定过渡带宽度确定过渡带宽度Bt,估算频域采样点数(即滤波器长度),估算频域采样点数(即滤波器长度)N。如果给定过渡带宽度。如果给定过渡带宽度Bt,则要求,则要求(m+1)2/NBt ,滤波器长,滤波器长度度N必须满足如下估算公式必须满足如下估算公式:(7.3.15) t2(1)NmB(3) 构造一个希望逼近的频率响应函数构造一个希望逼

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论