第6章 有限脉冲响应数字滤波器设计.pptVIP

第6章 有限脉冲响应数 字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.1 线性相位FIR数字滤波器的条件和特点 1. 线性相位条件 对于长度为N的h(n)，传输函数为 H(ejω)线性相位是指θ(ω)是ω的线性函数 第一类线性相位（严格线性相位） θ(ω)=—τω, τ为常数 (7.1.3) ? 第二类线性相位（广义线性相位） θ(ω)=θ0—τω,θ0是起始相位 (7.1.4) 均满足： 线性相位FIR数字滤波器也称线性相位FIR 系统。 线性相位FIR 系统广泛应用在数据通信、图像信号处理等领域，在实际工程中具有重要意义。 但并不是FIR 系统就具有线性相位，只有满足一定条件的FIR 系统才具有线性相位。 满足第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即 h(n)=h(N-n-1) (7.1.5) ? 满足第二类线性相位的条件是：h(n)是实序列且对(N-1)/2奇对称，即 h(n)=-h(N-n-1) (7.1.6) (1) 第一类线性相位条件证明： 按照上式可以将H(z)表示为 (2) 第二类线性相位条件证明： 3. 线性相位FIR滤波器零点分布特点 第一类和第二类线性相位的系统函数综合起来用下式表示： 7.2 窗函数设计法 分析矩形窗截断的影响 将H(ejω)写成下式： 窗函数对原理想低通Hd(ω)的影响： (1)在理想特性不连续点ω=ωc附近形成过渡带。过渡带的宽度，近似等于RN(ω)主瓣宽度，即4π/N。 (2)通带内增加了波动，最大的峰值（肩峰）在ωc－2π/N处，影响了滤波器通带的平稳性；阻带内产生了余振，最大的负峰在ωc+2π/N处，影响了阻带最小衰减，阻带内最大负肩峰为-8.95% ，对应的阻带最小衰减为-20lg0.089 ? 21dB。 (3) N增加最大的好处是过渡带减小,但N的改变不能改变主瓣与旁瓣的比例关系，即不能改变波动的幅度，只能使RN（ω）的绝对值幅值变大，频率轴变密。对于矩形窗最大相对肩峰为8.95%，不管N怎样改变，最大肩峰总是 8.95% 。 下面介绍几种常用的窗函数。设 h(n)=hd(n)w(n) ? 式中w(n)表示窗函数。 1. 矩形窗(Rectangle Window) wR(n)=RN(n) 前面已分析过，按照(7.2.5)式，其频率响应为 2. 三角形窗(Bartlett Window巴特利特) 3. 汉宁(Hanning)窗——升余弦窗 4. 哈明(Hamming)窗——改进的升余弦窗，海明窗 5. 布莱克曼(Blackman)窗 总结： 窗函数主瓣宽度要窄，以获得较陡的过渡带； 相对主瓣的幅度，旁瓣要尽可能小，使能量可以集中在主瓣，这样可以减小肩峰和余振，以提高阻带衰减和通带平稳性。 一般总是通过增加主瓣宽度来换取对旁瓣的抑制。 用窗函数设计FIR滤波器的步骤: (1)根据技术要求确定待求滤波器的单位取样响 应hd(n)。 给出待求滤波器的频响为Hd(ejω)，那么单位取样响应用下式求出： 例如，理想低通滤波器: (2)根据对阻带衰减的要求，选择窗函数的形式；并根据对过渡带的要求估计窗口长度N，滤波器的过渡带Δω近似等于窗函数主瓣宽度。 参见表7.2.2 原则：在保证阻带衰减的情况满足要求的情况下，尽量选择主瓣窄的窗函数。 (3) 计算滤波器的单位取样响应h(n) h(n)=hd(n)w(n) w(n)——选择好的窗函数 (4)验算技术指标是否满足要求。设计出的滤波器频率响应用下式计算： 例7.2.1 用矩形窗、汉宁窗和布莱克曼窗设计FIR低通滤波器，设N=11,ωc=0.2πrad。 解：用理想低通作为逼近滤波器，按照(7.2.2)式，有 用汉宁窗设计： 用矩形窗设计： 整个主瓣在积分区间内，一个大的旁瓣在积分区间外， 此时卷积值最大，响应出现峰值。 ?c ??c ? 0 WR(???) ③? =?c ?2? /N ? =?c ?2? /N H[?c?(2?/N) ] H(0) =1.0895 整个主瓣在积分区间