08-第八节 有限脉冲响应数字滤波器的设计.ppt

第八章 有限脉冲响应数字滤波器的设计 本章目录 利用窗函数法设计FIR滤波器 FIR数字滤波器很容易得到严格的线性相位。 FIR数字滤波器的单位脉冲响应是有限长的，因此总是稳定的。 FIR滤波器的设计方法： 窗函数法 频率取样法 等纹波逼近法 理想低通滤波器的频率响应 单位取样响应 将hd(n)截取长度为N的一段，构成h(n) 为了保证设计的滤波器具有线性相位，必须满足对称性要求 可以将h(n)看作是hd(n)与矩形窗w(n)相乘 窗函数设计法的频域解释 时域加窗 频域卷积 窗函数 时域表示 频域表示 幅度 相位 FIR滤波器的幅频特性 FIR滤波器的幅频特性 理想低通与矩形窗频谱函数卷积过程 (1）ω=0 时, 即H(0)等于WR(θ)在θ=-ωc到θ=+ωc一段的积分面积。通常ωc2π/N，H(0)实际上近似等于WR(θ)的全部积分（θ=-π到θ=+π）面积。 （2）ω=ωc时，Hd(θ)刚好与WR(ω-θ)的一半重叠，因此卷积值刚好是H(0)的一半，即H(ωc)/H(0)=1/2。 （3）ω=ωc-2π/N时，WR(ω-θ)的全部主瓣都在Hd(θ)的通带（|ω|≤ωc）之内。因此卷积结果有最大值，即H(ωc-2π/N)为最大值，频响出现正肩峰。 （4）ω=ωc+2π／Ｎ时，WR(ω-θ)的全部主瓣都在Hd(θ)的通带（|ω|≤ωc）之外，而通带内的旁瓣负的面积大于正的面积，因而卷积结果达到最负值，频响出现负肩峰。 （5）当ωωc+2π/N 当时，随着ω的继续增大，卷积值将随着WR(ω-θ)的旁瓣在Hd(θ)的通带内面积的变化而变化，H(ω)将围绕着零值波动。 (6)当ω由ωc-2π/N向通带内减小时，WR(ω-θ)的右旁瓣进入Hd(θ)的通带，使得H(ω)值围绕H(0)值而波动。 加窗对Hd(ω)的影响 在理想特性不连续点ωc附近形成过渡带。过滤带的宽度近似等于 W(θ)主瓣宽度，Δω=4π/N 。即正肩峰与负肩峰的间隔为 4π/N。窗函数的主瓣越宽，过渡带也越宽。 通带内增加了波动，最大的峰值在ωc- 2π/N 处。阻带内产生了余振，最大的负峰在ωc+2π/N处。通带与阻带中波动的情况与窗函数的幅度谱有关。 W(θ)波动愈快（加大时），通带与阻带内波动愈快， W(θ)旁瓣的大小直接影响波动的大小。 这些影响是对hd(n)加矩形窗引起的，称之为吉布斯效应。 减小吉布斯效应的方法 增加矩形窗口的宽度N不能减少吉布斯效应的影响。 N的改变只能改变ω坐标的比例和 的绝对大小，不能改变主瓣和旁瓣幅度相对值。加大N并不是减少吉布斯效应的有效方法。 寻找合适的窗函数形状，使其谱函数的主瓣包含更多的能量，相应旁瓣幅度就变小了；旁瓣的减少可使通带与阻带波动减少，从而加大阻带的衰减。但这样总是以加宽过渡带为代价的。 几种常用的窗函数 窗函数的要求： ? 窗谱函数主瓣尽可能窄，以获得较陡的过渡带； ? 尽量减少窗谱函数的最大旁瓣的相对幅度，使能量尽量 集中于主瓣，使肩峰和波纹减小，以增大阻带的衰减； ? 通常增加主瓣宽度以换取对旁瓣的抑制； 矩形窗 三角（Bartlett）窗 汉宁（Hanning）窗 汉明（Hamming）窗 布莱克曼（Blackman）窗 凯泽（Kaiser）窗 几种常用的窗函数 矩形窗 三角（Bartlett）窗 三角（Bartlett）窗 汉宁（Hanning）窗 汉宁（Hanning）窗 汉明（Hamming）窗 汉明（Hamming）窗 布莱克曼（Blackman）窗 布莱克曼（Blackman）窗 凯泽（Kaiser）窗 凯泽（Kaiser）窗 常用窗函数的波形 常用窗函数的频谱 用窗函数法设计FIR滤波器的步骤 给出希望设计的滤波器的频率响应函数Hd(ejω)。 根据允许的过渡带宽及阻带衰减，选定窗函数和N值。 计算hd(n) 如果Hd(ejω)不能用简单函数表示，可以用求和代替积分。 将hd(n)与窗函数相乘得FIR数字滤波器的冲激响应h(n) 计算FIR数字滤波器的频率响应，并验证是否达到所要求的指标 例 窗函数法设计FIR滤波器 例8.1 用窗函数法设计线性相位FIR低通滤波器，设N=11，ωc=0.2?rad 解：理想数字低通滤波器 单位取样响应 要求设计的FIR数字滤波器的单位取样响应 用不同窗函数设计的FIR滤波器 用矩形窗时过渡带最窄，而阻带衰减最小，布莱克曼窗过渡带最宽，但阻带衰减加大。为保证有同样的过渡带，必须加大窗口长度N 8.3 利用频率抽样法设计FIR滤波器 频率抽样法是在频率域对理想滤波器Hd(e