matlab数字信号处理fir,MATLAB数字信号处理工具箱的开发和应用——数字滤波器FIR的设计...

154 2010年2月刊算法语言 信息与电脑 China Computer&Communication 1. 引言 在数字信号处理中,数字滤波占有极其重要的地位,目前对数字滤波器的设计有多种方法,其中著名的MA TLAB软件包,功能强大、使用方便。相比传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难，影响了它的应用。利用matlab信号处理工具箱(signal processing toolbox)可以快速有效的设计由软件组成的常规数字滤波器的设计方法，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。 2. FIR滤波器简介 FIR(Finite Impulse Response)滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。 根据数字滤波器冲激响应的时域特征，可将数字滤波器分为2 种，即无限长冲激响应滤波器(IR DF)和有限长冲激响应滤波器((FR DF)。 FR DF具有突出的优点:系统总是稳定的、易于实现线性相位、允许设计多通带(或多阻带)滤波器。因此FR DF在数字信号处理中得到广泛的应用。但与IR D F相比，在满足同样的阻带衰减的情况下需要较高的阶数。滤波器阶数越高将占用更多的D SP运算时间。因此，对FR DF的设计日标是在满足指标要求的情况下尽量减少滤波器的阶数。数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。FR DF的冲激响应h (k)是有限长的M 阶FR DF系统函数可表示为 (1) 滤波器的输出为 (2) 它的设计问题实质上是确定能满足所要求的转移序列或脉冲响应的常数问题,设计方法主要有窗函数法、频率采样法和等波纹最佳逼近法等。 3. FIR滤波器设计步骤 数字滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。随着MATLAB软件尤其是MATLAB的信号处理工具箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。数字滤波器设计的基本步骤如下： 3.1确定指标 在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给出幅度和相位响应。幅度指 标主要以两种方式给出。第一种是绝对指标。它提供对幅度响应函数的要求，用于FIR滤波器的设计。第二种指标是相对指标。它以分贝值的形式给出，这种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中具有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点: 1) 只包含实数算法，不涉及复数运算; 2) 不存在延迟失真，只有固定数量的延迟; 3) 长度为N的滤波器〔阶数为N-1〕，计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位FIR滤波器的设计为例。 3.2逼近 确定了技术指标后，就可以建立一个目标的数字滤波器模型。通常采用理想的数字滤波器模型。之后，利用数字滤波器的设计方法，设计出一个