几种时频分析方法综述2——希尔伯特黄变换.doc

几种时频分析方法综述2希尔伯特黄变换夏巨伟(浙江大学空间结构研究中心)摘 要：希尔伯特黄变换由经验模态分解（empirical mode decomposition，简称EMD）和Hilbert谱分析两部分组成。经验模态分解方法是一种自适应的、高效的数据分解方法。由于这种分解是以局部时间尺度为基础，因此，它适应于非线性、非平稳过程。通过经验模型分解，任何复杂的数据集都可以被分解为个数有限的、而且常常是为数不多的几个固有模函数(intrinsic mode functions，简称IMF)的线性叠加。通过分解得到IMF后，就可以对每一个分量做希尔伯特变换，得到其瞬时频率和幅度。本文详细对Hilbert-Huang Transform的过程进行了阐述，并用算例分析指出了其优势所在。关键词：希尔伯特黄变换； 时频分析技术；1 希尔伯特黄变换（Hilbert-Huang Transform）1.1 希尔伯特变换与瞬时频率（Hilbert Transform and instantaneous frequency）对于任意一个时间序列X(t)，它的希尔伯特变换具有如下形式：其中，P积分的柯西主值；希尔伯特变换对于任何属于Lp 空间中的函数都存立，即上式中X(t)Lp（，+）。通过上述定义，X(t)和Y(t)成为一组复共轭对，同时能够构造一个实部和虚部分为X(t)和Y(t)的解析信号(Analytic Signal)Z(t)，Z(t)表示为：其中，理论上讲有无数种方式去定义虚部，但是希尔伯特变换是唯一能够得到解析信号结果的方法。X(t)的Hilbert变换实质上是将X(t)与函数1/t在时域上做卷积，这就决定了通过X(t)的Hilbert变换能够考察其局部特性。得到X(t)的瞬时相位函数后，其瞬时频率为：1.2 经验模态分解与固有模态函数（Empirical mode decomposition/EMD and Intrinsic mode function/IMF）固有模态函数需要满足两个条件：（1）极值与零点的数量必须相等或最多相差一个；（2）由局部极大值包络和局部极小值包络定义的平均包络曲线上任何一点的值为0；1.2.1 EMD筛选过程（Sifting process）1.3 Hilbert谱与Hilbert边际谱经过筛选过程后，X(t)可以表示为IMF与残差量的和：对X(t)的每一个IMF进行Hilbert变换可以得到X(t)的Hilbert谱：得到Hilbert谱后可以进一步定义Hilbert边际谱：1.4 算例分析1.4.1 一个有跳变的余弦信号图1:跳变信号及其分量图2:跳变信号EMD分量的瞬时相位与频率1.4.2 频率发生改变的余弦信号图3:频率改变余弦信号及其EMD分解分量图4:频率改变余弦信号IMF分量瞬时相位与瞬时频率1.4.3 余弦扫频信号图5：余弦扫频信号及其EMD分解分量图6：余弦扫频信号IMF分量瞬时相位与瞬时频率1.4.4 两个不同频率的正弦信号的叠加图7：两个不同频率叠加的正弦信号及其