Hilbert希尔伯特环变换.doc

《Hilbert希尔伯特环变换.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、.黄锷院士在《OnHolo-Hilbertspectralanalysis:afullinformationalspectralrepresentationfornonlinearandnon-stationarydata》中提出一种高维全息谱分析理论HHSA(Holo-Hilbertspectralanalysis)要理解HHSA方法，首先要了解希尔伯特变换、经验模态分解(EMD)、与希尔伯特-黄变换(HHT)。学术背景：在信号处理与频谱分析的目的是要描述信号的频谱含量在时间上变化，以便能在时间和频谱上同时表示信号的能量或者强度。傅里叶频谱并没有告诉我们哪些频率在什么时候

2、出现。因此傅里叶变换无法表现信号频率成分的时变性，因此学术界先后发展出了短时傅里叶变换、窗口傅里叶变换、小波等手段，近似的求信号某一时刻的瞬时频率。希尔伯特变换：希尔伯特变换是以著名数学家大卫·希尔伯特(DavidHilbert)来命名。通过希尔伯特变换，使得我们对短信号和复杂信号的瞬时参数的定义及计算成为可能，能够实现真正意义上的瞬时频率的提取，因而希尔伯特变换在信号处理上具有十分重要的地位，使得希尔伯特变换具有广泛的工程应用。但在进一步的工程应用中，希尔伯特变换具有以下缺陷：(1)希尔伯特变换只能近似应用于窄带信号。但实际应用中，存在许多非窄带信号，希尔伯特变换对这些信

3、号无能为力。即便是窄带信号，如果不能完全满足希尔伯特变换条件Word资料.，也会使结果发生错误。而实际信号中由于噪声的存在，会使很多原来满足希尔伯特变换条件的信号无法完全满足；(1)对于任意给定时刻，通过希尔伯特变换运算后的结果只能在一个频率值，即只能处理任何时刻为单一频率的信号；(2)对于一个非平稳的数据序列，希尔伯特变换得到的结果很大程度上失去了原有的物理意义。图1傅立叶、小波与希尔伯特-黄变换对瞬时频率的分辨率希尔伯特-黄变换：针对上述的三个问题，黄锷院士在1998年提出希尔伯特-黄变换(HHT)。其基本思想是：讲一个非稳态、非线性的信号分解为若干个稳态信号，在对分解

4、后的信号进行希尔伯特变换，分别求取对应的瞬时频率。Word资料.在这里将非稳态、非线性信号分解为多个稳态信号的算法成为经验模态分解(EMD)，EMD算法是希尔伯特-黄变换的核心，也是其能处理非稳态信号与非线性信号的关键。经验模态分解(EMD)EMD方法被认为是2000年来以傅立叶变换为基础的线性和稳态频谱分析的一个重大突破，该方法是依据数据自身的时间尺度特征来进行信号分解，无须预先设定任何基函数。这一点与建立在先验性的谐波基函数和小波基函数上的傅里叶分解与小波分解方法具有本质性的差别。正是由于这样的特点，EMD方法在理论上可以应用于任何类型的信号的分解，因而在处理非平稳及非

5、线性数据上，具有非常明显的优势,适合于分析非线性、非平稳信号序列，具有很高的信噪比。所以，EMD方法一经提出就在不同的工程领域得到了迅速有效的应用，例如用在海洋、大气、天体观测资料与地震记录分析、机械故障诊断、动力系统的阻尼识别以及大型土木工程结构的模态参数识别等方面。EMD能使复杂信号分解为有限个本征模函数（IntrinsicModeFunction，IMF），所分解出来的各IMF分量包含了原信号的不同时间尺度的局部特征信号。然后进行希尔伯特变换获得时频谱图，得到有物理意义的频率。本征模函数(IMF）在物理上，如果瞬时频率有意义，那么函数必须是对称的，局部均值为零，并且具

6、有相同的过零点和极值点数目。在此基础上，Word资料.NordneE.Huang等人提出了本征模函数（IntrinsicModeFunction，简称IMF）的概念。本征模函数任意一点的瞬时频率都是有意义的。Huang等人认为任何信号都是由若干本征模函数组成，任何时候，一个信号都可以包含若干个本征模函数，如果本征模函数之间相互重叠，便形成复合信号。EMD分解的目的就是为了获取本征模函数，然后再对各本征模函数进行希尔伯特变换，得到希尔伯特谱。Huang认为，一个本征模函数必须满足以下两个条件：(1)函数在整个时间范围内，局部极值点和过零点的数目必须相等，或最多相差一个；(2)

7、在任意时刻点，局部最大值的包络（上包络线）和局部最小值的包络（下包络线）平均必须为零。EMD分解为IMF过程：例如：某一信号如下所示：第一步，找出信号中的局部最大值，并使用三次样条拟合成一条包络线，如下图所示：Word资料.第二步，找出信号中的局部最最小值，并使用三次样条拟合成一条包络线，如下图所示：第三步，局部最大值包络线减去局部最小值包络线，如下图所示：第四步，原始信号减去均值包络线即使第一个分量IMF1，如下图所示Word资料.对IMF1进行同样操作可以获得第二个IMF2，如此反复，获得所有IMF，如下图所示