4.1FIR数字滤波器设计的窗函数法-数字信号处理试卷.pptVIP

第4章 有限冲激响应数字滤波器设计 窗函数法设计FIR数字滤波器 频率抽样法 最佳一致逼近法 FIR数字滤波器 FIR系统与IIR系统比较 数字信号处理研究的内容主要有以下三方面 1．FIR系统只有零点，要取得好的衰减特性，一般要求h（z）的阶次要高，也即M要大。优点：①系统总是稳定的。②易实现线性相位。③允许设计多通带（或多阻带）滤波器，目前，FIR DF的设计方法主要是建立在对理想滤波器的频率特性作某种近似的基础上的。 FIR系统与IIR系统比较 2．IIR系统零、极点全部包括，所以易取得比较好的通带与阻带衰减特性，对于IIR DF的设计方法，主要是借助与模拟滤波器的设计方法，这些面向极点的设计方法不适用于仅包含零点的FIR系统。 对于FIR DF的设计方法包括窗函数法、频率抽样法及最佳一致逼近法 . 4．1 FIR数字滤波器设计的窗函数法 对于一个理想低通数字滤波器，其频率特性为Hd（ejω），现假定其幅频特性 |Hd（ejω）|=1，相频特性Ф(ω)=0，那么，该滤波器的单位抽样响应 hd（n）是以hd（0）为对称的sinc函数， hd（0）=ωc/π，由于n的取值范围为　　　　　-∞~+∞，所以这样的系统是非因果的，因此是物理不可实现的。 4．1 FIR数字滤波器设计的窗函数法 如果我们将hd（n）截短，例如仅取hd（-M/2）…hd（0）…hd（M/2），并将截短后的hd（n）移位，得 那么hd（n）是因果的，且为有限长，长度为M+1，令 4．1 FIR数字滤波器设计的窗函数法 可得所设计的滤波器的转移函数，H(z)的频率响应将近似Hd(ejw)，且是线性相位的。 若 4．1 FIR数字滤波器设计的窗函数法 设计一FIR低通滤波器，所希望的频率响应H d（ejω）在　０≤｜ω｜≤0.25π之间为１，在0.25π＜｜ω｜≤π之间为０，分别取Ｍ＝10，20，40，观察其幅频响应的特点 由上式可得 当M=10时，求得 h（0）＝h（10）＝-０.045, h（1）＝h（9）＝０, h（2）＝h（8）＝０.075, h（3）＝h（7）＝０.1592, h（4）＝h（6）＝０.2251, h（5）=0.25 4．1 FIR数字滤波器设计的窗函数法 吉布斯现象的产生是由于对h d（n）突然截短的结果，将无穷长的h d（n）仅取长为0~M，等于在h d（n）上施加了长为(M+1)的矩形窗口，加窗的结果，等于H d（ejω）和矩形窗频谱的卷积，H d（ejω），矩形窗频谱D（ejω）为sinc(ω)函数，它有着较大的边瓣，正是这些边瓣在和H d（ejω）卷积时产生了吉布斯现象。 实际上，在对hd（n）自然截短时，H（ejω）是对H d（ejω）在最小平方意义上的逼近。 4．1 FIR数字滤波器设计的窗函数法 为了减少吉布斯现象，我们应选取边瓣较小的窗函数，如用汉宁窗或汉明窗来代替矩形窗。 由于本设计方法是选取某种好的窗函数对h d（n）予以截短来得到h（n）的，所以该设计方法称为“窗函数法”又由于H d（ejω）是以２π为周期的周期函数，h d（n）可看作H d（ejω）在频域展开为傅立叶级数时的系数，所以该方法又称“傅立叶级数法” 由上述设计过程可以看出，FIR DF设计的窗函数法不象IIR DF的设计那样能精确指定通带、阻带的边缘频率ωp,ωs，也不能精确地给定在这两个带内的衰减αp，αs，而是仅给出通带截止频率ωc（即ωp），其他几个参数是靠h（n）的长度M及所使用的窗函数的性能来决定的。选定了窗函数后，可以不断地变化M，以检查通带、阻带是否达到了希望的技术要求，直到满意时为止。这些工作在计算机上是极易实现的。 4．2 窗函数 在信号处理中不可避免地要遇到数据截短问题，例如，在应用DFT时，数据x（n）总是有限长，在FIR滤波器设计中遇到了对理想滤波器抽样响应h d（n）的截短问题，在功率谱估计中也要遇到对自相关函数的截短问题，总之，我们在实际工作中所能处理的离散序列总是有限长，把一个长序列变成有限长短序列不可避免地要用到窗函数，因此，窗函数本身的研究及应用是信号处理中的一个基本问题。 设 xn 为一长序列， 是长度为N的窗函数，则 吉伯斯现象 3dB带宽，最大边瓣峰值A（dB）,边瓣谱峰渐近衰减速度D（dB/oct） 一个理想的窗函数，应该具有最小的B和A及最大的D除以上三个指标外，我们对窗函数还有一些共同的要求ω(n)应是非负的实偶函数，且ω(n)从对称中心开始，应是非递增的。 几种常用的窗函数 矩形窗 汉宁窗（Hanning） 汉明窗（Hamming） 布莱克曼窗（Blackman） 凯泽窗（Kaiser） 常