有限脉冲响应数字滤波器设计.ppt

第八章 有限脉冲响应数字滤波器的设计 本章目录 利用窗函数法设计FIR滤波器 FIR数字滤波器很容易得到严格的线性相位。 FIR数字滤波器的单位脉冲响应是有限长的，因此总是稳定的。 FIR滤波器的设计方法： 窗函数法 频率取样法 等纹波逼近法 理想低通滤波器的频率响应 单位取样响应 将hd(n)截取长度为N的一段，构成h(n) 为了保证设计的滤波器具有线性相位，必须满足对称性要求 可以将h(n)看作是hd(n)与矩形窗w(n)相乘 窗函数设计法的频域解释 时域加窗 频域卷积 窗函数 时域表示 频域表示 幅度 相位 FIR滤波器的幅频特性 FIR滤波器的幅频特性 理想低通与矩形窗频谱函数卷积过程 加窗对Hd(ω)的影响 在理想特性不连续点ωc附近形成过渡带。过滤带的宽度近似等于 W(θ)主瓣宽度，Δω=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，ω=0.2?rad 解：理想数字低通滤波器 单位取样响应 要求设计的FIR数字滤波器的单位取样响应 用不同窗函数设计的FIR滤波器 用矩形窗时过渡带最窄，而阻带衰减最小，布莱克曼窗过渡带最宽，但阻带衰减加大。为保证有同样的过渡带，必须加大窗口长度N 8.3 利用频率抽样法设计FIR滤波器 频率抽样法是在频率域对理想滤波器Hd(ejω)采样，在采样点上设计的滤波器H(ejω)和理想滤波器Hd(ejω)幅度值相等，然后根据频率域的采样值求得实际设计的滤波器的频率特性H(ejω)。 对理想滤波器的频率特性Hd(ejω)在[0,2?]范围内等间隔地取样N个点 频率响应 Hd(k)选定原则 N为偶数时 滤波器性能的改进 采用频率抽样法设计的FIR数字滤波器在阻带内的衰减很小，在实际应用中往往达不到要求。 产生这种现象的原因是由于在通带边缘采样点的陡然变化而引起的起伏振荡。 增加阻带衰减的方法是在通带和阻带的边界处增加一些过渡的采样点，从而减小频带边缘的突变，也就减小了起伏振荡，增大了阻带最小衰减。 8.4 FIR数字滤波器的优化设计 FIR滤波器的优化设计采用“最大误差最小化”的优化准则，根据滤波器的设计指标，导出一组条件，要求在此条件下，在整个逼近的频带范围内使得逼近误差绝对值的最大值为最小，从而得到唯一的最佳解。 可以证明，采用最大误差最小化准则得到的最优滤波器，在通带和阻带内必然呈等纹波特性。 最大误差最小化准则也称为切比雪夫准则。 采用切比雪夫准则设计的滤波器，误差在整个频带均匀分布，对同样的技术指标，这种逼近法需要的滤波器阶数低，而对同样的滤波器阶数，这种逼近法的最大误差最小。 切比雪夫最佳一致逼近法准则 设所希望设计的滤波器幅度函数为 设计一个FIR滤波器，其幅度函数Hg(?)在通带和阻带内最佳地一致逼近Hd(?) 在滤波器的设计中，通带和