一种基于图论的轴承半监督故障诊断方法技术

本发明专利技术提供一种基于图论的轴承半监督故障诊断方法，该方法首先利用可视图算法将通过传感器得到的轴承原始振动加速度信号转换成复杂网络；再计算复杂网络的结构参数，提取度分布的均值和标准差及网络复杂指数；最后利用基于图的半监督学习处理无标签样本，实现轴承故障诊断。本发明专利技术基于少量的标签样本和无标签样本，本发明专利技术实现了变工况且样本类别不平衡下情况下的轴承故障诊断，故障识别准确率高，具有显著的使用价值。

【技术实现步骤摘要】
一种基于图论的轴承半监督故障诊断方法
本专利技术涉及一种基于图论的轴承半监督故障诊断方法，属于机械零部件故障诊断

技术介绍
滚动轴承是使用最频繁的机械零部件之一。据统计，在所有机械故障中，超过40％的故障由轴承引起。在工业现场中，标签样本的收集是一件艰难的工作，尤其是稀缺的故障样本，但却存在着大量的无标签样本。针对这种只有少量标签数据且类别不平衡的情况，传统的故障诊断算法往往表现很差。此外，轴承一般运行在不确定的变工况下，这也会对故障诊断造成干扰。因此，通过利用无标签数据来提高故障诊断准确率有着十分重要的意义。可视图算法是一种将时序数据转化成复杂网络的方法，该方法将样本中的每一个点当成图中的节点，利用物理可视的原理构建复杂网络。提取网络结构参数作为轴承故障特征，可视图特征对变工况具有天然的不敏感性。半监督学习是一种能利用无标签数据来提高故障诊断准确率的分类算法，该方法只需要少量的标签样本数据，通过求解满足最优平滑假设目标方程的最优解将少量标签传播到所有无标签数据上，实现无标签样本的状态辨识。
技术实现思路
本专利技术提出了一种基于图论的轴承半监督故障诊断方法，首先将信号通过可视图算法转换成复杂网络，提取复杂网络的结构特征作为轴承故障特征，作为半监督分类器的输入，然后利用基于图的半监督学习算法将标签传播到无标签数据上实现故障诊断。本专利技术提出了一种基于图论的轴承半监督故障诊断方法，具体包括以下步骤：(一)首先利用可视图算法将通过传感器得到的轴承原始振动加速度信号转换成复杂网络；(二)计算复杂网络的结构参数，提取度分布的均值和标准差及网络复杂指数。(三)最后利用基于图的半监督学习处理无标签样本，实现轴承故障诊断。优选地，步骤(一)采用可视图算法具体为：针对一个具有N个数据点的样本，首先通过如下方法将振动加速度信号转化成可视图：信号的任意两点(ta,xa)，(tb,xb)，对应可视图中的两个节点，设这两点之间任意一点为(tc,xc)，两个节点相连接的条件满足得到形容每个节点的连接情况邻接矩阵WD，矩阵WD为N*N维的对称矩阵，wij＝1代表节点i与节点j相连；对于任意一点，与其连接点的个数称为该点的度则度分布为DV＝[d1,…,dn]。优选地，步骤(二)复杂网络的结构特征参数具体计算方法如下：(1)度分布均值：(2)度分布标准差：(3)图复杂指数GIC定义如下：C＝4c(1-c)其中，λmax表示可视图邻接矩阵的最大特征值。优选地，步骤(三)实现轴承故障诊断的具体方法如下：(1)图的构建首先，通过核函数计算所有样本之间的相似度，计算得到点xi和xj之间邻接矩阵K：K∈Rn×n，Kij＝k(xi,xj)，通过将邻接矩阵乘以二值矩阵B∈Bn×n和距离矩阵H∈Rn×n完成稀疏化：稀疏化邻接矩阵后，通过高斯核加权算法计算重新计算两点之间的权重矩阵W，点xi和xj之间的权重计算公式如下：其中，d(xi,xj)代表xi与xj之间的欧式距离，δ表示数据点的标准差；(2)标签传播采用平滑项和标签匹配项作为最终优化目标：s.t.yij∈{0,1},yij＝1,forlabel(xi)＝j,j＝1,…,c.其中，基于度的权重矩阵Λ＝diag([λ1,…,λn])，Y＝{yij}∈Βn×c为数据的标签矩阵，ΛY为规范化标签变量；其中，d为某样本的度，p是类别数据的先验分布且A＝PTLP+μ(PT-I)(P-I)＝PTLP+μ(P-I)2其中，P为传播矩阵,L为归一化的拉普拉斯矩阵；通过贪婪梯度Max-Cut算法求解上式，最终得到无标签样本的标签，实现故障诊断。优选地，所述贪婪梯度Max-Cut算法流程为:(1)输入：图GA＝{X,A}标签样本Xl和标签Y；(2)初始化：通过标签样本Xl构建初始割集{Sj}:Sj＝{xi|yij＝1},j＝1,2,…,c；无标签样本集：Xu＝X\Xl；(3)以下步骤循环，直至满足条件对于z＝0至|Xu|做如下计算：计算连接性:将无标签样本加上最优标签，并将其放入标签割集{Sj}；将xi放入标签样本Xl:Xl←Xl+xi；从无标签样本中删除该样本Xu:Xu←Xu-xi；(4)输出Sj，j＝1,2,...,c。其中样本xi与标签样本子集S的连接性定义如下：Sj＝{xi|yij＝1},i＝1,2,…,n；j＝1,2,…,c。本专利技术的优点与积极效果在于：(1)可视图算法是一种将时序数列转换成复杂网络的算法，该方法将轴承振动信号映射到图中，图的结构特征对信号的水平和垂直方向的尺度变换(对应轴承的转速和负载)具有天然的不敏感性。(2)利用改进规范化标签变量的能平衡类别标签影响的基于图的半监督轴承故障诊断算法处理类别不平衡无标签数据。为了解决数据类别不平衡带来的影响，利用双变量的目标函数，其中规范化标签变量用来平衡类别的影响。(3)相比于其他故障诊断方法，半监督故障诊断算法只需要少量的标签样本即能获得高的故障诊断准确率。(4)本专利技术方法属于数据驱动方法，无需建立模型，即可实现故障的检测和识别，降低了专业要求，增加了工程应用性；附图说明图1是本专利技术故障诊断方法的整体步骤流程图。图2是可视图算法示意图。图3是某故障样本振动信号图。图4是对应故障样本度分布。图5是轻微故障不同类别比率下的分类准确率曲线图。图6是常规故障不同类别比率下的分类准确率曲线图。图7是严重故障不同类别比率下的分类准确率曲线图。具体实施方式下面将结合附图和实例对本专利技术作进一步的详细说明。本专利技术提出了一种可视图算法将信号转换成复杂网络，通过提取可视图特征作为半监督学习算法的输入，从而实现故障诊断的方法。如图1所示，具体步骤如下：步骤一、将轴承振动信号通过可视图算法转换成复杂网络。针对一个具有N个数据点的样本，首先通过如下方法将信号转化成可视图。如图2为转换示意图。数学计算方式如下：信号的任意两点(ta,xa)，(tb,xb)，对应可视图中的两个节点，两个节点相连接的条件如下：设这两点之间任意一点为(tc,xc)，满足由此可得到形容每个节点的连接情况邻接矩阵WD，矩阵为N*N维的对称矩阵，wij＝1代表节点i与节点j相连。对于任意一点，与其连接点的个数称为该点的度则度分布为DV＝[d1,…,dn]。步骤二、计算复杂网络的结构特征参数复杂网络的基本特征包括度分布均值及标准差和图复杂指数。计算方法如下：(1)度分布均值：(2)度分布标准差：(3)图复杂指数GIC定义如下：C＝4c(1-c)其中，λmax表示可视图邻接矩阵的最大特征值。步骤三、用基于图的半监督学习算法，实现故障识别与诊断基于图的半监督学习算法原理包括图的构建和标签传播两部分如下：数学符号定义：{(x1,y1),…,(xl,yl)}和{xl+1,…,xl+u}分别代表标签样本和无标签样本，Xl＝{x1,…,xl}和Xu＝{xl+1,…,xl+u}分别代表有标签和无标签输入数据。l和u表示标签样本和无标签样本的数量。Y＝{yij}∈Βn×c为数据的标签矩阵，其中yij＝1代表样本xi的标签为j，j∈{1,2,…,c}，每一个样本有且只有一个标签，即半监督学习的目标是识别无标签数据的标签{yl+1,…,yl+u}，其中l<<n(l+u＝n)。(3)图的构建首本文档来自技高网...

【技术保护点】
1.一种基于图论的轴承半监督故障诊断方法，其特征在于，具体包括以下步骤：(一)首先利用可视图算法将通过传感器得到的轴承原始振动加速度信号转换成复杂网络；(二)计算复杂网络的结构参数，提取度分布的均值和标准差及网络复杂指数；(三)最后利用基于图的半监督学习处理无标签样本，实现轴承故障诊断。

【技术特征摘要】
1.一种基于图论的轴承半监督故障诊断方法，其特征在于，具体包括以下步骤：(一)首先利用可视图算法将通过传感器得到的轴承原始振动加速度信号转换成复杂网络；(二)计算复杂网络的结构参数，提取度分布的均值和标准差及网络复杂指数；(三)最后利用基于图的半监督学习处理无标签样本，实现轴承故障诊断。2.如权利要求1所述的方法，其特征在于，步骤(一)采用可视图算法具体为：针对一个具有N个数据点的样本，首先通过如下方法将振动加速度信号转化成可视图：信号的任意两点(ta,xa)，(tb,xb)，对应可视图中的两个节点，设这两点之间任意一点为(tc,xc)，两个节点相连接的条件满足得到形容每个节点的连接情况邻接矩阵WD，矩阵WD为N*N维的对称矩阵，wij＝1代表节点i与节点j相连；对于任意一点，与其连接点的个数称为该点的度则度分布为DV＝[d1,…,dn]。3.如权利要求2所述的方法，其特征在于，步骤(二)复杂网络的结构特征参数具体计算方法如下：(1)度分布均值：(2)度分布标准差：(3)图复杂指数GIC定义如下：C＝4c(1-c)其中，λmax表示可视图邻接矩阵的最大特征值。4.如权利要求3所述的方法，其特征在于，步骤(三)实现轴承故障诊断的具体方法如下：(1)图的构建首先，通过核函数计算所有样本之间的相似度，计算得到点xi和xj之间邻接矩阵K：K∈Rn×n，Kij＝k(xi,xj)，通过将邻接矩阵乘以二值矩阵B∈Bn×n和距离矩阵H∈Rn×n完成稀疏化：稀疏化邻接矩阵后，通过高斯...

【专利技术属性】
技术研发人员：王志鹏，陈欣安，贾利民，张蛰，秦勇，王宁，耿毅轩，
申请(专利权)人：北京交通大学，
类型：发明
国别省市：北京,11