7有限长单位冲激响应FIR滤波器的设计方法试卷.ppt

图 7-12 理想低通滤波器加窗后的幅度响应（N=51）, A=20lg|H(ω)/H(0)| (a) 矩形窗; (b) 巴特利特窗（三角形窗）; (c) 汉宁窗; (d) 海明窗; (e) 布拉克曼窗 6. 凯塞（Kaiser）窗 这是一种适应性较强的窗，其窗函数的表示式为 0≤n≤N-1 式中，I0(x)是第一类变形零阶贝塞尔函数，β是一个可自由选择的参数。? 图7-13 凯塞窗函数 表7-2 凯塞窗的性能 表7-3 六种窗函数基本参数的比较 窗函数 窗谱性能指标 加窗后滤波器性能指标 旁瓣峰值 /dB 主瓣宽度 / 2π/N 过渡带宽 / 2π/N 阻带最小衰减 /dB 矩形窗 巴特列特 汉宁窗 海明窗 布拉克曼窗 凯泽窗 （β=7.865） -13 -25 -31 -41 -57 2 4 4 4 6 0.9 2.1 3.1 3.3 5.5 5 -21 -25 -44 -53 -74 -80 *最小阻带衰减只由窗形状决定，不受窗宽N的影响； 而过渡带的宽度既和窗形状有关，且随窗宽N的增加而减小。 下面将窗函数法的设计步骤归纳如下： （1） 给定希望逼近的频率响应函数Hd(ejω)。 （2） 利用式（7-36）求单位抽样响应hd(n)=IDTFT[Hd(ejω)] 。 如果Hd(ejω)很复杂或不能直接计算积分，则必须用求和代替积分，以便在计算机上计算，也就是要计算离散傅里叶反变换， 一般都采用FFT来计算。将积分限分成M段，也就是令采样频率为ωk=2πk/M，k=0, 1, 2, …, M-1，则有 频域的采样，造成时域序列的周期延拓，延拓周期是M， 即 由于hd(n)有可能是无限长的序列，因此严格说，必须当M→∞时，hＭ(n)才能等于hd(n)而不产生混叠现象，即 。实际上，由于hd(n)随n的增加衰减很快，一般只要M足够大，即MN，近似就足够了。 三、线性相位FIR滤波器的零点位置 线性相位FIR滤波器的系统函数为： H(z)=±z-(N-1)H(z-1) 因此，若z=zi是H(z)的零点，即H(zi)=0， 则z=1/zi=zi-1也一定是H(z)的零点，(H(zi-1)=±zi (N-1) H(zi)=0） 当h(n)是实数时，H(z)的零点必成共轭对出现， 所以 z=zi*及z=(z*i)-1也一定是H(z)的零点， 因而 线性相位FIR滤波器的零点必是互为倒数的共轭对。 这种互为倒数的共轭对有四种可能性： 图 7-5 线性相位FIR滤波器的零点位置图 由幅度响应的讨论可知， 第二种类型的线性相位滤波器 H(π)=0， 因此必然有单根 z= -1。 第四种类型的线性相位滤波器 H(0)=0, 因此必然有单根 z=1。 第三种类型的线性相位滤波器 H(0)=H(π）=0， 因此必然有两种单根 z=±1 。  了解了线性相位FIR滤波器的特点，便可根据实际需要选择合适类型的FIR滤波器，同时设计时需遵循有关的约束条件。下面讨论线性相位FIR滤波器的设计方法时，都要用到这些特点。 去逼近 , 如果希望得到的滤波器的理想频率响应为: 窗口设计法（时域逼近） 频率采样法（频域逼近） 最优化设计（等波纹逼近） 那么 FIR滤波器的设计就在于寻找一个传递函数 逼近方法有三种： 7.3 用窗函数法设计FIR滤波器 一、设计方法 窗函数法是设计FIR数字滤波器最简单的方法。这种方法一般是先给定所要求的理想滤波器的频率响应 ，要求设计一个FIR滤波器频率响应 , 去逼近理想的频率响应 。 因此，必须首先由理想频率响应 的傅里叶反变换推导出对应的单位抽样响应: 窗函数法设计FIR数字滤波器是在时域进行的, 从单位抽样响应序列着手，使h(n)逼近理想的单位抽样响应序列hd(n)。 （7-36） 由于许多理想化的系统均用分段恒定的或分段函数表示的频率响应来定义，因此hd(n)一定是无限长的序列，且是非因果的。而我们要设计的是FIR滤波器，其h(n)必定是有限长的，所以要用有限长的h(n)来逼近无限长的hd(n)，最简单且最有效的方法是截断hd (n) 0≤n≤N-1 其他