《数字信号处理教学课件》DSP-7.PPT

第7章 有限脉冲响应数字滤波器的设计 本章主要内容 线性相位FIR数字滤波器的条件和特点 利用窗函数法设计FIR滤波器 利用频率采样法设计FIR滤波器 IIR和FIR数字滤波器的比较 7.1 线性相位FIR数字滤波器的条件和特点 介绍FIR滤波器具有线性相位的条件、 FIR滤波器的幅度特性、零点分布特点和网络结构的特点。 1. 线性相位条件 对于长度为N的h(n)，传输函数为： H(ejω)线性相位： θ(ω)是ω的线性函数，即： θ(ω)= -τω, τ为常数 (7.1.3) ? 若θ(ω)满足下式： θ(ω)=θ0 -τω, θ0是起始相位 (7.1.4) 也称这种情况为线性相位。 以上两种情况都满足群时延是一个常数，即 3. 线性相位FIR滤波器零点分布特点 线性相位的系统函数满足： 相应的单位取样响应hd(n)为： 如何减少吉布斯效应的影响，设计一个满足要求的FIR滤波器呢？ 直观上，增加矩形窗口的宽度（即加大N）可以减少吉布斯效应的影响。N ? 时， 在主瓣附近， WRg(ω)近似为： 1. 矩形窗(Rectangle Window) 本章作业 3 5 20（实验） （7.2.7） 　　 为了描述方便，定义窗函数的几个参数: 　　 旁瓣峰值?n——窗函数的幅频函数|Wg(ω)|的最大旁瓣的最大值相对主瓣最大值的衰减值（dB）； 　　过渡带宽度Bg——用该窗函数设计的FIR数字滤波器（FIRDF）的过渡带宽度； 　　阻带最小衰减?s——用该窗函数设计的FIRDF的阻带最小衰减。 1． 矩形窗（Rectangle Window） wR(n)=RN(n) 其幅度函数为 图7.2.4 矩形窗函数长度的影响 2． 三角形窗（Bartlett Window） （7.2.8） 其频谱函数为 其幅度函数为 　　　　　　　　　　　　　　　　　 （7.2.10） 图7.2.5 三角窗的四种波形 主瓣能量更集中，但过渡带宽了很多 3． 汉宁（Hanning）窗——升余弦窗 （7.2.11） 当N1时， N－1≈N 汉宁窗的幅度函数WHng(ω)由三部分相加，旁瓣互相对消，使能量更集中在主瓣中。汉宁窗的四种波形如图7.2.6所示，参数为: ? n=－31 dB; Bg=8π/N; ? s=－44 dB。 图7.2.6 汉宁窗的四种波形 4． 哈明（Hamming）窗——改进的升余弦窗 　　　　　　　　　　　　　　　 （7.2.12）其频谱函数WHm(ejω)为 其幅度函数WHmg(ω)为 当N１时，其可近似表示为 图7.2.7 哈明窗的四种波形 这种改进的升余弦窗，能量更加集中在主瓣，但其主瓣宽度和汉宁窗的相同 （7.2.13） 5． 布莱克曼（Blackman）窗 　　　　　　　　　　　　　　　　　　　 其频谱函数为 其幅度函数为 幅度函数由五部分组成，它们都是移位不同，且幅度也不同的WRg(ω)函数，使旁瓣再进一步抵消。旁瓣峰值幅度进一步增加，其幅度谱主瓣宽度是矩形窗的3倍 图7.2.8 布莱克曼窗的四种波形 旁瓣进一步减小，但过渡带变宽 6． 凯塞—贝塞尔窗（Kaiser-Basel Window） 　　以上五种窗函数都称为参数固定窗函数，每种窗函数的旁瓣幅度都是固定的。凯塞—贝塞尔窗是一种参数可调的窗函数，是一种最优窗函数。 　　　　　　　　　　　　　　　 （7.2.15） 式中 I0(β)是零阶第一类修正贝塞尔函数，可用下面级数计算： 一般I0(β)取15~25项，便可以满足精度要求。 ?参数可以控制窗的形状。一般?加大，主瓣加宽，旁瓣幅度减小，典型数据为4 ? 9。当? =5.44时，窗函数接近哈明窗。 ? =7.865时，窗函数接近布莱克曼窗。在设计指标给定时，可以调整?值，使滤波器阶数最低，所以其性能最优。凯塞（Kaiser）给出的估算β和滤波器阶数N的公式如下: （7.2.17） 式中，Bt=|ωs－ωp|, 是数字滤波器过渡带宽度。应当注意，因为式（7.2.17）为阶数估算，所以必须对设计结果进行检验。另外，凯塞窗函数没有独立控制通带波纹幅度，实际中通带波纹幅度近似等于阻带波纹幅度。凯塞窗的幅度函数为 （7.2.16） （7.2.18） 对?的8种典型值，将凯塞窗函数的性能列于表7.2.1中，供设计者参考。由表可见, 当? =5.568时, 各项指标都好于哈明窗。 表7