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

图 7-12 理想低通滤波器加窗后的幅度响应（N51）, A20lg|Hω/H0| a 矩形窗; b 巴特利特窗（三角形窗）; c 汉宁窗; d 海明窗; e 布拉克曼窗 6. 凯塞（Kaiser）窗 这是一种适应性较强的窗，其窗函数的表示式为 0≤n≤N-1 式中，I0x是第一类变形零阶贝塞尔函数，β是一个可自由选择的参数。? 图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） 给定希望逼近的频率响应函数Hdejω。 （2） 利用式（7-36）求单位抽样响应hdnIDTFT[Hdejω] 。如果Hdejω很复杂或不能直接计算积分，则必须用求和代替积分，以便在计算机上计算，也就是要计算离散傅里叶反变换， 一般都采用FFT来计算。将积分限分成M段，也就是令采样频率为ωk2πk/M，k0, 1, 2, …, M-1，则有 频域的采样，造成时域序列的周期延拓，延拓周期是M， 即由于hdn有可能是无限长的序列，因此严格说，必须当M→∞时，hＭn才能等于hdn而不产生混叠现象，即。实际上，由于hdn随n的增加衰减很快，一般只要M足够大，即MN，近似就足够了。 三、线性相位FIR滤波器的零点位置线性相位FIR滤波器的系统函数为： Hz±z-N-1Hz-1因此，若zzi是Hz的零点，即Hzi0，则z1/zizi-1也一定是Hz的零点，Hzi-1±zi N-1 Hzi0）当hn是实数时，Hz的零点必成共轭对出现，所以 zzi*及zz*i-1也一定是Hz的零点，因而 线性相位FIR滤波器的零点必是互为倒数的共轭对。 这种互为倒数的共轭对有四种可能性： 图 7-5 线性相位FIR滤波器的零点位置图由幅度响应的讨论可知， 第二种类型的线性相位滤波器 Hπ0，因此必然有单根 z -1。 第四种类型的线性相位滤波器 H00,因此必然有单根 z1。 第三种类型的线性相位滤波器 H0Hπ）0，因此必然有两种单根 z±1 。 了解了线性相位FIR滤波器的特点，便可根据实际需要选择合适类型的FIR滤波器，同时设计时需遵循有关的约束条件。下面讨论线性相位FIR滤波器的设计方法时，都要用到这些特点。 去逼近,如果希望得到的滤波器的理想频率响应为:窗口设计法（时域逼近） 频率采样法（频域逼近） 最优化设计（等波纹逼近） 那么 FIR滤波器的设计就在于寻找一个传递函数 逼近方法有三种： 7.3 用窗函数法设计FIR滤波器 一、设计方法窗函数法是设计FIR数字滤波器最简单的方法。这种方法一般是先给定所要求的理想滤波器的频率响应，要求设计一个FIR滤波器频率响应 , 去逼近理想的频率响应。 因此，必须首先由理想频率响应的傅里叶反变换推导出对应的单位抽样响应:窗函数法设计FIR数字滤波器是在时域进行的, 从单位抽样响应序列着手，使hn逼近理想的单位抽样响应序列hdn。 （7-36）由于许多理想化的系统均用分段恒定的或分段函数表示的频率响应来定义，因此hdn一定是无限长的序列，且是非因果的。而我们要设计的是FIR滤波器，其hn必定是有限长的，所以要用有限长的hn来逼近无限长的hdn，最简单且最有效的方法是截断hd n 0≤n≤N-1 其他 式中如果采用简单截取，则窗函数为矩形窗。 0≤n≤N-1 其他 矩形窗通常，我们可以把hn表示为所需单位抽样响应与一个有限长的窗口函数序列wn的乘积，即 hnhdnwn 的波形如下图所示： 相应的单位抽样响应为： hdn是一个中心点在α的偶对称、无限长、非因果序列， 为了构造一个长度为N的线性相位滤波器，只有将hdn截取一段，并保证截取的一段对N-1/2 对称，故中心点a必须取aN-1/2 。 |ω|≤ωc ωc|ω|≤π 例如，要求设计一个线性相位FIR数字低通滤波器，假设理想低通滤波器的频率响应为： （7-39） 设截取的一段用hn表示，则 理想低通的单位抽样响应及矩形窗分析窗口函数法对频响产生的影响逼近程度根据复卷积定理，由可得 hn的频率特性为：Hejω能否逼近Hdejω取决于窗函数的频谱特性Wejω （7-42）