人民邮电数字信号处理.pptVIP

第五章 FIR数字滤波器的设计 5.1 线性相移FIR滤波器的特性 5.2 窗口法 5.3 频率取样法 5.4 FIR滤波器的优化设计 5.1 线性相移FIR滤波器的特性 1.线性相移FIR滤波器条件的推导 (1) 恒延时滤波 FIR数字滤波器的系统函数为 其频率响应为 其中 (2) 恒相延时与恒群延时同时成立 当恒相延时与恒群延时同时成立时，相频特性必须是线性的且过零点，即 θ(ω)=-τω 因此，由式(5-2)与式(5-3)得到 (3) 只要求恒群延时成立 当只要求恒群延时，即τg(ω)是一常数时，相频特性为 θ(ω)=θ0-τω θ0=±π/2 由式(5-2)与式(5-3)得到 2.线性相移FIR数字滤波器的频率响应 由以上讨论的对称条件，即式(5-7b)和式(5-10b)，可以导出线性相移FIR数字滤波器的特有形式，下边分4种不同情况来讨论。 (1) 偶对称，N为奇数 (2) 偶对称，N为偶数 (3) 奇对称，N为奇数 (4) 奇对称，N为偶数 3.线性相移FIR数字滤波器的零极点分布 FIR数字滤波器的系统函数即式(5-1)的表示，在z=0点有(N-1)个极点，有(N-1)个零点，而这(N-1)个零点分布的特点是互为倒数共轭对。 如果h(n)是偶对称，即h(n)=h(N-1-n)，则式(5-1)便可写为 5.2 窗口法 设计FIR滤波器最直接的方法就是把无限时宽冲激响应截短，得到有限长度的冲激响应。若Hd(ejω)是要求的理想频率响应，则 式中hd(n)是对应的冲激响应序列，即 1.基本思想 一般来讲，滤波器理想频率特性Hd(ejω)都是分段不连续的，因而，在时域内其单位取样响应hd(n)也是无限时宽的。 2.数学推导 (1) 非因果FIR低通滤波器 此时，理想低通滤波器的频率特性为 由傅氏变换关系有 (2) 因果FIR低通滤波器 理想低通滤波器的频率特性为 由傅氏变换关系有 由此可知，hd(n)是一个以a为中心的偶对称无限长序列。 通过以上讨论，得出结论，由于加了窗函数，将对理想特性产生以下3点影响。 (1) 使理想特性不连续的边沿加宽，形成一个过渡带，过渡带的宽度取决于窗口频谱主瓣的宽度。 (2) 在过渡带两旁产生肩峰和余振，它们取决于窗口频谱的旁瓣，旁瓣越多，余振也越多；旁瓣相对值越大，则肩峰越强。 (3) 增加选取长度N，只能缩小窗口频谱的主瓣宽度而不能改变旁瓣的相对值(但N太小则会影响旁瓣相对值)，旁瓣与主瓣的相对关系只决定于窗口函数的形状。因此，增加N只能减小过渡带宽度，而不能改变肩峰值。肩峰值的大小直接决定着通带内的平稳和阻带的衰减，与滤波器的性能有很大关系。例如矩形窗时肩峰达8.95%，致使阻带衰减只有21dB，这在工程上是不够的。要改善阻带的衰减特性只能改善窗口函数的形状。 改善窗口函数形状的标准是： (1) 尽量减少窗口频谱中的旁瓣，使能量尽量地集中在主瓣中，这样可以减少肩峰和余振，提高阻带的衰减； (2) 主瓣的宽度尽量窄，以获得较陡的过渡带。然而，这两个要求总不能兼得，往往需要用增加主瓣宽度来换取旁瓣的抑制。 3.几种常用窗函数 (1) 矩形窗 以上的讨论都是以矩形窗为例，现只将其表达式给出，即 (2) 升余弦(Hanning)窗一汉宁窗 表达式为 (3) 改进的升余弦(Hamming)窗一汉明窗 如果对余弦再作一点调整，可以得到旁瓣最小的效果，这样得到的窗函数为 (4) 二阶升余弦(Blackman)窗一布莱克曼窗 如果再要进一步抑制旁瓣的话，可以对升余弦再加一个二次谐波的余弦分量，这样得到的函数为 (5) 凯塞(Kaiser)窗 凯塞窗是利用零阶贝塞尔函数I0(x)构成的，即 4.Matlab实现FIR 窗函数设计法 在Matlab信号处理工具箱中，Matlab提供了几个子程序来实现上面的窗函数，同时还提供了两个基于窗函数法的FIR数字滤波器设计函数b=fir1(n，Wn，options)和b=fir2(n，f，m，options)。 5.3 频率取样法 从时域开始，以有限长单位取样响应h(n)近似理想的单位取样响应。然而，在第三章中曾经指出，一个有限时宽为N的序列其z变换可以用单位圆上的N个等间隔取样来表示，对于FIR滤波器来说式(3-46)意味着一个有限长序列，同样可以用N个频率取样值来惟一确定，其取样值表达z函数的内插公式为 这就为设计FIR滤波器的系统函数提供了另一途径，即直接从频域出发，对理想频响取样，以此来确定H(k)值，即令 这样可以使设计所得的系统函数H(z)去逼近理想的系统函数Hd(z)。至少在取样点的频率上，两者可以具有相同的频响，即使式(5-35)及