《数字滤波器.docVIP

数字滤波器 摘要 本文首先对数字滤波器的定义、作用及分类作简要介绍，然后重点对经典数字滤波器中的有限冲激响应（FIR）滤波器和无限冲激响应（IIR）滤波器进行说明。分别讨论了FIR数字滤波器和IIR数字滤波器的基本结构和设计方法，并在最后总结了这两种数字滤波器的优缺点。 1. 数字滤波器简介 数字滤波器是指输入和输出均为数字信号，且通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的软件或器件。数字滤波器和快速傅里叶变换一样，是数字信号处理学科的重要组成部分，应用非常广泛。数字滤波器通常是指一种算法或一种数字处理设备，它的功能是将一组输入的数字序列经过一定的运算后变换为另一组输出的数字序列。因此，数字滤波器既可以是用数字硬件装配成的一台完成给定运算的专用的数字计算机，也可以是根据所需要的运算而编写的程序，然后让计算机来执行。 对于数字滤波器而言，若系统函数为H(z)，其脉冲响应为h[n]，当输入为时间序列x[n]，它们在时域内的关系如下： y[n]=h[n]*x[n] 在Z域内，输入和输出关系如下： Y(z)=H(z)*X(z) 式中，X(z),Y(z)分别为x[n]和y[n]的Z变换。在频域内，输入和输出则存在如下关系： Y()=H()*X() 式中，H()是数字滤波器的频率特性，X()、Y()分别为x[n]和y[n]的频谱，而为数字角频率。 数字滤波器可以有很多分类方法，但总体上可分为两大类：一类为经典滤波器，即一般的滤波器，特点是输入信号中的有用成分和希望能滤除的成分占用不同的频带，通过合适的选频滤波器可以实现滤波。例如，若输入信号中有干扰，信号中的有用信号和干扰信号的频带互不重叠，即可滤除干扰得到想要的信号。另一类为现代滤波器，如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器，其输入信号中的有用信号和希望滤除的成分频带重叠，可按照随机信号内部的一些统计分布规律，从干扰中最佳地提取信号。 本文主要简介经典滤波器。该类滤波器，从功能上分类可以分为低通、高通、带通、带阻和全通滤波器；从实现的网络结构或单位冲激响应上分类，可以分为有限长单位冲激响应（FIR）数字滤波器和无限长单位冲激响应（IIR）数字滤波器。对于有限长冲击响应数字滤波器（FIR），其输出y[n]只取决于有限个过去的输入和现在的输入；而对于无线长冲击响应数字滤波器（IIR），它的输出不仅取决于过去和现在的输入，而且还取决于过去的输出。 2. 有限冲激响应数字滤波器 2.1 FIR数字滤波器介绍 FIR滤波器的单位冲激响应为有限长度，一般采用非递归形式实现。N阶因果有限冲激响应滤波器可以用如下传输函数H(z)来描述： （2.1） 该表达式是一个关于的N次多项式。在时域中，上述有限冲激响应滤波器的输入输出关系如下： （2.2） 其中，y[n]和x[n]分别是输出和输入序列。由于有限冲激响应滤波器可以设计成在整个频率范围内均可提供精确的线性相位，而且总是可以独立于滤波器系数保持BIBO稳定，因此在很多领域这种类型的滤波器是首选。 2.2 FIR数字滤波器的基本结构 FIR数字滤波器的基本结构有直接型、级联型、多相实现及线性相位结构。对于一个N阶FIR滤波器要用N+1个系数描述，用N+1个乘法器和N个两输入加法器来实现，而当乘法器的系数正好是传输函数的系数时，这样的结构称为直接型结构，它可以通过式（2.2）实现。FIR数字滤波器的级联结构的实现可以通过将高阶FIR传输函数进行因式分解为一阶或二阶传输函数的积即： 式中当N为偶数时K=N/2，当N为奇数时K=(N+1)/2且，其中的二阶部分可用直接型结构来实现。FIR滤波器的多相实现基于传输函数的多相位分解所得到的并联结构。一般情况下，式（2.1）中L支N阶多相分解的传输函数具有如下形式： 其中 当nN时，h[n]=0。由于多相分解较为抽象，现用一个长度为9的因果有限冲激响应的传输函数来说明，令H(z)为： 则其3支多相分解为： 其中 由于N阶线性相位FIR滤波器可以用对称冲激响应h[n]=h[N-n]或反对称冲激响应h[n]=-h[N-n]来描述，因此在传输函数的直接型实现中，利用线性FIR滤波器的对称、反对称性质可以减少近一半的乘法器。 2.3 FIR数字滤波器的设计 窗函数设计法：一般是先给定所要求的理想滤波器频率响应，由导出，由于理想滤波器的冲击响应是无限长的非因果序列，而我们要设计的是有限长的FIR滤波器，所以要用有限长序列来逼近无限长序列，设： 常用的方法是用有限长度的窗函数w[n]