Hilbert希尔伯特环变换

1、Hilbert 希尔伯特环变换黄锷院士在 On Holo-Hilbert spectral analysis: a full informational spectral representation for nonlinear and non-stationary data 中提出一种高维全 息谱分析理论 HHSA(Holo-Hilbert spectral analysis)要理解 HHSA 方法，首先要了解希尔伯特变换、 经验模态分解 (EMD) 、 与希尔伯特 -黄变换 (HHT) 。学术背景：在信号处理与频谱分析的目的是要描述信号的频谱含量在时间 上变化，以便能在时间和频谱上同时表示

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

3、近似应用于窄带信号。 但实际应用中，存在许多非窄带信号， 希尔伯特变换对这些信号无能为力。 即便是 窄带信号，如果不能完全满足希尔伯特变换条件， 也会使结果发生错误。而实际信号中由于噪声的存在，会使很多原来满足 希尔伯特变换条件的信号无法完全满足;(2) 对于任意给定时刻，通过希尔伯特变换运算后的结果只能在一个频率值，即只能处理任何时刻为单一频率的信号;(3) 对于一个非平稳的数据序列，希尔伯特变换得到的结果很大程度上失去了原有的物理意义。hrequeiicy (Hz)图1傅立叶、小波与希尔伯特-黄变换对瞬时频率的分辨率 希尔伯特-黄变换:针对上述的三个问题，黄铐院士在 1998年提出希尔伯特

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

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

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

7、部极值点和过零点的数目必须相等，或最多相差一个;在任意时刻点，局部最大值的包络（上包络线）和局部最小值的包络（下包络线） 平均必须为零。EMD分解为IMF过程:210203040 SO m 70'809010011 Wi 120第二步，找出信号中的局部最最小值，并使用三次样条拟合成一条包第三步，局部最大值包络线减去局部最小值包络线，如下图所示:1-2010203040500708090100110120第四步，原始信号减去均值包络线即使第一个分量IMF1，如下图所IMF，如下图所示:波与噪声信号乘积所产生的信号，虽然可以通过HHT变换识别出来,但不清楚具体细节。信号加上高斯噪声，与信号

8、乘上高斯噪声干 ill.;-ee-q J -正-I用土 ： L- - nil “ -utip 一 s-YFr- EF- r - 1二- 覚二 - L- -PL - 二 I 左-匸： 1 hl .pf * *1 Jd片 1 fl !. w-. V- - 吃，；，1JrvdTiencv iH/kX- t-tiwoTJefitv iH/l两种方式噪声引入信号后，所对应的频谱，其中绿色虚线为原正弦波 频谱，蓝色实线为噪声频谱，红色为信号频谱。可以发现正弦波乘上 噪声后的频谱中完全没有原信号的频率成分。而当前的HHT中可以观察到一部分的正弦波成分，但却缺失正弦波的关键信息。UmePHJi .4 昭鱼Mi

9、ne 如！； 罚臾弭出加；1 m舸：视 ！q恥 a & 仇 a a在这里黄铐认为，虽然在整个信号上的频率分布不清楚， 但在EMD分解后的IMF中，能量的分布与原信号相一致，即在IMF中的幅值是即关于频率也关于时间的函数。即黄铐最新提出的HHSA方法是一个同是体现时域-频域-能量域（幅值）的方法，如下图所示drhtaxi卜血/ejEu； ；inri j I'lcd? fugf Illi, 111 il)rr,L li iiyM/J /WHl-DiT-AM-卜M HilbertPM Hilhert打创j，"忤,!厂，I' ”，I叫IF臺£羣壬 臺枣暮至&

10、#163;g吓Uul Hi-Di HilhcMjf? '；' i" ill I ilcMPiiiKlf" MiliiiiL<,ir u'er/bCfiiiHiri'H主要应用 哈佛医学院用HHT来测量心率不整 约翰霍普金斯公共卫生学院用它来测量登革热的扩散 海军用它来探测潜艇 地震工程、地球物理探测、卫星资料分析 潜艇设计、结构损害侦测 潮汐、波浪场等各项研究 在脑机上的应用目前HHT作为经典的非线性系统分析方法，在脑机上早有大量 应用。上述的分析表明HHT与HHAS均是对信号进行 经验模态分解的 基础上进行的。而经验模态分解一个很重要

11、的特点是依靠信号自身的 特性分解为若干个IMF，通过对这些分解后的IMF分析，进而判断 信号的特性。因此从HHT与HHAS的应用场合上可以看出，其分析 的信号多是能反复出现的、有规律的、非线性的、非稳态的信号。在脑机上，HHT与H HAS的典型应用是针对脑信号经常稳定在某一状态上的场景, 大脑在这些状态下的脑电信号是可以稳定复现 的，即使是脑电信号本身非稳态、非线性，但EMD分解出的信号是 类似的。例如长期失眠、 抑郁状态、或高度睡眠、 注意力高度集中下， 其分解的信号是比较“纯净”的, 可以作为参考的。 而人脑在一般的 思维过程中产生的信号是很难分解出与有具体含义的成分的。 这是有EMD算法本身所决定的，因为EMD依靠信号本身特性分解， 没有 一个共同的“基”，因此一旦脑电信号解析出来的IMF没有对应参 考，就无法判断大脑状态。因此HHT与H HAS在神经学上检测睡眠、失眠、抑郁，即其他大脑长期处于某种状态的疾病，具有较好的 应用。此外课题组曾在2015年的暑假期间，使用HHT算法对14类 的运动想象信号进行识别分类，在离线数据集上达到7