第四章fir滤波器设计02.pptVIP

第四章 FIR数字滤波器的设计 §4.2 窗口设计法（时域）  1、设计原理 2、设计步骤 3、几种常用的窗函数 1、设计原理 去逼近 常用逼近方法：窗口设计法（时域逼近） 频率采样法（频域逼近） 最优化设计（等波纹逼近） 如果希望得到的滤波器的理想频率响应为 ，那么 FIR滤波器的设计就在于寻找一个传递函数 1、设计原理 时间窗口设计法：从单位脉冲响应序列着手 使h(n)逼近理想的单位脉冲响应序列hd(n)。而hd(n)可以从理想频响通过付氏反变换获得 由于Hd(e jω)是矩形频率特性，故hd(n)一定是无限长的序列，问题是怎样用有限长的h(n)来逼近无限长的hd(n) ？ 最简单办法：直接截取一段 hd(n) 代替 h(n) 即h(n)可用表达式为hd(n)和一个“窗函数”的乘积获得 h(n)=w(n) hd(n) 在这里窗口函数w(n)就是矩形脉冲函数RN(n) 为了改善设计滤波器的特性，窗函数w(n) 还可以有其它的形式，相当于在矩形窗内 对hd(n)作一定的加权处理 1）由定义 3）卷积 插值 2、设计步骤 以截止频率为ωc线性相位理想低通滤波器为例， 讨论FIR的设计问题 :低通滤波器的延时 理想特性的hd(n)和Hd(ω) 矩形窗口法 (a) 对于给定的理想低通滤波器 计算 则 Hd(ω) ω 0 ωc -ωc n hd(n) N-1 α=(N-1)/2 图 理想特性的hd(n)和Hd(ω) 这是一个以α为中心的偶对称的无限长非因果序列,如果截取一段n=0～N-1的hd(n)作为h(n), 则为保证所得到的是线性相位FIR滤波器,延时α 应为h(n)长度N的一半，即 其中 矩形窗口法 其中 （b）计算 矩形窗口法 （c）计算 矩形窗口法 设 为窗口函数的频谱： 用幅度函数和相位函数来表示，则有 其线性相位部分 则是表示延时一半长度 对频响起作用的是它的幅度函数 矩形窗口法 设 为窗口函数的频谱： 则实际FIR滤波器的幅度函数H(ω)为 1 0.5 窗口函数对理想特性的影响 ①改变了理想频响的边沿特性，形成过渡带 宽为4π/N，等于WR(ω)的主瓣宽度(决定于窗长) ②过渡带两旁产生肩峰和余振（带内、带外 起伏），取决于 WR(ω)的旁瓣，旁瓣多， 余振多（决定于窗口形状） ③N增加，过渡带宽减小，肩峰值不变。 因主瓣附近 窗口函数对理想特性的影响 其中x = Nω/2 所以 N 的改变不能改变主瓣与旁瓣比例关系 只能改变WR(ω)的绝对值大小和起伏的密度. 当N 增加时，幅值变大频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs）效应。 肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。 改变窗函数的形状,可改善滤波器的特性,窗函数有许多种，但要满足以下要求： ①窗谱主瓣宽度要窄，以获得较陡的过渡带 ②相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就可以减小肩峰和余振，以提高阻带衰减和通带平稳性 但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。 3、几种常用的窗函数 (1) 矩形窗 (2) 汉宁窗（升余弦窗） (3) 汉明窗（改进的升余弦窗） (4) 布莱克曼窗（二阶升余弦窗） (5) 凯塞窗(Kaiser) （1）矩形窗 矩形窗口函数为： 矩形窗的幅度函数表示为： 矩形窗滤波器的主瓣带宽约为4π/N 过渡带为1.8π/N （2） 汉宁窗（升余弦窗） 汉宁窗频率响应函数为： 三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，主瓣宽度增加1倍，为8π/N （2） 汉宁窗（升余弦窗） （3）汉明窗(改进的升余弦窗) 它是对汉宁窗的改进, 在主瓣宽度8π/N (对应第一零点宽度)相同的情况下,旁瓣进一步减小可使99.96%的能量集中在窗谱的主瓣内 其频率响应的幅度函数为： （4）布莱克曼窗(二阶升余弦窗) 它增加一个二次谐波余弦分量，可进一步 降低旁瓣, 但主瓣宽度进一步增加, 为12π/N 其频率响应的幅度函数为： N =51时的四种窗口函数谱的对比 四种窗函数的比较 窗口函数的频谱N=51，A=20lg|W(ω)/W(0)| 窗函数 主瓣宽度 过渡宽度 旁瓣峰值 衰减(dB) 阻带最小 衰