第7章 FIR数字滤波器的理论与设计

数字信号处理多媒体教学系统版权所有：yuning2003。3第2版第七章

FIR数字滤波器的理论与设计数字信号处理1

（1）以窗函数为基础的窗函数截取法：以传统的窗函数为基础，利用已有的窗函数特性曲线和设计数据进行FIR滤波器的设计。具有设计简单，工程实用价值高的优点。是本课程主要介绍的方法。（2）局部优化设计法：（等波纹逼近法）以理想滤波器特性为基础，设定一、二个过渡带逼近点，然后对FIR滤波器差分方程系数进行优化计算得H(z)。由于需要部分优化计算，所以计算量较大。局部优化设计法主要是针对过渡带进行优化，而通带波动，阻带特性等不一定很好。（3）最优化设计法：（计算机辅助设计）在某种最小化误差准则下，建立差分方程系数bi对理想特性的逼近方程，使用迭代方法解方程组得到最佳逼近系统。由于此方法计算量大，需要借助于计算机进行设计。

FIR滤波器的设计方法：27.1FIR滤波器的线性相位特性：

在一些实际应用中，要求滤波器具有线性相位特性。对于IIR滤波器其相位特性一般都是非线性的，并且很难于实现线性相位特性。而对于FIR滤波器就比较容易实现线性相位特性。7.1.1滤波器的线性相位特性1、系统的时延特性：32、FIR滤波器的线性相位（恒时延）条件4

可以证明，在其他几种对称情况下（冲击响应h(n)为偶/奇对称，无论N为偶数或奇数），FIR滤波器相位特性也都满足线性相位特性。（P149~150）57.1.2线性相位特性FIR滤波器的实现流图：6其余情况的实现流图结构类似。77.1.3线性相位特性FIR滤波器的零、极点：89窗函数截取法以传统的窗函数为基础，利用已有的窗函数特性曲线和设计数据进行FIR滤波器的设计。具有设计简单，工程实用价值高的优点。7.2FIR滤波器的窗函数截取设计方法：7.2.1理想滤波特性的傅立叶级数逼近：这里仍以低通滤波器为例来进行讨论。10理想滤波器冲击序列的加窗处理117.2.2截断特性的误差分析1213)()()(wwwjjdjeWeHeH*=|Hd(ejω)|N-1|W(ejω)|nhd(n)0h(n)=hd(n-N/2).w(n)N-1w(n)14吉布斯效应及改善：

这种由于截断而产生的理想滤波器的幅度特性的波动现象称为吉布斯效应。它使得截断后产生的FIR滤波器特性与理想特性之间有误差，分析误差产生的原因和影响误差的因素可以设计出特性更加好的FIR滤波器。N=8N=16N=3215结论：为了改善逼近效果，需要选用主瓣狭窄，旁瓣相对低的窗函数。但这两项要求往往不能同时满足。因为过渡带宽可以用提高阶数N来满足设计要求，所以工程上主要是选取旁瓣相对低的窗函数来对理想冲击响应序列进行截断。讨论：167.2.3几种常用的窗函数1718矩形窗N=8汉宁窗N=8哈明窗N=83阶Blackman窗N=819矩形窗N=8哈明窗N=82阶Blackman窗N=8汉宁窗N=8窗函数幅度频率特性20矩形窗N=64汉宁窗N=64哈明窗N=642阶Blackman窗N=6421对理想低通特性的逼近矩形窗N=8汉宁窗N=8哈明窗N=82阶Blackman窗N=822矩形窗N=16汉宁窗N=16哈明窗N=162阶Blackman窗N=1623矩形窗N=64汉宁窗N=64哈明窗N=642阶Blackman窗N=6424矩形窗N=8汉宁窗N=8哈明窗N=82阶Blackman窗N=8不同窗函数的阻带特性25矩形窗N=64汉宁窗N=64哈明窗N=642阶Blackman窗N=642627矩形窗N=64汉宁窗N=64哈明窗N=643阶Blackman窗N=64不同的窗函数（N相同）的阻带特性28阶数N对窗函数阻带特性的影响哈明窗N=8哈明窗N=16哈明窗N=64哈明窗N=32297.2.4FIR滤波器的窗函数设计法设计步骤：303132333435（6）序列向右平移（N-1)/2个时间点，成为因果序列：36373839（6）序列向右平移（N-1)/2个时间点，成为因果序列：4041说明：*窗函数设计方法设计的FIR滤波器一般可以满足实际应用需要。如果系统设计指标要求较高，则可能需要复核技术指标，通过调整窗函数和阶数N进一步优化系统设计。或考虑采用计算机辅助设计等其他优化设计方法达到高指标要求。*FIR滤波器的窗函数设计方法不仅可以设计低通滤波器，而且可以设计非低通滤波器。甚至其它非基本类型的滤波特性的系统，只要其频率特性可以进行分段积分，就可以直接使用窗函数设计方法得到相应的FIR滤波系统。*为了得到因果线性相位特性的滤波系统，在窗函数设计方法中一个是使系统频率特性为偶对称，另一个是对理想系统冲击响应序列进行移位再截断。如果是先移位再截断，使用0≤n≤N-1范围的窗函数；如果是先截断再移位，使用（N-1）/2≤n≤（N-1）/2范围的窗函数。*FIR滤波系统始终是稳定系统。实现中没有误差积累，但有较大的时延（NT）。427.3FIR滤波器的频率取样设计法7.3.1设计原理：4344N=64,无过渡点N=64,过渡点=0.5N=64,过渡点=0.6N=64,过渡点=0.390445N=64,无过渡点N=512,无过渡点N=64,过渡点=0.6N=64,过渡点=0.390446结论：频率取样法所得FIR滤波器通带和阻带波动主要是由于过渡带的突变引起的。通过在过渡带优化过渡样点可以得到较好的通带和阻带特性，缺点是过渡带加宽，但可以通过增加样点数N来克服。一般一个过渡样点可以使最小阻带衰耗达-40dB，而2个过渡样点可以使最小阻带衰耗达-60dB。所以频率取样法设计FIR滤波器主要是确定样点数N和过渡样点。N=64,过渡点1=0.5886,过渡点2=0.1065477.3.2频率采样法设计FIR滤波器的步骤4849无过渡点，N=64过渡点H1=0.5，N=64过渡点H1=0.3904，N=64过渡点H1=0.5898,H2=0.1065，N=645051计算结果如下：52537.4IIR滤波器和FIR滤波器的比较1、滤波器特性方面：

（1）IIR滤波器可以通过调整另、极点位置实现较陡峭的幅频特性。所以一般逼近函数阶数较低，但难以实现线性相位特性。

FIR滤波器极点在z=0，无法调整。实现同样的特性需要较高的阶数。容易实现线性相位特性。（2）IIR滤波器用以逼近比较有规则的特性。对于复杂特性，难以确定逼近函数。

FIR滤波器可用分段积分或特性样本点来逼近复杂特性。（3）IIR滤波器时延较小，特别是采取并联实现方案时，容易满足实时系统要求。

FIR滤波器时延为NT，一般比较大。主要用于非实时系统。（4）对系数误差的灵敏度（由于系数误差引起系统特性变化，或不稳定），IIR滤波器较高，特别当极点位置靠近单位圆时。对计算中的有限字长效应敏感。

FIR不会产生稳定问题。542、在实现方法上：（1）IIR滤波器用递归算法、有反馈系统实现，有积累误差产生。

FIR滤波器用非递归算法、无反馈系统实现，没有积累误差。（2）IIR滤波器用差分方程实现（可以用低阶差分方程级联或并联实现）。

FIR滤波器为有限冲击响应，可以用直接卷积实现；也可以用FFT方法实现（分段卷积），有利于大量连续数据的滤波处理。（3）在满足相同的逼近要求下，

IIR滤波器阶数较低，便于用集成电路数字硬件实现。

FIR滤波器阶数高，一般需要用信号处理机和软件实现。55（1）IIR滤波器的设计可以借助于成熟的模拟滤波器设计技术，简化设计过程。

FIR滤波器设计的窗函数设计法中，对通带波动和阻带衰耗不容易控制，一般对特性要求较高的滤波器设计需要反复试算，设计计算量大。3、在设计方法上：（2）IIR滤波器的设计借助于模拟滤波器设计技术，对于非低通滤波特性需要进行频率坐标变换，增加了计算量。对于其他复杂特性要求的滤波器（非基本类型）用IIR形式不容易实现。

FIR滤波器设计的窗函数设计法中，窗函数是对时域的冲击响应进行处理。所以，对于非低通滤波特性或其他复杂的滤波特性只要可以计算分段积分，求取冲击响应序列h(n),都可以直接使用FIR滤波器的窗函数设计法。（3）两种滤波器都可以使用最优化设计方法进行设计。但对于IIR滤波器因为有稳定性问题，所以在最优化设计后需要进行稳定性校验。而FIR滤波器就没有这一问题。56结论：