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

新疆师范大学物理与电子工程学院 FIR数字滤波器最大的优点是容易设计成线性相位。在数字信号传输与处理及图像信号处理中，要求系统具有线性相位特性。但由于FIR滤波器为全零点系统，所以对同一幅频特性要求，用FIR滤波器实现要比用IIR滤波器实现阶数高较多。所以，一般在要求线性相位时，选用FIR滤波器。 本章先总结线性相位FIR数字滤波器的特点及条件，然后介绍了几种窗函数的设计方法及仿真实现方法。 本节主要介绍FIR滤波器具有线性相位的条件及幅度特性以及零点、网络结构的特点。 1. 线性相位 对于长度为N的h(n)，传输函数为 H(ejω)线性相位是指θ(ω)是ω的线性函数， 即： 如果θ(ω)满足： 以上两种情况都满足群时延是一个常数， 即: ,τ为常数。 h(n)是实序列且对(N-1)/2偶对称，即h(n)=h(N-n-1) 当h(n)=h(N-n-1) 时 h(n)是实序列且对(N-1)/2奇对称，即h(n)=-h(N-n-1) 满足第二类线性相位的系统，在ω=0处有一个-π/2的截距，即信号通过系统时不仅有一个τ的延时，还有一个π/2的相移，这种网络称为正交变换网络（如希尔伯特变换）。 3.线性相位FIR滤波器幅度特性Hg(ω)的特点 分以下四种情况讨论其幅度特性的特点： sin项对ω=0、π、2π皆为奇对称，且在这三点sinωn=0，因此幅度特性Hg(ω) 对ω=0、π、2π奇对称，且在这三点Hg(ω) =0。 仅能够用来设计BP滤波器。 （4）h(n)=-h(N-n-1) ，N=偶数 由于N为偶数，则序列正好分成对称的两组： 结论：了解了线性相位FIR滤波器的各种特性，便可根据实际需要选择合适的FIR滤波器，同时设计时要遵循有关约束条件。 （1）第3、4种情况，对于任何频率都有一固定的相移，一般微分器及90°相移器采用这两种情况，而选频性滤波器则用第1、2种情况。 （2）仅从幅度考虑： LPDF：第1、2种情况均可； HPDF：第1、4种情况均可； BSDF：仅可选第1种情况； BPDF：4种情况均可； （1）设希望设计的滤波器传输函数为Hd(ejω) hd(n)是与其对应的单位脉冲响应，因此 相应的单位脉冲响应hd(n)为： 2、设计步骤 (1)先由Hd(ejw)求傅里叶反变换hd(n)。 (2)砍头去尾。因为要设计FIR滤波器h(n)是： a、具有因果性的，t0时，为0，砍头 b、由于要求所设计滤波器是线性相位，所以要求其偶/奇对称。由于砍头，所以必须去尾，让它们中心对称。即用有限长的h(n)去逼近无限长的hd(n)。 (3)利用卷积过程。即h(n)=W(n)? hd(n) 因此，窗函数序列的形状及长度的选择是设计关键。 3、频率特性 Hd(ejω)是一个以2π为周期的函数，可以展为傅氏级数，即 7.2.2 几种常用的窗函数 （1）矩形窗(Rectangle Window)： 频率响应为： 幅度函数为： 主瓣宽度： 第一副瓣比主瓣低13dB；阻带最小衰减－21dB 一般I0(x)取15～25项，便可以满足精度要求。 ?参数可以控制窗的形状，典型数据为4?9。 ?增大，主瓣加宽，旁瓣幅度减小。 当?=5.44时，窗函数接近哈明窗。 ?=7.865时，窗函数接近布莱克曼窗。 凯塞窗的幅度函数为： MATLAB信号处理工具箱提供了14种窗函数的产生函数，下面列出6种及其调用格式。 7.2.3 用窗函数法设计FIR滤波器的步骤 (1)根据技术要求确定待求滤波器的单位取样响应hd(n)。 如果给出待求滤波器的频响为Hd(ejω)，那么单位脉冲响应为： 若Hd(ejω)较复杂，或不能用封闭公式表示时，可在[0，2π]上采样M个点，得： (2)根据对过渡带（Δω）及阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。 过渡带Δω近似等于窗函数主瓣宽度。 选窗原则：在保证阻带衰减的要求下，选择主瓣宽度较小的。 (3) 计算滤波器的单位取样响应h(n)=hd(n)w(n) 若要求线性相位，则hd(n)与窗函数w(n)均应关于(N-1)/2点对称。 前面所提到的几种窗函数都关于(N-1)/2偶对称。 (4) 计算： 若不满足设计要求，则重复2、3、4步。 【例7.2.1】 用矩形窗、汉宁窗和布莱克曼窗设计FIR低通滤波器，设N=11，ωc=0.2πrad。 解：用理想低通作为逼近滤波器，有 延时： 则： 调用工具箱fir1实现窗函数法设计步骤。 实现线性相位FIR数字滤波器的标准窗函数法设计，“标准”是指在设计低通、高通、带通和带阻FIR滤波