[发明专利]一种经验模态分解筛选迭代过程终止准则的方法

说明书

本发明公开了一种经验模态分解筛选迭代过程的终止准则的方法，确定数据的气候态，用数据x(t)的线性拟合来代表其气候态；迭代求解，每次迭代得到IMF的雏形后，计算余量相对方差；找出相对方差最小时对应的迭代次数，此时的IMF雏形就是本次迭代过程的终点。进入下一个迭代循环过程，获取下一个IMF，直至余量xsubgt;i+1/subgt;(t)为趋势项为止。本发明的有益效果是本发明能够使EMD的分解过程更加规范、客观，提高EMD分解的可操作性和结果的一致性，弱化其“经验”特征。

技术领域

本发明属于数据分析技术领域，涉及筛选迭代余量的相对方差作为经验模态分解方法筛选迭代过程的终止准则。

背景技术

一维的Fourier分解、小波分析，二维的PCA/EOF等方法，都是从低频开始分解，获取不同的模态/特征函数，且低频部分通常取得最大的方差贡献率。经验模态分解(EMD)是一种自适应的数据分析方法，没有固定的基底函数。EMD将一个时间序列分解为一系列从高频到低频的IMF，以及一个非线性趋势项。分解得到的各本征模态函数(IMF)通常并不正交。因而，各IMF在数据中所占的方差百分比之和，可能大于1，不能解释为各IMF的方差贡献率。但是，仍然可以根据低频部分的相对方差来确定IMF分解结果是否合理。EMD分解过程中每个IMF的求解都是通过一个称为“筛选(sifting)”的迭代过程得到。目前EMD筛选迭代的终止准则有如下几种：

(1)Cauchy型准则(Huang et al.,1998)：相邻两次迭代的IMF雏形的相对均方根误差小于某个预先给定的充分小量；

(2)平均值准则(Flandrin，2004)：相邻两次迭代的IMF雏形，在每一个时间点上的偏差都小于预先给定的充分小量；

(3)S准则(Huang et al.,2003)：连续S次迭代的IMF雏形，其极值点及跨零点的个数保持不变，或者最多相差1；

(4)固定迭代次数(Wu and Huang，2010)：不考虑数据的特征，也不考虑将要求解的IMF为第几模态，所有筛选的迭代次数取固定值。这个固定值通常取为8～12。

以上筛选迭代的终止准则中，前3种准则通常有较大的主观性，比如Cauchy性准则和平均值准则中的“充分小量”如何选取？S准则中的S如何选取？一般来说，小的“充分小量”和大的“S”，需要更多的迭代次数才能达到迭代终止的条件。Wang et al.(2010)证明，筛选过程中迭代次数的增加将导致分解得到的IMF振幅趋于常数，即IMF接近线性函数。这与EMD追求IMF“非线性非稳定”特征是背离的。第4种准则(即固定次数的迭代终止准则)，从其迭代次数的选择上，也容易引起困惑。筛选中迭代8次或10次，得到的IMF雏形可能有较大的差异，更不能保证满足前3种准则中的任意一个。

综上，经验模态分解(EMD)通过一系列迭代筛选过程将数据分解为多个不同频率段的本征模态函数(IMF)和一个非线性趋势。每个IMF的产生，都要经过一个筛选迭代过程。目前采用的4种迭代终止准则，缺乏客观的判定标准，迭代次数的选择要靠经验给出，可能存在迭代不足或者迭代过度的情况。

发明内容

本发明的目的在于提供一种经验模态分解筛选迭代过程的终止准则的方法，本发明的有益效果是本发明能够使EMD的分解过程更加规范、客观，提高EMD分解的可操作性和结果的一致性，淡化其“经验”特征。

本发明所采用的技术方案按照以下步骤进行：

步骤1：确定数据的气候态。用数据x(t)的线性拟合来代表其气候态；

步骤2：迭代求解。每次迭代得到IMF的雏形后，计算余量的相对方差；

步骤3：找出相对方差最小时对应的迭代次数，此时的IMF雏形就是本次迭代过程的终点；