CN104748961A - 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 - Google Patents
基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 Download PDFInfo
- Publication number
- CN104748961A CN104748961A CN201510146250.0A CN201510146250A CN104748961A CN 104748961 A CN104748961 A CN 104748961A CN 201510146250 A CN201510146250 A CN 201510146250A CN 104748961 A CN104748961 A CN 104748961A
- Authority
- CN
- China
- Prior art keywords
- signal
- noise
- gear
- noise reduction
- sigma
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种基于SVD分解降噪和相关性EEMD熵特征的齿轮故障诊断方法,利用加速度振动传感器采集实验台齿轮振动信号,得到的信号包括齿轮正常、齿轮断齿、齿轮少齿、齿轮磨损四种故障类型;利用通过相关性分析和信噪比优化的SVD分解降噪方法对包含高斯白噪声的模拟强噪声背景的四种齿轮状态信号进行降噪处理;利用EEMD分解方法分别对降噪之后的四类信号进行分解,根据相关系数选取有效的IMF分量;将得到的每组有效IMF分量进行样本熵计算,并构建由IMF样本熵组成的特征向量;利用PNN神经网络识别四种不同的齿轮故障。本发明能够在强噪声背景下有效地识别齿轮故障类型,是一种有效的齿轮故障诊断方法。
Description
技术领域
本发明属于故障诊断技术领域,涉及一种基于SVD分解降噪和相关性EEMD熵特征的齿轮故障诊断方法。
背景技术
齿轮作为旋转机械设备中重要的零部件,其发生故障会导致机器工作效率降低,甚至造成重大经济损失。因此,研究齿轮的状态监测与故障诊断技术对于提高机械设备的运行效率及维修效能,避免人员财产损失具有重要的现实意义。
基于振动信号分析的齿轮故障诊断技术是一种行之有效的诊断方式,具有较高的精度。常见的诊断方法如:时频特征参数法、倒频谱法、EMD分解等。但是在进行信号处理的过程中,对于齿轮传动产生的非线性、非平稳信号,这些方法都有各自的局限性。时频特征参数法中的不同指标只对特定的齿轮缺陷判别较为有效;针对强噪声背景等复杂状况,倒频谱法难以发现齿轮的缺陷频率;传统的经验模式分解(Empirical Mode Decomposition,EMD)方法于1998年由N.E.Huang提出,因其适用于非线性、非平稳信号的研究,近年来得到了广泛的应用,但是EMD方法仍存在诸多问题,包括端点效应、模态混叠、迭代循环次数多等。EEMD方法在EMD的基础上加以改进,利用了高斯白噪声具有频率均匀分布的统计特性,当信号加入白噪声之后,将使信号在不同尺度上具有连续性,以减小模态混叠的程度。
发明内容
本发明的目的是提供一种基于SVD分解降噪和相关性EEMD熵特征的齿轮故障诊断方法,针对强噪声背景对齿轮运转状态进行监测,以发现并判断齿轮故障,避免机械设备出现较为严重的故障。
本发明所采用的技术方案是,一种基于SVD分解降噪和相关性EEMD熵特征的齿轮故障诊断方法,包括以下步骤:
步骤1,利用加速度振动传感器采集实验台齿轮振动信号,得到的信号包括齿轮正常、齿轮断齿、齿轮少齿、齿轮磨损四种故障类型;
步骤2,利用通过相关性分析和信噪比优化的SVD分解降噪方法对包含高斯白噪声的模拟强噪声背景的四种齿轮状态信号进行降噪处理;
步骤3,利用EEMD分解方法分别对降噪之后的四类信号进行分解,根据相关系数选取有效的IMF分量;
步骤4,将得到的每组有效IMF分量进行样本熵计算,并构建由IMF样本熵组成的特征向量;
步骤5,利用PNN神经网络识别四种不同的齿轮故障。
本发明的特征还在于,
步骤2中SVD分解降噪方法对包含高斯白噪声的模拟强噪声背景的四种齿轮状态信号进行降噪处理的过程包括以下步骤:
对于含有噪声的齿轮故障振动信号y(k)(k=1,2,…,N),根据相空间重构理论,将其映射到m×n(m<n)维相空间内,得到满足m+n+1=N的Hankle矩阵Bm,对Bm进行奇异值分解,求取矩阵Bm的奇异值,保留前k个的奇异值而置零其他位置的奇异值,利用奇异值分解的逆过程得到B'm,B'm即为Bm的一个最佳逼近,这样就达到了降噪的效果,再对B'm中的反对角元素取平均即完成信号降噪过程;
轨道矩阵Bm的重构阶次的选择,通过信号的信噪比和相关系数确定有效重构阶次,其中,
1)相关系数的计算公式:
其中:yk无噪声信号的第k个数据点;mk为降噪后的信号的第k个数据点,n为数据长度;
2)信噪比的计算公式:
其中,y(k)为无噪声信号的第k个数据点,为含噪声信号的第k个数据点,N为信号长度;
在信号与噪声的共同作用下,降噪后信号与原始信号的相关系数和降噪后信号的信噪比会随着重构阶次的增加而快速增大,当阶次达到一定值时,相关系数和信噪比的增长速度逐步变缓并趋于稳定,这时重构信号包含的有效信息趋于饱和,所以根据相关系数和信噪比达到最大的阶次作为重构阶次,可以有效地保留有含噪信号的有用信息。
步骤3中,EEMD分解方法包括以下步骤:
1)在信号y(t)中加入白噪声mj(t),其中幅值均值为0,标准差为原信号标准差的0.3倍则:
yi(t)=y(t)+mj(t) (14)
式中,i为yi(t)分解的次数;
2)对yi(t)进行EMD分解,得到若干IMF分量djk(t)与余项ej(t);其中djk(t)为第j次加入白噪声后所得的第k个IMF分量;
3)重复步骤1和步骤2N次,得到消除模态混叠的IMF为:
信号EEMD分解的最终结果为:
步骤3中,根据相关系数选取有效的IMF分量包括以下步骤:
相关系数的计算公式为:
其中:yk为EEMD分解前信号的第k个数据点;fk为IMF分量的第k个数据点,n为数据长度。
步骤4中,样本熵的计算其计算步骤如下:
1)对于一个拥有Nt个点的数据序列,y(1),y(2),…,y(Nt)可以组成一组m维的矢量:
y(i)=[y(i),y(i+1),…,y(i+m-1)]
i=1,2,…,Nt-m+1 (18)
2)定义两个m维的矢量Y(i)和Y(j)之间的最大距离为:
k=0,1,…,m-1 (19)
3)对于给定的阀值r,从计算d(i,j)<r的数目除以Nt-m+1的值,记为Bi m(r),即:
4)求Bi m(r)平均值:
5)根据维数m,重复上面的步骤1~4可以得到
6)计算样本熵Se(m,r):
步骤5中PNN神经网络识别过程为:根据输入特征向量构建合适的PNN网络,初始化网络后,利用训练样本对网络进行训练,训练结束后将测试样本输入到网络进行诊断识别并输出结果。
本发明的有益效果是,该方法针对传统的齿轮故障诊断方法在强噪声背景下提取有效故障特征困难的缺陷,提出了优化的SVD分解降噪方法和相关性EEMD熵特征相结合的齿轮故障特征提取方法,并利用PNN神经网络对齿轮故障类型进行识别,一定程度上丰富和完善了故障诊断方法。利用信噪比和相关系数确定重构阶次的SVD分解降噪方法对含噪信号降噪,能够有效提高信号的信噪比。优化的SVD分解降噪方法和相关性EEMD熵特征相结合可以准确地提取出齿轮故障特征信息,通过PNN神经网络可以有效地识别齿轮的故障类型。该方法可有效应用于工业部门中机械设备状态监测及诊断,特别适用于强噪声背景下工作的煤矿机械。
附图说明
图1是本发明的齿轮故障诊断方法的流程图。
图2是降噪后信号与原始信号的相关系数随重构阶次的变化图。
图3是降噪后信号的信噪比随重构阶次的变化图。
图4a是具有断齿故障的齿轮振动信号时域图。
图4b是具有断齿故障的齿轮加噪后的振动信号时域图。
图4c是本发明所述方法降噪处理后的断齿齿轮振动信号时域图。
图5a是降噪之前的断齿齿轮振动信号FFT频谱图。
图5b是加噪之后的断齿齿轮振动信号FFT频谱图。
图5c是降噪之后的断齿齿轮振动信号FFT频谱图。
图6是降噪之后的断齿齿轮振动信号EEMD分解图。
图7是本发明的PNN神经网络诊断识别流程图。
图8是四种齿轮状态的PNN神经网络识别结果图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
一种基于SVD分解降噪和相关性EEMD熵特征的齿轮故障诊断方法,流程如图1所示,综合应用了相关性分析、信噪比、SVD分解方法、EEMD分解方法、样本熵理论和概率神经网络理论,
具体来说,包括以下步骤:
1)利用加速度振动传感器采集实验台齿轮振动信号,得到的信号包括齿轮正常、齿轮断齿、齿轮少齿、齿轮磨损四种故障类型;
2)利用通过相关性分析和信噪比优化的SVD分解降噪方法对包含高斯白噪声的模拟强噪声背景的四种齿轮状态信号进行降噪处理;
3)利用EEMD分解方法分别对降噪之后的四类信号进行分解,根据相关系数选取有效的IMF分量;
4)将得到的每组有效IMF分量进行样本熵计算,并构建由IMF样本熵组成的特征向量;
5)利用PNN神经网络识别四种不同的齿轮故障。
其中,步骤2)中SVD分解降噪方法过程为:对于含有噪声的齿轮故障振动信号y(k)(k=1,2,…,N),根据相空间重构理论,将其映射到m×n(m<n)维相空间内,得到满足m+n+1=N的Hankle矩阵Bm,对Bm进行奇异值分解,求取矩阵Bm的奇异值,根据Frobenious范数意义下矩阵最佳逼近定理,保留前k个较大的奇异值而置零其他位置的奇异值,利用奇异值分解的逆过程得到B'm,B'm即为Bm的一个最佳逼近,这样就达到了降噪的效果。再对B'm中的反对角元素取平均即完成信号降噪过程。
根据SVD分解降噪的过程可以发现,能否对含噪信号进行有效降噪关键在于准确确定轨道矩阵Bm的有效重构阶次k。本发明通过信号的信噪比和相关系数来确定有效重构阶次。
SVD奇异值分解降噪方法还包括轨道矩阵Bm的重构阶次的选择,其选择方法依据求取的相关系数和信噪比:
1)相关系数的计算公式:
其中:yk无噪声信号的第k个数据点;mk为降噪后的信号的第k个数据点,n为数据长度。
2)信噪比的计算公式:
其中,y(k)为无噪声信号的第k个数据点,为含噪声信号的第k个数据点,N为信号长度。
在信号与噪声的共同作用下,降噪后信号与原始信号的相关系数和降噪后信号的信噪比会随着重构阶次的增加而快速增大,当阶次达到一定值时,相关系数和信噪比的增长速度逐步变缓并趋于稳定,这时重构信号包含的有效信息趋于饱和,所以根据相关系数和信噪比达到最大的阶次作为重构阶次,可以有效地保留有含噪信号的有用信息。在实际应用中发现,利用相关系数和信噪比可以很好的判断重构阶次的大小。
步骤3)中,EEMD分解方法包括以下步骤:
1)在信号y(t)中加入白噪声mj(t),其中幅值均值为0,标准差为原信号标准差的0.3倍则:
yi(t)=y(t)+mj(t) (25)
式中,i为yi(t)分解的次数。
2)对yi(t)进行EMD分解,得到若干IMF分量djk(t)与余项ej(t)。其中djk(t)为第j次加入白噪声后所得的第k个IMF分量。
3)重复步骤1和步骤2N次。得到消除模态混叠的IMF为:
信号EEMD分解的最终结果为:
步骤3)中的EEMD分解方法还包括有效IMF分量的选取,是根
据相关系数决定的,相关系数的计算公式为:
其中:yk为EEMD分解前信号的第k个数据点;fk为IMF分量的第k个数据点,n为数据长度。
步骤4)中的特征向量的构建,根据相关系数选取的有效IMF分量进行样本熵计算并构成特征向量。
其中,样本熵的计算,其计算步骤如下:
1)对于一个拥有Nt个点的数据序列,y(1),y(2),…,y(Nt)可以组成一组m维的矢量:
y(i)=[y(i),y(i+1),…,y(i+m-1)]
i=1,2,…,Nt-m+1 (29)
2)定义两个m维的矢量Y(i)和Y(j)之间的最大距离为:
3)对于给定的阀值r,从计算d(i,j)<r的数目除以Nt-m+1的值,记为Bi m(r),即:
4)求Bi m(r)平均值:
5)根据维数m,重复上面的步骤1~4可以得到
6)计算样本熵Se(m,r):
步骤5)中,PNN神经网络识别过程为:根据输入特征向量构建合适的PNN网络,初始化网络后,利用训练样本对网络进行训练,训练结束后将测试样本输入到网络进行诊断识别并输出结果。
下面予以分述:
本实验在美国Spectra Quest公司的机械故障综合模拟实验台上进行,该实验台配备了齿面磨损与断齿等故障。本发明主要对齿轮正常、齿轮断齿、齿轮少齿和齿面磨损四种状态进行实验,每种齿轮状态采集50个数据样本,共200个样本。采样频率设置为10kHZ,样本长度为5000,电机转速20HZ。下面以齿轮断齿故障为例,实现信号预处理降噪、特征提取及故障诊断过程。
本实验利用加速度振动传感器采集实验台齿轮振动信号,得到的信号包括齿轮正常、齿轮断齿、齿轮少齿、和齿轮磨损四种类型。为了模拟强噪声背景的工作环境,在采集的齿轮振动信号中加入高斯白噪声,使得信噪比为0.25,此外在采集齿轮振动信号时,往往受到设备振动、其它零部件振动以及外界因素的影响,产生一定的噪声。因此,在对齿轮振动信号进行分析之前需要对信号进行降噪处理,以消除故障信号的噪声背景提高信噪比。利用SVD分解降噪是一种有效的降噪处理方法,已经被应用于许多工程领域。
在步骤2)中的SVD分解降噪方法过程为:
对于一个含有噪声的齿轮故障振动信号y(k)(k=1,2,…,N),根据相空间重构理论,将原始信号映射到m×n(m<n)维相空间内,得到满足m+n+1=N的Hankle矩阵Bm
对于含噪齿轮故障振动信号而言,Bm可以表示成Bm=B+W。其中,B、W分别为有用信号的轨道矩阵和噪声信号的轨道矩阵,对原始信号降噪就是寻找B的最佳逼近问题,降噪效果随着逼近程度变好而明显变好。
对Bm进行奇异值分解,可以得到Bm=UΣVH,其中:矩阵Σ的主对角线元素λi(1,2,3,…,m)为矩阵Bm的奇异值。如果保留前k个较大的奇异值而置零其它较小的奇异值,利用奇异值分解的逆过程就可以得到Bm',Bm'即为Bm的一个最佳逼近,这样即达到了降噪的效果。为了得到降噪后的信号y'(n),需要对Bm'中的反对角元素取平均,即
其中:α=max(1,k-m+1);β=min(n,k)。
根据SVD分解降噪的过程可以发现,能否对含噪信号进行有效降噪关键在于准确确定轨道矩阵Bm的有效重构阶次k。本发明通过信号的信噪比和相关系数来确定有效重构阶次。
1)相关系数的计算公式:
其中:yk无噪声信号的第k个数据点;mk为降噪后的信号的第k个数据点,n为数据长度。
2)信噪比的计算公式:
其中,y(k)为无噪声信号的第k个数据点,为含噪声信号的第k个数据点,N为信号长度。
在信号与噪声的共同作用下,降噪后信号与原始信号的相关系数和降噪后信号的信噪比会随着重构阶次的增加而快速增大,当阶次达到一定值时,相关系数和信噪比的增长速度逐步变缓并趋于稳定,这时重构信号包含的有效信息趋于饱和,所以根据相关系数和信噪比达到最大的阶次作为重构阶次,可以有效地保留有含噪信号的有用信息。在实际应用中发现,利用相关系数和信噪比可以很好的判断重构阶次的大小。
本发明根据信噪比和相关系数判断的重构阶次为10。其中信噪比和相关系数随重构阶次的变化如图2、图3所示。由图可以看出,随着重构阶次的增加,信噪比和相关系数快速增大到最大值,之后出现微小波动,最后趋于稳定;另一方面,通过信噪比和相关系数的大小很容易发现SVD分解重构起到了很好的降噪效果。
断齿故障信号的时域波形如图4a所示,图4b为加入高斯白噪声之后的断齿齿轮故障信号时域图,图4c是采用本发明所述的SVD分解降噪方法处理后的断齿信号时域图,对比降噪前后的断齿故障时域图可以发现,经过SVD分解降噪处理后的信号,其所含噪声得到了有效的消除。对降噪前后的断齿故障信号进行FFT变换,得到FFT频谱如图5a、5b、5c所示,表明本发明采用的降噪方法可以有效地保留齿轮故障信息,为进一步处理做准备。
在步骤3)中的利用EEMD分解方法对降噪之后的信号y(t)进行分解的
过程如下:
1)在信号y(t)中加入白噪声mj(t),其中幅值均值为0,标准差为原信号标准差的0.3倍则:
yi(t)=y(t)+mj(t) (38)
式中,i为yi(t)分解的次数。
2)对yi(t)进行EMD分解,得到若干IMF分量djk(t)与余项ej(t)。其中djk(t)为第j次加入白噪声后所得的第k个IMF分量。
3)重复步骤1和步骤2N次。得到消除模态混叠的IMF为:
信号EEMD分解的最终结果为:
本发明中确定噪声与故障信号幅值标准差比值为0.3,平均运算次数为100次,得到分解结果如图6所示,由图可见断齿信号被分解成频率从高到低的9个IMF分量和一个残项。
根据式(28)计算断齿齿轮EEMD分解得到的各IMF分量的相关系数,如表1所示:
表1断齿齿轮EEMD分解得到的各IMF分量相关系数
相关系数可以说明各IMF分量与原始信号的相关程度,相关系数越大,说明相关程度越大,包含的故障特征信息越多。由表1的相关系数可以发现,前3个IMF分量的相关系数相比于其它系数更大,包含了齿轮故障的主要特征信息,所以可以对前3个IMF分量进行近一步的分析。
在步骤4)中的对得到的有效IMF分量进行样本熵计算过程为:
1)对于一个拥有Nt个点的数据序列,y(1),y(2),…,y(Nt)可以组成一组m维的矢量:
y(i)=[y(i),y(i+1),…,y(i+m-1)]
i=1,2,…,Nt-m+1 (41)
2)定义两个m维的矢量Y(i)和Y(j)之间的最大距离为:
3)对于给定的阀值r,从计算d(i,j)<r的数目除以Nt-m+1的值,记为Bi m(r),即:
4)求Bi m(r)平均值:
5)根据维数m,重复上面的步骤1~4可以得到
6)计算样本熵Se(m,r):
得到的样本熵如表2所示:
表2四种齿轮状态下的样本熵
由表2可以发现,正常齿轮EEMD样本熵较小,其它故障样本上均大于正常齿轮,说明齿轮发生故障时,样本熵值有明显的变化,区别于正常状态。但仅根据样本熵还不能得到准确的故障类型,故需要对信号进行进一步分析。
在步骤5)中的利用概率神经网络识别不同齿轮故障的过程为:
本发明采用PNN神经网络对不同故障类型进行识别。PNN神经网络是基于Bayes分类规则与Parzen窗的概率密度函数估计方法发展而来的一种并行算法,在解决分类问题方面得到了广泛应用。在实际应用中,根据输入特征向量构建合适的PNN网络,初始化网络后,利用训练样本对网络进行训练,训练结束后将测试样本输入到网络进行诊断识别并输出结果。其诊断识别流程如图7所示。
齿轮故障识别中,取每种齿轮故障类型的30个数据样本,共120个组成训练样本,再各取20个数据样本,共80个组成测试样本,用以验证识别模式的有效性。诊断识别结果如图8所示。观察图8可以发现,在类别1中,即齿轮正常状态下,识别率为100%;在类别2中,即齿轮断齿状态下,有两个错误识别,识别率为90%;在类别3中,即齿轮少齿状态下,识别率为100%;在类别4中,即齿轮磨损状态下,识别率为100%。由此可见识别率最低为齿轮断齿故障,识别率为90%,但总体识别率达到97.5%,表明本发明所述方法可以有效识别齿轮故障类型。
Claims (6)
1.一种基于SVD分解降噪和相关性EEMD熵特征的齿轮故障诊断方法,其特征在于,包括以下步骤:
步骤1,利用加速度振动传感器采集实验台齿轮振动信号,得到的信号包括齿轮正常、齿轮断齿、齿轮少齿、齿轮磨损四种故障类型;
步骤2,利用通过相关性分析和信噪比优化的SVD分解降噪方法对包含高斯白噪声的模拟强噪声背景的四种齿轮状态信号进行降噪处理;
步骤3,利用EEMD分解方法分别对降噪之后的四类信号进行分解,根据相关系数选取有效的IMF分量;
步骤4,将得到的每组有效IMF分量进行样本熵计算,并构建由IMF样本熵组成的特征向量;
步骤5,利用PNN神经网络识别四种不同的齿轮故障。
2.根据权利要求1所述的齿轮故障诊断方法,其特征在于,所述步骤2中SVD分解降噪方法对包含高斯白噪声的模拟强噪声背景的四种齿轮状态信号进行降噪处理的过程包括以下步骤:
对于含有噪声的齿轮故障振动信号y(k)(k=1,2,…,N),根据相空间重构理论,将其映射到m×n(m<n)维相空间内,得到满足m+n+1=N的Hankle矩阵Bm,对Bm进行奇异值分解,求取矩阵Bm的奇异值,保留前k个的奇异值而置零其他位置的奇异值,利用奇异值分解的逆过程得到B'm,B'm即为Bm的一个最佳逼近,这样就达到了降噪的效果,再对B'm中的反对角元素取平均即完成信号降噪过程;
轨道矩阵Bm的重构阶次的选择,通过信号的信噪比和相关系数确定有效重构阶次,其中,
1)相关系数的计算公式:
其中:yk无噪声信号的第k个数据点;mk为降噪后的信号的第k个数据点,n为数据长度;
2)信噪比的计算公式:
其中,y(k)为无噪声信号的第k个数据点,为含噪声信号的第k个数据点,N为信号长度;
在信号与噪声的共同作用下,降噪后信号与原始信号的相关系数和降噪后信号的信噪比会随着重构阶次的增加而快速增大,当阶次达到一定值时,相关系数和信噪比的增长速度逐步变缓并趋于稳定,这时重构信号包含的有效信息趋于饱和,所以根据相关系数和信噪比达到最大的阶次作为重构阶次,可以有效地保留有含噪信号的有用信息。
3.根据权利要求1所述的齿轮故障诊断方法,其特征在于,所述步骤3中,EEMD分解方法包括以下步骤:
1)在信号y(t)中加入白噪声mj(t),其中幅值均值为0,标准差为原信号标准差的0.3倍则:
yi(t)=y(t)+mj(t) (3)式中,i为yi(t)分解的次数;
2)对yi(t)进行EMD分解,得到若干IMF分量djk(t)与余项ej(t);其中djk(t)为第j次加入白噪声后所得的第k个IMF分量;
3)重复步骤1和步骤2N次,得到消除模态混叠的IMF为:
信号EEMD分解的最终结果为:
4.根据权利要求1所述的齿轮故障诊断方法,其特征在于,所述步骤3中,根据相关系数选取有效的IMF分量包括以下步骤:
相关系数的计算公式为:
其中:yk为EEMD分解前信号的第k个数据点;fk为IMF分量的第k个数据点,n为数据长度。
5.根据权利要求1所述的齿轮故障诊断方法,其特征在于,所述步骤4中,样本熵的计算其计算步骤如下:
1)对于一个拥有Nt个点的数据序列,y(1),y(2),…,y(Nt)可以组成一组m维的矢量:
y(i)=[y(i),y(i+1),…,y(i+m-1)]
i=1,2,…,Nt-m+1 (7)
2)定义两个m维的矢量Y(i)和Y(j)之间的最大距离为:
k=0,1,…,m-1 (8)
3)对于给定的阀值r,从计算d(i,j)<r的数目除以Nt-m+1的值,记为Bi m(r),即:
4)求Bi m(r)平均值:
5)根据维数m,重复上面的步骤1~4可以得到
6)计算样本熵Se(m,r):
6.根据权利要求1所述的齿轮故障诊断方法,其特征在于,所述步骤5中PNN神经网络识别过程为:根据输入特征向量构建合适的PNN网络,初始化网络后,利用训练样本对网络进行训练,训练结束后将测试样本输入到网络进行诊断识别并输出结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510146250.0A CN104748961A (zh) | 2015-03-30 | 2015-03-30 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510146250.0A CN104748961A (zh) | 2015-03-30 | 2015-03-30 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104748961A true CN104748961A (zh) | 2015-07-01 |
Family
ID=53588977
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510146250.0A Pending CN104748961A (zh) | 2015-03-30 | 2015-03-30 | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104748961A (zh) |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105092023A (zh) * | 2015-08-11 | 2015-11-25 | 西安科技大学 | 基于白噪声统计特征的锚杆振动信号矫正方法 |
CN105258942A (zh) * | 2015-09-30 | 2016-01-20 | 南阳理工学院 | 一种采煤机摇臂齿轮故障诊断系统及方法 |
CN105758634A (zh) * | 2016-02-26 | 2016-07-13 | 中国航空工业集团公司上海航空测控技术研究所 | 一种应用于直升机传动系统振动信号的故障检测方法 |
CN105928701A (zh) * | 2016-04-29 | 2016-09-07 | 石家庄铁道大学 | 基于相关分析的emd过程中有效imf的判定方法 |
CN106197999A (zh) * | 2016-07-08 | 2016-12-07 | 安徽德衍智控科技有限公司 | 一种行星齿轮故障诊断方法 |
CN106226077A (zh) * | 2016-07-01 | 2016-12-14 | 中国科学技术大学 | 一种基于时变奇异值分解的周期性暂态信号的检测方法 |
CN106289780A (zh) * | 2016-09-20 | 2017-01-04 | 北京工业大学 | 一种基于Sparsogram和样本熵的滚动轴承损伤程度识别方法 |
CN106383028A (zh) * | 2016-09-05 | 2017-02-08 | 北京航空航天大学 | 一种齿轮箱故障的诊断方法 |
CN106442368A (zh) * | 2016-07-22 | 2017-02-22 | 天津理工大学 | 基于emd的二氧化硫紫外差分吸收光谱降噪重构方法 |
CN106524962A (zh) * | 2016-09-30 | 2017-03-22 | 中国矿业大学 | 一种采煤机行走轮磨损量检测装置和磨损量检测预警方法 |
CN106644484A (zh) * | 2016-09-14 | 2017-05-10 | 西安工业大学 | Eemd与邻域粗糙集结合的涡桨发动机转子系统故障诊断方法 |
CN107480619A (zh) * | 2017-08-03 | 2017-12-15 | 中国地质大学(武汉) | 基于eemd和排列熵的探地雷达b扫描图像的降噪方法及系统 |
CN107907327A (zh) * | 2017-11-14 | 2018-04-13 | 上海电力学院 | 一种风电机组行星齿轮箱故障诊断方法 |
CN108551412A (zh) * | 2018-05-03 | 2018-09-18 | 网宿科技股份有限公司 | 监控数据降噪处理方法和装置 |
CN108593286A (zh) * | 2018-04-12 | 2018-09-28 | 中国神华能源股份有限公司 | 旋转机械的故障诊断方法及旋转机械的故障诊断装置 |
CN108937847A (zh) * | 2017-05-22 | 2018-12-07 | 北京大学 | 一种评价人体运动协调性的方法 |
CN109858109A (zh) * | 2019-01-14 | 2019-06-07 | 北京工业大学 | 一种基于相关性的emd与形态奇异值分解相结合的齿轮信号降噪方法 |
CN109975013A (zh) * | 2019-04-30 | 2019-07-05 | 南京工程学院 | 基于ivmd-se的风电机组齿轮箱故障特征提取方法 |
CN110133396A (zh) * | 2019-04-09 | 2019-08-16 | 上海电力学院 | 一种基于二分递推svd分解的居民电器开关事件检测方法 |
CN110501603A (zh) * | 2019-07-24 | 2019-11-26 | 昆明理工大学 | 利用emd与神经网络的永富直流换相失败故障诊断方法 |
CN110657985A (zh) * | 2019-10-11 | 2020-01-07 | 重庆邮电大学 | 基于奇异值谱流形分析的齿轮箱故障诊断方法及系统 |
CN110677125A (zh) * | 2019-10-11 | 2020-01-10 | 国网冀北电力有限公司秦皇岛供电公司 | 电弧故障检测的方法及装置 |
CN112578019A (zh) * | 2019-09-27 | 2021-03-30 | 北京化工大学 | 一种基于表面漏磁信号的高速长钢轨内部温度应力集中早期预警方法 |
CN112817056A (zh) * | 2020-12-30 | 2021-05-18 | 吉林大学 | 一种大地电磁信号去噪方法及系统 |
US11188065B2 (en) | 2017-09-23 | 2021-11-30 | Nanoprecise Sci Corp. | System and method for automated fault diagnosis and prognosis for rotating equipment |
CN116929753A (zh) * | 2023-09-18 | 2023-10-24 | 杭州景业智能科技股份有限公司 | 传动齿轮状态检测方法、装置、计算机设备以及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0288938A (ja) * | 1988-09-27 | 1990-03-29 | Kobe Steel Ltd | 減速機歯車の故障診断方法 |
CN103048137A (zh) * | 2012-12-20 | 2013-04-17 | 北京航空航天大学 | 一种变工况下的滚动轴承故障诊断方法 |
CN103226649A (zh) * | 2013-03-25 | 2013-07-31 | 西安交通大学 | 机械早期及复合故障的集成噪声重构经验模式分解方法 |
CN104330255A (zh) * | 2014-11-04 | 2015-02-04 | 徐州隆安光电科技有限公司 | 一种基于多传感器信息融合的齿轮故障诊断方法 |
CN104390781A (zh) * | 2014-11-26 | 2015-03-04 | 中国矿业大学 | 一种基于lmd和bp神经网络的齿轮故障诊断方法 |
-
2015
- 2015-03-30 CN CN201510146250.0A patent/CN104748961A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0288938A (ja) * | 1988-09-27 | 1990-03-29 | Kobe Steel Ltd | 減速機歯車の故障診断方法 |
CN103048137A (zh) * | 2012-12-20 | 2013-04-17 | 北京航空航天大学 | 一种变工况下的滚动轴承故障诊断方法 |
CN103226649A (zh) * | 2013-03-25 | 2013-07-31 | 西安交通大学 | 机械早期及复合故障的集成噪声重构经验模式分解方法 |
CN104330255A (zh) * | 2014-11-04 | 2015-02-04 | 徐州隆安光电科技有限公司 | 一种基于多传感器信息融合的齿轮故障诊断方法 |
CN104390781A (zh) * | 2014-11-26 | 2015-03-04 | 中国矿业大学 | 一种基于lmd和bp神经网络的齿轮故障诊断方法 |
Non-Patent Citations (3)
Title |
---|
宋耀文等: "基于级联奇异值分解和Hilbert变换的滚动轴承故障诊断", 《煤矿机械》 * |
楼军伟等: "EEMD样本熵在轴承故障SVM识别中的研究", 《机械传动》 * |
陈仁祥等: "基于相关系数的EEMD转子振动信号降噪方法", 《振动、测试与诊断》 * |
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105092023A (zh) * | 2015-08-11 | 2015-11-25 | 西安科技大学 | 基于白噪声统计特征的锚杆振动信号矫正方法 |
CN105258942A (zh) * | 2015-09-30 | 2016-01-20 | 南阳理工学院 | 一种采煤机摇臂齿轮故障诊断系统及方法 |
CN105758634A (zh) * | 2016-02-26 | 2016-07-13 | 中国航空工业集团公司上海航空测控技术研究所 | 一种应用于直升机传动系统振动信号的故障检测方法 |
CN105928701A (zh) * | 2016-04-29 | 2016-09-07 | 石家庄铁道大学 | 基于相关分析的emd过程中有效imf的判定方法 |
CN106226077A (zh) * | 2016-07-01 | 2016-12-14 | 中国科学技术大学 | 一种基于时变奇异值分解的周期性暂态信号的检测方法 |
CN106226077B (zh) * | 2016-07-01 | 2018-08-21 | 中国科学技术大学 | 一种基于时变奇异值分解的周期性暂态信号的检测方法 |
CN106197999B (zh) * | 2016-07-08 | 2019-07-09 | 山东中衡光电科技有限公司 | 一种行星齿轮故障诊断方法 |
CN106197999A (zh) * | 2016-07-08 | 2016-12-07 | 安徽德衍智控科技有限公司 | 一种行星齿轮故障诊断方法 |
CN106442368A (zh) * | 2016-07-22 | 2017-02-22 | 天津理工大学 | 基于emd的二氧化硫紫外差分吸收光谱降噪重构方法 |
CN106383028A (zh) * | 2016-09-05 | 2017-02-08 | 北京航空航天大学 | 一种齿轮箱故障的诊断方法 |
CN106644484A (zh) * | 2016-09-14 | 2017-05-10 | 西安工业大学 | Eemd与邻域粗糙集结合的涡桨发动机转子系统故障诊断方法 |
CN106289780A (zh) * | 2016-09-20 | 2017-01-04 | 北京工业大学 | 一种基于Sparsogram和样本熵的滚动轴承损伤程度识别方法 |
CN106524962B (zh) * | 2016-09-30 | 2019-03-19 | 中国矿业大学 | 一种采煤机行走轮磨损量检测装置和磨损量检测预警方法 |
CN106524962A (zh) * | 2016-09-30 | 2017-03-22 | 中国矿业大学 | 一种采煤机行走轮磨损量检测装置和磨损量检测预警方法 |
CN108937847A (zh) * | 2017-05-22 | 2018-12-07 | 北京大学 | 一种评价人体运动协调性的方法 |
CN107480619A (zh) * | 2017-08-03 | 2017-12-15 | 中国地质大学(武汉) | 基于eemd和排列熵的探地雷达b扫描图像的降噪方法及系统 |
US11188065B2 (en) | 2017-09-23 | 2021-11-30 | Nanoprecise Sci Corp. | System and method for automated fault diagnosis and prognosis for rotating equipment |
CN107907327A (zh) * | 2017-11-14 | 2018-04-13 | 上海电力学院 | 一种风电机组行星齿轮箱故障诊断方法 |
CN107907327B (zh) * | 2017-11-14 | 2019-06-11 | 上海电力学院 | 一种风电机组行星齿轮箱故障诊断方法 |
CN108593286A (zh) * | 2018-04-12 | 2018-09-28 | 中国神华能源股份有限公司 | 旋转机械的故障诊断方法及旋转机械的故障诊断装置 |
CN108551412A (zh) * | 2018-05-03 | 2018-09-18 | 网宿科技股份有限公司 | 监控数据降噪处理方法和装置 |
CN109858109A (zh) * | 2019-01-14 | 2019-06-07 | 北京工业大学 | 一种基于相关性的emd与形态奇异值分解相结合的齿轮信号降噪方法 |
CN110133396B (zh) * | 2019-04-09 | 2021-07-20 | 上海电力学院 | 一种基于二分递推svd分解的居民电器开关事件检测方法 |
CN110133396A (zh) * | 2019-04-09 | 2019-08-16 | 上海电力学院 | 一种基于二分递推svd分解的居民电器开关事件检测方法 |
CN109975013A (zh) * | 2019-04-30 | 2019-07-05 | 南京工程学院 | 基于ivmd-se的风电机组齿轮箱故障特征提取方法 |
CN110501603A (zh) * | 2019-07-24 | 2019-11-26 | 昆明理工大学 | 利用emd与神经网络的永富直流换相失败故障诊断方法 |
CN112578019A (zh) * | 2019-09-27 | 2021-03-30 | 北京化工大学 | 一种基于表面漏磁信号的高速长钢轨内部温度应力集中早期预警方法 |
CN110677125A (zh) * | 2019-10-11 | 2020-01-10 | 国网冀北电力有限公司秦皇岛供电公司 | 电弧故障检测的方法及装置 |
CN110657985B (zh) * | 2019-10-11 | 2021-07-06 | 重庆邮电大学 | 基于奇异值谱流形分析的齿轮箱故障诊断方法及系统 |
CN110657985A (zh) * | 2019-10-11 | 2020-01-07 | 重庆邮电大学 | 基于奇异值谱流形分析的齿轮箱故障诊断方法及系统 |
CN110677125B (zh) * | 2019-10-11 | 2022-07-01 | 国网冀北电力有限公司秦皇岛供电公司 | 电弧故障检测的方法及装置 |
CN112817056A (zh) * | 2020-12-30 | 2021-05-18 | 吉林大学 | 一种大地电磁信号去噪方法及系统 |
CN112817056B (zh) * | 2020-12-30 | 2021-12-07 | 吉林大学 | 一种大地电磁信号去噪方法及系统 |
CN116929753A (zh) * | 2023-09-18 | 2023-10-24 | 杭州景业智能科技股份有限公司 | 传动齿轮状态检测方法、装置、计算机设备以及存储介质 |
CN116929753B (zh) * | 2023-09-18 | 2024-01-09 | 杭州景业智能科技股份有限公司 | 传动齿轮状态检测方法、装置、计算机设备以及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104748961A (zh) | 基于svd分解降噪和相关性eemd熵特征的齿轮故障诊断方法 | |
CN103575523B (zh) | 基于FastICA-谱峭度-包络谱分析的旋转机械故障诊断方法 | |
CN101900789B (zh) | 基于小波变换与分形维数的容差模拟电路故障诊断方法 | |
CN104390781A (zh) | 一种基于lmd和bp神经网络的齿轮故障诊断方法 | |
CN105760839A (zh) | 基于多特征流形学习与支持向量机的轴承故障诊断方法 | |
CN105758644A (zh) | 基于变分模态分解和排列熵的滚动轴承故障诊断方法 | |
Huang et al. | Periodic feature oriented adapted dictionary free OMP for rolling element bearing incipient fault diagnosis | |
CN102156042A (zh) | 一种基于信号多特征匹配的齿轮故障诊断方法 | |
CN105928702B (zh) | 基于形态分量分析的变工况齿轮箱轴承故障诊断方法 | |
CN105241666A (zh) | 一种基于信号稀疏表示理论的滚动轴承故障特征提取方法 | |
CN102998118B (zh) | 一种基于形态学滤波和复杂度测度的轴承定量诊断方法 | |
CN103226649B (zh) | 机械早期及复合故障的集成噪声重构经验模式分解方法 | |
CN107817106A (zh) | 基于贝叶斯残余变换‑奇异值分解‑高斯混合隐马尔科夫模型框架的滚动轴承故障诊断方法 | |
Sheng et al. | Applications in bearing fault diagnosis of an improved Kurtogram algorithm based on flexible frequency slice wavelet transform filter bank | |
CN108731945B (zh) | 一种航空发动机转子系统故障信号特征信息的提取方法 | |
CN104408302A (zh) | 一种基于lmd-svd和极限学习机的轴承变工况故障诊断方法 | |
CN104200065A (zh) | 基于随机共振与核主元分析相结合的机械振动信号特征提取方法 | |
CN106404386A (zh) | 一种用于采集、提取及诊断齿轮箱早期故障特征信号的方法 | |
CN104504403A (zh) | 一种基于散射变换的旋转机械故障预测方法 | |
CN104330258A (zh) | 一种基于特征参量的滚动轴承故障灰色关联度辨识方法 | |
Meng et al. | Loose parts detection method combining blind deconvolution with support vector machine | |
CN103926097B (zh) | 一种用于采集和提取低速重载设备故障特征信息的方法 | |
CN103335840A (zh) | 一种矿用钻机变速箱故障智能诊断方法 | |
Huijie et al. | Structured sparsity assisted online convolution sparse coding and its application on weak signature detection | |
CN104089778B (zh) | 一种水轮机振动故障诊断方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150701 |
|
RJ01 | Rejection of invention patent application after publication |