频谱分析中如何选择合适的窗函数.pdf

频谱分析中如何选择合适的窗函数 本文章转自扬州晶明 加窗是为了减小泄漏 ! 1、信号截断及能量泄漏效应 数字信号处理的主要数学工具是傅里叶变换。应注意到，傅里叶变 换是研究整个时间域和频率域的关系。然而，当运用计算机实现工程测试 信号处理时，不可能对无限长的信号进行测量和运算，而是取其有限的时 间片段进行分析。做法是从信号中截取一个时间片段，然后用观察的信号 时间片段进行周期延拓处理，得到虚拟的无限长的信号，然后就可以对信 号进行傅里叶变换、相关分析等数学处理。 周期延拓后的信号与真实信号是不同的，下面从数学的角度来看这 种处理带来的误差情况。设有余弦信号 x （t ）在时域分布为无限长（ - ∞， ∞），将截断信号的谱 XT( ω)与原始信号的谱 X （ω）相比，它已不是原来 的两条谱线，而是两段振荡的连续谱。这表明原来的信号被截断以后，其 频谱发生了畸变， 原来集中在 f0 处的能量被分散到两个较宽的频带中去了， 这种现象称之为频谱能量泄漏（ Leakage ）。 信号截断以后产生的能量泄漏现象是必然的， 因为窗函数 w(t) 是一个 频带无限的函数，所以即使原信号 x(t) 是限带宽信号，而在截断以后也必 然成为无限带宽的函数，即信号在频域的能量与分布被扩展了。又从采样 定理可知，无论采样频率多高，只要信号一经截断，就不可避免地引起混 叠，因此信号截断必然导致一些误差，这是信号分析中不容忽视的问题。 如果增大截断长度 T ，即矩形窗口加宽，则窗谱 W （ω）将被压缩变 窄（ π／T 减小）。虽然理论上讲，其频谱范围仍为无限宽，但实际上中 心频率以外的频率分量衰减较快，因而泄漏误差将减小。当窗口宽度 T 趋 于无穷大时，则谱窗 W （ω）将变为 δ（ω）函数，而 δ（ω）与 X （ω） 的卷积仍为 H （ω），这说明，如果窗口无限宽，即不截断，就不存在泄 漏误差。 为了减少频谱能量泄漏，可采用不同的截取函数对信号进行截断， 截断函数称为窗函数，简称为窗。泄漏与窗函数频谱的两侧旁瓣有关，如 果两侧 p 旁瓣的高度趋于零，而使能量相对集中在主瓣，就可以较为接近 于真实的频谱，为此，在时间域中可采用不同的窗函数来截断信号。 2 、 常用窗函数 实际应用的窗函数，可分为以下主要类型： 幂窗：采用时间变量某种幂次的函数，如矩形、三角形、梯形或其它时间 函数 x(t) 的高次幂； 三角函数窗：应用三角函数，即正弦或余弦函数等组合成复合函数，例如 汉宁窗、海明窗等； 指数窗。：采用指数时间函数，如 e-st 形式，例如高斯窗等。 下面介绍几种常用窗函数的性质和特点。 (l) 矩形窗 矩形窗使用最多，习惯上不加窗就是使信号通过了矩形窗。这种窗的优点 是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了高频 干扰和泄漏，甚至出现负谱现象。 (2) 三角窗 三角窗亦称费杰（ Fejer ）窗，是幂窗的一次方形式， 三角窗与矩形窗比较，主瓣宽约等于矩形窗的两倍，但旁瓣小，而且无负 旁瓣 (3) 汉宁窗 汉宁（ Hanning ）窗又称升余弦窗，汉宁窗可以看作是 3 个矩形时间窗的 频谱之和，它可以使用旁瓣互相抵消，消去高频干扰和漏能。 汉宁窗与矩形窗的谱图对比，可以看出，汉宁窗主瓣加宽（第一个零点在 2 π/T 处）并降低，旁瓣则显著减小。第一个旁瓣衰减一 32dB ，而矩形窗