CN108647707B - 概率神经网络创建方法、故障诊断方法及装置、存储介质 - Google Patents
概率神经网络创建方法、故障诊断方法及装置、存储介质 Download PDFInfo
- Publication number
- CN108647707B CN108647707B CN201810379172.2A CN201810379172A CN108647707B CN 108647707 B CN108647707 B CN 108647707B CN 201810379172 A CN201810379172 A CN 201810379172A CN 108647707 B CN108647707 B CN 108647707B
- Authority
- CN
- China
- Prior art keywords
- matrix
- data
- input
- training
- probability
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
- G01M13/04—Bearings
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Evolutionary Biology (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本申请公开了一种概率神经网络创建方法、故障诊断方法及装置、存储介质,所述概率神经网络创建方法包括:输入训练数据以及相对应的训练标签矩阵,生成概率神经网络模型的训练参数,所述训练参数为训练数据的属性概率矩阵;计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵,根据联合概率矩阵和训练标签矩阵计算各类别的求和概率;根据计算得出的各类别的求和概率,对输入测试数据进行分类判别,并获得所创建的概率神经网络的分类准确度。本申请创建的的概率神经网络,加快了学习速度,有效避免了陷入局部最优。
Description
技术领域
本发明涉及计算机技术领域,具体涉及一种概率神经网络创建方法、故障诊断方法及装置、存储介质。
背景技术
故障预测与健康管理(Prognostic and Health Management,PHM)是利用各类先进传感器实时监测装备运行状态参数及特征信号,借助智能算法和模型来评估装备健康状态,预测剩余使用寿命,诊断故障类型并在故障发生前提供一系列的故障维修决策的系统。PHM技术是由先进的诊断技术、测试技术装备维修管理理论相结合的产物。装备使用人员与维修人员借助该项技术的故障诊断能力来识别故障类型,从而采取有效的维修方式,有效降低故障风险,节约装备资源,减少因装备故障误诊带来的经济损失。
目前,PHM技术对于实施装备故障诊断的方法主要有:基于专家系统故障诊断方法,基于实例的诊断方法,基于故障树的诊断方法,基于人工神经网络故障诊断方法。基于专家系统故障诊断方法的缺陷在于:该方法存在知识经验获取的瓶颈,对于复杂系统所对应的症状之间的联系难以给予清楚的解释,此方法通常适用于单一故障诊断,难以用于多故障并存情况下的诊断;基于实例的诊断方法的缺陷在于:它是一种根据过去经验实例指导解决新问题的一种方法,能够搜索到的问题解决实例是有限的,不能覆盖搜索整个有解空间,可能会漏掉最优解。基于故障树的诊断方法的缺陷在于:不能诊断不可预知的故障,严重依赖于故障树信息的完整度。相比于上述几种方法,基于人工神经网络的故障诊断方法的优点在于:它有很强的非线性问题处理能力,以及分布式信息储存,并行处理,自适应学习等优点。
然而,目前已经实现的人工神经网络故障诊断缺陷在于:学习速度慢,简单的问题要学习上百次甚至上千次才能够收敛;容易陷入局部最优化,对于推断能力十分有限,网络模型参数易选取不当。
发明内容
为了解决上述技术问题,本发明提供了一种概率神经网络创建方法、故障诊断方法及装置、存储介质,能够加快学习速度并有效避免陷入局部最优。
为了达到本发明目的,本发明实施例的技术方案是这样实现的:
本发明实施例提供了一种概率神经网络创建方法,包括:
输入训练数据以及相对应的训练标签矩阵,生成概率神经网络模型的训练参数,所述训练参数为训练数据的属性概率矩阵;
计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵,根据联合概率矩阵和训练标签矩阵计算各类别的求和概率;
根据计算得出的各类别的求和概率,对输入测试数据进行分类判别,并获得所创建的概率神经网络的分类准确度。
进一步地,所述生成概率神经网络模型的训练参数,包括:
根据所述输入训练数据traindata生成所述输入训练数据的属性概率矩阵wtrain:其中,输入训练数据traindata为M*N的矩阵,M为输入训练数据的样本数目,N为特征属性个数,IM为M*1的单位矩阵;wtrain为N*M的概率矩阵;
根据所述输入训练数据traindata生成所述输入训练数据的训练标签矩阵B,其中B为M*C的矩阵,C为训练标签类别数目,所述标记训练标签矩阵B中每个训练样本序号所对应的标签类别值被标记为1。
进一步地,所述根据联合概率矩阵和训练标签矩阵计算各类别的求和概率,包括:
根据所述输入训练数据的训练标签矩阵BMj生成每类联合概率对应的标签标记矩阵d:d=BMj*I1,其中,d为M*P的单位矩阵,P为输入测试数据的样本数目,BMj为第j类类别的训练标签矩阵,j为1至C之间的自然数,I1为1*P的单位矩阵;
本发明实施例还提供了一种存储介质,其特征在于,所述存储介质上存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如以上任一项所述的概率神经网络创建方法的步骤。
本发明实施例还提供了一种故障诊断方法,包括:
获取各故障类别的传感器采集信号样本,对所述样本经过信号处理得到特征数据,将得到的特征数据划分为输入训练数据和输入测试数据;
根据输入训练数据以及相对应的训练标签矩阵,生成概率神经网络模型的训练参数,所述训练参数为训练数据的属性概率矩阵;
计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵,根据联合概率矩阵和训练标签矩阵计算各故障类别的求和概率;
根据计算出的各故障类别的求和概率,对输入测试数据进行分类判别,并获得所创建的概率神经网络的分类准确度;
利用创建的概率神经网络对待诊断的传感器采集信号的特征数据进行故障诊断。
进一步地,在将得到的特征数据划分为输入训练数据和输入测试数据之前,所述方法还包括:
对所述得到的特征数据进行加入白化处理的主要成分分析。
本发明实施例还提供了一种存储介质,所述存储介质上存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如以上任一项所述的故障诊断方法的步骤。
本发明实施例还提供了一种概率神经网络创建装置,包括输入单元、样本单元、求和单元和竞争单元,其中:
输入单元,用于输入训练数据以及相对应的训练标签矩阵,并输入测试数据;
样本单元,用于根据输入训练数据以及相对应的训练标签矩阵生成概率神经网络模型的训练参数,所述训练参数为输入训练数据的属性概率矩阵;计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵;
求和单元,用于根据联合概率矩阵和训练标签矩阵计算各类别的求和概率;
竞争单元,用于根据计算出的各类别的求和概率,得到输入测试数据被判别出的分类类别。
本发明实施例还提供了一种包括如上所述的概率神经网络创建装置的故障诊断装置,所述故障诊断装置还包括特征提取单元、划分单元、准确度评估单元和故障诊断单元,其中:
特征提取单元,用于获取各故障类别的传感器采集信号样本,对所述样本经过信号处理得到特征数据;
划分单元,用于将得到的特征数据划分为输入所述概率神经网络创建装置的输入训练数据和输入测试数据;
准确度评估单元,用于计算所创建的概率神经网络的分类准确度;
故障诊断单元,用于利用所述概率神经网络创建装置创建的概率神经网络对待诊断的传感器采集信号的特征数据进行故障诊断。
本发明的技术方案,具有如下有益效果:
本发明提供的概率神经网络创建方法、故障诊断方法及装置、存储介质,通过以输入训练数据的属性概率矩阵为训练参数,加快了学习速度,有效避免了陷入局部最优,并有较强的推断能力,有效避免了网络参数选取不当的问题,对于装备的故障诊断有很高的检测准确率。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为本发明实施例的一种创建概率神经网络的方法的流程示意图;
图2为本发明实施例创建的概率神经网络的分层结构示意图;
图3为本发明实施例的一种故障诊断方法的流程示意图;
图4为本发明实施例的对得到的特征数据进行加入白化处理的主要成分分析的流程示意图;
图5为本发明实施例的一种创建概率神经网络的装置的结构示意图;
图6为本发明实施例的一种故障诊断装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下文中将结合附图对本发明的实施例进行详细说明。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
如图1所示,根据本发明的一种创建概率神经网络的方法,包括如下步骤:
步骤101:输入训练数据以及相对应的训练标签矩阵,生成概率神经网络模型的训练参数,所述训练参数为训练数据的属性概率矩阵;
需要说明的是,如图2所示,所创建的概率神经网络分为四层,包括:输入层,样本层,求和层,竞争层。其中,各层的作用分别为:
输入层,用于接受输入数据;
样本层,用于数据的收集整理;
求和层,用于判断数据的类别;
竞争层,用于类别判断的结果数据输出。
进一步地,所述生成概率神经网络模型的训练参数,包括:
根据所述输入训练数据traindata生成所述输入训练数据的属性概率矩阵wtrain:其中,输入训练数据traindata为M*N的矩阵,M为输入训练数据的样本数目,N为特征属性个数,IM为M*1的单位矩阵;wtrain为N*M的概率矩阵;
根据所述输入训练数据traindata生成所述输入训练数据的训练标签矩阵B,其中B为M*C的矩阵,C为训练标签类别数目,所述标记训练标签矩阵B中每个训练样本序号所对应的标签类别值被标记为1。
步骤102:计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵,根据联合概率矩阵和训练标签矩阵计算各类别的求和概率;
具体地,所述计算输入测试数据的属性概率矩阵wtest为:
其中,输入测试数据testdata为P*N的矩阵,P为输入测试数据的样本数目,N为特征属性个数;IM是大小为M*1的单位矩阵,wtest为概率神经网络生成的输入测试数据的N*P的属性概率矩阵,wtest中的每个值代表的是每个输入样本中的每个属性的概率值。
将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵:
wnew=(wtrain)T*wtest;
其中wnew是大小为M*P的矩阵,其表示为输入训练数据各样本属性概率与输入测试数据各样本属性概率的联合概率。
进一步地,所述根据联合概率矩阵和训练标签矩阵计算各类别的求和概率,包括:
根据所述输入训练数据的训练标签矩阵BMj生成每类联合概率对应的标签标记矩阵d:d=BMj*I1,其中,d为M*P的单位矩阵,P为输入测试数据的样本数目,BMj为第j类类别的训练标签矩阵,j为1至C之间的自然数,I1为1*P的单位矩阵;
步骤103:根据计算得出的各类别的求和概率,得到输入测试数据被判别出的分类类别,并获得所创建的概率神经网络的分类准确度。
具体地,比较计算出的求和概率值并找出其中的最大值,其最大值所对应的分类类别即为最终该测试数据被判别出的分类类别。
进一步地,利用决策矩阵法获得所创建的概率神经网络的分类准确度。
例如,假设选用轴承加速度的信号采集数据作为原始数据,不同轴承状态包括四种:正常轴承,轴承内环故障、轴承外环故障、轴承滚珠故障。每一个轴承故障类别数据集取相同的样本量,经过概率神经网络的分类得到样本数据所对应的分类结果,假设判别结果矩阵为A,其中,Aij表示实际故障类型为第j种,判别故障类型为第i种的个数,在本实例中,i和j均为1至4之间的自然数,则各类故障诊断精确度以及整体检测精确度计算公式如下:
在本例中,数据特征参数选取为时域7个不同的特征参数,以及时-频域6个特征参数共计13个特征参数。在本例中,数据集不采用特征选择处理,直接将数据集合按照样本数量4:1的比例将其划分为训练集与测试集,即训练样本数量为320,测试样本数量为80。各种轴承故障类型的样本数量分别在训练集、测试集中均等。因此,本例的训练数据traindata大小为320*13的矩阵,测试数据testdata大小为80*13的矩阵。
在本发明一实施例中,假设平滑参数σ取值为0.01,测试分类结果数量为80。计算各个故障种类的分类结果,实验结果如表1所示:
表1
在本例中,各类故障诊断精确度如下:
正常检测准确度:97.6%
内环故障检测准确度:97.4%
外环故障检测准确度:100%
滚珠故障检测准确度:100%
整体检测准确度:98.75%
经检测该,PHM概率神经网络算法对于本实施例中的故障数据诊断正确率为98.75%。
各模型在相同样本下的诊断结果统计和比较,如表2:
表2
根据各类模型在相同数据下轴承故障诊断检测准确度的比较下,概率神经网络模型相比于另外两种模型在故障诊断分类效果更佳。
本发明还公开了一种存储介质,所述存储介质上存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如以上任一项所述的概率神经网络创建方法的步骤。
如图3所示,根据本发明的一种故障诊断方法,包括如下步骤:
步骤301:获取各故障类别的传感器采集信号样本,对所述样本经过信号处理得到特征数据(这里使用信号处理中的特定的特征提取方法得到均方根,波峰因子,脉冲因子等特征组成的特征数据),将得到的特征数据划分为输入训练数据和输入测试数据;
在本发明一实施例中,选用轴承加速度的传感器采集信号作为原始数据,将不同的轴承状态分类为:正常轴承,轴承内环故障、轴承外环故障、轴承滚珠故障。因此,将获取的原始数据分为正常轴承数据集,轴承内环故障数据集、轴承外环故障数据集、轴承滚珠故障数据集,(原始数据集中含有上述四种数据集)每一个轴承故障类别数据集取相同的样本量。
以上四种轴承数据集合经过信号处理后的特征提取得到时域,时-频域特征参数数据。其中时域特征参数包括:均方根、歪度因子、峭度因子、波峰因子、裕度因子、波形因子、脉冲因子。时-频域特征参数为:信号经过经验模态分解得出的IMF本征模函数值。
根据上述特征参数数据组建成用于诊断的数据集。
进一步地,在将得到的特征数据划分为输入训练数据和输入测试数据之前,所述方法还包括:
对所述得到的特征数据进行加入白化处理的主要成分分析。
具体地,如图4所述,所述对所述得到的特征数据进行加入白化处理的主要成分分析,具体包括:
计算得到的特征数据矩阵X各列的平均值,将特征数据矩阵X中的每个数值减去该数所在列的列平均值;
计算该特征数据矩阵X的特征协方差矩阵C,采用公式:C=X*X′;
对计算出的特征协方差矩阵C进行奇异值分解,采用公式:C=U*D0*U′其中U为酉矩阵,D0为特征值对角矩阵;
计算白化后的数据矩阵Z,采用公式:Z=M*X;
计算该数据矩阵Z的特征协方差矩阵C,并求出特征协方差的特征值D与特征向量V;
将特征值的对角矩阵D转换为列向量Dn,并将其进行降序排列;
计算列向量特征值总和,并计算累加每次列向量总和与列向量特征值总和的比值,进而与预先设置的信息贡献度h比较大小;
若第K次累加列向量总和与列向量特征值总和大于预先设置的信息贡献度h,则将其所对应的K个特征向量作为新的特征向量;反之,继续增加累加次数,直到累加列向量总和与列向量特征值总和的比值大于预先设置的信息贡献度h;
将输入数据矩阵Z投影到新的特征向量上,产生降低维数后的新数据矩阵Xnew。
进一步地,所述在将得到的特征数据划分为输入训练数据和输入测试数据时,所述输入训练数据和输入测试数据的分配比例为4:1到2:1之间。
需要说明的是,本发明采用留出法将整个数据集合划分为两个互斥的集合,其中一个集合作为输入训练数据,另外一个作为输入测试数据。在输入训练数据训练出模型参数后,再用输入测试数据来评估其测试误差。若输入训练数据过大接近整个数据集,输入测试数据过小评估结果会不够准确。反之,降低了评估结果的保真性。
步骤302:根据输入训练数据以及相对应的训练标签矩阵,生成概率神经网络模型的训练参数,所述训练参数为训练数据的属性概率矩阵;
需要说明的是,如图2所示,所创建的概率神经网络分为四层,包括:输入层,样本层,求和层,竞争层。其中,各层的作用分别为:
输入层,用于接受输入数据;
样本层,用于数据的收集整理;
求和层,用于判断数据的类别;
竞争层,用于类别判断的结果数据输出。
进一步地,所述生成概率神经网络模型的训练参数,包括:
根据所述输入训练数据traindata生成所述输入训练数据的属性概率矩阵wtrain:其中,输入训练数据traindata为M*N的矩阵,M为输入训练数据的样本数目,N为特征属性个数,IM为M*1的单位矩阵;wtrain为N*M的概率矩阵;
根据所述输入训练数据traindata生成所述输入训练数据的训练标签矩阵B,其中B为M*C的矩阵,C为训练标签类别数目,所述标记训练标签矩阵B中每个训练样本序号所对应的标签类别值被标记为1。
步骤303:计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵,根据联合概率矩阵和训练标签矩阵计算各故障类别的求和概率;
具体地,所述计算输入测试数据的属性概率矩阵wtest为:
其中,输入测试数据testdata为P*N的矩阵,P为输入测试数据的样本数目,N为特征属性个数;IM是大小为M*1的单位矩阵,wtest为概率神经网络生成的输入测试数据的N*P的属性概率矩阵,wtest中的每个值代表的是每个输入样本中的每个属性的概率值。
将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵:
wnew=(wtrain)T*wtest;
其中wnew是大小为M*P的矩阵,其表示为输入训练数据各样本属性概率与输入测试数据各样本属性概率的联合概率。
进一步地,所述根据联合概率矩阵和训练标签矩阵计算各类别的求和概率,包括:
根据所述输入训练数据的训练标签矩阵BMj生成每类联合概率对应的标签标记矩阵d:d=BMj*I1,其中,d为M*P的单位矩阵,P为输入测试数据的样本数目,BMj为第j类类别的训练标签矩阵,j为1至C之间的自然数,I1为1*P的单位矩阵;
步骤304:根据计算出的各类别的求和概率,对输入测试数据进行分类判别,并获得所创建的概率神经网络的分类准确度;
具体地,比较计算出的求和概率值并找出其中的最大值,其最大值所对应的分类类别即为最终该测试数据被判别出的分类类别。
进一步地,利用决策矩阵法获得所创建的概率神经网络的分类准确度。
例如,假设选用轴承加速度的信号采集数据作为原始数据,不同轴承状态包括四种:正常轴承,轴承内环故障、轴承外环故障、轴承滚珠故障。每一个轴承故障类别数据集取相同的样本量,假设判别结果矩阵为A,其中,Aij表示实际故障类型为第j种,判别故障类型为第i种的个数,在本实例中,i和j均为1至4之间的自然数,则各类故障诊断精确度以及整体检测精确度计算公式如下:
步骤305:利用创建的概率神经网络对待诊断的传感器采集信号的特征数据进行故障诊断。
需要说明的是,至于如何利用创建的概率神经网络对待诊断的传感器采集信号的特征数据进行故障诊断,其具体过程与对所述输入测试数据进行分类类别判别的过程原理相同,此处不再赘述。
在本发明另一实施例中,采用与上述表1实施例相同的实验数据,再划分训练集与测试集前,先将数据进行特征选择。数据集划分比例与实验1相同。平滑参数设置也相同。实验结果如表3:
表3
在本例中,各类故障诊断精确度如下:
正常检测准确度:97.4%
内环故障检测准确度:95%
外环故障检测准确度:100%
滚珠故障检测准确度:97.6%
整体检测准确度:97.5%
经检测该PHM概率神经网络算法对于本例中的故障数据诊断正确率为97.5%。
本发明还公开了一种存储介质,所述存储介质上存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如以上任一项所述的故障诊断方法的步骤。
如图5所示,本发明还公开了一种概率神经网络创建装置,包括输入单元501、样本单元502、求和单元503和竞争单元504,其中:
输入单元501,用于输入训练数据以及相对应的训练标签矩阵,并输入测试数据;
样本单元502,用于根据输入训练数据以及相对应的训练标签矩阵生成概率神经网络模型的训练参数,所述训练参数为输入训练数据的属性概率矩阵阵;计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵;
求和单元503,用于根据联合概率矩阵和训练标签矩阵计算各类别的求和概率;
竞争单元504,用于根据计算出的各类别的求和概率,得到输入测试数据被判别出的分类类别。
进一步地,所述样本单元502的生成概率神经网络模型的训练参数,包括:
根据所述输入训练数据traindata生成所述输入训练数据的属性概率矩阵wtrain:其中,输入训练数据traindata为M*N的矩阵,M为输入训练数据的样本数目,N为特征属性个数,IM为M*1的单位矩阵;wtrain为N*M的概率矩阵;
根据所述输入训练数据traindata生成所述输入训练数据的训练标签矩阵B,其中B为M*C的矩阵,C为训练标签类别数目,所述标记训练标签矩阵B中每个训练样本序号所对应的标签类别值被标记为1。
所述样本单元502的计算输入测试数据的属性概率矩阵wtest为:
其中,输入测试数据testdata为P*N的矩阵,P为输入测试数据的样本数目,N为特征属性个数;IM是大小为M*1的单位矩阵,wtest为概率神经网络生成的输入测试数据的N*P的属性概率矩阵,wtest中的每个值代表的是每个输入样本中的每个属性的概率值。
将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵:
wnew=(wtrain)T*wtest;
其中wnew是大小为M*P的矩阵,其表示为输入训练数据各样本属性概率与输入测试数据各样本属性概率的联合概率。
进一步地,所述求和单元503的根据联合概率矩阵和训练标签矩阵计算各类别的求和概率,包括:
根据所述输入训练数据的训练标签矩阵BMj生成每类联合概率对应的标签标记矩阵d:d=BMj*I1,其中,d为M*P的单位矩阵,P为输入测试数据的样本数目,BMj为第j类类别的训练标签矩阵,j为1至C之间的自然数,I1为1*P的单位矩阵;
需要说明的是,所述竞争单元504比较计算出的求和概率值并找出其中的最大值,其最大值所对应的分类类别即为最终该测试数据被判别出的分类类别。
进一步地,在竞争单元504获得判别分类的结果后,可以利用决策矩阵法获得所创建的概率神经网络的分类准确度。
例如,假设选用轴承加速度的信号采集数据作为原始数据,不同轴承状态包括四种:正常轴承,轴承内环故障、轴承外环故障、轴承滚珠故障。每一个轴承故障类别数据集取相同的样本量,假设判别结果矩阵为A,其中,Aij表示实际故障类型为第j种,判别故障类型为第i种的个数,在本实例中,i和j均为1至4之间的自然数,则各类故障诊断精确度以及整体检测精确度计算公式如下:
如图6所示,本发明还公开了一种故障诊断装置,包括特征提取单元601、划分单元602、概率神经网络创建单元603、准确度评估单元604和故障诊断单元605,其中:
特征提取单元601,用于获取各故障类别的传感器采集信号样本,对所述样本经过信号处理得到特征数据(这里使用信号处理中的特定的特征提取方法得到均方根,波峰因子,脉冲因子等特征组成的特征数据);
划分单元602,用于将得到的特征数据划分为输入训练数据和输入测试数据;
概率神经网络创建单元603,用于根据输入训练数据以及相对应的训练标签矩阵,生成概率神经网络模型的训练参数,所述训练参数为训练数据的属性概率矩阵;计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵,根据联合概率矩阵和训练标签矩阵计算各故障类别的求和概率;根据计算出的求和概率,对输入测试数据进行分类判别;
准确度评估单元604,用于计算所创建的概率神经网络的分类准确度;
故障诊断单元605,用于输入测试数据,利用创建的概率神经网络进行故障诊断。
在本发明一实施例中,选用轴承加速度的传感器采集信号作为原始数据,将不同的轴承状态分类为:正常轴承,轴承内环故障、轴承外环故障、轴承滚珠故障。因此,将获取的原始数据分为正常轴承数据集,轴承内环故障数据集、轴承外环故障数据集、轴承滚珠故障数据集,每一个轴承故障类别数据集取相同的样本量。
以上四种轴承数据集合经过信号处理后的特征提取得到时域,时-频域特征参数数据。其中时域特征参数包括:均方根、歪度因子、峭度因子、波峰因子、裕度因子、波形因子、脉冲因子。时-频域特征参数为:信号经过经验模态分解得出的IMF本征模函数值。
根据上述特征参数数据组建成用于诊断的数据集。
进一步地,所述特征提取单元601还用于:
对所述得到的特征数据进行加入白化处理的主要成分分析。
具体地,所述特征提取单元601的对所述得到的特征数据进行加入白化处理的主要成分分析,包括:
计算得到的特征数据矩阵X各列的平均值,将特征数据矩阵X中的每个数值减去该数所在列的列平均值;
计算该特征数据矩阵X的特征协方差矩阵C,采用公式:C=X*X′;
对计算出的特征协方差矩阵C进行奇异值分解,采用公式:C=U*D0*U′其中U为酉矩阵,D0为特征值对角矩阵;
计算白化后的数据矩阵Z,采用公式:Z=M*X;
计算该数据矩阵Z的特征协方差矩阵C,并求出特征协方差的特征值D与特征向量V;
将特征值的对角矩阵D转换为列向量Dn,并将其进行降序排列;
计算列向量特征值总和,并计算累加每次列向量总和与列向量特征值总和的比值,进而与预先设置的信息贡献度h比较大小;
若第K次累加列向量总和与列向量特征值总和大于预先设置的信息贡献度h,则将其所对应的K个特征向量作为新的特征向量;反之,继续增加累加次数,直到累加列向量总和与列向量特征值总和的比值大于预先设置的信息贡献度h;
将输入数据矩阵Z投影到新的特征向量上,产生降低维数后的新数据矩阵Xnew。
进一步地,所述划分单元602在将得到的特征数据划分为输入训练数据和输入测试数据时,所述输入训练数据和输入测试数据的分配比例为4:1到2:1之间。
需要说明的是,本发明采用留出法将整个数据集合划分为两个互斥的集合,其中一个集合作为输入训练数据,另外一个作为输入测试数据。在输入训练数据训练出模型参数后,再用输入测试数据来评估其测试误差。若输入训练数据过大接近整个数据集,输入测试数据过小评估结果会不够准确。反之,降低了评估结果的保真性。
进一步地,所述概率神经网络创建单元603的生成概率神经网络模型的训练参数,包括:
根据所述输入训练数据traindata生成所述输入训练数据的属性概率矩阵wtrain:其中,输入训练数据traindata为M*N的矩阵,M为输入训练数据的样本数目,N为特征属性个数,IM为M*1的单位矩阵;wtrain为N*M的概率矩阵;
根据所述输入训练数据traindata生成所述输入训练数据的训练标签矩阵B,其中B为M*C的矩阵,C为训练标签类别数目,所述标记训练标签矩阵B中每个训练样本序号所对应的标签类别值被标记为1。
所述概率神经网络创建单元603的计算输入测试数据的属性概率矩阵wtest为:
其中,输入测试数据testdata为P*N的矩阵,P为输入测试数据的样本数目,N为特征属性个数;IM是大小为M*1的单位矩阵,wtest为概率神经网络生成的输入测试数据的N*P的属性概率矩阵,wtest中的每个值代表的是每个输入样本中的每个属性的概率值。
所述概率神经网络创建单元603将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵:
wnew=(wtrain)T*wtest;
其中wnew是大小为M*P的矩阵,其表示为输入训练数据各样本属性概率与输入测试数据各样本属性概率的联合概率。
进一步地,所述概率神经网络创建单元603的根据联合概率矩阵和训练标签矩阵计算各类别的求和概率,包括:
根据所述输入训练数据的训练标签矩阵BMj生成每类联合概率对应的标签标记矩阵d:d=BMj*I1,其中,d为M*P的单位矩阵,P为输入测试数据的样本数目,BMj为第j类类别的训练标签矩阵,j为1至C之间的自然数,I1为1*P的单位矩阵;
所述概率神经网络创建单元603比较计算出的求和概率值并找出其中的最大值,其最大值所对应的分类类别即为最终该测试数据被判别出的分类类别。
进一步地,所述准确度评估单元604利用决策矩阵法获得所创建的概率神经网络的分类准确度。
需要说明的是,至于故障诊断单元605如何利用创建的概率神经网络对待诊断的传感器采集信号的特征数据进行故障诊断,其具体过程与所述概率神经网络创建单元603对所述输入测试数据进行分类类别判别的过程原理相同,此处不再赘述。
本领域普通技术人员可以理解上述方法中的全部或部分步骤可通过程序来指令相关硬件完成,所述程序可以存储于计算机可读存储介质中,如只读存储器、磁盘或光盘等。可选地,上述实施例的全部或部分步骤也可以使用一个或多个集成电路来实现,相应地,上述实施例中的各模块/单元可以采用硬件的形式实现,也可以采用软件功能模块的形式实现。本发明不限制于任何特定形式的硬件和软件的结合。
尽管上文对本发明进行了详细说明,但是以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种故障诊断方法,其特征在于,包括:
获取各故障类别的传感器采集信号样本,对所述样本经过信号处理得到特征数据,将得到的特征数据划分为输入训练数据和输入测试数据;
根据输入训练数据以及相对应的训练标签矩阵,生成概率神经网络模型的训练参数,所述训练参数为训练数据的属性概率矩阵;
计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵,根据联合概率矩阵和训练标签矩阵计算各类别的求和概率;
根据计算得出的各类别的求和概率,对输入测试数据进行分类判别,并获得所创建的概率神经网络的分类准确度;
利用创建的概率神经网络对待诊断的传感器采集信号的特征数据进行故障诊断;
其中,所述生成概率神经网络模型的训练参数,包括:
根据所述输入训练数据traindata生成所述输入训练数据的属性概率矩阵wtrain:其中,输入训练数据traindata为M*N的矩阵,M为输入训练数据的样本数目,N为特征属性个数,IM为M*1的单位矩阵;wtrain为N*M的概率矩阵;
根据所述输入训练数据traindata生成所述输入训练数据的训练标签矩阵B,其中B为M*C的矩阵,C为训练标签类别数目,所述训练标签矩阵B中每个训练样本序号所对应的标签类别值被标记为1。
3.根据权利要求1所述的故障诊断方法,其特征在于,在所述将得到的特征数据划分为输入训练数据和输入测试数据之前,所述方法还包括:
对所述得到的特征数据进行加入白化处理的主要成分分析。
4.一种存储介质,其特征在于,所述存储介质上存储有一个或者多个程序,所述一个或者多个程序可被一个或者多个处理器执行,以实现如权利要求1至3中任一项所述的故障诊断方法的步骤。
5.一种故障诊断装置,其特征在于,包括输入单元、样本单元、求和单元、竞争单元、特征提取单元、划分单元、准确度评估单元和故障诊断单元,其中:
输入单元,用于输入训练数据以及相对应的训练标签矩阵,并输入测试数据;
样本单元,用于根据输入训练数据以及相对应的训练标签矩阵生成概率神经网络模型的训练参数,所述训练参数为输入训练数据的属性概率矩阵;计算输入测试数据的属性概率矩阵,将输入训练数据和输入测试数据的属性概率矩阵相乘,得到输入训练数据和输入测试数据的联合概率矩阵;
求和单元,用于根据联合概率矩阵和训练标签矩阵计算各类别的求和概率;
竞争单元,用于根据计算出的各类别的求和概率,得到输入测试数据被判别出的分类类别;
特征提取单元,用于获取各故障类别的传感器采集信号样本,对所述样本经过信号处理得到特征数据;
划分单元,用于将得到的特征数据划分为输入所述概率神经网络创建装置的输入训练数据和输入测试数据;
准确度评估单元,用于计算所创建的概率神经网络的分类准确度;
故障诊断单元,用于利用所述概率神经网络创建装置创建的概率神经网络对待诊断的传感器采集信号的特征数据进行故障诊断;
其中,所述生成概率神经网络模型的训练参数,包括:
根据所述输入训练数据traindata生成所述输入训练数据的属性概率矩阵wtrain:其中,输入训练数据traindata为M*N的矩阵,M为输入训练数据的样本数目,N为特征属性个数,IM为M*1的单位矩阵;wtrain为N*M的概率矩阵;
根据所述输入训练数据traindata生成所述输入训练数据的训练标签矩阵B,其中B为M*C的矩阵,C为训练标签类别数目,所述训练标签矩阵B中每个训练样本序号所对应的标签类别值被标记为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810379172.2A CN108647707B (zh) | 2018-04-25 | 2018-04-25 | 概率神经网络创建方法、故障诊断方法及装置、存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810379172.2A CN108647707B (zh) | 2018-04-25 | 2018-04-25 | 概率神经网络创建方法、故障诊断方法及装置、存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108647707A CN108647707A (zh) | 2018-10-12 |
CN108647707B true CN108647707B (zh) | 2022-09-09 |
Family
ID=63747539
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810379172.2A Active CN108647707B (zh) | 2018-04-25 | 2018-04-25 | 概率神经网络创建方法、故障诊断方法及装置、存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108647707B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111160454B (zh) * | 2019-12-27 | 2023-12-22 | 北京旋极信息技术股份有限公司 | 一种速变信号检测方法和装置 |
CN112182214B (zh) * | 2020-09-27 | 2024-03-19 | 中国建设银行股份有限公司 | 一种数据分类方法、装置、设备及介质 |
CN112417766B (zh) * | 2020-12-03 | 2024-05-14 | 深制科技(苏州)有限公司 | 一种以无故障数据为主的故障诊断方法 |
CN112766327B (zh) * | 2021-01-05 | 2024-05-24 | 格力电器(武汉)有限公司 | 空调故障预测方法、电子设备和储存介质 |
CN113240000B (zh) * | 2021-05-10 | 2022-08-23 | 北京航空航天大学 | 机器状态监测方法、可读存储介质及电子设备 |
CN113420813B (zh) * | 2021-06-23 | 2023-11-28 | 北京市机械工业局技术开发研究所 | 一种车辆尾气检测设备颗粒物过滤棉状态的诊断方法 |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5566092A (en) * | 1993-12-30 | 1996-10-15 | Caterpillar Inc. | Machine fault diagnostics system and method |
US6076083A (en) * | 1995-08-20 | 2000-06-13 | Baker; Michelle | Diagnostic system utilizing a Bayesian network model having link weights updated experimentally |
CN101984340A (zh) * | 2010-10-28 | 2011-03-09 | 浙江工业大学 | 一种基于人工智能的软性磨粒两相湍流流型识别方法 |
CN102254193A (zh) * | 2011-07-16 | 2011-11-23 | 西安电子科技大学 | 基于相关向量机的多类数据分类方法 |
CN102622510A (zh) * | 2012-01-31 | 2012-08-01 | 龚波 | 一种软件缺陷量化管理系统和方法 |
CN103822793A (zh) * | 2014-01-20 | 2014-05-28 | 北京邮电大学 | 一种复杂设备声学故障识别定位方法 |
CN104200077A (zh) * | 2014-08-22 | 2014-12-10 | 广西师范大学 | 基于子空间学习的嵌入式属性选择方法及其应用 |
CN104268381A (zh) * | 2014-09-16 | 2015-01-07 | 哈尔滨工业大学 | 一种基于AdaBoost算法的卫星故障诊断方法 |
CN104712542A (zh) * | 2015-01-12 | 2015-06-17 | 北京博华信智科技股份有限公司 | 一种基于物联网的往复压缩机敏感特征提取与故障诊断方法 |
CN106649696A (zh) * | 2016-12-19 | 2017-05-10 | 北京云知声信息技术有限公司 | 信息分类方法及装置 |
CN106951197A (zh) * | 2017-01-19 | 2017-07-14 | 北京印刷学院 | 一种基于信号测试的印刷设备远程故障诊断方法及系统 |
CN107256546A (zh) * | 2017-05-23 | 2017-10-17 | 上海海事大学 | 基于pca卷积池化softmax的海流机叶片附着物故障诊断方法 |
CN107609063A (zh) * | 2017-08-29 | 2018-01-19 | 重庆邮电大学 | 一种多标签分类的手机应用推荐系统及其方法 |
CN107831013A (zh) * | 2017-10-11 | 2018-03-23 | 温州大学 | 一种利用概率主分量分析增强循环双谱的轴承故障诊断方法 |
CN107886123A (zh) * | 2017-11-08 | 2018-04-06 | 电子科技大学 | 一种基于辅助判决更新学习的合成孔径雷达目标识别方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050010444A1 (en) * | 2003-06-06 | 2005-01-13 | Iliff Edwin C. | System and method for assisting medical diagnosis using an anatomic system and cause matrix |
WO2007055222A1 (ja) * | 2005-11-08 | 2007-05-18 | Tohoku University | ネットワーク異常検知方法およびネットワーク異常検知システム |
US8943011B2 (en) * | 2011-06-28 | 2015-01-27 | Salesforce.Com, Inc. | Methods and systems for using map-reduce for large-scale analysis of graph-based data |
CN103777123A (zh) * | 2014-01-27 | 2014-05-07 | 国家电网公司 | 用于gis设备的局部放电故障综合诊断方法 |
CN103886405B (zh) * | 2014-02-20 | 2017-02-22 | 东南大学 | 一种基于信息熵特征和概率神经网络的锅炉燃烧工况识别方法 |
CN106570513B (zh) * | 2015-10-13 | 2019-09-13 | 华为技术有限公司 | 大数据网络系统的故障诊断方法和装置 |
CN105426839A (zh) * | 2015-11-18 | 2016-03-23 | 清华大学 | 基于稀疏自动编码器的电力系统过电压分类方法 |
CN106408687B (zh) * | 2016-11-24 | 2019-04-05 | 沈阳航空航天大学 | 一种基于机器学习方法的汽车发动机故障预警方法 |
-
2018
- 2018-04-25 CN CN201810379172.2A patent/CN108647707B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5566092A (en) * | 1993-12-30 | 1996-10-15 | Caterpillar Inc. | Machine fault diagnostics system and method |
US6076083A (en) * | 1995-08-20 | 2000-06-13 | Baker; Michelle | Diagnostic system utilizing a Bayesian network model having link weights updated experimentally |
CN101984340A (zh) * | 2010-10-28 | 2011-03-09 | 浙江工业大学 | 一种基于人工智能的软性磨粒两相湍流流型识别方法 |
CN102254193A (zh) * | 2011-07-16 | 2011-11-23 | 西安电子科技大学 | 基于相关向量机的多类数据分类方法 |
CN102622510A (zh) * | 2012-01-31 | 2012-08-01 | 龚波 | 一种软件缺陷量化管理系统和方法 |
CN103822793A (zh) * | 2014-01-20 | 2014-05-28 | 北京邮电大学 | 一种复杂设备声学故障识别定位方法 |
CN104200077A (zh) * | 2014-08-22 | 2014-12-10 | 广西师范大学 | 基于子空间学习的嵌入式属性选择方法及其应用 |
CN104268381A (zh) * | 2014-09-16 | 2015-01-07 | 哈尔滨工业大学 | 一种基于AdaBoost算法的卫星故障诊断方法 |
CN104712542A (zh) * | 2015-01-12 | 2015-06-17 | 北京博华信智科技股份有限公司 | 一种基于物联网的往复压缩机敏感特征提取与故障诊断方法 |
CN106649696A (zh) * | 2016-12-19 | 2017-05-10 | 北京云知声信息技术有限公司 | 信息分类方法及装置 |
CN106951197A (zh) * | 2017-01-19 | 2017-07-14 | 北京印刷学院 | 一种基于信号测试的印刷设备远程故障诊断方法及系统 |
CN107256546A (zh) * | 2017-05-23 | 2017-10-17 | 上海海事大学 | 基于pca卷积池化softmax的海流机叶片附着物故障诊断方法 |
CN107609063A (zh) * | 2017-08-29 | 2018-01-19 | 重庆邮电大学 | 一种多标签分类的手机应用推荐系统及其方法 |
CN107831013A (zh) * | 2017-10-11 | 2018-03-23 | 温州大学 | 一种利用概率主分量分析增强循环双谱的轴承故障诊断方法 |
CN107886123A (zh) * | 2017-11-08 | 2018-04-06 | 电子科技大学 | 一种基于辅助判决更新学习的合成孔径雷达目标识别方法 |
Non-Patent Citations (6)
Title |
---|
基于LSSVM和PNN的车轮状态安全域估计及故障诊断;冯坚强等;《机械制造与自动化》;20170220(第01期);第146-150页 * |
基于双层长短时记忆网络的齿轮故障诊断方法;王维锋等;《装甲兵工程学院学报》;20180415(第02期);第86-90页 * |
基于支持向量机的旋转机械振动故障诊断研究;于波等;《电子设计工程》;20160320(第06期);第110-113页 * |
基于改进多分类概率SVM模型的变压器故障诊断;彭刚等;《机械与电子》;20180424(第04期);第44-49页 * |
基于概率神经网络的设备故障诊断及仿真分析;姬东朝等;《火力与指挥控制》;20090115(第01期);第84-87页 * |
概率神经网络在往复压缩机多故障诊断技术研究中的应用;彭琦等;《机械设计与制造》;20160408(第04期);第246-249页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108647707A (zh) | 2018-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108647707B (zh) | 概率神经网络创建方法、故障诊断方法及装置、存储介质 | |
CN107066759B (zh) | 一种汽轮机转子振动故障诊断方法及装置 | |
CN111858231B (zh) | 一种基于运维监控的单指标异常检测方法 | |
CN106570513B (zh) | 大数据网络系统的故障诊断方法和装置 | |
KR101964412B1 (ko) | 이동통신데이터 처리시스템의 이상로그 발생을 진단하는 방법 및 그 시스템 | |
CN109948647A (zh) | 一种基于深度残差网络的心电图分类方法及系统 | |
CN112859822A (zh) | 基于人工智能的设备健康分析及故障诊断的方法及系统 | |
CN113255848A (zh) | 基于大数据学习的水轮机空化声信号辨识方法 | |
CN112906764B (zh) | 基于改进bp神经网络的通信安全设备智能诊断方法及系统 | |
CN111126820A (zh) | 反窃电方法及系统 | |
Tang et al. | Simulation-based diagnostic model for automatic testability analysis of analog circuits | |
Cheng et al. | Online power system event detection via bidirectional generative adversarial networks | |
CN111259949A (zh) | 飞机环控系统故障识别模型构建方法、模型及识别方法 | |
CN111368911A (zh) | 一种图像分类方法、装置和计算机可读存储介质 | |
JP2004152306A (ja) | 双峰性データを判別する方法 | |
CN113177357B (zh) | 一种电力系统暂态稳定评估方法 | |
CN115712871A (zh) | 一种结合重采样和集成学习的电力电子系统故障诊断方法 | |
Tan et al. | Fault diagnosis method for rolling bearing based on VMD and improved SVM optimized by METLBO | |
CN113935413A (zh) | 一种基于卷积神经网的配网录波文件波形识别方法 | |
CN115878992A (zh) | 综合管廊供电系统的监控方法和监控系统 | |
CN116776245A (zh) | 一种基于机器学习的三相逆变器设备故障诊断方法 | |
Sim et al. | Attention masking for improved near out-of-distribution image detection | |
Razavi-Far et al. | Imputation of missing data for diagnosing sensor faults in a wind turbine | |
CN114528906A (zh) | 一种旋转机械的故障诊断方法、装置、设备和介质 | |
CN114565051A (zh) | 基于神经元影响程度的产品分类模型的测试方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |