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

第七章 FIR 有限脉冲响应数字滤波器的设计 本章目录 利用窗函数法设计FIR滤波器（吉布斯现象） 引言 IIR数字滤波器： FIR滤波器的设计方法： 7.1 线性相位FIR滤波器的条件和特点 FIR滤波器的单位冲激响应： 1、线性相位条件 h(n)为实序列时，其频率响应： 4、线性相位FIR滤波器频率响应的特点 频率响应： 频率响应： 结论： 当h(n)为实序列 且当h(n)奇或偶对称 3）、幅度函数的特点 A、h(n)偶对称，N为奇数 说明 B、h(n)偶对称，N为偶数 C、h(n)奇对称，N为奇数 小结： 说明 5、线性相位FIR DF的零点分布特点 即共轭成对且镜像成对 设N为偶数，有 小结： 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 采取的措施： 将hd(n)截取长度为N的一段，构成h(n) 为了保证设计的滤波器具有线性相位，必须满足对称性要求 可以将h(n)看作是hd(n)与矩形窗RN(n)相乘 二、改进的设计思路 三、时域加窗对频域的影响 时域加窗处理 频域卷积 窗函数 时域表示 频域表示 幅度 相位 FIR滤波器的幅频特性 FIR滤波器的幅度特性 理想低通与矩形窗频谱函数卷积 加窗对H(ω)的影响:Gibbs效应 在理想特性不连续点ωc附近形成过渡带。过滤带的宽度近似等于RN(θ)主瓣宽度，Δω=4π/N ，过渡带越小越好。 通带内增加了波动，最大的峰值在ωc- 2π/N 处。阻带内产生了余振，最大的负峰在ωc+2π/N处。 通带与阻带中波动的情况与窗函数的幅度谱有关。 RN(θ)波动愈快（加大时），通带与阻带内波动愈快， RN(θ)旁瓣的大小直接影响波动的大小。 这些影响是对hd(n)加矩形窗引起的，称之为吉布斯效应。 结论一 时域截断，频域产生吉布斯现象 频域截断，时域产生吉布斯现象 结论二 四、减小吉布斯效应的方法 增加矩形窗口的宽度N，可以减少过渡带的宽度， Δω=4π/N 。 但加大N不能减少吉布斯效应的影响。 N加大时，主、旁瓣幅度同时加高，相对值不变，总是最大肩峰比H(0)高8.95%,最大负峰比0值超过8.95%。 不能改变主瓣和旁瓣幅度相对值，即不能加大阻带衰减。---见p206页图 选择合适的窗函数形状，使其谱函数的主瓣包含更多的能量，相应旁瓣幅度就变小了；旁瓣的减少可使通带与阻带波动减少，从而加大阻带的衰减。但这样总是以加宽过渡带为代价的。 五、几种常用的窗函数 矩形 (Rectangle)窗 三角（Bartlett）窗 汉宁（Hanning）窗 汉明（Hamming）窗 布莱克曼（Blackman）窗 凯泽（Kaiser）窗 三角形（Bartlett）窗 汉宁（Hanning）窗 （升余弦窗） 哈明（Hamming）窗 （改进的升余弦窗） 布莱克曼（Blackman）窗 （二阶升余弦窗） 常用窗函数的频谱 理想低通加窗后的幅度特性N=51,ωc=0.5π 凯泽（Kaiser）窗 六种窗函数的基本参数 注意事项： 六、用窗函数法设计FIR滤波器的步骤 给出希望设计的滤波器的频率响应函数Hd(ejω)。 计算hd(n) 如果Hd(ejω)不能用简单函数表示，可以用求和 代替积分。 根据允许的过渡带宽及阻带衰减，选定窗函数和N值，见书表7-2, 7-3 5. 将hd(n)与窗函数相乘得FIR数字滤波器的冲激响应h(n)。 6. 计算FIR数字滤波器的频率响应，并验证是否达到所要求的指标，若不能满足要求，根据情况重复4、5、6步。 例 窗函数法设计FIR滤波器 例7.2.1 用窗函数法设计线性相位FIR低通滤波器，设N=11，ω=0.2?rad 解：理想数字低通滤波器 单位取样响应 要求设计的FIR数字滤波器的单位取样响应 用不同窗函数设计的FIR滤波器 用矩形窗时过渡带最窄，而阻带衰减最小，布莱克曼窗过渡带最宽，但阻带衰减加大。为保证有同样的过渡带，必须加大窗口长度N Matlab函数fir1的调用格式 函数fir1的调用格式为 b= fir1(n, Wc, ftype, Window) 其中，n为滤波器阶数，Wc为截止频率 ftype决定滤波器类型，ftype= high，设计高通FIR滤波器，ftype=bandpass、 stop，设计带通、带阻FIR滤波器。 Window指定窗函数类型，默认为Hamming窗；可选Hanning、Hamming、Blackman、bartlett和boxcar窗，每种窗都可以由Matlab的相应函数生成。 MATLAB实现 2）求hd(n) 5）确定FIR滤波器的h(n) 例7.2.3 线性相位FIR高通滤波器的设