基于相关分析的EMD过程中有效IMF的判定方法技术

本发明专利技术公开了一种基于相关分析的EMD过程中有效IMF的判定方法，涉及专门适用于特定应用的数字计算或数据处理的设备或方法技术领域。所述方法包括如下步骤：采集故障振动信号，将该振动信号作为原始信号；将原始信号使用EMD方法进行分解得到若干个IMF分量；计算每一个IMF分量与原始信号的相关系数，画出相关系数分布图；利用噪声判断准则确定噪声IMF分量；利用伪分量判断准则确定伪IMF分量；将判定为噪声的IMF分量和伪IMF分量去除，得到有效的IMF分量。所述方法有效的解决了EMD过程产生的IMF分量中的噪声分量和伪分量的判定问题，进而得到有效分量，提高EMD的分析效果。

【技术实现步骤摘要】

本专利技术涉及专门适用于特定应用的数字计算或数据处理的设备或方法
，尤其涉及一种基于相关分析的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，使得IMF1～IMFk是以噪声为主导模态的IMF分量；伪分量主要集中在阶数较大的IMF分量中，因此存在一个数字h，使得从第h个IMF往后的所有分量IMFh～IMFN都是伪IMF分量；剩下的IMFk+1～IMFh-1是以有用信号为主导模态的分量。进一步的技术方案在于：所述的计算每一个IMF分量与原始信号的相关系数的方法包括如下步骤：相关系数ρxy的计算公式为： ρ x y = Σ n = 0 ∞ x ( n ) y ( n ) [ Σ n = 0 ∞ x 2 ( n ) Σ n = 0 ∞ y 2 ( n ) ] 1 / 2 ]]>式中x(n)表示EMD分解前的原始信号，y(n)表示IMF分量。进一步的技术方案在于：所述的根据相关系数分布图并利用噪声判断准则确定噪声IMF分量的方法包括如下步骤：计算各个IMF分量与原信号的相关系数，画出相关系数分布图，在所述相关系数分布图中找到第一个相关系数发生突变的点K，该点对应的IMF分量为IMFk，将IMF1～IMFk分量视为噪声IMF分量。进一步的技术方案在于：将各个IMF分量与原信号的相关系数点依次连接构成相关系数曲线，曲线方向发生反转的相关系数对应的点为相关系数发生突变的点。进一步的技术方案在于：所述的根据相关系数分布图并利用伪分量判断准则确定伪IMF分量的方法包括如下步骤：计算各个IMF分量与原信号的相关系数，画出相关系数分布图，从图中找到第一个相关系数小于A的点h，该点对应的IMF分量为IMFh，将IMFh+1～IMFN分量视为伪分量。进一步的技术方案在于：所述A的值为0.1。采用上述技术方案所产生的有益效果在于：假设原始信号经EMD方法分解得到N个频率从高到低排列的IMF分量，噪声成分主要集中在信号的高频段，即EMD方法分解得到的阶数较小的IMF当中，因此存在一个数字k，使得IMF1～IMFk是以噪声为主导模态的IMF分量；伪分量主要集中在阶数较大的IMF分量中，因此存在一个数字h，使得从第h个IMF往后的所有分量IMFh～IMFN都是伪IMF分量；剩下的IMFk+1～IMFh-1是以有用信号为主导模态的分量。因此，所述方法有效的解决了EMD过程产生的IMF分量中的噪声分量和伪分量的判定问题，进而得到有效分量，提高EMD的分析效果。附图说明图1是本专利技术所述方法的流程图；图2齿轮磨损故障振动信号时域图；图3a-3b齿轮磨损故障振动信号EMD分解得到的若干个IMF图；图4是若干个IMF分量的相关系数分布图；图5重构信号的时域图。具体实施方式下面结合本专利技术实施例中本文档来自技高网...

【技术保护点】
一种基于相关分析的EMD过程中有效IMF的判定方法，其特征在于所述方法包括如下步骤：采集故障振动信号，将该振动信号作为原始信号；将原始信号使用EMD方法进行分解得到若干个IMF分量；计算每一个IMF分量与原始信号的相关系数，画出相关系数分布图；根据相关系数分布图并利用噪声判断准则确定噪声IMF分量；根据相关系数分布图并利用伪分量判断准则确定伪IMF分量；将判定为噪声的IMF分量和伪IMF分量去除，得到有效的IMF分量。

【技术特征摘要】
1.一种基于相关分析的EMD过程中有效IMF的判定方法，其特征在于所述方法包括如下步骤：采集故障振动信号，将该振动信号作为原始信号；将原始信号使用EMD方法进行分解得到若干个IMF分量；计算每一个IMF分量与原始信号的相关系数，画出相关系数分布图；根据相关系数分布图并利用噪声判断准则确定噪声IMF分量；根据相关系数分布图并利用伪分量判断准则确定伪IMF分量；将判定为噪声的IMF分量和伪IMF分量去除，得到有效的IMF分量。2.如权利要求1所述的基于相关分析的EMD过程中有效IMF的判定方法，其特征在于：采用压电加速度传感器采集故障振动信号。3.如权利要求1所述的基于相关分析的EMD过程中有效IMF的判定方法，其特征在于：假设原始信号经EMD方法分解得到N个频率从高到低排列的IMF分量，噪声成分主要集中在信号的高频段，即EMD方法分解得到的阶数较小的IMF当中，因此存在一个数字k，使得IMF1～IMFk是以噪声为主导模态的IMF分量；伪分量主要集中在阶数较大的IMF分量中，因此存在一个数字h，使得从第h个IMF往后的所有分量IMFh～IMFN都是伪IMF分量；剩下的IMFk+1～IMFh-1是以有用信号为主导模态的分量。4.如权利要求1所述的基于相关分析的EMD过程中有效IMF的判定方法，其特征在于，所述的计算每一个IMF分量与原始信号的相关系数的方法包括如下步骤：相关系数ρxy的计算公式为： ρ x y = Σ n = 0 ∞ x ( n ) y ( n ) [ Σ n ...

【专利技术属性】
技术研发人员：郝如江，王胜江，牛志雷，杨红娜，李辉，
申请(专利权)人：石家庄铁道大学，
类型：发明
国别省市：河北;13