DSP“数字信号处理”07第七章lzh.pptVIP

七、线性相位FIR带阻滤波器的设计 一个理想的FIR带阻滤波器可通过一个全通系统减去两个FIR低通滤波器的差实现(三者有相同的相位响应)。 理想的FIR带阻滤波器的频率响应为 理想的FIR带阻滤波器的单位冲击响应为 结束 return 幅频特性 幅度特性 相位特性 返回 return return return 偶对称 return 奇对称 不能用来设计高通、带阻滤波器 return 奇对称 不能用来设计低通、高通、带阻滤波器 return 不能用来设计低通滤波器 偶对称 jIm[z] Re[z] 0 return 四个一组 复习 return return return return return return return return return return return return return return return return return return * (digital signal processing\chap7\linear_phase_FIR\magVSamp.m) Piecewise：分段地 线性相位 复 习 线性相位FIR： h(n)是实数 h(n)=±h(N-1-n) 线性相位FIR滤波器的零点位置 偶对称 奇对称 偶数 奇数 FIR分四种类型 P330表7-1 7.3 窗函数设计法 一、 FIR滤波器设计方法 如果希望得到的滤波器的理想频率响应为 那么FIR滤波器的设计就是寻找一个频率响应 逼近方法有三种： 去逼近 窗函数设计法(时域逼近) 频率抽样设计法(频域逼近) 最优化设计法(等波纹逼近) 窗函数设计法就是使所要设计的滤波器的单位冲激响应 逼近理想滤波器的单位冲激响应 一般来说，理想频率响应 是分段恒定的，在边界频率处有突变点，所以，得到的理想单位冲激响应 往往是无限长非因果序列。 但FIR的 是因果有限长的。 最简单的办法是直接截取一段 得到 这种截取可以想象为 是通过一个“窗口”所看到的一段 ，因此， 可表达为 和一个“窗函数”的乘积 问题是怎样用一个有限长的 去近似无限长的 设计思路： 在这里窗口函数就是矩形脉冲函数 ，以后还可看到，为了改善所设计滤波器的性能，窗函数还可以有其它的形式，相当于在矩形窗内对 作一定的加权处理。 理想的滤波器 IDTFT 以线性相位的理想低通滤波器为例： 图 以α为中心的偶对称无限长非因果序列。 可以看出：hd(n)偶对称 图 按照第一类线性相位条件 窗函数截断的影响 幅度函数 卷积图 加窗处理对理想频率响应的影响： 1、改变了理想频率响应的边沿特性，形成了过渡带，过渡带的宽度等于窗函数的主瓣宽度。 2、过渡带两侧产生肩峰和余振，即带内、带外产生起伏振荡，这些取决于窗函数的旁瓣，旁瓣多则振荡多，旁瓣相对值大则肩峰强度大，与N无关。 3、改变N，只能改变窗函数主瓣的宽度，窗函数频谱的绝对值大小和窗函数频谱的起伏密度，不能改变主瓣与旁瓣的比例关系，也就是不能改变肩峰的相对值。增大N，其最大肩峰总是8.95%，这种现象称为吉布斯效应。 改变窗函数的形状，可以改善滤波器的特性，窗函数有许多种，选取时希望能满足以下两点要求： 但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。窗函数不仅起截断作用，还能起平滑作用。 肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。 1、窗函数的主瓣要尽量窄，以获得较陡的过渡带； 2、相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就可以减小肩峰和余振，以提高阻带的衰减和通带的平稳性。 二、各种窗函数 矩形窗 最简单的窗函数；从阻带衰减特性来说，性能最差；主瓣宽度为 巴特列特窗（三角形窗） 主瓣宽度 幅度函数 窗函数 汉宁窗（升余弦窗） 主瓣宽度 幅度函数 窗函数 海明窗（改进的升余弦窗） 主瓣宽度 幅度函数 窗函数 布莱克曼窗（二阶升余弦窗） 主瓣宽度 幅度函数 窗函数 凯泽窗 是第一类变形零阶贝塞尔函数, 是一个可以自由选侧的参数，用于调整主瓣宽度和旁瓣的幅度。 矩形窗函数 接近汉明窗 接近布莱克曼 窗函数基本参数的比较(p342) p486 窗函数 窗性能指标 加窗后滤波器性能 旁瓣峰值 主瓣宽度 过渡带 宽度 阻带最小 衰减 Rectangular -13 dB -21 dB Bartlett -25 dB -25 dB Hanning -31 dB -44 dB Hamming -41 dB -53 dB Blackman -57 dB -74 dB 阻带衰减只与窗函数形状有关；过