第七章FIR滤波器的窗函数设计.doc

第7章 FIR滤波器设计 第六章我们介绍了无限冲激响应（IIR）滤波器的设计方法。其中最常用的由模拟滤波器转换为数字滤波器的方法为双线性变换法，因为这种方法无混叠效应，效果较好。但通过前面的例子我们看到，IIR数字滤波器相位特性不好（非线性，如图 6-11、图6-13、图6-15 ），也不易控制。然而在现代信号处理中，例如图像处理、数据传输、雷达接收以及一些要求较高的系统中对相位特性要求较为严格，这种滤波器就无能为力了。改善相位特性的方法是采用有限冲激响应滤波器。本章首先对FIR滤波器原理及其使用函数作基本介绍，然后重点介绍窗函数法设计FIR滤波器，并对最优滤波器设计函数进行介绍。 7.1 FIR滤波器原理概述及滤波函数 7.1.1 FIR滤波器原理及设计方法分类 根据第 6 章对数字滤波器的介绍，我们知道FIR滤波器的传递函数为： (7-1) 可得FIR滤波器的系统差分方程为： 因此，FIR滤波器又称为卷积滤波器。根据第 4 章中所描述的系统频率响应，FIR滤波器的频率响应表达式为： （7-2） 信号通过FIR滤波器不失真条件与（6-6）式所描述的相同，即滤波器在通带内具有恒定的幅频特性和线性相位特性。理论上可以证明（这里从略）：当FIR滤波器的系数满足下列中心对称条件： （7-3） 时，滤波器设计在逼近平直幅频特性的同时，还能获得严格的线性相位特性。线性相位FIR滤波器的相位滞后和群延迟在整个频带上是相等且不变的。对于一个 N 阶的线性相位FIR滤波器，群延迟为常数，即滤波后的信号简单地延迟常数个时间步长。这一特性使通带频率内信号通过滤波器后仍保持原有波形形状而无相位失真。 本章主要介绍的FIR数字滤波器设计方法及 MATLAB 信号处理工具箱提供的了FIR数字滤波器设计函数，见表7-1。由于篇幅所限，本章我们主要介绍窗函数法和最优化设计方法。 表7-1 FIR滤波器设计的主要方法 函数设计方法 说明 工具函数 窗函数法 理想滤波器加窗处理 fir1(单频带) ， fir2(多频带) ， kaiserord 最优化设计 平方误差最小化逼近理想幅频响应或Park-McClellan 算法产生等波纹滤波器 firls ， remez,remezord 约束最小二乘逼近 在满足最大误差限制条件下使整个频带平方误差最小化 fircls,fircls1 升余弦函数 具有光滑、正弦过渡带的低通滤波器设计 Fircos 7.1.2 FIR数字滤波器滤波函数 相对于IIR 滤波器的滤波函数，FIR数字滤波器滤波函数除了dimpulse和dstep仅适用于IIR滤波器外，其他各种函数可直接应用于FIR滤波器，只是输入的分母多项式向量a=1。另外，MATLAB还提供了一个函数fftfilt,该函数利用效率高的基于FFT算法实现对数据的滤波，该函数只适用于FIR滤波器，调用形式为： y=fftfilt(b,x[,n]) 式中，b为FIR滤波器的系数向量；x为输入数据；n为FFT长度，缺省时，函数选用最佳的FFT长度，y为滤波器的输出。该函数执行下面的操作： n=length(x); y=ifft(fft(x).*fft(b,n)./fft(a,n)); 应注意，y=fftfilt(b,x)等价于y=filter(b,a,x)。 7.2 FIR滤波器的窗函数设计 7.2.1 窗函数的基本原理 FIR滤波器设计的主要任务是根据给定的性能指标确定滤波器的系数b，即系统单位脉冲序列h(n)，它是一个有限长序列。 FIR滤波器的理想频率响应，可写成复数形式的Fourier级数形式： (7-4) 式中，hd(n)是对应的单位脉冲响应序列。这说明滤波器的频率响应和单位脉冲响应互为Fourier变换对。因此其单位脉冲响应可由下式求得， (7-5) 求得序列后，通过z变换，可得到 (7-6) 注意，这里为无限长序列，因此是物理上不可实现的。如何变成物理上可实现呢？一个自然的想法是只取其中的某些项，即只截取中的一部分，比如n=0,…,N-1，N为正整数。这种处理相当于将,n=-∞~∞与函数w(n)相乘，w(n)具有下列形式： w(n)相当于一个矩形，我们称之为