第07章信号的自适应分解方法要点.ppt

机械工程与应用电子技术学院 北京市先进制造技术重点实验室 Key Laboratory of Advanced Manufacturing Technology 1 胥永刚/张建宇 现代测试信号分析与处理 (Advanced Signal Analysis and Processing) 2 期末大作业 1、时间：2015.12.23之前，由张志新收齐交给胥老师。 2、内容：选用一种或多种现代信号处理方法，完成数据分析，以论文形式提交（不少于3页）。信号类型不限定，由自己课题组提供。 3、格式：按照《机械工程学报》的模板。 3 讲授提纲 调制信号的解调分析方法 5 非平稳信号的时频分析 6 信号的自适应分解方法 7 本章内容： 课程回顾 何为EMD？ EMD的实现过程 Hilbert-Huang变换 EMD方法的主要问题 EEMD的提出 EMD的其他衍生算法 非平稳信号的自适应分解 4 §7.1 课程回顾 5 §7.1 课程回顾 6 §7.1 课程回顾 7 8 §7.2 何为EMD？ 1937年12月生于湖北，1960年毕业于台湾大学土木工程系，1967年获约翰霍普金斯大学流体力学博士学位。1975年起进入太空总署工作，是该署海洋科学首席科学家。在美国NASA工作超过三十年，为NASA 海洋科学首席科学家。曾当选中央研究院院士（2004）、NASA年度发明家（2003），以及美国国家工程研究院院士（2000）。 (1) EMD的基本思想 HHT发明人：黄锷 为了研究瞬态与非平稳现象，频率与能量必须是时间的函数，因此我们需要瞬时频率与能量的定义。信号的瞬时能量与瞬时包络的概念被广泛接受，而瞬时频率的概念在Hilbert变换方法产生之前，却一直具有争议性。大多数观点认为其不存在或认为只在特定的条件下存在，比如单分量信号。 接受瞬时频率这一概念主要有两个基本的困难： (1)首先是受到了傅立叶分析的根深蒂固的影响； (2)其次是定义瞬时频率的方法不统一。 9 §7.2 何为EMD？ (2) 瞬时频率的概念 对任意的时间序列x(t)，可得到它的Hilbert变换y(t)为： 构造解析函数 其中幅值函数 相位函数 相位函数的导数即为瞬时频率 10 §7.2 何为EMD？ (2) 瞬时频率的概念 然而在某些情况下，按上述定义求解的瞬时频率可能会出现没有意义的负频率。 11 (2) 瞬时频率的概念 §7.2 何为EMD？ 考虑如下信号 为了简单起见，假设信号幅值A1和A2是恒定的，ω1和ω2是正的。信号的频谱应由两个在ω1和ω2的δ函数组成，即 12 (2) 瞬时频率的概念 §7.2 何为EMD？ 这个信号是解析的，按前式可以求解其相位和幅值，得到 取相位的导数，得到其瞬时频率，有 可以发现，按照前述公式计算得到的瞬时频率，不仅和信号原有频率相关，还和信号中不同频率信号的幅值有关。 13 (2) 瞬时频率的概念 §7.2 何为EMD？ A1=0.2 A2=1 瞬时频率不仅是连续的，而且还出现了负值，而我们已知信号的频率是离散的和正的。可见，对任一信号做简单的Hilbert变换可能会出现无法解释的、缺乏实际物理意义的频率成分。 A1=-1.2 A2=1 14 (2) 瞬时频率的概念 §7.2 何为EMD？ 基本模式分量是为了得到有意义的瞬时频率而提出的。基本模式分量需要满足的两个条件为： (b) 在任一时间点ti 上，信号局部极大值确定的上包络线fmax(t)和局部极小值确定的下包络线fmin(t)的均值为零。即 15 (3) 基本模式分量的定义 §7.2 何为EMD？ 16 (3) 基本模式分量的定义 §7.2 何为EMD？ 典型的基本模式分量（IMF） 满足以上两个条件的基本模式分量，只包括一个基本模式的振荡，没有复杂的叠加波存在。需要注意的是，如此定义的基本模式分量并不被限定为窄带信号，可以是具有一定带宽的非平稳信号，例如纯粹的频率和幅度调制函数。 17 (3) 基本模式分量的定义 §7.2 何为EMD？ 对满足基本模式分量两个限定条件的信号可以通过Hilbert变换求出其瞬时频率。 但不幸的是，大多数信号或数据并不是基本模式分量，任何时刻，信号中可能包括多个振荡模式，这就是为什么简单的Hilbert变换不能给出一般信号的完全的频率内容的原因。 我们必须把复杂的非平稳信号按一定的规则提取出所包含的基本模式分量。 基于此，Norden E. Huang等人创造