经验模态分解算法

经验模态分解摘要——黄提出了经验模态分解（EMD）的数据处理方法，也对这种技术应用的有效性进行了讨论。

许多变种算法（新的停止准则，即时版本的算法）也产生出来。

数值模拟用来作经验性的评估执行单元运用于语音识别和分离方面，得出的实验结果认为这种方法是根据自适应的常数Q的滤波器组提出的。

1．介绍近来，一种被称为EMD的新的非线性方法被黄等人提出，这种方法能够自适应的把非平稳信号分解成一系列零均值的AMFM信号(调频调幅) 的总和。

尽管这种方法经常有着显著的效果，但是这个方法在算法方面的定义是困难的，因此这种方法没有作为一种分析方法得到承认，一般一种分析方法是需要有理论分析和性能评估。

因此本文的目的是用实验的方式使得该算法更容易理解，并且提出了基于原算法的各种各样的改进的算法。

设置实验性能评估的许多初始条件是为了获取一种有效的分解并且使得该算法更容易理解。

2．EMD基础EMD的出发点是把信号内的震荡看作是局部的。

实际上，如果我们要看评估信号x(t)的2个相邻极值点之间的变化（2个极小值，分别在t-和t+处），我们需要定义一个（局部）高频成分{d(t),t-<=t<=t+}(局部细节),这个高频成分与震荡相对应，震荡在2个极小值之间并且通过了极大值（肯定出现在2极小值之间）。

为了完整这个图形，我们还需要定义一个（局部）低频成分m(t)（局部趋势），这样x(t)=m(t)+d(t),(t-<=t<=t+)。

对于整个信号的所有震动成分，如果我们能够找到合适的方法进行此类分解，这个过程可以应用于所有的局部趋势的残余成分，因此一个信号的构成成分能够通过迭代的方式被抽离出来。

对于一个给定的信号x(t),进行有效的EMD分解步骤如下：1）找出想x(t)的所有极值点2）用插值法对极小值点形成下包络emint(t)，对极大值形成上包络emax(t)3）计算均值m(t)=(emint(t)+emax(t))/24）抽离细节d(t)=x(t)-m(t)5）对残余的m(t)重复上诉步骤在实际中，上述过程需要通过一个筛选过程进行重定义，筛选过程的第一个迭代步骤是对细节信号d(t)重复从1-4步，直到d(t)的均值是0，或者满足某种停止准则才停止迭代。

一旦满足停止准则，此时的细节信号d(t)就被称为IMF，d(t)对应残量信号用第5步计算。

通过以上过程，极值点的数量伴随着残量信号的产生而越来越少，整个分解过程会产生有限个模函数（IMF）。

模函数和残量信号可以进行谱分析，但是这个谱分析不能从狭隘的角度来看。

首先,需要强调一下，即使是谐振荡，应用上述方法产生的高频和低频也只是局部的，没办法产生一个预设的频带过滤（例如小波变换）进行辨识。

选择的模函数对应了一个自适应（依赖于信号自身的）的时变滤波器。

一个这方面的例子：一个信号由3个部分组成（这3个部分是时间频率上都明显叠加的信号），用上述方法成功的分解了。

分解如图1所示。

这个例子的程序是emd_fmsin2.m另外一个例子（emd_sawtooth.m）强调了EMD潜在的非谐振性质如图2所示。

在这些例子中，线性的非线性的震荡都能被有效的识别和分离。

因而，任何谐振分析（傅里叶，小波，…）可能结束在同类文章中，更少的紧凑和更少的实际意义的分解。

3．算法的改进正如第二部分所定义的，EMD算法依赖于一系列的选项，这些选项需要用户控制，并且需要专业的知识。

在此我们的目的找出更准确的选项，并且给予原来的算法进行改进。

3.1采样率，插值方法和边缘效应EMD的基础操作是估计出上包络和下包络作为极值点之间的插值曲线。

选择的插值法的性能是非常重要的参数。

我们的实验是要确定三次样条插值法作为首选。

其他的插值法（线性的或者多项式的）会增加筛选的迭代次数，并且会产生过分解信号，这些过分解信号散布在临近的模函数内。

其次，自从这个算法运用到实际的离散时间信号中，需要注意的是极值点必须能够正确的找出，其必要条件是要求大量的过采样（关于这点将会在第4段进行进一步的研究）最后，边界状态也要进行考虑，因为极小的的误差会扩散到有限的测量数据长度内。

作为这方面的考量，通过镜像沿拓加入极值点到边界附近能够得到较好的结果。

3.2筛选的停止准则当筛选过程结束时，抽取的模函数会是满意的。

在这方面有2个必备的条件：第一个是极值点数和过零点数最多相差1个，第二个是上包络和下包络的均值必须近似等于零或者满足某种准则。

均值的幅度有多小需要通过计算与对应的模函数的幅值，但是用一个过低的阈值来终止迭代过程会导致如前面提到的问题（过迭代产生过分解）。

作为一个改进的准则，我们目前所认识的，我们提出（在emd.m）一个新的准备建立在2个阈值上的。

目的是保证整体的小的波动在整体大的过程中有意义。

这个总量引进了模幅值a(t)=(emax(t)-emin(t))/2，和估计函数o(t)=|m(t)/a(t)|，为了让筛选迭代停止，当满足o(t)< θ1对于规定的部分（1-α）在整个段，当o(t)< θ2对于剩余的部分。

一旦设定α≈0.05，θ1≈0.05，θ2≈0.5（这些在emd.m是默认值）3.3局部EMD经典EMD应用中，筛选迭代适合于整个数据长度，并且追求只要有局部区域存在（出现在包络均值处），就不认为足够小。

然而，正如文中已经提到的，因为好的局部逼近会有对其他部分信号进行数据污染的缺点，所以整个信号会出现过迭代。

这种情况特别容易发生在在一样的幅值和在过分解通过扩散到其他临近的模函数。

原始算法的分层的和非线性的并不能保证序列信号的EMD会是每个EMD的序列。

因此，本文提出的第一个改进算法是基于原始的EMD基础上的。

这个新的算法称其为局部EMD(local_emd.m)，这个新算法在赛选过程中引入了一个中间步骤：对有异常的较大的数据的局部区域进行识别和隔离，这个附加的步骤只适合处理这些异常区域。

通过引进一个重量函数w(t)来实现上述功能。

当满足o(t)> θ1，的时间范围时，w(t)=1,否则w(t)缓慢的衰减至0。

第二部分所描述的原EMD算法的第4步就改成了d(t)=x(t)-w(t)m(t)。

3.4在线EMD第二个改进算法是基于这样一种事实提出的：筛选过程依赖于对极值点的插值，因此对给定的数据长度的数据点只需要使用有限个极值点（三次样条插值5个极小值和5个极大值）。

因此模函数的抽取不需要了解整个信号（或者前一次残余信号）。

这就为在线运行，并且将之应用于数据流（emd_online.m）的改进的EMD算法铺平了道路。

模函数的抽取的必要条件同样适用于筛选过程，这样做能够防止出现可能的不连续的点，这样就需要了解整个信号了，筛选运行次数是与先验知识相符的，并且少量的迭代次数（少于10次，特别是是4次）能够产生足量的有意义IMF。

为了更有效应用在线EMD算法，我们建议对局部算法上再加一个平滑窗。

新的有效数据进入窗的前边沿，当满足停止准则时，整块数据进入后边沿。

基于这样一种远离，从而IMF和其相应的残量成分可以计算出来。

因此整个算法能够应用于残量成分，继而从剩余的模函数中抽离有效成分。

这个新的算法怎样工作例子的可以通过运行ex_omline.m这个程序观看到，这个例子中的解析信号是由图1中的3个周期信号组成的。

运用在线EMD算法，对2000个数据点组成的基础数据块进行分解，最终分解得到16000个数据点。

除了将在线算法的基本使用在数据流的分解中，在计算负担上(处理较长的数据，计算负担会变的更重)，有人指出它较标准算法有更好的优越性。

4． 性能基础自从EMD 本质上定义为一种算法，但是并没有作为一种分析方法所承认，它的性能估计是困难的，并且需要大量的模拟实验。

我们将从两方面报告他的原理，一方面是从先前将EMD 应用到实数信号所表现出的非凡性，另一方面是从对分解的理解。

4.1分量和采样率当我们分析一个单分量信号，EMD 作为分辨单元理想状态下只分离出1个模函数（这个模函数正好识别这个单分量）且不包含残余分量。

然而，即时出去边缘效应的影响，也很难得到这样理想化的分解结果，这是因为无法避免采样率对EMD 分解的影响。

如果一个周期内采样的点很少，那么在这样一个采样率下，采样得到的极值点也就不准确了。

图3（emd_sam;ling.m ）通过固定频率f 的方程221/21()(([][])/[])f f n ne f x n d n x n =-∑∑产生了一个信号进行EMD 分解验证以上由于采样率造成的错误。

1[]d n 代表了从频率为f 的分量[]f x n 抽离出来的第一个EMD 模函数。

结果是，即时是这样一个分量估计，也完全取决于频率f ：当分量的周期是采样率的倍数时，错误发生在极小值时，我们会观察到2()e f Cf ≤。

4.2分量的分离在由2个分量构成的信号1122[]cos2cos2x n a f n a f n ππ=+中，其中211/2f f <<,理想状态下EMD 抽离出2个模函数，尽管采用了适当的采样率，但是第一个模函数要求采样率为1f ,采样率为2f 时对模函数1的影响比或函数2要大。

发生在抽离过程中的错误可以通过准则（1）的扩展重量得到量化。

1f 对应于模函数1，2f 对应于模函数2（会出现较少的错误点）（emd_separation.m ）。

上面的结果用图4描述的，结果表明对于复杂结构的信号，对整个区域进行分量分离是困难的，特别是当1f >1/4时。

观察模型依赖于赋值比12:/a a ρ=，但是在第一个逼近出都显示了同样一个特性：许多错误都包含在三角区域内，这些三角区域被2条通过这个区域的直线所限制。

换句话说，对于一个给定的频率1f ，对于每一个幅值比必然存在一个适当的1ρα<使得111():[,]f f f ρβα=（121,()f f f β∈）不能被分离。

这个发现认为EMD 可以作为是一组常数Q 的滤波器组来理解，这个结论与文献[1,4,7]提到的包含宽带噪声的随机过程的结论是一致的5． 小结EMD 是新的有前景的非平稳非线性处理方法，但是仍需要更好的解释。

这篇文章讨论了算法现存的问题，目的是使该算法能够得到更有效的应用，并且提出了许多具有初始的性能估计。

这篇文章在使用EMD 方面提供了一个新的视角和做了一些改进的算法，但这些工作都是经验性的，需要对该算法做进一步的理论研究。

参考资料[1] K.T. Coughlin and K.K. Tung, “11-year solar cycle in the stratosphere extracted by the empirical mode decomposition method,” Adv. Space Res.,Nov. 2002 (submitted).[2] N.E. Huang, Z. Shen, S.R. Long, M.L. Wu, H.H. Shih, Q. Zheng, N.C. Yen, C.C. Tung andH.H. Liu, “The empirical mode decomposition and Hilbert spectrum for nonlinear and nonstationary time series analysis,” Proc. Roy. Soc. London A,V ol. 454, pp. 903–995, 1998.[3] P. Flandrin, Time-Frequency/Time-Scale Analysis, Academic Press, 1999.[4] P. Flandrin, G. Rilling and P. Gon¸ calv` es, “Empirical Mode Decomposition as a filter bank,” IEEE Sig. Proc. Lett., 2003 (in press).[5] R. Fournier, “Analyse stochastique modale du signal stabilom´ etrique. Application ` al’´ etude del’´ equilibre chez l’Homme,” Th` ese de Doctorat, Univ. Paris XII Val de Marne, 2002.[6] E.P. Souza Neto, M.A. Custaud, C.J. Cejka,P. Abry, J. Frutoso, C. Gharib and P. Flandrin, “Assessment of cardiovascular autonomic control by the Empirical Mode Decomposition,”4th Int. Workshop on Biosignal Interpretation,Como (I), pp. 123-126, 2002.[7] Z. Wu and N.E. Huang, “A study of the characteristics of white noise using the EmpiricalMode Decomposition method,” Proc. Roy. Soc. London A, Dec. 2002 (submitted).[8] Z. Wu, E.K. Schneider, Z.Z. Hu and L. Cao,“The impact of global warming on ENSO varia bility in climate records,” COLA Technical Report, CTR 110, Oct. 2001.[9] www.ens-lyon.fr/~flandrin/software.html。