(数字信号处理)第6章有限长单位脉冲响应.ppt

第一类线性相位滤波器幅度函数H(ω)关于ω=0, π, 2π为偶对称，即 （6-60） 如果采样值H(k)=H(ej2πk/N)也用幅值Hk（纯标量）与相角θk表示， 即 （6-61） 并在ω=0～2π之间等间隔采样N点 k=0, 1, 2, …, N-1 将ω=ωk代入式（6-59）与式（6-60)中，并写成k的函数, 有： （6-62） （6-63） 由式（6-63）可知，Hk满足偶对称要求。 3. 汉宁（Hanning）窗 汉宁窗又称升余弦窗。 (6-42) 利用傅里叶变换特性，可得 (6-43) 当N1 时，N-1≈N, 所以窗函数的幅频函数为 （6-44） 这三部分之和，使旁瓣互相抵消，能量更集中在主瓣，它的最大旁瓣值比主瓣值约低31dB。但是代价是主瓣宽度比矩形窗的主瓣宽度增加一倍， 即为 8π/N。 4. 海明（Hamming）窗 海明窗又称改进的升余弦窗。 把升余弦窗加以改进， 可以得到旁瓣更小的效果， 窗形式为 (6-45) w(n)的频率响应的幅度特性为 (6-46) 与汉宁窗相比，主瓣宽度相同，为 8π/N，但旁瓣又被进一步压低， 结果可将99.963%的能量集中在窗谱的主瓣内，它的最大旁瓣值比主瓣值约低41dB。 5. 布拉克曼（Blackman）窗 布拉克曼窗又称二阶升余弦窗。 为了进一步抑制旁瓣，对升余弦窗函数再加上一个二次谐波的余弦分量， 变成布拉克曼窗，故又称二阶升余弦窗。 (6-47) w(n)的频率响应的幅度特性为 (6-48) 图 6-10 五种常用的窗函数 图 6-11 图 6-10 的各种窗函数的傅里叶变换（N=51），A=20 lg|W(ω)/W(θ)| (a) 矩形窗; (b) 巴特利特窗（三角形窗）; (c) 汉宁窗; (d) 海明窗; (e) 布拉克曼窗 图 6-12 理想低通滤波器加窗后的幅度响应（N=51）, A=20lg|H(ω)/H(0)| (a) 矩形窗; (b) 巴特利特窗（三角形窗）; (c) 汉宁窗; (d) 海明窗; (e) 布拉克曼窗 6. 凯塞（Kaiser）窗 这是一种适应性较强的窗，其窗函数的表示式为 0≤n≤N-1 （6-49） 式中，I0(x)是第一类变形零阶贝塞尔函数，β是一个可自由选择的参数。 ? 零阶贝塞尔函数的曲线如图6-13所示。 图 6-13 零阶贝塞尔函数 图6-14 凯塞窗函数 表6-2 凯塞窗的性能 虽然凯塞窗看上去没有初等函数的解析表达式。但是，在设计凯塞窗时，对零阶变形贝塞尔函数可采用无穷级数来表达 （6-50） 这个无穷级数，可用有限项级数去近似，项数多少由要求的精度来确定。 因而采用计算机是很容易求解的。 表6-3 六种窗函数基本参数的比较 下面将窗函数法的设计步骤归纳如下： （1） 给定希望逼近的频率响应函数Hd(ejω)。 （2） 根据式（6-24）求单位脉冲响应hd(n)。 如果Hd(ejω)很复杂或不能直接计算积分，则必须用求和代替积分，以便在计算机上计算，也就是要计算离散傅里叶反变换， 一般都采用FFT来计算。将积分限分成M段，也就是令采样频率为ωk=2πk/M，k=0, 1, 2, …, M-1，则有 （6-51） 频域的采样，造成时域序列的周期延拓，延拓周期是M， 即 （6-52） 由于hd(n)有可能是无限长的序列，因此严格说，必须当M→∞时，hＭ(n)才能等于hd(n)而不产生混叠现象，即 。实际上，由于hd(n)随n的增加衰减很快，一般只要M足够大，即MN，近似就足够了。 （3） 由过渡带宽及阻带最小衰减的要求，可选定窗形状， 并估计窗口长度N。设待求滤波器的过渡带用Δω表示，它近似等于窗函数主瓣宽度。因过渡带Δω近似与窗口长度成反比， N≈A/Δω，A决定于窗口形式。例如，矩形窗A=4π，海明窗A=8π等，A参数选择参考表6-3。按照过渡带及阻带衰减情况，选择窗函数形式。原则是在保证阻带衰减满足要求的情况下， 尽量选择主瓣窄的窗函数。 （4） 最后，计算所设计的FIR滤波器的单位脉冲响应。 h(n)=hd(n)w(n) 0≤n≤N-1 （5）由h(n)求FIR滤波器的系统函数H(z) 通常整个设计过程可利用计算机编程来实现，可多选择几种窗函数来试探，从而设计出性能良好的FIR滤波器。 例 6-6 用矩形窗设计一个线性相位带通滤波器 -ωc≤ω-ω