[发明专利]基于相关分析的EMD过程中有效IMF的判定方法

说明书

技术领域

本发明涉及专门适用于特定应用的数字计算或数据处理的设备或方法技术领域，尤其涉及一种基于相关分析的EMD过程中有效IMF的判定方法。

背景技术

经验模态分解(Empirical Mode Decomposition,EMD)方法由Huang在1998年提出，至今已在多个领域得到了十分广泛的应用。信号经过EMD的“筛分”过程，将信号进行分解得到若干个不同频率的本征模态分量(IMF)，频率从高到低排列，其频率成分和带宽随信号的不同而不同，体现了该方法的自适应性，因此EMD方法可以作为一个自适应滤波器。由于EMD方法的这一特点，我们可以根据所要研究信号的特点和要求，选择性地将适当的IMF相加，进行信号重构，把我们感兴趣的某一频率范围内的信号特征凸显出来，从而达到滤波器滤波的效果。信号经EMD分解后得到有限个频率从高到低的IMF分量，其中分解得到的前几个阶数较小的IMF对应于信号的高频成分，通常认为它们主要包含了原信号中的噪声成分；分解得到的后几个阶数较大的IMF对应于信号的低频成分，通常认为它们包含的噪声成分很少。

EMD去噪方法的主要思想是把原信号经过层层分解成为不同特征时间尺度的IMF分量，有选择性地将某些IMF分量重新组合起来，得到一个新的信号，就可以组成低通滤波器、高通滤波器、带通滤波器或带阻滤波器，从而达到去噪的目的。大多数基于EMD的去噪方法都是将IMF分量中阶次较低的高频分量作为噪声直接去除，但是具体应该去除几个、去除哪几个分量目前还没有特定的准则，很多情况下都要依靠人的经验，不能排除人的主观因素对其造成的影响。

发明内容

本发明所要解决的技术问题是提供一种基于相关分析的EMD过程中有效IMF的判定方法，所述方法有效的解决了EMD过程产生的IMF分量中的噪声分量和伪分量的判定问题，进而得到有效分量，提高EMD的分析效果。

为解决上述技术问题，本发明所采取的技术方案是：一种基于相关分析的EMD过程中有效IMF的判定方法，其特征在于所述方法包括如下步骤：

采集故障振动信号，将该振动信号作为原始信号；

将原始信号使用EMD方法进行分解得到若干个IMF分量；

计算每一个IMF分量与原始信号的相关系数，画出相关系数分布图；

根据相关系数分布图并利用噪声判断准则确定噪声IMF分量；

根据相关系数分布图并利用伪分量判断准则确定伪IMF分量；

将判定为噪声的IMF分量和伪IMF分量去除，得到有效的IMF分量。

进一步的技术方案在于：采用压电加速度传感器采集故障振动信号。

进一步的技术方案在于：假设原始信号经EMD方法分解得到N个频率从高到低排列的IMF分量，噪声成分主要集中在信号的高频段，即EMD方法分解得到的阶数较小的IMF当中，因此存在一个数字k，使得IMF₁～IMF_k是以噪声为主导模态的IMF分量；伪分量主要集中在阶数较大的IMF分量中，因此存在一个数字h，使得从第h个IMF往后的所有分量IMF_h～IMF_N都是伪IMF分量；剩下的IMF_k+1～IMF_h-1是以有用信号为主导模态的分量。

进一步的技术方案在于：所述的计算每一个IMF分量与原始信号的相关系数的方法包括如下步骤：

相关系数ρ_xy的计算公式为：

式中x(n)表示EMD分解前的原始信号，y(n)表示IMF分量。

进一步的技术方案在于：所述的根据相关系数分布图并利用噪声判断准则确定噪声IMF分量的方法包括如下步骤：

计算各个IMF分量与原信号的相关系数，画出相关系数分布图，在所述相关系数分布图中找到第一个相关系数发生突变的点K，该点对应的IMF分量为IMF_k，将IMF₁～IMF_k分量视为噪声IMF分量。

进一步的技术方案在于：将各个IMF分量与原信号的相关系数点依次连接构成相关系数曲线，曲线方向发生反转的相关系数对应的点为相关系数发生突变的点。

进一步的技术方案在于：所述的根据相关系数分布图并利用伪分量判断准则确定伪IMF分量的方法包括如下步骤：