FIR数字滤波器的理论和设计.pptVIP

数字信号处理基础 第7章 FIR数字滤波器的理论和设计 主讲教师：常 军 联系方式：thinkbank@ 武汉大学计算机学院 计算机应用系 FIR数字滤波器的特点 单位抽样响应h(n)是有限长的，因此FIR数字滤波器一定是稳定的。 经延时，h(n)总可变成因果序列，所以FIR数字滤波器总可以由因果系统实现。 h(n)为有限长，可以用FFT实现FIR数字滤波器。 FIR的系统函数是Z-1的多项式，故IIR的方法不适用。 FIR的相位特性可以是线性的，因此，它有更广泛的应用。 7.1 FIR滤波器的线性相位特性： 7.1.1 滤波器的线性相位特性 1、系统的时延特性： 7.1.3 线性相位特性FIR 滤波器的零、极点： FIR 数字滤波器的系统函数只在 Z=0 处有N-1 阶极点。在Z平面有 N-1 个零点，如系统具有线性相位特性，则系统零点有一些规律。 具有线性相位特性的FIR数字滤波器的系统函数零点具有对称和共轭对称特性(互为倒数的共轭对)。 7.2 FIR滤波器的窗函数截取设计方法： 7.2.1 理想滤波特性的傅立叶级数逼近： 截断后的h(n)作为FIR滤波器的冲击响应(也是差分方程系数),与理想滤波器的差别主要有： (1)存在误差； (2)非因果特性； (3)如何保证线性相位特性(恒时延特性)。 7.2.2 截断特性误差分析 吉布斯效应及改善：这种由于截断而产生的理想滤波器的幅度特性的波动现象称为吉布斯效应。它使得截断后产生的FIR滤波器特性与理想特性之间有误差，分析误差产生的原因和影响误差的因素可以设计出特性更加好的FIR滤波器。 说明： 要取得好的逼近效果，除阶数N外，选择适合的窗函数也是必要的。 通带的波动幅度与窗函数的旁瓣有关，所以选用哈明窗、布拉克曼窗等窗函数可以有较平直的通带特性。 阻带衰耗特性主要取决于窗函数的旁瓣特性，所以为了得到较大的阻带衰耗应选择旁瓣幅度低的窗函数，如哈明窗、布拉克曼窗等。这一点与通带平直特性要求是一致的。 因为特性比较好的窗函数：哈明窗、布拉克曼窗等过渡带都比较宽。所以要采用较高的阶数N来满足过渡带宽要求。 其他还有一些较复杂的窗函数，如凯塞窗、高斯窗等，可以满足较高的设计指标要求(主要是主瓣能量集中，旁瓣幅度低)。但因为计算较复杂，所以一般只用于特殊要求的情况。 说明： 以上方法设计的FIR滤波器一般可以满足实际应用需要。如果系统设计指标要求较高，则可能需要复核技术指标，通过调整窗函数和阶数N进一步优化系统设计。或考虑采用计算机辅助设计等优化设计方法达到高指标要求。 FIR滤波器的窗函数设计方法不仅可以设计低通滤波器，而且可以设计非低通滤波器。甚至其它非基本类型的滤波特性的系统，只要其频率特性可以进行分段积分，也可以直接使用窗函数设计方法得到相应的FIR滤波系统。 为了得到因果线性相位特性的滤波系统，在窗函数设计方法中一个是使系统频率特性为偶对称，另一个是对理想系统冲击响应序列进行移位再截断。(FIR滤波系统始终是稳定系统) 例： 某数字信号处理系统采样频率为1KHz。设计一FIR滤波器其截止频率为125Hz。要求滤波器的冲击响应延时不大于20ms，阻带峰值小于－50dB。 （3）选定滤波器阶数N： 根据冲击响应延时要求 t=(N-1)T，取N-1=20，N=21 例： （3）选定滤波器阶数N：根据要求取N=51。 7.3 FIR滤波器的频率取样设计法 7.3.1 设计原理：通过有限的频率特性取样值去逼近理想滤波特性，然后由有限的频率特性取样值（如系统冲击响应的DFT）取得系统函数 。 7.4 IIR滤波器和FIR滤波器的比较 1、滤波器特性方面： 7.2.4 FIR 滤波器的窗函数设计法设计步骤： （6）序列向右平移（N-1)/2 个时间点，成为因果序列： （6）序列向右平移（N-1)/2 个时间点，成为因果序列： N=64,无过渡点 N=64,过渡点=0.5 N=64,过渡点=0.6 N=64,过渡点=0.3904 N=64,无过渡点 N=512,无过渡点 N=64,过渡点=0.6 N=64,过渡点=0.3904 说明： 频率取样法所得FIR滤波器通带和阻带波动主要是由于过渡带的突变引起的。通过在过渡带优化过渡样点可以得到较好的通带和阻带特性，缺点是过渡带加宽，但可以通过增加样点数N来克服。一般一个过渡样点可以使最小阻带衰耗达-40dB，而2个过渡样点可以使最小阻带衰耗达-60dB。所以频率取样法设计FIR滤波器主要是确定样点数N和过渡样点。 N=64,过渡点1=0.5886,过渡点2=0.1065 7.3.2 频率采样法设计FIR滤波器的步骤 例 : 无过渡点，N=64 过渡点H1=0.5，N=64 过渡点H1=0.3904