基于matlab的FIR数字滤波器设计

1、基于matlab的FIR数字滤波器设计摘 要本文针对窗函数法、频率采样法和最优化设计法三种设计方法，采用MATLAB 进行FIR滤波器的设计与仿真，并比较了三种方法的特点。结果表明，在同样的设计指标下，利用等波纹切比雪夫逼近法则的设计可以获得最佳的频率特性和衰耗特性，具有通带和阻带平坦，过渡带窄等优点。 1. 引言数字滤波器设计，简单地说，就是要找到一组能满足特定滤波要求的系数向量a和b。而滤波器设计完成后还需要进一步考虑如何将其实现，即选择什么样的滤波器结构来完成滤波运算。FIR数字滤波器的设计方法很多，其中较为常用的是窗函数设计法、频率采样设计法和最优化设计法。本文讨论利用窗函数法、频率采

2、样法和等波纹切比雪夫逼近法（调用remez函数）来分别实现各种FIR滤波器的设计。2. FIR滤波器设计方法FIR滤波器设计的可以是选择有限长度的，使传输函数满足一定的幅度特性和线性相位要求。由于FIR 滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。设计过程一般包括以下三个基本问题： （1） 根据实际要求确定数字滤波器性能指标；（2）用一个因果稳定的系统函数去逼近这个理想性能指标；（3）用一个有限精度的运算去实现这个传输函数。 2.1 窗函数法设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数法。FIR数字滤波器

3、的设计首先给出要求的理想滤波器的频率响应,设计一个FIR数字滤波器频率响应，去逼近理想的滤波响应。然而，窗函数法设计FIR数字滤波器是在时域进行的，因而必须由理想的频率响应推导出对应的单位取样响应，再设计一个FIR数字滤波器的单位取样响应去逼近。设计过程如下： 加窗的作用是通过把理想滤波器的无限长脉冲响应乘以窗函数来产生一个被截断的脉冲响应，即并且对频率响应进行平滑。MATLAB工具箱提供的窗函数有：矩形窗(Rectangularwindow)、三角窗(Triangular window)、布拉克曼窗(Blackman window)、汉宁窗(Hanningwindow)、海明窗(Hammin

4、g window)、凯塞窗(Kaiser window)、切比雪夫窗(Chebyshev window)。窗函数主要用来减少序列因截断而产生的Gibbs效应。但当这个窗函数为矩形时，得到的FIR滤波器幅频响应会有明显的Gibbs效应，并且任意增加窗函数的长度（即FIR滤波器的抽头数）Gibbs效应也不能得到改善。为了克服这种现象,窗函数应该使设计的滤波器：（1） 频率特性的主瓣宽度应尽量窄,且尽可能将能量集中在主瓣内；（2） 窗函数频率特性的旁瓣趋于 的过程中，其能量迅速减小为零。2.2频率采样法频率采样法是从频域出发,根据频域采样定理，对给定的理想滤波器的频率响应加以等间隔的抽样 ，得到：

5、k=0,1,N-1 再利用可求得FIR滤波器的系统函数及频率响应。 而在各采样点间的频率响应则是其的加权内插函数延伸叠加的结果。但对于一个无限长的序列，用频率采样法必然有一定的逼近误差,误差的大小取决于理想频响曲线的形状, 理想频响特性变换越平缓, 则内插函数值越接近理想值,误差越小。为了提高逼近的质量，可以通过在频率相应的过渡带内插入比较连续的采样点，扩展过渡带使其比较连续，从而使得通带和阻带之间变换比较缓慢，以达到减少逼近误差的目的。 选取w0,2内N个采样点的约束条件为： 增大阻带衰减三种方法：1）加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。2）过渡带的优化设计利用线性最优化的方法确定

6、过渡带采样点的值，得到要求的滤波器的最佳逼近（而不是盲目地设定一个过渡带值）。3）增大N。如果要进一步增加阻带衰减，但又不增加过渡带宽，可增加采样点数N。代价是滤波器阶数增加，运算量增加。直接从频域进行设计，物理概念清楚，直观方便；适合于窄带滤波器设计，这时频率响应只有少数几个非零值，但是截止频率难以控制。典型应用：用一串窄带滤波器组成多卜勒雷达接收机，覆盖不同的频段，多卜勒频偏可反映被测目标的运动速度。2.3最优化设计最优化设计方法是指采用最优化准则来设计的方法。在 FIR DF的最优化设计中 ,最优化准则有均方误差最小化准则和等波纹切比雪夫逼近(也称最大误差最小化)准则两种。实际设计中 ,

7、只有采用窗函数法中的矩形窗 才能满足前一种最优化准则 ,但由于吉布斯 (Gibbs )效应的存在，使其根本不能满足设计的要求。为了满足设计的要求 ,可以采用其它的窗函数来消除吉布斯效应 ,但此时的设计已经不能满足该最优化准则了。因此 ,要完成 FIR DF的最优化设计 ,只能采用后一种优化准则来实现。3.仿真函数 利用数字信号处理工具箱中的 remezord 和 remez函数可以实现 FIRDF的最优化设计。在此先介绍这两个函数:(1)n ,fo ,ao ,weights =remezordf ,a ,dev功能:利用 remezord 函数可以通过估算得到滤波器的近似阶数 n ,归一化频率

8、带边界fo ,频带内幅值ao 及各个频带内的加权系数weights。输入参数f为频带边缘频率 ,a 为各个频带所期望的幅度值 ,dev是各个频带允许的最大波动。(2)h =remez(n ,fo ,ao ,weights, ftype)功能:利用 remez 函数可以得到最优化设计的FIR DF的系数 ,输入参数 n 是滤波器的阶数 ,fo ,ao ,weights参数含义说明同 (1)。ftype 是所设计的滤波器类型 ,它除了可以设计普通的滤波器外 ,它还可以设计数字希尔钞特变换器以及数字微分器。实际设计中 ,由于 remezord 函数可跑高估或低估滤波器的阶数 n ,因此在得到滤波器的

9、系数后 ,必须检查其阻带最小衰减是否满足设计要求。如果此时的技术指标不能满足设计要求 ,则必须提高滤波器的阶数到 n +1 ,n +2等。故等波纹切比雪夫逼近法设计FIR数字滤波器的步骤是： 给出所需的频率响应，加权函数和滤波器的单位取样响应的长度N。 由中给定的参数来形成所需的、和的表达式。 根据Remez算法,求解逼近问题。 利用傅立叶逆变换计算出单位取样响应 。4 .结果分析（1）窗函数法仿真结果图1 窗函数设计的FIR低通滤波器频率响应从仿真结果可以看出，采用特殊的窗函数如Hamming窗，可以减小Gibbs效应，但同时也会使滤波器的过度带变宽。波动幅度取决于窗函数幅度频谱旁瓣的相对幅

10、度, 而波纹的多少取决于窗函数旁瓣的多少。以上两点是就是窗函数直接截断Hd(n)引起的截断效应在频域的反映, 截断效应直接影响滤波器的性能, 因为通带内的波动会影响滤波器痛带中的平稳性, 阻带内的波动则影响阻带最小衰减,因此,减少截断效应也是FIR数字滤波器设计的关键之一。（2）频率采样法仿真结果图 2 FIR的单位取样响应图3 FIR的低通衰减幅频特性图2为在间断点处增加一个过渡点后的情况。从图3,中可以看出滤波器的带外衰减指标有了明显的改善!但这同时增加了滤波器的过渡带宽。所以，在带外衰减和过渡带宽这两个指标之间需要有一个折衷。因频率取样点都局限在2/N的整数倍点上，所以在指定通带和阻带截

11、止频率时，这种方法受到限制，比较死板。充分加大N，可以接近任何给定的频率，但计算量和复杂性增加。频率采样法偏离设计指标明显，阻带衰减最小，只有适当选取过渡带样点值，才会取得较好的衰耗特性。（3）最优化设计图 4 滤波器输出的幅频及相频响应特性在设计中 ,如果该滤波器的特性不满足要求 ,那么 ,原有参数必须作适当调整。这在程序中很容易实现 ,只需对参数进行重新设定 ,就可以得到新条件下滤波器的特性。采用最优化设计方法时大大减小了滤波器的阶数,从而减小了滤波器的体积,并最终降低了滤波器的成本。这样使得设计出来的滤波器更为简单经济。因而在实际的滤波器设计中,这种最优化方法是完全可行的。在实际应用中

12、,如果需要对某一信号源进行特定的滤波 ,并要检验滤波效果 ,应用传统方法实施起来比较繁琐。在Matlab环境下,可先用软件模拟产生信号源 ,再设计滤波器对其进行滤波 。同样是设计一个FIR低通数字滤波器，综合分析可以看出： （1）窗函数法在阶数较低时，阻带特性不满足设计要求，只有当滤波器阶数较高时，使用海明窗和凯塞窗基本可以达到阻带衰耗要求； （2）频率采样法偏离设计指标最明显，阻带衰减最小，而且设计比采用窗函数法复杂。只有适当选取过渡带样点值，才会取得较好的衰耗特性； （3）利用等波纹切比雪夫逼近法则的设计可以获得最佳的频率特性和衰耗特性，具有通带和阻带平坦，过渡带窄等优点。 5.结论综上所述，FIR 滤波器很容易实现具有严格线性相位的系统, 使信号经过处理后不产生相位失真,舍入误差小,而且稳定，因此越来越受到广泛的重视。MATLAB软件的诞生, 使数字信号处理系统的分析与设计得简单,它已经成为电子工程师必备的一个工具软件。参考文献1 张亚妮. 基于 MATLAB 的数字滤波器设计J. 辽宁工程技术大学学报, 2005, 24 (5): 716-718. 2 饶志强