7 有限长单位冲激响应fir滤波器的设计方法_ppt课件

1、第7章 有限长单位冲激响应FIR滤波器的设计方法1第第7 7章章 有限长单位冲激呼应有限长单位冲激呼应FIRFIR滤波器滤波器的设计方法的设计方法7.1 7.1 引言引言7.2 7.2 线性相位线性相位FIRFIR滤波器的特点滤波器的特点 7.3 7.3 用窗函数法设计用窗函数法设计FIRFIR滤波器滤波器 7.4 7.4 用频率采样法设计用频率采样法设计FIRFIR滤波器滤波器 7.5 FIR7.5 FIR滤波器和滤波器和IIRIIR滤波器的比较滤波器的比较第7章 有限长单位冲激响应FIR滤波器的设计方法210( )()Niiy nbx ni7.1 7.1 引言引言 10( )NiiiH z

2、b z 10)()()(Niinxihny10( )( )( )iNiibh iH zh i z对应的系统函数对应的系统函数: : 由于它是一种线性时不变系统，可用卷积和方式表示由于它是一种线性时不变系统，可用卷积和方式表示 比较、得：比较、得：FIRFIR数字滤波器的差分方程描画数字滤波器的差分方程描画: :10( )( )NnnH zh n z第7章 有限长单位冲激响应FIR滤波器的设计方法3FIRFIR数字滤波器的单位抽样呼应数字滤波器的单位抽样呼应h(n)h(n)是有限长的。任何是有限长的。任何一个非因果的有限长序列，总可以经过一定的延时，一个非因果的有限长序列，总可以经过一定的延时，

3、转变为因果序列，所以因果性总是满足；转变为因果序列，所以因果性总是满足；极点全部在原点永远稳定，无稳定性问题；极点全部在原点永远稳定，无稳定性问题；无反响运算，运算误差小，构造普通是非递归的。无反响运算，运算误差小，构造普通是非递归的。很容易获得严厉的线性相位，防止被处置的信号产生很容易获得严厉的线性相位，防止被处置的信号产生相位失真，这一特点在宽频带信号处置、阵列信号处相位失真，这一特点在宽频带信号处置、阵列信号处置、数据传输等系统中非常重要；置、数据传输等系统中非常重要；FIR数字滤波器的特点：数字滤波器的特点：第7章 有限长单位冲激响应FIR滤波器的设计方法47.2 7.2 线性相位线性

4、相位FIRFIR滤波器的特点滤波器的特点 假设假设FIRFIR数字滤波器的单位抽样呼应数字滤波器的单位抽样呼应h(n)h(n)是实数序列，是实数序列， 而而且满足偶对称或奇对称的条件，即且满足偶对称或奇对称的条件，即( )(1)h nh Nn 那么滤波器就具有严厉的线性相位特点。那么滤波器就具有严厉的线性相位特点。第7章 有限长单位冲激响应FIR滤波器的设计方法5一、线性相位特性一、线性相位特性 (1) h(n) (1) h(n)偶对称的情况偶对称的情况: : h(n)=h(N-1-n) h(n)=h(N-1-n) 0nN-1 0nN-1 其系统函数为：其系统函数为： 1010)1()()(N

5、nnNnnznNhznhzH 将将m=N-1-n代入代入 1010)1()1()()()(NmmNmNmNzmhzzmhzH第7章 有限长单位冲激响应FIR滤波器的设计方法6即即 )()(1)1(zHzzHN上式进一步写成：上式进一步写成： (1)11(1)0111221201( )( )()21( )2( )2NNnNnnNNnnNNnH zH zzH zh n zzzzzzh n第7章 有限长单位冲激响应FIR滤波器的设计方法7滤波器的频率呼应为滤波器的频率呼应为 1120()( )|1( )cos2jjz eNNjnH eH zNeh nn可以看到，上式的可以看到，上式的以内全部是标量，

6、假设将频率响运用以内全部是标量，假设将频率响运用相位函数相位函数()()及幅度函数及幅度函数H()H()表示表示 )()()(jjeHeH第7章 有限长单位冲激响应FIR滤波器的设计方法8那么有：那么有： 101( )( )cos21( )2NnNHh nnN 幅度函数幅度函数H()H()是标量函数，可以包括正值、负值和零，是标量函数，可以包括正值、负值和零， 而且是而且是的偶对称函数和周期函数的偶对称函数和周期函数; ; 而而|H(ej)|H(ej)|取值大于等于零取值大于等于零, , 两两者在某些者在某些值上相位相差值上相位相差。相位函数相位函数()()具有严厉的线性相位，如图具有严厉的线

7、性相位，如图7-37-3所示。所示。 第7章 有限长单位冲激响应FIR滤波器的设计方法9图图7-3. h(n)7-3. h(n)偶对称时的线性相位特性偶对称时的线性相位特性 第7章 有限长单位冲激响应FIR滤波器的设计方法10数字滤波器的群延迟数字滤波器的群延迟()()定义为定义为 )()()(ddeHgrdj式中，式中，grd(groupgrd(groupdelay)delay)为群延迟函数。由上式可知，当为群延迟函数。由上式可知，当h(n)h(n)满足偶对称时，满足偶对称时，FIRFIR数字滤波器具有数字滤波器具有(N-1)/2(N-1)/2个采样的延时，个采样的延时， 它它等于单位抽样呼

8、应等于单位抽样呼应h(n)h(n)长度的一半。也就是说，长度的一半。也就是说，FIRFIR数字滤波器数字滤波器的输出呼应整体相对于输入延时了的输出呼应整体相对于输入延时了(N-1)/2(N-1)/2个采样周期。个采样周期。第7章 有限长单位冲激响应FIR滤波器的设计方法11其系统函数为其系统函数为 11001(1)01(1)0( )( )(1)( )( )NNnnnnNNmmNNmmH zh n zh Nn zh m zzh m z 因此因此 H(z)=-z-(N-1)H(z-1) h(n)=-h(N-1-n) 0nN-1 h(n)h(n)奇对称的情况奇对称的情况: : 第7章 有限长单位冲激

9、响应FIR滤波器的设计方法12同样可以改写成 (1)1N-1(1)n 0111221201( )( )()21( )2( )2NnNnNNnnNNnH zH zzH zh nzzzzzzh n第7章 有限长单位冲激响应FIR滤波器的设计方法13其频率呼应为其频率呼应为 112011/220()( )|1( )sin21( )sin2jjz eNNjnNNjjnH eH zNjeh nnNeh nn所以有所以有: : 101( )( )sin21( )22NnNHh nnN 第7章 有限长单位冲激响应FIR滤波器的设计方法14 幅度函数H()可以包括正值、负值和零，而且是的奇对称函数和周期函数。

10、相位函数既是线性相位，又包括/2的相移，如图7-4所示。可以看出，当h(n)为奇对称时，FIR滤波器不仅有(N-1)/2 个采样的延时， 还产生一个90的相移。这种使一切频率的相移皆为90的网络，称为移相器，或称正交变换网络。它和理想低通滤波器、理想微分器一样，有着极重要的实际和实践意义。 当h(n)为奇对称时，FIR滤波器将是一个具有准确的线性相位的正交变换网络。 第7章 有限长单位冲激响应FIR滤波器的设计方法15图图7-4 h(n)7-4 h(n)奇对称时的奇对称时的90o90o线性相位特性线性相位特性 2o()223 N第7章 有限长单位冲激响应FIR滤波器的设计方法16二、二、 幅度

11、呼应特性幅度呼应特性 1. 第一种类型：第一种类型： h(n)为偶对称，为偶对称，N为奇数为奇数 h(n)偶对称的幅度函数式为：偶对称的幅度函数式为： 101( )( )cos2NnNHh nn可以看出，不但可以看出，不但h(n)对于对于(N-1)/2 呈偶对称，而且呈偶对称，而且 也对也对(N-1)/2 呈偶对称，即呈偶对称，即: 1cos2NnnNnNnNNnNhnh21cos21cos)1(21cos)1()(第7章 有限长单位冲激响应FIR滤波器的设计方法17将将内两两相等的项合并，幅度函数就可以表示为内两两相等的项合并，幅度函数就可以表示为 2/ )3(021cos)(221)(Nn

12、nNnhNhH令令 ，那么上式可改写为：，那么上式可改写为： mNn212/ )1(1)cos(21221)(NmmmNhNhH可表示为可表示为 2/ )1(0)cos()()(NnnnaH第7章 有限长单位冲激响应FIR滤波器的设计方法18式中式中: : 21)0(NhanNhna212)(n=1,2,3,(N-1)/2 由于由于cos(n)项对于项对于=0,2皆为偶对称，因此幅度函数皆为偶对称，因此幅度函数H()对于对于=0, ，2也呈偶对称。也呈偶对称。 第7章 有限长单位冲激响应FIR滤波器的设计方法19 2. 第二种类型：第二种类型：h(n)为偶对称，为偶对称，N为偶数为偶数 nNn

13、hHNn21cos)(2)(12/0令令，代入上式可得，代入上式可得 mNn221cos22)(2/1mmNhHNm101( )( )cos2NnNHh nn由于由于N为偶数，因此式中无单独项，全部可以两两合并得为偶数，因此式中无单独项，全部可以两两合并得第7章 有限长单位冲激响应FIR滤波器的设计方法20式中式中: nNhnb22)(n=1,2, 3, , N/2 因此因此 21cos)()(2/1nnbHNn第7章 有限长单位冲激响应FIR滤波器的设计方法21假设数字滤波器在假设数字滤波器在=处不为零，例如高通滤波器、带阻滤波处不为零，例如高通滤波器、带阻滤波器，那么不能用这类数字滤波器来

14、设计。器，那么不能用这类数字滤波器来设计。 当当=时，时， ，余弦项对，余弦项对=呈奇对称，呈奇对称，因此因此H()=0，即，即H(z)在在z=ej=-1 处必然有一个零点，而且处必然有一个零点，而且H()对对=呈奇对称。呈奇对称。 当当=0或或2时，时， 或或-1，余弦项对，余弦项对=0, 2为偶对称，为偶对称，幅度函数幅度函数H()对于对于=0, 2也呈偶对称。也呈偶对称。 121cosn021cosn第7章 有限长单位冲激响应FIR滤波器的设计方法223. 第三种类型：第三种类型： h(n)为奇对称，为奇对称，N为奇数为奇数 h(n)奇对称的幅度函数式如下：奇对称的幅度函数式如下： 10

15、21sin)()(NnnNnhH 由于由于h(n)对于对于(N-1)/2 呈奇对称，即呈奇对称，即h(n)=-h(N-1-n)，当，当n=(N-1)/2时，时， 2121121NhNNhNh因此，因此，, 即即h(n)奇对称时，中间项一定为零。此外，奇对称时，中间项一定为零。此外，式中，式中， 也对也对(N-1)/2 呈奇对称。呈奇对称。 021NhnN21sin第7章 有限长单位冲激响应FIR滤波器的设计方法23nNnNnNN21sin21sin)1(21sin因此，在因此，在中第中第n项和第项和第(N-1-n)项是相等的，将这两两相等的项项是相等的，将这两两相等的项合并，即合并，即 nNn

16、hHNn21sin)(2)(2/ )3(0第7章 有限长单位冲激响应FIR滤波器的设计方法24令令 ， 那么上式可改写为那么上式可改写为 mNn21)sin(21)(2)(2/ )1(1mmNnhHNm即即 )sin()()(2/ )1(1nncHNn式中式中: nNhnc212)(n=1, 2, 3, , (N-1)/2 第7章 有限长单位冲激响应FIR滤波器的设计方法25 由于由于sin(n)在在=0, , 2处都为零，并对这些点呈奇对称，处都为零，并对这些点呈奇对称，因此幅度函数因此幅度函数H()在在=0,2处为零，即处为零，即H(z)在在z=1上都有零上都有零点，且点，且H()对于对于

17、=0,2也呈奇对称。也呈奇对称。 假设数字滤波器在假设数字滤波器在=0, , 2处不为零，例如低通滤波器、处不为零，例如低通滤波器、 高通高通滤波器、带阻滤波器，那么不能用这类数字滤波器来设计，滤波器、带阻滤波器，那么不能用这类数字滤波器来设计， 除除非不思索这些频率点上的值。非不思索这些频率点上的值。 第7章 有限长单位冲激响应FIR滤波器的设计方法26 4. 第四种类型：第四种类型：h(n)为奇对称，为奇对称，N为偶数为偶数12/01021sin)(221sin)(2)(NnNnnNnhnNnhH令令， 那么有那么有 mNn221sin22)(2/1mmNhHNm由于由于N为偶数，因此式中

18、无单独项，全部可以两两合并得为偶数，因此式中无单独项，全部可以两两合并得1021sin)()(NnnNnhH第7章 有限长单位冲激响应FIR滤波器的设计方法27因此因此 21sin)()(2/1nndHNm式中式中: 2/, 3 , 2 , 122)(NnnNhnd当当=0, 2时，时， ，且对，且对=0, 2呈奇对称，呈奇对称，因此因此H()在在=0, 2处为零，即处为零，即H(z)在在z=1处有一个零点，且处有一个零点，且H()对对=0, 2也呈奇对称。也呈奇对称。 021sinn第7章 有限长单位冲激响应FIR滤波器的设计方法28 当当=时，时， 或或1，那么，那么 对对=呈偶对称，幅度

19、函数呈偶对称，幅度函数H()对于对于=也呈偶对称。也呈偶对称。 121sinn21sinn上述四种线性相位上述四种线性相位FIR滤波器的特性示于表滤波器的特性示于表7-1中。中。假设数字滤波器在假设数字滤波器在=0, 2处不为零，例如低通滤波器、处不为零，例如低通滤波器、 带阻滤波器，那么不能用这类数字滤波器来设计。带阻滤波器，那么不能用这类数字滤波器来设计。第7章 有限长单位冲激响应FIR滤波器的设计方法29表表7-1 四种线性相位四种线性相位FIR滤波器特性滤波器特性第7章 有限长单位冲激响应FIR滤波器的设计方法30表表7-1 四种线性相位四种线性相位FIR滤波器特性滤波器特性第7章 有

20、限长单位冲激响应FIR滤波器的设计方法31三、线性相位三、线性相位FIR滤波器的零点位置滤波器的零点位置 线性相位线性相位FIR滤波器的系统函数为：滤波器的系统函数为： H(z)=z-(N-1)H(z-1) 因此，假设因此，假设z=zi是是H(z)的零点，即的零点，即H(zi)=0， 那么那么z=1/zi=zi-1也一定是也一定是H(z)的零点，的零点，(H(zi-1)=zi (N-1) H(zi)=0 当当h(n)是实数时，是实数时，H(z)的零点必成共轭对出现，的零点必成共轭对出现， 所以所以 z=zi*及及z=(z*i)-1也一定是也一定是H(z)的零点，的零点， 因此因此 线性相位线性

21、相位FIR滤波器的零点必是互为倒数的共滤波器的零点必是互为倒数的共轭对。轭对。这种互为倒数的共轭对有四种能够性：这种互为倒数的共轭对有四种能够性： 第7章 有限长单位冲激响应FIR滤波器的设计方法32图图 7-5 线性相位线性相位FIR滤波器的零点位置图滤波器的零点位置图 第7章 有限长单位冲激响应FIR滤波器的设计方法33 由幅度呼应的讨论可知，由幅度呼应的讨论可知，第二种类型的线性相位滤波器第二种类型的线性相位滤波器 H()=0， 因此必然有单根因此必然有单根 z= -1。第四种类型的线性相位滤波器第四种类型的线性相位滤波器 H(0)=0, 因此必然有单根因此必然有单根 z=1。第三种类型

22、的线性相位滤波器第三种类型的线性相位滤波器 H(0)=H(=0， 因此必然有两种单根因此必然有两种单根 z=1 。 了解了线性相位了解了线性相位FIR滤波器的特点，便可根据实践需求滤波器的特点，便可根据实践需求选择适宜类型的选择适宜类型的FIR滤波器，同时设计时需遵照有关的约束滤波器，同时设计时需遵照有关的约束条件。下面讨论线性相位条件。下面讨论线性相位FIR滤波器的设计方法时，都要用滤波器的设计方法时，都要用到这些特点。到这些特点。 第7章 有限长单位冲激响应FIR滤波器的设计方法34去逼近去逼近 , 假设希望得到的滤波器的理想频率呼应为假设希望得到的滤波器的理想频率呼应为: 窗口设计法时域

23、逼近窗口设计法时域逼近 频率采样法频域逼近频率采样法频域逼近 最优化设计等波纹逼近最优化设计等波纹逼近那么那么 FIR滤波器的设计就在于寻觅一个传送函数滤波器的设计就在于寻觅一个传送函数逼近方法有三种：逼近方法有三种：10()( )NjjnddnHeh n e10( )NjjnnH eh n e()jdHe第7章 有限长单位冲激响应FIR滤波器的设计方法357.3 用窗函数法设计用窗函数法设计FIR滤波器滤波器 一、设计方法一、设计方法 窗函数法是设计窗函数法是设计FIR数字滤波器最简单的方法。这数字滤波器最简单的方法。这种方法普通是先给定所要求的理想滤波器的频率呼种方法普通是先给定所要求的理

24、想滤波器的频率呼应应 ，要求设计一个，要求设计一个FIR滤波器频率呼应滤波器频率呼应, 去逼近理想的频率呼应去逼近理想的频率呼应 。10)()(NnnjjenheHdeeHnhnjjdd)(21)(因此，必需首先由理想频率呼应因此，必需首先由理想频率呼应 的傅里叶反变换推导出的傅里叶反变换推导出对应的单位抽样呼应对应的单位抽样呼应: 窗函数法设计窗函数法设计FIR数字滤波器是在时域进展的数字滤波器是在时域进展的, 从单位抽从单位抽样呼应序列着手，使样呼应序列着手，使h(n)逼近理想的单位抽样呼应序列逼近理想的单位抽样呼应序列hd(n)。()jdHe()jdHe()jdHe7-36第7章 有限长

25、单位冲激响应FIR滤波器的设计方法36 由于许多理想化的系统均用分段恒定的或分段函数表示由于许多理想化的系统均用分段恒定的或分段函数表示的频率呼应来定义，因此的频率呼应来定义，因此hd(n)一定是无限长的序列，且是一定是无限长的序列，且是非因果的。而我们要设计的是非因果的。而我们要设计的是FIR滤波器，其滤波器，其h(n)必定是有必定是有限长的，所以要用有限长的限长的，所以要用有限长的h(n)来逼近无限长的来逼近无限长的hd(n)，最，最简单且最有效的方法是截断简单且最有效的方法是截断hd (n) 0)()(nhnhd0nN-1 其他其他 第7章 有限长单位冲激响应FIR滤波器的设计方法37式

26、中假设采用简单截取，那么窗函数为矩形窗。式中假设采用简单截取，那么窗函数为矩形窗。 01)()(nRnwN0nN-1 其他 矩形窗矩形窗 通常，我们可以把通常，我们可以把h(n)表示为所需单位抽样呼应与一个有限长的表示为所需单位抽样呼应与一个有限长的窗口函数序列窗口函数序列w(n)的乘积，即的乘积，即 h(n)=hd(n)w(n) 的波形如以下图所示：的波形如以下图所示：hd( n )RN( n )h ( n )nnnooo( N 1) / 2N 1( N 1) / 2N 1( a )( b )( c )1第7章 有限长单位冲激响应FIR滤波器的设计方法38相应的单位抽样呼应为：相应的单位抽样

27、呼应为： 1( )2sin()()ccjj ndch neednnhd(n)是一个中心点在是一个中心点在的偶对称、无限长、非因果序列，的偶对称、无限长、非因果序列，为了构造一个长度为为了构造一个长度为N的线性相位滤波器，只需将的线性相位滤波器，只需将hd(n)截截取一段，并保证截取的一段对取一段，并保证截取的一段对(N-1)/2 对称，故中心点对称，故中心点a必必需取需取a=(N-1)/2 。0)(jjdeeH|c c2/Nc2/N，H(0)H(0)实践上近似等于实践上近似等于WR()WR()的全部积分的全部积分=-=-到到=+=+面积。面积。 dWHHRd)()(21)(第7章 有限长单位冲

28、激响应FIR滤波器的设计方法462=c时的呼应时的呼应H(c)，Hd()刚好与刚好与WR(-)的一半重的一半重叠，如图叠，如图(c) 。因此卷积值刚好是。因此卷积值刚好是H(0)的一半，即的一半，即H(c)/H(0)=1/2，如图，如图f。第7章 有限长单位冲激响应FIR滤波器的设计方法47 (4)当当 时，时， 主瓣全部在通带外都在主瓣全部在通带外都在Hd()的的 通带通带|c之外，而通带内的旁瓣负的面积大于正之外，而通带内的旁瓣负的面积大于正 的面积，因此卷积结果到达最负值，频响出现负肩峰。的面积，因此卷积结果到达最负值，频响出现负肩峰。2 /cN(3)当当 时，时， 的主瓣全部在的主瓣全

29、部在 的通带的通带内，这时应出现正的肩峰。内，这时应出现正的肩峰。 2cN()RW( )dH第7章 有限长单位冲激响应FIR滤波器的设计方法48(6)当当 时，时， 的右边旁瓣将进入的右边旁瓣将进入 的通的通带，右边旁瓣的起伏呵斥带，右边旁瓣的起伏呵斥 值围绕值围绕 值而动摇。值而动摇。(5)当当 时，随时，随 添加，添加， 左边旁瓣的左边旁瓣的起伏部分扫过通带，卷积起伏部分扫过通带，卷积 也随着也随着 的旁瓣在通的旁瓣在通带内的面积变化而变化，故带内的面积变化而变化，故 将围绕着零值而动摇。将围绕着零值而动摇。2cN()RW( )H( )H()RW2cN()RW( )dH( )H(0)H第7

30、章 有限长单位冲激响应FIR滤波器的设计方法49 综上所述，综上所述， 加窗函数处置后，加窗函数处置后， 对理想频率呼应产生以下几对理想频率呼应产生以下几点影响：点影响： 1H()将将Hd()在截止频率处的延续点变成了延续曲线，在截止频率处的延续点变成了延续曲线， 使理想频率特性不延续点处边沿加宽，构成一个过渡带，过渡使理想频率特性不延续点处边沿加宽，构成一个过渡带，过渡带的宽度等于窗的频率呼应带的宽度等于窗的频率呼应WR()的主瓣宽度的主瓣宽度=4/N，即正，即正肩峰与负肩峰的间隔为肩峰与负肩峰的间隔为 4/N。窗函数的主瓣越宽，过渡带也越。窗函数的主瓣越宽，过渡带也越宽。宽。 2在截止频率

31、在截止频率c的两边即的两边即=c(2/N)的地方，的地方，H()出出现最大的肩峰值，肩峰的两侧构成起伏振荡，其振荡幅度取决现最大的肩峰值，肩峰的两侧构成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的多少，那么取决于旁瓣的多少。于旁瓣的相对幅度，而振荡的多少，那么取决于旁瓣的多少。 第7章 有限长单位冲激响应FIR滤波器的设计方法50 3改动改动N，只能改动窗谱函数的主瓣宽度，改动，只能改动窗谱函数的主瓣宽度，改动的坐标的坐标比例以及改动比例以及改动WR()的绝对值大小。例如，在矩形窗情况下，的绝对值大小。例如，在矩形窗情况下， xxNNNWRsin2/)2/sin()2/sin()2/si

32、n()(式中，式中，x=N/2。 当截取长度当截取长度N添加时，只会减小过渡带宽度添加时，只会减小过渡带宽度4/N，但不，但不能改动主瓣与旁瓣幅值的相对比例能改动主瓣与旁瓣幅值的相对比例; 同样，也不会改动肩峰的相同样，也不会改动肩峰的相对值。这个相对比例是由窗函数外形决议的，与对值。这个相对比例是由窗函数外形决议的，与N无关。换句话无关。换句话说，添加截取窗函数的长度说，添加截取窗函数的长度N只能相应的减少过渡带，而不能改只能相应的减少过渡带，而不能改动肩峰值。动肩峰值。 第7章 有限长单位冲激响应FIR滤波器的设计方法51 由于肩峰值的大小直接影响通带特性和阻带衰减，所以对由于肩峰值的大小

33、直接影响通带特性和阻带衰减，所以对滤波器的性能影响较大。例如滤波器的性能影响较大。例如, 在矩形窗情况下，最大相对肩在矩形窗情况下，最大相对肩峰值为峰值为8.95%，N添加时，添加时，2/N减小，起伏振荡变密，减小，起伏振荡变密， 最大相最大相对肩峰值那么总是对肩峰值那么总是8.95%，这种景象称为吉布斯效应。，这种景象称为吉布斯效应。 第7章 有限长单位冲激响应FIR滤波器的设计方法52二、二、 各种窗函数各种窗函数 矩形窗截断呵斥的肩峰值为矩形窗截断呵斥的肩峰值为8.95%，那么阻带最小衰减为，那么阻带最小衰减为20 lg(8.95%)=-21 dB, 这个衰减量在工程上经常是不够大的。这

34、个衰减量在工程上经常是不够大的。 为了为了加大阻带衰减，加大阻带衰减， 只能改动窗函数的外形。只需当窗谱逼近冲激函只能改动窗函数的外形。只需当窗谱逼近冲激函数时，也就是绝大部分能量集中于频谱中点时，数时，也就是绝大部分能量集中于频谱中点时，H()才会逼近才会逼近Hd()。这相当于窗的宽度为无限长，等于不加窗口截断，这没。这相当于窗的宽度为无限长，等于不加窗口截断，这没有实践意义。有实践意义。 从以上讨论中看出，窗函数序列的外形及长度的选择很关键，从以上讨论中看出，窗函数序列的外形及长度的选择很关键， 普通希望窗函数满足两项要求：普通希望窗函数满足两项要求： 第7章 有限长单位冲激响应FIR滤波

35、器的设计方法53 1窗谱主瓣尽能够地窄，以获取较陡的过渡带。窗谱主瓣尽能够地窄，以获取较陡的过渡带。 2尽量减少窗谱的最大旁瓣的相对幅度。也就是能量尽量尽量减少窗谱的最大旁瓣的相对幅度。也就是能量尽量集中于主瓣，集中于主瓣， 这样使肩峰和波纹减小，就可增大阻带的衰减。这样使肩峰和波纹减小，就可增大阻带的衰减。 但是这两项要求是不能同时都满足的。中选用主瓣宽度较但是这两项要求是不能同时都满足的。中选用主瓣宽度较窄时，虽然得到较陡的过渡带，但通带和阻带的动摇明显添加；窄时，虽然得到较陡的过渡带，但通带和阻带的动摇明显添加； 中选用最小的旁瓣幅度时，虽能得到平坦的幅度呼应和较小的阻中选用最小的旁瓣幅

36、度时，虽能得到平坦的幅度呼应和较小的阻带波纹，但过渡带加宽，也即主瓣会加宽。因此带波纹，但过渡带加宽，也即主瓣会加宽。因此, 实践所选用的实践所选用的窗函数往往是它们的折衷。在保证主瓣宽度到达一定要求的前提窗函数往往是它们的折衷。在保证主瓣宽度到达一定要求的前提下，适当牺牲主瓣宽度以换取相对旁瓣的抑制。下，适当牺牲主瓣宽度以换取相对旁瓣的抑制。 以上是从幅频特性的改善对窗函数提出的要求。实践上设计以上是从幅频特性的改善对窗函数提出的要求。实践上设计的的FIR滤波器往往要求具有线性相位滤波器往往要求具有线性相位:第7章 有限长单位冲激响应FIR滤波器的设计方法54 因此，除了要求因此，除了要求h

37、d(n)满足线性相位条件外，对满足线性相位条件外，对w(n)也要求也要求长度长度N有限，且以有限，且以(N-1)/2为其对称中心，即为其对称中心，即综上所述，窗函数不仅起截断作用，还能起平滑作用，在很多领综上所述，窗函数不仅起截断作用，还能起平滑作用，在很多领域都得到广泛运用。因此域都得到广泛运用。因此, 设计一个特性良好的窗函数有着重要设计一个特性良好的窗函数有着重要的实践意义。的实践意义。 设计设计FIR滤波器常用的窗函数有：滤波器常用的窗函数有： ( )( ) ( )dh nh n w n( )(-1- )w nw Nn第7章 有限长单位冲激响应FIR滤波器的设计方法551. 矩形窗矩形

38、窗 01)()(nRnwN0nN-1 其他其他 )2/sin()2/sin()()()(21NWeWeWRNjRjR第7章 有限长单位冲激响应FIR滤波器的设计方法562. 三角形三角形Bartlett窗窗 12112221012)(NnNNnNnNnnww(n)的傅里叶变换为的傅里叶变换为 212212)2/sin()4/sin(2)2/sin(41sin12)(NjNjjeNNeNNeW 近似结果在近似结果在N1 时成立。时成立。 此时，主瓣宽度为此时，主瓣宽度为8/N, 比矩形窗比矩形窗主瓣宽度添加一倍，主瓣宽度添加一倍， 但旁瓣却小很多。但旁瓣却小很多。 第7章 有限长单位冲激响应FI

39、R滤波器的设计方法573. 汉宁汉宁Hanning窗窗汉宁窗又称升余弦窗。汉宁窗又称升余弦窗。 12( )1 cos( )21Nnw nRnN利用傅里叶变换特性，可得利用傅里叶变换特性，可得 2121)(121225. 0)(5 . 0)(NjNjRRRjeWeNWNWWeW第7章 有限长单位冲激响应FIR滤波器的设计方法58NWNWWWRRR2225. 0)(5 . 0)(当当N1 时，时，N-1N, 所以窗函数的幅频函数为所以窗函数的幅频函数为 这三部分之和，使旁瓣相互抵消，能量更集中在主瓣，但这三部分之和，使旁瓣相互抵消，能量更集中在主瓣，但是代价是主瓣宽度比矩形窗的主瓣宽度添加一倍，是

40、代价是主瓣宽度比矩形窗的主瓣宽度添加一倍， 即为即为 8/N。 第7章 有限长单位冲激响应FIR滤波器的设计方法594. 海明海明Hamming窗窗海明窗又称改良的升余弦窗。海明窗又称改良的升余弦窗。把升余弦窗加以改良，把升余弦窗加以改良， 可以得到旁瓣更小的效果，可以得到旁瓣更小的效果， 窗方式为窗方式为 )(12cos46. 054. 0)(nRNnnwNw(n)的频率呼应的幅度特性为的频率呼应的幅度特性为 NWNWWNWNWWWRRRRRR2223. 0)(54. 0121223. 0)(54. 0)( 与汉宁窗相比，主瓣宽度一样，为与汉宁窗相比，主瓣宽度一样，为 8/N，但旁瓣又被进一

41、步，但旁瓣又被进一步压低，压低， 结果可将结果可将99.963%的能量集中在窗谱的主瓣内。的能量集中在窗谱的主瓣内。 第7章 有限长单位冲激响应FIR滤波器的设计方法60 5. 布拉克曼布拉克曼Blackman窗窗 布拉克曼窗又称二阶升余弦窗。布拉克曼窗又称二阶升余弦窗。 为了进一步抑制旁瓣，对升余弦窗函数再加上一个二次谐为了进一步抑制旁瓣，对升余弦窗函数再加上一个二次谐波的余弦分量，波的余弦分量， 变成布拉克曼窗，故又称二阶升余弦窗。变成布拉克曼窗，故又称二阶升余弦窗。 )(14cos08. 012cos5 . 042. 0)(nRNnNnnwNw(n)的频率呼应的幅度特性为的频率呼应的幅度

42、特性为 22()0.42()0.2511440.0411RRRRRWWWWNNWWNN主瓣宽度是矩形窗的主瓣宽度的主瓣宽度是矩形窗的主瓣宽度的3 3倍倍12/N12/N第7章 有限长单位冲激响应FIR滤波器的设计方法61图图 7-10 五种常用的窗函数五种常用的窗函数 w(n)0.20矩形窗三角窗布拉克曼窗海宁窗海明窗(N1) / 2N1n第7章 有限长单位冲激响应FIR滤波器的设计方法62图图 7-11 图图 7-10 的各种窗函数的傅里叶变换的各种窗函数的傅里叶变换N=51，A=20 lg|W()/W(0)|(a) 矩形窗矩形窗; (b) 巴特利特窗三角形窗巴特利特窗三

43、角形窗; (c) 汉宁窗汉宁窗; (d) 海明窗海明窗; (e) 布拉克曼窗布拉克曼窗 20 40 60 80 10000A / dB 40 80 200 40 80 120 160 2000 20 40 60 80 100 40 80 120 160 2000(d)(e)0 120 160A / dBA / dBA / dBA / dB(a)(b)(c)第7章 有限长单位冲激响应FIR滤波器的设计方法63图图 7-12 理想低通滤波器加窗后的幅度呼应理想低通滤波器加窗后的幅度呼应N=51, A=20lg|H()/H(0)|(a) 矩形窗矩形窗; (b) 巴特利特窗三角形窗巴特利特窗三角形窗;

44、 (c) 汉宁窗汉宁窗; (d) 海明窗海明窗; (e) 布拉克曼窗布拉克曼窗 0306090120150A / dB(a)0306090120150A / dB(b)0306090120150(c)0306090120130A / dB(d )A / dB(e)030609012013000A / dBccccc第7章 有限长单位冲激响应FIR滤波器的设计方法646. 凯塞凯塞Kaiser窗窗这是一种顺应性较强的窗，其窗函数的表示式为这是一种顺应性较强的窗，其窗函数的表示式为 )()1/(21 1()(020INnInw0nN-1 式中，式中，I0(x)是第一类变形零阶贝塞尔函数，是第一类变

45、形零阶贝塞尔函数，是一个可自在选择是一个可自在选择的参数。的参数。 图图7-13 7-13 凯塞窗函数凯塞窗函数 第7章 有限长单位冲激响应FIR滤波器的设计方法65表表7-2 凯塞窗的性能凯塞窗的性能第7章 有限长单位冲激响应FIR滤波器的设计方法66表表7-3 7-3 六种窗函数根本参数的比较六种窗函数根本参数的比较* *最小阻带衰减只由窗外形决议，不受窗宽最小阻带衰减只由窗外形决议，不受窗宽N N的影响；的影响；而过渡带的宽度既和窗外形有关，且随窗宽而过渡带的宽度既和窗外形有关，且随窗宽N N的添加而减小。的添加而减小。第7章 有限长单位冲激响应FIR滤波器的设计方法67下面将窗函数法的

46、设计步骤归纳如下：下面将窗函数法的设计步骤归纳如下：1 给定希望逼近的频率呼应函数给定希望逼近的频率呼应函数Hd(ej)。2 利用式利用式7-36求单位抽样呼应求单位抽样呼应hd(n)=IDTFTHd(ej) 。 deeHnhnjjdd)(21)( 假设假设Hd(ej)很复杂或不能直接计算积分，那么必需用求和很复杂或不能直接计算积分，那么必需用求和替代积分，以便在计算机上计算，也就是要计算离散傅里叶反变替代积分，以便在计算机上计算，也就是要计算离散傅里叶反变换，换， 普通都采用普通都采用FFT来计算。将积分限分成来计算。将积分限分成M段，也就是令采段，也就是令采样频率为样频率为k=2k/M，k

47、=0, 1, 2, , M-1，那么有，那么有 第7章 有限长单位冲激响应FIR滤波器的设计方法6822101( )()kknMjjMMMdkhnHeeM频域的采样，呵斥时域序列的周期延拓，延拓周期是频域的采样，呵斥时域序列的周期延拓，延拓周期是M， 即即 rdMrMnhnh)()( 由于由于hd(n)有能够是无限长的序列，因此严厉说，必需当有能够是无限长的序列，因此严厉说，必需当M时，时，h(n)才干等于才干等于hd(n)而不产生混叠景象，即而不产生混叠景象，即 。实践上，由于。实践上，由于hd(n)随随n的添加衰减很的添加衰减很快，普通只需快，普通只需M足够大，即足够大，即MN，近似就足够

48、了。，近似就足够了。 )(lim)(nhnhMMd第7章 有限长单位冲激响应FIR滤波器的设计方法69 3 由阻带最小衰减及过渡带宽的要求，可选定窗外形，由阻带最小衰减及过渡带宽的要求，可选定窗外形， 并并估计窗口长度估计窗口长度N。设待求滤波器的过渡带用。设待求滤波器的过渡带用表示，它近似等于表示，它近似等于窗函数主瓣宽度。因过渡带窗函数主瓣宽度。因过渡带近似与窗口长度成反比，近似与窗口长度成反比， NA/，A决议于窗口方式。例如，矩形窗决议于窗口方式。例如，矩形窗A=4，海明窗，海明窗A=8等，等，A参数参数选择参考表选择参考表7-3。按照过渡带及阻带衰减情况，选择窗函数方式。按照过渡带及

49、阻带衰减情况，选择窗函数方式。原那么是在保证阻带衰减满足要求的情况下，原那么是在保证阻带衰减满足要求的情况下， 尽量选择主瓣窄的尽量选择主瓣窄的窗函数。窗函数。 4 求得所设计的求得所设计的FIR滤波器的单位抽样呼应。滤波器的单位抽样呼应。 h(n)=hd(n)w(n) 0nN-1 第7章 有限长单位冲激响应FIR滤波器的设计方法705由由h(n)求求FIR滤波器的系统函数滤波器的系统函数H(z) 或或H (ej)=DTFTh(n) 检查能否满足设计要求。检查能否满足设计要求。 通常整个设计过程可利用计算机编程来实现，可多项选择通常整个设计过程可利用计算机编程来实现，可多项选择择几种窗函数来试

50、探，从而设计出性能良好的择几种窗函数来试探，从而设计出性能良好的FIR滤波器。滤波器。 第7章 有限长单位冲激响应FIR滤波器的设计方法71【例【例7-1】根据以下技术目的，设计一个线性相位】根据以下技术目的，设计一个线性相位FIR低低 通滤波器。通滤波器。 抽样频率为抽样频率为 s=2*1.5*104(rad/sec) 通带截止频率为通带截止频率为 p=2*1.5*103(rad/sec) 阻带起始频率为阻带起始频率为st=2*3*103(rad/sec) 阻带衰减不小于阻带衰减不小于50dB。解解: (1) 求对应的数字频率求对应的数字频率通带截止频率：通带截止频率：p=pT= p / f

51、s=2p /s=0.2阻带截止频率：阻带截止频率：st=stT= st / fs=2st /s=0.4阻带最小衰减：阻带最小衰减：2=50dB 第7章 有限长单位冲激响应FIR滤波器的设计方法72(2) 求求hd(n)。设。设Hd(ejw)为理想线性相位低通滤波器为理想线性相位低通滤波器0)(jjdeeH|c c| 0.32spc- )11( )22ccjj njndh needed（理想低通滤波器通带截止频率为理想低通滤波器通带截止频率为 : 由式由式7-39可知，理想低通滤波器的单位抽样呼应为可知，理想低通滤波器的单位抽样呼应为 )()(sin)(nnnhcd21N线性相位要求：线性相位要

52、求： 第7章 有限长单位冲激响应FIR滤波器的设计方法73(3) 求窗函数。求窗函数。 由阻带最小衰减由阻带最小衰减2确定窗外形，由过渡带确定窗外形，由过渡带宽度确定宽度确定N。2=50dB海明窗海明窗stp6.6N过渡带宽度为：过渡带宽度为：由于海明窗过渡带宽度满足：由于海明窗过渡带宽度满足：所以所以6.66.6330.2N1162N第7章 有限长单位冲激响应FIR滤波器的设计方法74(4) 求求h(n)。海明窗为海明窗为 )(12cos46. 054. 0)(nRNnnwN那么所设计的滤波器的单位抽样呼应为那么所设计的滤波器的单位抽样呼应为 ( )( ) ( )sin()

53、20.540.46cos( )()1sin0.3 (16)0.540.46cos( )(16)16dcNNh nh n w nnnRnnNnnRnnsin()( ) ()cdnh nn第7章 有限长单位冲激响应FIR滤波器的设计方法75所设计的滤波器的频率呼应为所设计的滤波器的频率呼应为 10)()(NnnjjenheH 设计结果如设计结果如P345. 图图7-15 所示，满足要求所示，满足要求 。5由由h(n)求求FIR滤波器的滤波器的H (ej)=DTFTh(n)。检查能。检查能否满足设计要求。否满足设计要求。 如不满足要求，那么要改动如不满足要求，那么要改动N，或改动窗外形，或两者，或改

54、动窗外形，或两者都改动，然后重新计算都改动，然后重新计算 。第7章 有限长单位冲激响应FIR滤波器的设计方法76 窗口法设计的主要优点是简单，运用方便。窗口函数大多窗口法设计的主要优点是简单，运用方便。窗口函数大多有封锁的公式可循，性能、参数都已有表格、资料可供参考，有封锁的公式可循，性能、参数都已有表格、资料可供参考， 计算程序简便，计算程序简便， 所以很适用。缺陷是通带和阻带的截止频率不所以很适用。缺陷是通带和阻带的截止频率不易控制。易控制。 第7章 有限长单位冲激响应FIR滤波器的设计方法77 工程上，常给定频域上的技术目的，所以采用频域设计更直接。根本思想 使所设计的FIR数字滤波器的

55、频率特性在某些离散频率点上的值准确地等于所需理想滤波器在这些频率点处的值，在其它频率处的特性那么要有较好的逼近。内插公式内插公式7.4 频率采样设计法频率采样设计法确定确定kNkHeHNkj2)()(102/ )1()2/sin()2/sin()(NjeNN内插函数内插函数:2( )( )( )dkjIDFTjjNddNNhnHeHeHkh nH e 频率取样点点不同于第7章 有限长单位冲激响应FIR滤波器的设计方法78)(ejHN2H(k )(a )(b )oo)(ejdH)(ejH)(ejH图图 7-16 7-16 频率采样的呼应频率采样的呼应 在各频率采样点上，滤波器的实践频率呼应是严厉

56、地和理在各频率采样点上，滤波器的实践频率呼应是严厉地和理想频率呼应数值相等的。但是在采样点之间的频响那么是想频率呼应数值相等的。但是在采样点之间的频响那么是由各采样点的加权内插函数的延伸叠加而成的由各采样点的加权内插函数的延伸叠加而成的, , 因此有一因此有一定的逼近误差，定的逼近误差， 误差大小取决于理想频率呼应曲线外形。误差大小取决于理想频率呼应曲线外形。第7章 有限长单位冲激响应FIR滤波器的设计方法79一、一、 线性相位的约束线性相位的约束 设计线性相位的设计线性相位的FIRFIR滤波器，那么其采样滤波器，那么其采样值值H(k)H(k)的幅度和相位一定要满足前面所讨论的四类的幅度和相位

57、一定要满足前面所讨论的四类线性相位滤波器的约束条件。线性相位滤波器的约束条件。 1 1 对于第一类线性相位滤波器，对于第一类线性相位滤波器， 即即h(n)h(n)偶对称，偶对称， 长度长度N N为奇数时，为奇数时， 12()( )NjjH eHe7-91第7章 有限长单位冲激响应FIR滤波器的设计方法80 第一类线性相位滤波器幅度函数第一类线性相位滤波器幅度函数H()关于关于=0, , 2为偶为偶对称，即对称，即 )2()( HH假设采样值假设采样值H(k)=H(ej2k/N)也用幅值也用幅值Hk纯标量与相角纯标量与相角k表示，表示， 即即 kjkNkjeHeHkH)()(/2并在并在=02之

58、间等间隔采样之间等间隔采样N点点 kNk2k=0, 1, 2, , N-1 7-92第7章 有限长单位冲激响应FIR滤波器的设计方法81由式由式7-91可知，必需有：可知，必需有： 21112kNkkNN 由式由式7-92可知，可知，Hk满足偶对称要求：满足偶对称要求： kN kHH第7章 有限长单位冲激响应FIR滤波器的设计方法82 2 对于第二类线性相位对于第二类线性相位FIR滤波器，即滤波器，即h(n)偶对称，偶对称，N为为偶数，那么其偶数，那么其H(ej)的表达式仍为的表达式仍为: 21)()()()(NeHeHjj 其幅度函数其幅度函数H()关于关于=是奇对称的，关于是奇对称的，关于

59、=0, 2为偶为偶对称，对称， H()=-H(2-) Hk也应满足奇对称要求：也应满足奇对称要求： Hk= - HN-k 21112kNkkNN 第7章 有限长单位冲激响应FIR滤波器的设计方法83 3对于第三类线性相位对于第三类线性相位FIR滤波器，即滤波器，即h(n)奇对称，奇对称，N为为奇数时，奇数时， H(ej)=H()ej()式中：式中： 221)(N 第三类线性相位滤波器幅度函数第三类线性相位滤波器幅度函数H()关于关于=0, , 2为奇对为奇对称，即称，即 )2()(HH2112212NkNkNkkNkHH这样有：这样有：第7章 有限长单位冲激响应FIR滤波器的设计方法84 4对

60、于第四类线性相位对于第四类线性相位FIR滤波器，即滤波器，即h(n)奇对称，奇对称，N为为偶数，那么其偶数，那么其H(ej)的表达式仍为的表达式仍为: 221)()()()(NeHeHjj 其幅度函数其幅度函数H()关于关于=是偶对称的，关于是偶对称的，关于=0, 2为奇对称，为奇对称， 即即 )2()( HH所以，这时的所以，这时的Hk也应满足偶对称要求也应满足偶对称要求 kNkHH而而k那么与前面那么与前面3中的一样。中的一样。 第7章 有限长单位冲激响应FIR滤波器的设计方法85例：设计一个例：设计一个FIR数字低通滤波器，其理想特性为数字低通滤波器，其理想特性为 采样点数采样点数 N=