数字信号处理及应用工业和信息化普通高等教育“十二五”规划教材立项项目新教学课件作者卢光跃黄庆东包志强第7章节有限脉冲响应数字滤波器的设计.pptVIP

第7章 有限长单位脉冲响应（FIR）数字滤波器的设计方法 第7章 有限脉冲响应数字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.3 利用频率采样法设计FIR滤波器 7.4 利用切比雪夫逼近法设计FIR滤波器 7.5 IIR和FIR数字滤波器的比较 IIR数字滤波器： (1) 第一类线性相位条件证明： 按照上式可以将H(z)表示为 (2) 第二类线性相位条件证明： 2、幅度函数的特点 1）h(n)偶对称，N为奇数 3）h(n)奇对称，N为奇数 4）h(n)奇对称，N为偶数 3 零点位置 §7.2 窗函数设计法 1、设计原理（频率响应逼近理想滤波器特性） 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 取矩形窗： 加窗处理后对频率响应的影响： 时域乘积相当于频域卷积 理想滤波器的频率响应： 2、各种窗函数 窗函数的要求： 三角形（Bartlett）窗 汉宁（Hanning）窗 （升余弦窗） 海明（Hamming）窗 （改进的升余弦窗） 布莱克曼（Blackman）窗 （二阶升余弦窗） 凯泽（Kaiser）窗 3、窗函数法的设计步骤 2）求hd(n) 5）确定FIR滤波器的h(n) 7.5 IIR和FIR数字滤波器的比较 首先，从性能上来说，IIR滤波器传输函数的极点可位于单位圆内的任何地方，因此可用较低的阶数获得高的选择性，所用的存贮单元少，所以经济而效率高。但是这个高效率是以相位的非线性为代价的。 从结构上看，IIR滤波器必须采用递归结构，极点位置必须在单位圆内，否则系统将不稳定。 从设计工具看，IIR滤波器可以借助于模拟滤波器的成果，因此一般都有有效的封闭形式的设计公式可供准确计算，计算工作量比较小，对计算工具的要求不高。 而矩形窗的频率响应： 其幅度函数： 则FIR滤波器的频率响应： 不连续点处边沿加宽形成过渡带，其宽度（两肩峰之间的宽度）等于窗函数频率响应的主瓣宽度,即 。 加窗函数的影响： 在 处出现肩峰值，两侧形成起伏振 荡，振荡的幅度和多少取决于旁瓣的幅度和多少 改变N只能改变窗谱的主瓣宽度，但不能改变主瓣与旁瓣的相对比例。其相对比例由窗函数形状决定，称为Gibbs效应 1 窗谱主瓣尽可能窄以获得较陡的过渡带 2 尽量减少窗谱最大旁瓣的相对幅度以减小肩峰和波纹 矩形窗 主瓣宽度最窄： 旁瓣幅度大 窗谱： 幅度函数： 主瓣宽度宽： 旁瓣幅度较小 窗谱： 幅度函数： 主瓣宽度宽： 旁瓣幅度小 幅度函数： 主瓣宽度宽： 旁瓣幅度更小 幅度函数： 主瓣宽度最宽： 旁瓣幅度最小 幅度函数： ：第一类变形零阶 贝塞尔函数 -21 -25 -44 -53 -74 -80 0.9 2.1 3.1 3.3 5.5 5 2 4 4 4 6 -13 -25 -31 -41 -57 -57 矩形窗 三角形窗 汉宁窗 海明窗 布拉克曼窗 凯泽窗 阻带最小衰减 /dB 过渡带宽 主瓣宽度 旁瓣峰值 /dB 加窗后滤波器性能指标 窗谱性能指标 窗函数 阻带最小衰减只由窗形状决定 过渡带宽则与窗形状和窗宽N都有关 给定理想的频率响应函数 及技术指标 求出理想的单位抽样响应 根据阻带衰减选择窗函数 计算频率响应 ，验算指标是否满足要求 根据过渡带宽度确定N值 求所设计的FIR滤波器的单位抽样响应 公式法： IFFT法： 计算其IFFT，得： 对 M点等间隔抽样： 4、线性相位FIR低通滤波器的设计 通带截止频率为 ， 阻带起始频率为 ， 解： 1） 例：设计一个线性相位FIR低通滤波器， 阻带衰减不小于-50dB， 4）确定N 值 3）选择窗函数：由 确定海明窗（-53dB） 6）求 ，验证 若不满足，则改变N或窗形状重新设计 离散时间信号—序列 尚辅网 / 通信与信息工程学院 数字信号处理教学团队 可以利用模拟滤波器设计 但相位非线性 FIR数字滤波器： 可以严格线性相位，又可任意幅度特性 因果稳定系统 1、线性相位条件 h(n)为单位取样响应，其频率响应： 第一类线性相位： 线