7有限长单位冲激响应FIR平安滤波器的设计方法.pptVIP

6. 凯塞（Kaiser）窗 这是一种适应性较强的窗，其窗函数的表示式为 0≤n≤N-1 式中，I0(x)是第一类变形零阶贝塞尔函数，β是一个可自由选择的参数。? 图7-13 凯塞窗函数 表7-2 凯塞窗的性能 表7-3 六种窗函数基本参数的比较 下面将窗函数法的设计步骤归纳如下： （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，近似就足够了。 （3） 由过渡带宽及阻带最小衰减的要求，可选定窗形状， 并估计窗口长度N。设待求滤波器的过渡带用Δω表示，它近似等于窗函数主瓣宽度。因过渡带Δω近似与窗口长度成反比， N≈A/Δω，A决定于窗口形式。例如，矩形窗A=4π，海明窗A=8π等，A参数选择参考表7-3。按照过渡带及阻带衰减情况，选择窗函数形式。原则是在保证阻带衰减满足要求的情况下， 尽量选择主瓣窄的窗函数。 （4） 求得所设计的FIR滤波器的单位脉冲响应。 h(n)=hd(n)w(n) 0≤n≤N-1 （5）由h(n)求FIR滤波器的系统函数H(z) 或H (ejω)=DTFT[h(n)] 检查是否满足设计要求。 通常整个设计过程可利用计算机编程来实现，可多选择几种窗函数来试探，从而设计出性能良好的FIR滤波器。 例7-1 ：根据下列技术指标，设计一个FIR低通滤波器。  通带截止频率ωp=0.2π，通带允许波动δ1=0.25dB；  阻带截止频率ωst=0.3π，阻带衰减δ2=50dB。  解 查表7-3可知，海明窗和布拉克曼窗均可提供大于 50 dB的衰减。但海明窗具有较小的过渡带从而具有较小的长度N。 根据题意，所要设计的滤波器的过渡带为 由表7-3可知，利用海明窗设计的滤波器的过渡带宽Δω=8π/N，所以低通滤波器单位脉冲响应的长度为 3 dB通带截止频率为 由式（7-39）可知，理想低通滤波器的单位脉冲响应为 海明窗为 则所设计的滤波器的单位脉冲响应为 N=80 所设计的滤波器的频率响应为 利用计算机编程实现，结果如下图所示。 (a)是理想低通滤波器的单位脉冲响应hd(n); (b)是海明窗函数; (c) 是实际低通滤波器的单位脉冲响应h(n); (d)是实际低通滤波器的幅频特性|H(ejω)|，以dB为单位。滤波器长N=80，实际阻带衰减为As=53dB，通带波动为Ap=0.0316 dB，均满足设计要求。 例7-A中低通滤波器设计结果 图 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滤波器频