专利名称：一种机械故障预测的特征提取方法
技术领域：
本发明涉及一种机械故障预测方法，特别是关于一种机械故障预测的特征提取方法。
背景技术：
故障预测是保障机电设备长期安全、满负荷运行的关键技术，是机电故障诊断学研究的重点之一。在机电设备故障预测中，特征提取是一个重要环节，也是故障预测中的难 题。复杂机电系统是非线性系统，其运行状态具有非线性特点，设备运行的数据虽然提供了 设备运行状况的极其丰富、详细的信息，但给故障敏感特征提取带来很大困难。流形学习(Manifold Learning)算法是近年来发展起来的非线性降维机器学习算法，以保持数据局部结构的方式将高维输入投影到低维空间，发现隐藏在数据中的内在几 何结构与规律性。随着流形学习理论研究的深入，其在医学、图像处理、语音信号处理、天文 数据处理等领域应用广泛。目前，有将流形学习方法用于冲击故障特征提取和喘振监测技 术、轴承故障分类和设备状态趋势分析中的应用，但这些方法均为流形学习局部方法，一方 面是不能很好的保持数据的整体性，在从高维到低维的映射过程中，很可能会将相隔较远 的点映射到近邻点的位置，因此不能有效提取低维特征；另一方面是流形学习所应用的对 象大多为轴承故障和压缩机，所应用的领域较具体，因此没有通用性。
发明内容
针对上述问题，本发明的目的是提供一种针对性和准确性较高、能保持数据整体 性、有效提取低维特征、且应用较为广泛的机械故障预测的特征提取方法。为实现上述目的，本发明采取以下技术方案一种机械故障预测的特征提取方法， 其步骤如下(1)通过远程监测诊断中心获取代表设备运行状况的一段连续振动信号，进 行分析；(2)根据连续振动信号，定义一个判断设备运行状态是否异常的综合评判函数f =max(fi; f2，f3)，当出现异常时，进行特征提取；否则不提取特征；其中，f” f2和f3分别 表示振动信号绝对量异常与否的函数、振动信号变化量正常与否的函数和振动信号变化率 正常与否的函数；(3)采用提升小波方法对异常振动信号中所包含的噪声进行降噪；(4) 利用时间延迟法对降噪后的振动信号进行相空间重构到一个m维的相空间中，进而获知 原始动力学系统完整的信息；(5)根据相空间重构后的两样本点之间的欧氏距离判断两该 样本点是否为近邻点；当两样本点为近邻点，则两点之间有边连接；当两样本点为非近邻 点，则没有边连接；(6)根据样本点的边连接距离确定邻近图，边连接的最短路径为De，其 中包含了临近图中任意两个点之间的最短路径距离；(7)规范化变换矩阵S，进而确定矩 阵ι (De)，矩阵t (De)表示流形局部性质；(8)利用平移不变的方法，根据矩阵t (De)构 造一中间计算矩阵K1和一测地距离核矩阵K，并根据中间计算矩阵K1的最大特征值b，使 得测地距离核矩阵K满足正定性；(9)利用谱分解方法，计算正定测地距离核矩阵K的特征值和特征向量；(10)根据亥维塞德函数cn( ε)确定拓扑维数d为乂 = limlofC"(g)^<formula>formula see original document page 5</formula>
中o)=^r:Sl>“-|h— I) ’*9⑴H1坐…；m)根据拓扑维数d进行低
维流形特征提取，得到用于故障预测的低维流形特征τ为:r二论其中， Xli^dXd为正定测地距离核矩阵κ的d个最大特征值；υ = [Ul，...，ud]，为d个最大特 征值对应的特征向量。
<formula>formula see original document page 5</formula>所述综合评判函数f中的函数f\、f2和f3分别为y; = Λ，所述步骤(5)中，所述近邻点的判定方法为当欧氏距离IlXi-XjII彡ε时，则 样本点Xi和Xj为近邻点，两点之间有边连接，且边长为dx(i，j) = Il Xi-Xj Il ；当欧氏距 离IlXi-XjII > ε时，则样本点Xi和Xj为非近邻点，即两点之间没有边连接；其中ε为一 个无穷小量。所述步骤(7)中，所述变换矩阵S和矩阵t (De)分别为Sij = (JlU)， ,(Dg) = -|瓜丑，其中，H 为中心化矩阵，H = I — eN6rNlN, eN = [1，· · ·，1]T e Rn ；
dG(i, j)为两样本点Xi和Xj的最短路径距离，dG(i, j) = dx(i，j) = Il Xi-Xj Il。所述
<formula>formula see original document page 5</formula>
步骤(8)中，所述中间计算矩阵&和测地距离核矩阵K分别为r，
<formula>formula see original document page 5</formula> 为单位矩阵。所述步骤(8)中，所述测地距离核矩阵K为正定矩阵的条件为所述测地距离核矩 阵K中的参数a要满足a > b > 0，其中参数b为所述中间计算矩阵K1的最大特征值。本发明由于采取以上技术方案，其具有以下优点1、本发明由于以故障预测为目 的，通过判断振动信号的异常与否，进而确定是否进行特征提取，因此提高了以故障预测为 目的的特征提取的针对性和准确性。2、本发明由于是基于流形的内在几何结构通过拓扑映 射进行的，是非线性方法，因此更能体现异常数据的本质，实现了特征提取与能量的解耦和 分离。3、本发明由于采用非参数方法，因此不需要对流形做很多参数假设，只需转化为求解 特征值问题，不需要迭代算法，因此大大简化了求解过程，避免了局部极值问题。4、本发明 由于采用全局优化算法，因此能够很好的保持数据的整体性，有效提取流形距离的特征。5、 本发明由于采用变换矩阵满足正定条件，有利于增强特征提取的泛化能力，可以将测试数 据投影到相关的低维流形。本发明可广泛应用于各种机电设备的故障预测特征提取中。
图1是本发明的整体流程示意图。
具体实施例方式下面结合附图和实施例对本发明进行详细的描述。本发明是为了有针对性的进行以故障预测为目的的特征提取，判断设备运行状态 是否异常基础上，进而对运行状态异常的信号进行特征提取，利用提升小波对设备运行的 振动数据进行滤波降低噪声，将降噪后的信号序列通过相空间重构方法重构到高维，引入 流形学习的全局方法，对机电设备运行的异常数据重构进行非线性拓扑映射处理，从而找 到数据集中最本质的非线性结构，找到数据中流形分布变化，提取能用作故障预测、能揭示 设备运行异常的几何结构与规律性的特征。其具体步骤如下(1)通过现有远程监测诊断中心获取能代表设备运行状况的数据，由于振动信号 能够反映设备的机械动特性(即能代表机械设备运行状态性能的指标)，因此选取对故障 敏感的某测点传感器输出的一段连续振动信号Ix1, ...，％}进行分析；(2)根据步骤⑴中的连续振动信号Ix1, . . .，％}，定义一个能判断设备状态异常 与否的综合评判函数f = max(fi; f2, f3)，进而判断设备运行状态是否异常，当设备运行状 态出现异常时，进行特征提取；否则不进行特征提取；函数f = max(fi; f2, f3)中的函数f” f2和f3分别为
<formula>formula see original document page 6</formula>
<formula>formula see original document page 6</formula>(2)
<formula>formula see original document page 6</formula>(3)上述各公式中，f\、f2和f3分别表示振动信号绝对量异常与否的函数、振动信号变 化量正常与否的函数和振动信号变化率正常与否的函数；A、B和C分别为函数f\、f2和函数 f3中定义的相应阈值，当函数值为1时表示信号异常，函数值为0时表示信号正常；μ为振
动信号振动量的均值…；(3)采用提升小波方法对异常振动信号Ix1,. . .，％}中所包含的噪声进行降噪，因 为基于非线性流形的特征提取方法本身对噪声很敏感；(4)利用时间延迟法对降噪后的振动信号Ix1,. . .，％}进行相空间重构，将异常振 动信号{X1，...，XN}重构到一个m维的相空间中，由于信号重构前后在拓扑上是等价的，因 此可以通过重构的动力学系统获知原始动力学系统完整的信息；
(5)根据相空间重构后的两样本点Xi和Xj之间的欧氏距离Il Xi-Xj Il判断这两样 本点Xi和X^是否为近邻点，当Il Xi-Xj Il彡ε时，则样本点Xi和&为近邻点，即两点之间 有边连接，且边长为dx(i，j) = Il Xi-Xj Il ；当IlXi-XjII > ε时，则样本点Xi和Xj为非近 邻点，即两点之间没有边连接；其中ε为一个无穷小量；(6)根据样本点的边连接距离确定邻近图，若两样本点Xi和Xj之间有边连接，则设置边连接的最短路径距离为de(i，j) =dx(i，j)；反之dc(i，j) =C ，计算‘⑴j)= min{dG(i, j), dc(i，l)+dc(l，j)}，其中，1 为 1，2，...，n，η 为样本点数；则最短路径 Dc = {dG(i, j)}中包含了临近图中任意两个点之间的最短路径距离；(7)规范化变换矩阵S，进而确定矩阵t (De)，矩阵t (De)表示流形局部性质；其 中，变换矩阵S和t (De)分别为Sij =d2G(i,j),(4)I(Dc) = HSH,(5)上述公式(5)中，H为中心化矩阵，// = /- </7V，eN= [1，· · ·，1]T e Rn ；(8)利用平移不变的方法，根据矩阵t (De)构造一中间计算矩阵K1和一测地距离 核矩阵K，并根据中间计算矩阵K1的最大特征值b，使得测地距离核矩阵K满足正定性；中 间计算矩阵K1和测地距离核矩阵K分别为
“0 2i(Dr) ,、K1 =G/ ,(6)
L-I -MDg)]K = i(DG) + 2a(DG) + ^a2H,(7)其中，I为单位矩阵；当参数a满足a彡b > 0时，则测地距离核矩阵K为正定；(9)利用谱分解方法，计算正定测地距离核矩阵K的特征值和特征向量；(10)根据亥维塞德(Heaviside)函数Cn( ε )确定拓扑维数d为^ = Iim1T^'(8)
1 η π[β 当；f < O其中，C^ipHhiMtj1当…·’(11)根据拓扑维数d进行低维流形特征提取，得到用于故障预测的低维流形特征 T = diag{^,...^)UT，其中，A1,..., Ad为正定测地距离核矩阵K的d个最大特征值；U =[U1,...，ud]，为d个最大特征值对应的特征向量。上述各实施例仅是本发明的优选实施方式，在本技术领域内，凡是基于本发明技 术方案上的变化和改进，不应排除在本发明的保护范围之外。
权利要求
一种机械故障预测的特征提取方法，其步骤如下(1)通过远程监测诊断中心获取代表设备运行状况的一段连续振动信号，进行分析；(2)根据连续振动信号，定义一个判断设备运行状态是否异常的综合评判函数f＝max(f1，f2，f3)，当出现异常时，进行特征提取；否则不提取特征；其中，f1、f2和f3分别表示振动信号绝对量异常与否的函数、振动信号变化量正常与否的函数和振动信号变化率正常与否的函数；(3)采用提升小波方法对异常振动信号中所包含的噪声进行降噪；(4)利用时间延迟法对降噪后的振动信号进行相空间重构到一个m维的相空间中，进而获知原始动力学系统完整的信息；(5)根据相空间重构后的两样本点之间的欧氏距离判断两该样本点是否为近邻点；当两样本点为近邻点，则两点之间有边连接；当两样本点为非近邻点，则没有边连接；(6)根据样本点的边连接距离确定邻近图，边连接的最短路径为DG，其中包含了临近图中任意两个点之间的最短路径距离；(7)规范化变换矩阵S，进而确定矩阵ι(DG)，矩阵ι(DG)表示流形局部性质；(8)利用平移不变的方法，根据矩阵ι(DG)构造一中间计算矩阵K1和一测地距离核矩阵K，并根据中间计算矩阵K1的最大特征值b，使得测地距离核矩阵K满足正定性；(9)利用谱分解方法，计算正定测地距离核矩阵K的特征值和特征向量；(10)根据亥维塞德函数Cn(ε)确定拓扑维数d为 <mrow><mi>d</mi><mo>=</mo><munder> <mi>lim</mi> <mrow><mi>ϵ</mi><mo>→</mo><mn>0</mn> </mrow></munder><mfrac> <mrow><mi>log</mi><msub> <mi>C</mi> <mi>n</mi></msub><mrow> <mo>(</mo> <mi>ϵ</mi> <mo>)</mo></mrow> </mrow> <mrow><mi>log</mi><mi>ϵ</mi> </mrow></mfrac> </mrow>其中，(11)根据拓扑维数d进行低维流形特征提取，得到用于故障预测的低维流形特征T为 <mrow><mi>T</mi><mo>=</mo><mi>diag</mi><mrow> <mo>(</mo> <msubsup><mi>λ</mi><mn>1</mn><mfrac> <mn>1</mn> <mn>2</mn></mfrac> </msubsup> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mo>,</mo> <msubsup><mi>λ</mi><mi>d</mi><mfrac> <mn>1</mn> <mn>2</mn></mfrac> </msubsup> <mo>)</mo></mrow><msup> <mi>U</mi> <mi>T</mi></msup><mo>,</mo> </mrow>其中，λ1，…，λd为正定测地距离核矩阵K的d个最大特征值；U＝[u1，…，ud]，为d个最大特征值对应的特征向量。FSA00000007027600012.tif
2.如权利要求1所述的一种机械故障预测的特征提取方法，其特征在于所述步骤(2) 中，所述综合评判函数f中的函数f\、f2和f3分别为<formula>formula see original document page 2</formula>其中，Α、Β和C分别为函数f\、f2和函数f3中定义的相应阈值，当函数值为1时表示信 号异常，函数值为O时表示信号正常；μ为振动信号振动量的均值，<formula>formula see original document page 2</formula>
3.如权利要求1所述的一种机械故障预测的特征提取方法，其特征在于所述步骤(5) 中，所述近邻点的判定方法为当欧氏距离IlXi-XjII彡ε时，则样本点Xi和Xj为近邻点， 两点之间有边连接，且边长*dx(i，j) = Il Xi-Xj Il ；当欧氏距离IlXi-XjII > ε时，则样本 点Xi和\为非近邻点，即两点之间没有边连接；其中ε为一个无穷小量。
4.如权利要求1所述的一种机械故障预测的特征提取方法，其特征在于所述步骤(7) 中，所述变换矩阵S和矩阵t (De)分别为<formula>formula see original document page 3</formula>其中，H为中心化矩阵<formula>formula see original document page 3</formula>)为两样本点Xi*Xj 的最短路径距离，<formula>formula see original document page 3</formula>
5.如权利要求1所述的一种机械故障预测的特征提取方法，其特征在于所述步骤(8) 中，所述中间计算矩阵K1和测地距离核矩阵K分别为<formula>formula see original document page 3</formula>其中，I为单位矩阵。
6.如权利要求1所述的一种机械故障预测的特征提取方法，其特征在于所述步骤(8) 中，所述测地距离核矩阵K为正定矩阵的条件为所述测地距离核矩阵K中的参数a要满足 a^b> 0，其中参数b为所述中间计算矩阵K1的最大特征值。
全文摘要
本发明涉及一种机械故障预测的特征提取方法，(1)获取代表设备运行状况的连续振动信号进行分析；(2)定义判断设备运行状态是否异常的综合评判函数；(3)对振动信号中包含的噪声降噪；(4)用时间延迟法对振动信号进行相空间重构；(5)根据欧氏距离判断两样本点是否为近邻点；(6)根据样本点的边连接距离确定邻近图；(7)规范化变换矩阵S，确定矩阵ι(DG)；(8)利用平移不变的方法，根据矩阵ι(DG)构造中间计算矩阵K1和测地距离核矩阵K，使矩阵K满足正定性；(9)计算矩阵K的特征值和特征向量；(10)根据亥维塞德函数确定拓扑维数；(11)根据拓扑维数进行低维流形特征提取。本发明通过判断振动信号的异常与否，提高了特征提取的针对性和准确性。本发明可广泛应用于各种机械设备的故障预测特征提取中。
文档编号G01M99/00GK101799366SQ20101010132
公开日2010年8月11日 申请日期2010年1月27日 优先权日2010年1月27日
发明者徐小力, 王少红, 王红军, 陈涛 申请人:北京信息科技大学