基于相关分析的emd过程中有效imf的判定方法
【专利摘要】本发明公开了一种基于相关分析的EMD过程中有效IMF的判定方法，涉及专门适用于特定应用的数字计算或数据处理的设备或方法技术领域。所述方法包括如下步骤：采集故障振动信号，将该振动信号作为原始信号；将原始信号使用EMD方法进行分解得到若干个IMF分量；计算每一个IMF分量与原始信号的相关系数，画出相关系数分布图；利用噪声判断准则确定噪声IMF分量；利用伪分量判断准则确定伪IMF分量；将判定为噪声的IMF分量和伪IMF分量去除，得到有效的IMF分量。所述方法有效的解决了EMD过程产生的IMF分量中的噪声分量和伪分量的判定问题，进而得到有效分量，提高EMD的分析效果。
【专利说明】
基于相关分析的EMD过程中有效IMF的判定方法
技术领域
[0001] 本发明涉及专门适用于特定应用的数字计算或数据处理的设备或方法技术领域， 尤其涉及一种基于相关分析的EMD过程中有效頂F的判定方法。
【背景技术】
[0002] 经验模态分解（Empirical Mode Decomposition,EMD)方法由Huang在 1998年提 出，至今已在多个领域得到了十分广泛的应用。信号经过EMD的"筛分"过程，将信号进行分 解得到若干个不同频率的本征模态分量(IMF)，频率从高到低排列，其频率成分和带宽随信 号的不同而不同，体现了该方法的自适应性，因此EMD方法可以作为一个自适应滤波器。由 于EMD方法的这一特点，我们可以根据所要研究信号的特点和要求，选择性地将适当的MF 相加，进行信号重构，把我们感兴趣的某一频率范围内的信号特征凸显出来，从而达到滤波 器滤波的效果。信号经EMD分解后得到有限个频率从高到低的MF分量，其中分解得到的前 几个阶数较小的IMF对应于信号的高频成分，通常认为它们主要包含了原信号中的噪声成 分;分解得到的后几个阶数较大的頂F对应于信号的低频成分，通常认为它们包含的噪声成 分很少。
[0003] EMD去噪方法的主要思想是把原信号经过层层分解成为不同特征时间尺度的IMF 分量，有选择性地将某些MF分量重新组合起来，得到一个新的信号，就可以组成低通滤波 器、高通滤波器、带通滤波器或带阻滤波器，从而达到去噪的目的。大多数基于EMD的去噪方 法都是将IMF分量中阶次较低的高频分量作为噪声直接去除，但是具体应该去除几个、去除 哪几个分量目前还没有特定的准则，很多情况下都要依靠人的经验，不能排除人的主观因 素对其造成的影响。
【发明内容】
[0004] 本发明所要解决的技术问题是提供一种基于相关分析的EMD过程中有效MF的判 定方法，所述方法有效的解决了EMD过程产生的MF分量中的噪声分量和伪分量的判定问 题，进而得到有效分量，提高EMD的分析效果。
[0005] 为解决上述技术问题，本发明所采取的技术方案是:一种基于相关分析的EMD过程 中有效IMF的判定方法，其特征在于所述方法包括如下步骤：
[0006] 采集故障振动信号，将该振动信号作为原始信号；
[0007] 将原始信号使用EMD方法进行分解得到若干个頂F分量；
[0008]计算每一个頂F分量与原始信号的相关系数，画出相关系数分布图；
[0009] 根据相关系数分布图并利用噪声判断准则确定噪声IMF分量；
[0010] 根据相关系数分布图并利用伪分量判断准则确定伪IMF分量；
[0011]将判定为噪声的IMF分量和伪IMF分量去除，得到有效的IMF分量。
[0012] 进一步的技术方案在于:采用压电加速度传感器采集故障振动信号。
[0013] 进一步的技术方案在于:假设原始信号经EMD方法分解得到N个频率从高到低排列 的IMF分量，噪声成分主要集中在信号的高频段，即EMD方法分解得到的阶数较小的IMF当 中，因此存在一个数字k，使得頂Fi-MFk是以噪声为主导模态的MF分量;伪分量主要集中 在阶数较大的頂F分量中，因此存在一个数字h，使得从第h个頂F往后的所有分量頂Fh~頂F N 都是伪MF分量;剩下的頂Fk+1~頂Fh-i是以有用信号为主导模态的分量。
[0014] 进一步的技术方案在于:所述的计算每一个IMF分量与原始信号的相关系数的方 法包括如下步骤：
[0015] 相关系数pxy的计算公式为：
[0017]式中x(n)表示EMD分解前的原始信号，y(n)表示IMF分量。
[0018] 进一步的技术方案在于:所述的根据相关系数分布图并利用噪声判断准则确定噪 声頂F分量的方法包括如下步骤：
[0019] 计算各个IMF分量与原信号的相关系数，画出相关系数分布图，在所述相关系数分 布图中找到第一个相关系数发生突变的点K，该点对应的頂F分量为頂Fk，将頂Fi~頂Fk分量 视为噪声頂F分量。
[0020] 进一步的技术方案在于:将各个IMF分量与原信号的相关系数点依次连接构成相 关系数曲线，曲线方向发生反转的相关系数对应的点为相关系数发生突变的点。
[0021] 进一步的技术方案在于:所述的根据相关系数分布图并利用伪分量判断准则确定 伪IMF分量的方法包括如下步骤：
[0022]计算各个MF分量与原信号的相关系数，画出相关系数分布图，从图中找到第一个 相关系数小于A的点h，该点对应的IMF分量为IMFh，将IMFh+1~IMFN分量视为伪分量。
[0023]进一步的技术方案在于:所述A的值为0.1。
[0024]采用上述技术方案所产生的有益效果在于:假设原始信号经EMD方法分解得到N个 频率从高到低排列的頂F分量，噪声成分主要集中在信号的高频段，即EMD方法分解得到的 阶数较小的頂F当中，因此存在一个数字k，使得頂Fi~頂Fk是以噪声为主导模态的頂F分量； 伪分量主要集中在阶数较大的IMF分量中，因此存在一个数字h，使得从第h个IMF往后的所 有分量頂Fh~頂Fn都是伪MF分量;剩下的MF k+1~頂Fh-i是以有用信号为主导模态的分量。 因此，所述方法有效的解决了 EMD过程产生的頂F分量中的噪声分量和伪分量的判定问题， 进而得到有效分量，提高EMD的分析效果。
【附图说明】
[0025]图1是本发明所述方法的流程图；
[0026]图2齿轮磨损故障振动信号时域图；
[0027]图3a_3b齿轮磨损故障振动信号EMD分解得到的若干个頂F图；
[0028] 图4是若干个IMF分量的相关系数分布图；
[0029] 图5重构信号的时域图。
【具体实施方式】
[0030] 下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整 地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0031] 在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以 采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的 情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。
[0032]如图1所示，本发明公开了一种基于相关分析的EMD过程中有效頂F的判定方法，包 括如下步骤：
[0033] 采用压电加速度传感器采集故障振动信号，将该振动信号作为原始信号；
[0034] 将原始信号使用EMD方法进行分解得到若干个頂F分量。
[0035] 计算每一个IMF分量与原始信号的相关系数，画出相关系数分布图，相关系数pxy的 计算公式为：
[0037]式中x (n)表示EMD分解前的原始信号，y (n)表示頂F分量。
[0038]根据相关系数分布图并利用噪声判断准则确定噪声IMF分量，确定噪声IMF分量的 方法如下:计算各个IMF分量与原信号的相关系数，画出相关系数分布图，在所述相关系数 分布图中找到第一个相关系数发生突变的点K，该点对应的頂F分量为頂Fk，将頂Fi~頂Fk分 量视为噪声IMF分量。将各个IMF分量与原信号的相关系数点依次连接构成相关系数曲线， 曲线方向发生反转的相关系数对应的点为相关系数发生突变的点。
[0039]根据相关系数分布图并利用伪分量判断准则确定伪IMF分量，确定伪IMF分量的方 法如下:计算各个IMF分量与原信号的相关系数，画出相关系数分布图，从图中找到第一个 相关系数小于〇. 1的点h，该点对应的IMF分量为IMFh，将IMFh+1~IMFN分量视为伪分量。
[0040]将判定为噪声的頂F分量和伪頂F分量去除，得到有效的頂F分量。
[00411设原始信号经EMD方法分解得到N个频率从高到低排列的MF分量，噪声成分主要 集中在信号的高频段，即EMD方法分解得到的阶数较小的頂F当中，因此存在一个数字k，使 得頂Fi~頂Fk是以噪声为主导模态的頂F分量;伪分量主要集中在阶数较大的頂F分量中，因 此存在一个数字h，使得从第h个IMF往后的所有分量IMFh~IMFn都是伪IMF分量；剩下的 頂是以有用信号为主导模态的分量。因此，所述方法有效的解决了 EMD过程产生 的IMF分量中的噪声分量和伪分量的判定问题，进而得到有效分量，提高EMD的分析效果。 [0042]实验分析：
[0043]采集齿轮磨损故障振动信号，信号时域波形如图2所示。将故障信号进行EMD分解， 得到12个頂F和一个趋势项，趋势项对本发明没有影响，故将其忽略不计，各个頂F的时域图 如图3a-3b所示。
[0044]计算每一个IMF与齿轮磨损故障振动信号的相关系数，如下表所示：
[0045]表1齿轮磨损信号EMD分解得到的IMFs的相关系数
[0047]画出各个MF的相关系数分布图，如图4所示。相关系数分布突变的第一个分量是 IMF3,根据噪声分量判断准则，噪声和有效分量的分界点是k = 3,贝iJIMFdMFs就是以噪声 为主导模态的分量;根据伪分量判断准则，相关系数小于0.1的第一个分量是IMF8,有效分 量和伪分量的分界点是h = 8,则IMFs~IMF12就是伪分量。将以噪声为主导模态的分量和伪 分量全部去除，剩下的所有分量即为有效的IMF，将其相加得到重构信号。重构信号时域图 如图5所示。与图3相比，峰值之间具有明显的周期性，计算得出周期为16.67Hz，与齿轮磨损 故障频率16.97Hz很接近。因此，当本方法应用到齿轮故障判断中，可以判断出该齿轮是否 发生故障。
【主权项】
1. 一种基于相关分析的EMD过程中有效MF的判定方法，其特征在于所述方法包括如下 步骤： 采集故障振动信号，将该振动信号作为原始信号； 将原始信号使用EMD方法进行分解得到若干个頂F分量； 计算每一个IMF分量与原始信号的相关系数，画出相关系数分布图； 根据相关系数分布图并利用噪声判断准则确定噪声IMF分量； 根据相关系数分布图并利用伪分量判断准则确定伪IMF分量； 将判定为噪声的IMF分量和伪IMF分量去除，得到有效的IMF分量。2. 如权利要求1所述的基于相关分析的EMD过程中有效頂F的判定方法，其特征在于： 采用压电加速度传感器采集故障振动信号。3. 如权利要求1所述的基于相关分析的EMD过程中有效頂F的判定方法，其特征在于： 假设原始信号经EMD方法分解得到N个频率从高到低排列的IMF分量，噪声成分主要集 中在信号的高频段，即EMD方法分解得到的阶数较小的頂F当中，因此存在一个数字k，使得 ~頂Fk是以噪声为主导模态的頂F分量;伪分量主要集中在阶数较大的頂F分量中，因此 存在一个数字h，使得从第h个IMF往后的所有分量IMFh~IMFn都是伪IMF分量;剩下的IMFk+i ~nffVi是以有用信号为主导模态的分量。4. 如权利要求1所述的基于相关分析的EMD过程中有效MF的判定方法，其特征在于，所 述的计算每一个IMF分量与原始信号的相关系数的方法包括如下步骤： 相关系数pxy的计算公式为：式中X (η)表示EMD分解前的原始信号，y (η)表示IMF分量。5. 如权利要求1所述的基于相关分析的EMD过程中有效MF的判定方法，其特征在于，所 述的根据相关系数分布图并利用噪声判断准则确定噪声MF分量的方法包括如下步骤： 计算各个IMF分量与原信号的相关系数，画出相关系数分布图，在所述相关系数分布图 中找到第一个相关系数发生突变的点K，该点对应的頂F分量为頂Fk，将頂F!~頂Fk分量视为 噪声頂F分量。6. 如权利要求5所述的基于相关分析的EMD过程中有效MF的判定方法，其特征在于:将 各个IMF分量与原信号的相关系数点依次连接构成相关系数曲线，曲线方向发生反转的相 关系数对应的点为相关系数发生突变的点。7. 如权利要求1所述的基于相关分析的EMD过程中有效MF的判定方法，其特征在于，所 述的根据相关系数分布图并利用伪分量判断准则确定伪MF分量的方法包括如下步骤： 计算各个IMF分量与原信号的相关系数，画出相关系数分布图，从图中找到第一个相关 系数小于A的点h，该点对应的IMF分量为IMFh，将IMFh+1~IMFN分量视为伪分量。8. 如权利要求7所述的基于相关分析的EMD过程中有效MF的判定方法，其特征在于，所 述A的值为0.1。
【文档编号】G01M13/02GK105928701SQ201610279173
【公开日】2016年9月7日
【申请日】2016年4月29日
【发明人】郝如江, 王胜江, 牛志雷, 杨红娜, 李辉
【申请人】石家庄铁道大学