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

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

频谱分析中如何选择合适的窗函数
加窗是为了减小泄漏!1、信号截断及能量泄漏效应\x0d数字信号处理的主要数学工具是傅里叶变换.应注意到,傅里叶变换是研究整个时间域和频率域的关系.然而,当运用计算机实现工程测试信号处理时,不可能对无限长的信号进行测量和运算,而是取其有限的时间片段进行分析.做法是从信号中截取一个时间片段,然后用观察的信号时间片段进行周期延拓处理,得到虚拟的无限长的信号,然后就可以对信号进行傅里叶变换、相关分析等数学处理.\x0d周期延拓后的信号与真实信号是不同的,下面从数学的角度来看这种处理带来的误差情况.设有余弦信号x（t）在时域分布为无限长（- ∞,∞）,将截断信号的谱XT(ω)与原始信号的谱X（ω）相比,它已不是原来的两条谱线,而是两段振荡的连续谱.这表明原来的信号被截断以后,其频谱发生了畸变,原来集中在f0处的能量被分散到两个较宽的频带中去了,这种现象称之为频谱能量泄漏（Leakage）.\x0d信号截断以后产生的能量泄漏现象是必然的,因为窗函数w(t)是一个频带无限的函数,所以即使原信号x(t)是限带宽信号,而在截断以后也必然成为无限带宽的函数,即信号在频域的能量与分布被扩展了.又从采样定理可知,无论采样频率多高,只要信号一经截断,就不可避免地引起混叠,因此信号截断必然导致一些误差,这是信号分析中不容忽视的问题.\x0d如果增大截断长度T,即矩形窗口加宽,则窗谱W（ω）将被压缩变窄（π／T减小）.虽然理论上讲,其频谱范围仍为无限宽,但实际上中心频率以外的频率分量衰减较快,因而泄漏误差将减小.当窗口宽度T趋于无穷大时,则谱窗W（ω）将变为δ（ω）函数,而δ（ω）与X（ω）的卷积仍为H（ω）,这说明,如果窗口无限宽,即不截断,就不存在泄漏误差.\x0d为了减少频谱能量泄漏,可采用不同的截取函数对信号进行截断,截断函数称为窗函数,简称为窗.泄漏与窗函数频谱的两侧旁瓣有关,如果两侧p旁瓣的高度趋于零,而使能量相对集中在主瓣,就可以较为接近于真实的频谱,为此,在时间域中可采用不同的窗函数来截断信号.2、 常用窗函数\x0d实际应用的窗函数,可分为以下主要类型：\x0d幂窗：采用时间变量某种幂次的函数,如矩形、三角形、梯形或其它时间函数x(t)的高次幂；\x0d三角函数窗：应用三角函数,即正弦或余弦函数等组合成复合函数,例如汉宁窗、海明窗等；\x0d指数窗.：采用指数时间函数,如e-st形式,例如高斯窗等. 下面介绍几种常用窗函数的性质和特点.\x0d(l) 矩形窗矩形窗使用最多,习惯上不加窗就是使信号通过了矩形窗.这种窗的优点是主瓣比较集中,缺点是旁瓣较高,并有负旁瓣,导致变换中带进了高频干扰和泄漏,甚至出现负谱现象.\x0d(2) 三角窗三角窗亦称费杰（Fejer）窗,是幂窗的一次方形式,\x0d三角窗与矩形窗比较,主瓣宽约等于矩形窗的两倍,但旁瓣小,而且无负旁瓣(3) 汉宁窗汉宁（Hanning）窗又称升余弦窗,汉宁窗可以看作是3个矩形时间窗的频谱之和,它可以使用旁瓣互相抵消,消去高频干扰和漏能.\x0d汉宁窗与矩形窗的谱图对比,可以看出,汉宁窗主瓣加宽（第一个零点在2π/T处）并降低,旁瓣则显著减小.第一个旁瓣衰减一32dB,而矩形窗第一个旁瓣衰减-13dB.此外,汉宁窗的旁瓣衰减速度也较快,约为60dB/（10oct）,而矩形窗为20dB/（10oct）.由以上比较可知,从减小泄漏观点出发,汉宁窗优于矩形窗.但汉宁窗主瓣加宽,相当于分析带宽加宽,频率分辨力下降.\x0d(4) 海明窗海明（Hamming）窗也是余弦窗的一种,又称改进的升余弦窗,海明窗与汉宁窗都是余弦窗,只是加权系数不同.海明窗加权的系数能使旁瓣达到更小.分析表明,海明窗的第一旁瓣衰减为-42dB.海明窗的频谱也是由 3个矩形时窗的频谱合成,但其旁瓣衰减速度为20dB/（10oct）,这比汉宁窗衰减速度慢.海明窗与汉宁窗都是很有用的窗函数.\x0d(5) 高斯窗是一种指数窗,高斯窗谱无负的旁瓣,第一旁瓣衰减达一55dB.高斯窗谱的主瓣较宽,故而频率分辨力低.高斯窗函数常被用来截断一些非周期信号,如指数衰减信号等.\x0d除了以上几种常用窗函数以外,尚有多种窗函数,如平顶窗、帕仁（Parzen）窗、布拉克曼（Blackman）窗、凯塞（kaiser）窗等.\x0d对于窗函数的选择,应考虑被分析信号的性质与处理要求.如果仅要求精确读出主瓣频率,而不考虑幅值精度,则可选用主瓣宽度比较窄而便于分辨的矩形窗,例如测量物体的自振频率等；如果分析窄带信号,且有较强的干扰噪声,则应选用旁瓣幅度小的窗函数,如汉宁窗、三角窗等；对于随时间按指数衰减的函数,可采用指数窗来提高信噪比.3、窗函数选择指南如果在测试中可以保证不会有泄露的发生,则不需要用任何的窗函数（在软件中可选择uniform）.但是如同刚刚讨论的那样,这种情况只是发生在时间足够长的瞬态捕捉和一帧数据中正好包含信号整周期的情况.如果测试信号有多个频率分量,频谱表现的十分复杂,且测试的目的更多关注频率点而非能量的大小.在这种情况下,需要选择一个主畔够窄的窗函数,汉宁窗是一个很好的选择.如果测试的目的更多的关注某周期信号频率点的能量值,比如,更关心其EUpeak,EUpeak-peak,EUrms或者EUrms2,那么其幅度的准确性则更加的重要,可以选择一个主畔稍宽的窗,flattop窗在这样的情况下经常被使用.对冲击实验的数据进行分析时,因为在数据帧开始段的一些重要信息会被一般的窗函数所衰减,因此可以使用force/exponential窗.Force窗一移去了数据帧末端的噪声,对激励信号有用.而exponential窗则确保响应信号在末端的振动衰减为零值.