FIR数字滤波器的基本理及设计方法.doc

第一章 FIR数字滤波器的基本原理及设计方法 有限长单位脉冲响应数字滤波器（FIRDF，Finite Impulse Response Digital Filter）的最大优点是可以实现线性相位滤波。而IIRDF主要对幅频特性进行逼近，相频特性会存在不同程度非线性。我们知道，无失真传输与滤波处理的条件是，在信号的有效频谱范围内系统幅频响应为常数，相频响应具有线性相位。在数字通信和图像处理与传输等应用场合都要求滤波器具有线性相位特性。另外FIRDF是全零点滤波器，硬件和软件实现结构简单，不用考虑稳定性问题。所以，FIRDF是一种很重要的滤波器，在数字信号处理领域得到广泛应用。当幅频特性指标相同时，FIRDF的阶数比IIRDF高的多，但是同时考虑幅频特性指标和线性相位要求时，IIRDF要附加复杂的相位校正网络，而且难以实现严格线性相位特性。所以，在要求线性相位滤波的应用场合，一般都用FIRDF。 FIRDF的设计方法主要有两类：第一类是基于逼近理想滤波器特性的方法，包括窗函数法，频率采样法和等波纹最佳逼近法。第二类是最优设计法，我们主要讨论第一类设计法，侧重与滤波器的设计方法和相应的MATLAB工具箱函数的介绍。 FIR数字滤波器的设计方法有窗函数法、频率采样法和基于firls函数和remez函数的最优化方法。MATLAB语言中的数字信号处理工具箱，提供了一些滤波器的函数，使FIR滤波器的运算更加方便和快捷。在MATLAB中提供的滤波函数有fir1(),此函数以经典的方法实现加窗线性相位FIR数字滤波器设计，可以设计出低通、高通、带通和带阻滤波器；fir2函数设计的FIR滤波器，其滤波的频率特性由矢量f和m决定，f和m分别为滤波器的期望幅频响应的频率相量和幅值相量。Firls()和remez()的基本格式用于设计I型和II型线性相位FIR滤波器，I型和II型的区别是偶函数还是奇函数。freqz()用于求数字滤波器的频率响应。并且提供了各种窗函数的函数，比如，hamming()是海明窗函数，hanning()是汉宁窗函数，kaiser()是凯泽窗函数，使在设计的过程中，不用自己重新设计窗函数。 1.1 窗函数法设计FIR数字滤波器 设我们所要设计的FIR滤波器的传输函数是(e),(n)是与其对应的单位脉冲响应，因此 （1-1） （1-2) 如果我们能够在已知的情况下，求出，经过Z变换可得到滤波器的系统函数。通常情况下理想数字滤波器的单位脉冲相应是无限长的，且是非因果序列。获得有限脉冲响应滤波器的一种可能方法是对截取一段来近似代替，可是这样会改变原来的滤波器指标，出现吉布斯效应误差。 窗函数法就是用被称为窗函数的有限加权序列w(n)来修正式(1)的傅里叶基数以求得要求的有限脉冲响应序列，即 （1-3） w(n)是有限长序列，当n0或nN-1时，w(n)=0。 这种方法的重点在于选择某种合适的窗函数。要求窗函数主瓣宽度尽可能窄，以获得最小的过渡带；旁瓣相对值尽可能小，以使得通带波纹小，并且阻带衰减大。 下面介绍几种常用的窗函数： 1.矩形窗(Rectangle Window) （3-4） 其频率为： （3-5） 2.三角形窗(Bartlett Window) （3-6） 其频率为： （3-7) 3.汉宁(Hanning)窗，又称升余弦窗 （3-8） 利用傅里叶变换得到频率为： （3-9） 当时，，所以窗函数的幅度函数为 (3-10) 4.汉明(Hamming)窗改进的升余弦窗 （3-11） 其幅度为： （3-12） 5.布莱克曼(Blankman)窗，又称二阶升余弦窗 （3-13) 其幅度为： (3-14） 6.凯泽(Kaiser)窗 （3-15） 其中：β是一个可选参数，I0()是第一类修正零阶贝塞尔函数[/N -21 三角形窗 -26 8/N -25 汉宁窗 -31 8/N -44 汉明窗 -40 8/N -53 布莱克曼窗 -57 12/N -74 凯泽窗 -57 10/N -80 窗函数法设计滤波器的步骤： 1）根据技术要求确定待求滤波器的单位取样响应。 2）根据对过渡带和阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。 3）计算滤波器的单位取样响应h(n)： （3-16) 式中,是前面所选择