HHT变换讲义EMD讲解.pdfVIP

HHT变换讲义 1.1简介 传统的信号处理方法，如傅立叶分析是一种纯频域的分析方法。它用频率不       同的各复正弦分量的叠加来拟合原函数，也即用F  来分辨f  。而F  在有   限频域上的信息不足以确定在任意小范围内的函数f  ，特别是非平稳信号在 时间轴上的任何突变，其频谱将散布在整个频率轴上。而且，非平稳动态信号的 统计特性与时间有关，对非平稳信号的处理需要进行时频分析，希望得到时域和 频域中非平稳信号的全貌和局域化结果。在傅立叶变换中，人们若想得到信号的 时域信息，就得不到频域信息。反之亦然。后来出现的小波（Wavelet ）变换通 过一种可伸缩和平移的小波对信号变换，从而达到时频局域化分析的目的。但这 种变换实际上没有完全摆脱傅立叶变换的局限，它是一种窗口可调的傅立叶变 换，其窗内的信号必须是平稳的。另外，小波变换是非适应性的，小波基一旦选 定，在整个信号分析过程中就只能使用这一个小波基了。 HHT （Hilbert-Huang Transform ）技术是(1998年由NASA 的Norden E Huang 等提出的新的信号处理方法。该方法适用于非线性非平稳的信号分析，被认为是 近年来对以傅立叶变换为基础的线性和稳态谱分析的一个重大突破。目前HHT 技术已用于地球物理学和生物医学等领域的研究，并取得了较好的结果。 存在的问题 尽管HHT技术在处理非线性、非稳态信号方面有很大的优势，但是这个方法 本身还是有许多的问题有待进一步研究。 正如Huang 在文章中指出的那样，对于这种新的信号处理方法，其基的完 备性还需要严密的证明。另外，在做Hilbert变换时出现的边界效应也需要更好的 方法来解决。但是，HHT技术中最严重，也是现今研究的最多的是EMD 分解中 的包络过程。从对EMD分解方法的介绍可以看出，包络线的构造影响着整个分 解的结果，也决定了后面的Hilbert变换。Huang 采用的三次样条插值来拟和包络 线。在实际应用中，发现这样做会产生严重的边界效应，污染了原始数据。特别 是对短数据而言，这种影响使分析所得的结果失去了原有的意义。对此，Huang 等采用的是根据信号端点处的振幅和频率，分别增加两组特征波的方式进行数据 延拓。Huang的这种延拓方法已经向NASA申请了专利。除此之外，还有人提出 了其它方法进行端点延拓。比如国家海洋局的黄大吉等提出的镜像闭合延拓法和 极值延拓法。其中镜像闭合延拓法是根据信号的分布特性，把镜子放在具有对称 性的极值位置，通过镜像法把镜内信号映射成一个周期性的环形信号，不存在端 点，从根本上避免了端点效应。而极值延拓法具有镜像闭合延拓法相当的效果， 而不增加信号序列本身的长度，计算较快，尤其在处理非对称波形信号时比较优 越。还有青岛海洋大学的邓拥军等，利用神经网络分析方法来延拓数据端。这些 方法对边界效应的抑制都有一定的作用，但是也都需要更深一步的研究。 1.2 瞬时频率 频率是个极其重要的物理量，定义为信号周期倒数，其物理含义显而易见。 对于正弦信号，它的频率为恒值。但是对于大部分信号，它的频率是随时间变化 的函数，故提出瞬时频率概念。瞬时频率即表征信号在局部时间点上瞬态频率特 性，整个持续期上的瞬时频率反映了信号频率的时变规律。 [39] 对于随机时间序列X(t)，对其进行Hilbert 变换 ，可以得到Y(t)如下： 1 X       Y t  PV( d ) (1-1)   t