发明名称 基于EEMD和小波包变换的早期故障敏感特征提取方法
摘要 本发明涉及一种基于EEMD和小波包变换的早期故障敏感特征提取方法,其包括以下步骤:1)将采集到的机电设备原始振动信号进行EEMD分解,并加入白噪声,分解得到IMF分量;2)选择与故障密切相关的敏感IMF分量,忽略其它不相关的IMF分量;3)对经步骤2)选择的敏感IMF分量进行正交小波包分解,得到各个节点的小波系数;4)对获得的小波包系数采用Hilbert变换和傅里叶变换提取其包络,并计算其功率谱,获得各个节点小波包系数所对应的功率谱作为早期故障的敏感特征,实现敏感特征的自动获取。本发明能实现自适应信号分解,便于实现敏感特征的自动化获取,提高诊断精度和速度,实现机电系统的快准稳诊断。本发明可以广泛在机电设备故障诊断领域中应用。
申请公布号 CN103091096A 申请公布日期 2013.05.08
申请号 CN201310025831.X 申请日期 2013.01.23
申请人 北京信息科技大学 发明人 王红军;徐小力
分类号 G01M13/00(2006.01)I;G01M13/02(2006.01)I;G01M13/04(2006.01)I 主分类号 G01M13/00(2006.01)I
代理机构 北京远大卓悦知识产权代理事务所(普通合伙) 11369 代理人 贺持缓
主权项 1.一种基于EEMD和小波包变换的早期故障敏感特征提取方法,其包括以下步骤:1)将采集到的机电设备原始振动信号y(t)进行EEMD分解,并加入白噪声,分解得到IMF分量c<sub>j</sub>(t);2)选择与故障密切相关的敏感IMF分量c<sub>j</sub>(t),忽略其它不相关的IMF分量,采用将IMF分量和原始振动信号y(t)之间的相关系数作为判断指标,选择与故障密切相关的敏感IMF分量的步骤如下:①将步骤1)中得到的所有IMF分量和原始振动信号y(t)进行归一化处理,得到原始振动信号归一化处理后的信号X<sub>i</sub>,归一化处理后的IMF分量值X<sub>i</sub>′,原始振动信号归一化处理后的均值<img file="FDA00002769632800011.GIF" wi="45" he="53" />和归一化化处理后得到的IMF分量的均值<img file="FDA00002769632800012.GIF" wi="82" he="64" />②计算归一化处理后的所有IMF分量与原始振动信号y(t)的相关系数μ<sub>i</sub>,i=1,2,...,n:<maths num="0001"><![CDATA[<math><mrow><msub><mi>&mu;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>X</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>-</mo><mover><msup><mi>X</mi><mo>&prime;</mo></msup><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mrow><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></msqrt><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>X</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>-</mo><mover><msup><mi>X</mi><mo>&prime;</mo></msup><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>式中,M为振动信号的采样点数;③当相关系数μ<sub>i</sub>≥λ时,保留第i个IMF分量c<sub>i</sub>(t);当μ<sub>i</sub>≤λ,剔除第i个IMF分量c<sub>i</sub>(t),并且令r<sub>n</sub>=r<sub>n</sub>+c<sub>i</sub>,其中r<sub>n</sub>是分解余项,n表示余项阶数;λ为一个固定阈值,取最大相关系数的一个比值:λ=max(μ<sub>i</sub>)/η,i=1,2,...,n;式中,η是一个大于1.0的比例系数,取η=10.0;3)对经步骤2)选择的敏感IMF分量进行正交小波包分解,得到各个节点的小波系数;4)对获得的小波包系数采用Hilbert变换和傅里叶变换提取其包络,并计算其功率谱,获得各个节点小波包系数所对应的功率谱作为早期故障的敏感特征,实现敏感特征的自动获取。
地址 100192 北京市海淀区清河小营东路12号166信箱