CN104751229A - 一种改进bp神经网络估值的缺失数据轴承故障诊断方法 - Google Patents
一种改进bp神经网络估值的缺失数据轴承故障诊断方法 Download PDFInfo
- Publication number
- CN104751229A CN104751229A CN201510172600.0A CN201510172600A CN104751229A CN 104751229 A CN104751229 A CN 104751229A CN 201510172600 A CN201510172600 A CN 201510172600A CN 104751229 A CN104751229 A CN 104751229A
- Authority
- CN
- China
- Prior art keywords
- data
- sample
- hidden layer
- attribute
- layer
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 24
- 238000003745 diagnosis Methods 0.000 title claims abstract description 15
- 238000012549 training Methods 0.000 claims abstract description 54
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 230000008034 disappearance Effects 0.000 claims description 44
- 239000011159 matrix material Substances 0.000 claims description 35
- 238000005096 rolling process Methods 0.000 claims description 22
- 230000006870 function Effects 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 12
- 230000005284 excitation Effects 0.000 claims description 9
- 238000013480 data collection Methods 0.000 claims description 7
- 210000002569 neuron Anatomy 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 238000007621 cluster analysis Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 230000036541 health Effects 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000007781 pre-processing Methods 0.000 abstract 1
- 238000001574 biopsy Methods 0.000 description 4
- 210000004072 lung Anatomy 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Landscapes
- Complex Calculations (AREA)
Abstract
一种改进BP神经网络估值的缺失数据轴承故障诊断方法,步骤如下:1)轴承数据预处理;2)确定和优化训练样本;3)初始化网络;4)改进后基于缺失数据的BP网络的训练;5)对缺失属性进行估值;6)对数据集进行聚类分析,运用本发明方法能够对存在缺失的轴承数据进行处理,采用模糊C均值聚类算法对恢复后的完整数据集进行聚类,从而对轴承健康度进行评价。
Description
技术领域
本发明涉及一种改进BP神经网络估值的缺失数据轴承故障诊断方法。
背景技术
现代生产中,滚动轴承在旋转机械中被广泛使用,滚动轴承的健康状况是对整个机械运行的重要影响之一。滚动轴承需要具有较高的可靠性,机械运行中轴承故障的发生可能会导致致命的机械故障。因此,对滚动轴承健康度的评价技术是极其重要的。
近年来,健康度评价技术迅速发展,研究成果不断出现,采用的方法也是多种多样,应用最为广泛的是模糊C均值聚类算法。在工业实际生产中,由于采集设备的精度限制,噪声的影响或数据漏读等多方面原因造成轴承采集数据缺失而产生不完整数据集。但是模糊C均值聚类算法不能直接对不完整数据进行聚类。目前,还没有轴承缺失数据集健康度评价方法。
发明内容
为了解决上述存在的技术问题,本发明提供种改进BP神经网络估值的缺失数据轴承故障诊断方法,将不完整的缺失数据通过局部距离公式处理得到训练样本集,在针对得到的训练样本集使用改进后的BP网络进行训练,从而得到权值和阈值,再利用得到的权值和阈值对每个缺失属性进行估值,进而将缺失数据恢复完整。再通过模糊C均值算法对恢复完整后的滚动轴承数据进行聚类,得到轴承数据的故障分类结果。
本发明的目的是通过下述技术方案实现的:一种改进BP神经网络估值的缺失数据轴承故障诊断方法,其步骤如下:
1)轴承数据预处理:将采集到的滚动轴承的原始数据进行特征提取,选取其中的9个特征,并将确定的滚动轴承数据进行人工随机缺失处理,得到缺失样本;
2)确定和优化训练样本:采用局部距离公式(1)计算每个缺失样本与其他所有样本的相似度,将得到的相似度从大到小排列,为每个缺失属性选取相似度最大的样本作为相应的预训练样本集,再针对每个样本缺失属性的位置对其预训练样本集中数据的相应属性做缺失处理,将处理后的数据集作为训练样本集,训练样本集作为网络的输入,同时每个输入值也作为网络的期望输出Y;局部距离公式如下:
其中,对于缺失数据集 与均是中的数据样本,xia和xib分别是和的第i个属性,s表示样本属性的数目,N表示数据集中样本的总数;
3)初始化网络:根据选取的缺失轴承数据训练样本集来确定网络输入层节点数n=9,隐含层节点数l=14,输出层节点数m=9;初始化输入层,隐含层和输出层神经元之间的连接权值wij,wjk,初始化隐含层阈值a和输出层阈值b,给定学习速率和神经元激励函数,确定最大训练次数M,误差精度ε1;
4)改进后基于缺失数据的BP网络的训练:运用每个缺失属性的训练样本集对改进的BP网络进行训练,得到针对每个缺失属性训练的神经网络,获得相应的权值wij,wjk和阈值a,b:
4a)隐含层输出计算:根据训练样本输入向量输入层和隐含层之间的连接权值wij以及隐含层阈值a,计算隐含层输出H;
且
式中,n表示输入层节点数,l表示隐含层节点数,是训练样本集中的一个数据样本,xi表示数据样本的第i个属性,为输入层节点数目恢复系数,f为隐含层激励函数,激励函数为:
4b)输出层计算输出:根据隐含层输出H,连接权值wjk和阈值b,计算输出层的预测输出O;
式中,l表示隐含层节点数,m表示输出层节点数。
4c)误差计算:根据网络预测输出O和期望输出Y,计算网络预测误差e;
且
式中,m表示输出层节点数,Yk表示期望输出数据样本的第k个属性,Ok表示预测输出数据的第k个属性,为输出层节点数目恢复系数;
4d)权值更新:根据网络预测误差e更新连接权值wij和wjk;
wjk=wjk+ηHjek,j=1,2,...,l;k=1,2,...m; (13)
式中,n表示输入层节点数,l表示隐含层节点数,m表示输出层节点数,x(i)表示数据样本的第i个属性,η=0.1为学习速率;
4e)阈值更新:根据网络预测误差e更新网络节点阈值a,b;
bk=bk+ek,k=1,2,...m (15)
4f)算法终止条件判定:当e<ε1或者训练次数大于最大训练次数M时,得到相应的权值和阈值,转到步骤4g);否则,则返回步骤4a);
4g)利用步骤4f)得到的针对相应缺失属性训练的网络的权值和阈值,对隐含层和输出层节点之间的连接权值wij和wjk,隐含层阈值a和输出层阈值b进行赋值;
5)对缺失属性进行估值:利用训练好的改进BP网络对每个缺失属性进行估值,进而将缺失数据恢复完整,最后得到恢复完整后的滚动轴承数据集:
5a)将得到的wij和隐含层阈值a带入公式(3)计算隐含层输出H;
5b)将得到的隐含层输出H、权值wjk和输出层阈值b带入公式(7),从得到的输出层输出值中获取相应缺失属性的估计值,将整个缺失数据集填补成完整数据集。
6)对数据集进行聚类分析:利用模糊C均值算法对恢复完整后的滚动轴承数据集进行聚类,最终得到轴承数据的故障分类结果。
所述的步骤6)的具体步骤为:
6a)初始化参数:设定聚类中心数目c=4,即缺失轴承数据集故障分类数目,迭代最大次数为G;确定模糊加权系数m和迭代终止阈值ε,m的值一般取2,ε的值一般取0.001至0.01之间的数;初始化隶属度矩阵U(0);
6b)计算聚类中心矩阵V:当迭代到第l(l=1,2,...)次时,根据U(l-1),利用式(16)计算聚类中心矩阵V(l);
式中,vi表示聚类中心矩阵V中的第i个中心,uik表示隶属度矩阵中第k个数据样本隶属于第i类的程度,xk表示轴承数据集中的第k个样本;
6c)计算隶属度矩阵:根据V(l),利用式(17)计算隶属度矩阵U(l);
式中,m为模糊加权系数,n为轴承数据集的样本数目;
6d)迭代终止阈值判定:对于给定的阈值ε,如果max|U(l+1)-U(l)|≤ε,或者迭代次数l>G,则迭代终止,否则l=l+1,转到步骤6b),最终,得到聚类中心矩阵V和隶属度矩阵U;
6f)通过判断每个数据样本隶属于每个故障类别的程度,来对缺失的轴承数据进行故障诊断。
本发明的有益效果:本发明采用上述方法,能够充分利用数据样本及属性之间的关联性和完整数据样本及缺失数据样本的分布信息,得到合理的属性估值,对不完整的数据进行聚类,从而得到缺失数据轴承故障的诊断结果,解决了现有技术中存在的由于数据的缺失而导致的不能应用模糊C均值聚类算法直接对不完整数据进行聚类的技术问题。
附图说明
图1是BP神经网络的拓补结构图。
图2是改进BP神经网络的拓补结构图。
具体实施方式
1、BP神经网络
BP神经网络通常是由输入层、隐含层和输出层这三层组成,层与层之间全互连,但每层节点之间不相连。如图1所示,为一个具有单隐含层的BP神经网络模型。该BP网络由输入层,隐含层和输出层构成。每个圆圈表示一个节点,每层各包含n,l,m个节点。节点之间的链接用箭头来表示,每个箭头表示一个权重。wij表示输入层和隐含层之间的连接权值,wjk表示隐含层和输出层得连接权值。数据的处理和计算将有隐含层和输出层的每个节点执行,隐含层节点的具体数目将在实验中确定。
2、改进BP神经网络
基本的BP神经网络的训练样本都必须是完整的数据样本,而本方法中对每个缺失属性的训练样本都是缺失样本集。因而,基本BP神经网络是不能直接使用的,需要对其进行改进。如图2所示,假设一个输入样本中的第三个属性是缺失的,用“?”来表示。则在计算隐含层输出时缺失属性不参与计算隐含层的值,而是有其余完整属性进行计算。训练过程中需要通过误差的反向传播更新网络的权值和阈值,因此在计算网络预测输出和期望输出的误差时,缺失属性不参与误差的计算,以免影响网络的权值和阈值的更新。经过多次迭代学习,可以得到网络的权值和阈值,即训练好的网络。
3、模糊C均值(FCM)聚类算法
模糊C均值聚类算法(Bezdek,1981)将s维的数据集中的数据样本分为c类,且c(2≤c≤n),聚类中心为V=[v1,v2,...,vc],第j类的聚类中心用vj∈Rs表示。它的基本思想是:建立基于隶属度和聚类中心的目标函数,通过对隶属度矩阵以及聚类中心的迭代优化,达到目标函数极小化的目的,从而实现对样本的聚类。聚类结果用隶属度矩阵U(c×n)表示,其中,n表示聚类数据集中的样本数目,c表示聚类中心数目,隶属度矩阵中的元素uij表示第j个数据样本隶属于第i类的程度,并且满足以下条件:
uik∈[0,1],i=1,2,...,c;k=1,2,...,n; (18)
目标函数定义如下:
其中,xk=[x1k,x2k,...,xsk]T是第k个数据样本,xjk是xk的第j个属性值;vi是第i个聚类中心;m(m≥1)为影响隶属度矩阵模糊化程度的指数权重;||·||2表示欧氏距离。
聚类中心和隶属度的迭代更新公式如下:
在公式(19)的约束下,交替迭代U和V,使公式(21)达到极小值。从而可以获得隶属度矩阵,并确定数据集聚类结果。
将本发明一种改进BP神经网络估值的缺失数据轴承故障诊断方法用于机械轴承的故障诊断,具体实施步骤如下:
1)采集原始信号:滚动轴承数据时来自美国Case Western ReserveUniversity电气工程实验室。在不同负载(0,1,2,3hp)和不同故障深度(7,14,21mil)下有四种状态,分别是正常、内圈故障、外圈故障、滚动体故障。数据的采样频率为12K和48K。每个状态有100组样本,一共是400组。
2)对原始信号进行特征提取:对信号进行预处理,提取特征,震动信号特征值得计算有很多种方法,选择峰峰值,均值,绝对平均值,均方值,均方根值,方差,标准偏差,偏度,峰值这9中特征值来处理原始信号。
①峰峰值是指信号的变化范围。公式为:
max(xi)-min(xi) (24)
②平均值是信号的平均值
③绝对平均值是信号幅值绝对值的算术平均值
④均方值不仅反应信号的平均值还反应信号的波动和离散程度
⑤均方根植反应信号振动强度和能量的大小
⑥方差描述信号偏离中心趋势的波动强度,公式为
⑦标准偏差公式;标准偏差是一种量度数据分布的分散程度之标准
μ为均值
⑧偏度是指振动信号的偏态是指数据分布的偏斜方向和程度
s为标准偏差μ为均值
⑨峰值是指指数据分布的振动信号尖峭程度或峰凸程度
选择这九个特征作为轴承数据样本的特征属性。
3)轴承数据预处理:将采集到的滚动轴承的原始数据进行特征提取,选取其中的9个特征,并将确定的滚动轴承数据进行人工随机缺失处理,得到缺失样本;
4)确定和优化训练样本:采用局部距离公式(1)计算每个缺失样本与其他所有样本的相似度,将得到的相似度从大到小排列,为每个缺失属性选取相似度最大的样本作为相应的预训练样本集,再针对每个样本缺失属性的位置对其预训练样本集中数据的相应属性做缺失处理,将处理后的数据集作为训练样本集,训练样本集作为网络的输入,同时每个输入值也作为网络的期望输出Y;局部距离公式如下:
其中,对于缺失数据集 与均是中的数据样本,xia和xib分别是和的第i个属性,s表示样本属性的数目,N表示数据集中样本的总数;
5)初始化网络:根据选取的缺失轴承数据训练样本集来确定网络输入层节点数n=9,隐含层节点数l=14,输出层节点数m=9;初始化输入层,隐含层和输出层神经元之间的连接权值wij,wjk,初始化隐含层阈值a和输出层阈值b,给定学习速率和神经元激励函数,确定最大训练次数M,误差精度ε1;
6)改进后基于缺失数据的BP网络的训练:运用每个缺失属性的训练样本集对改进的BP网络进行训练,得到针对每个缺失属性训练的神经网络,获得相应的权值wij,wjk和阈值a,b:
6a)隐含层输出计算:根据训练样本输入向量输入层和隐含层之间的连接权值wij以及隐含层阈值a,计算隐含层输出H;
且
式中,n表示输入层节点数,l表示隐含层节点数,是训练样本集中的一个数据样本,xi表示数据样本的第i个属性,为输入层节点数目恢复系数,f为隐含层激励函数,激励函数为:
6b)输出层计算输出:根据隐含层输出H,连接权值wjk和阈值b,计算输出层的预测输出O;
式中,l表示隐含层节点数,m表示输出层节点数。
6c)误差计算:根据网络预测输出O和期望输出Y,计算网络预测误差e。
且
式中,m表示输出层节点数,Yk表示期望输出数据样本的第k个属性,Ok表示预测输出数据的第k个属性,为输出层节点数目恢复系数。
6d)权值更新:根据网络预测误差e更新连接权值wij和wjk。
wjk=wjk+ηHjek,j=1,2,...,l;k=1,2,...m; (13)
式中,n表示输入层节点数,l表示隐含层节点数,m表示输出层节点数,x(i)表示数据样本的第i个属性,η=0.1为学习速率。
6e)阈值更新:根据网络预测误差e更新网络节点阈值a,b;
bk=bk+ek,k=1,2,...m (15)
6f)算法终止条件判定:当e<ε1或者训练次数大于最大训练次数M时,得到相应的权值和阈值,转到步骤6g);否则,则返回步骤6a);
6g)利用步骤6f)得到的针对相应缺失属性训练的网络的权值和阈值,对隐含层和输出层节点之间的连接权值wij和wjk,隐含层阈值a和输出层阈值b进行赋值;
7)对缺失属性进行估值:利用训练好的改进BP网络对每个缺失属性进行估值,进而将缺失数据恢复完整,最后得到恢复完整后的滚动轴承数据集:
7a)将得到的wij和隐含层阈值a带入公式(3)计算隐含层输出H;
7b)将得到的隐含层输出H、权值wjk和输出层阈值b带入公式(7),从得到的输出层输出值中获取相应缺失属性的估计值,将整个缺失数据集填补成完整数据集。
8)对数据集进行聚类分析:利用模糊C均值算法对恢复完整后的滚动轴承数据集进行聚类,最终得到轴承数据的故障分类结果,具体步骤为:
8a)初始化参数:设定聚类中心数目c=4,即缺失轴承数据集故障分类数目,迭代最大次数为G;确定模糊加权系数m和迭代终止阈值ε,m的值一般取2,ε的值一般取0.001至0.01之间的数;初始化隶属度矩阵U(0);
8b)计算聚类中心矩阵V:当迭代到第l(l=1,2,...)次时,根据U(l-1),利用式(16)计算聚类中心矩阵V(l);
式中,vi表示聚类中心矩阵V中的第i个中心,uik表示隶属度矩阵中第k个数据样本隶属于第i类的程度,xk表示轴承数据集中的第k个样本;
8c)计算隶属度矩阵:根据V(l),利用式(17)计算隶属度矩阵U(l);
式中,m为模糊加权系数,n为轴承数据集的样本数目;
8d)迭代终止阈值判定:对于给定的阈值ε,如果max|U(l+1)-U(l)|≤ε,或者迭代次数l>G,则迭代终止,否则l=l+1,转到步骤8b),最终,得到聚类中心矩阵V和隶属度矩阵U;
8f)通过判断每个数据样本隶属于每个故障类别的程度,来对缺失的轴承数据进行故障诊断。
实验结果分析:将特征提取后的轴承数据通过人工处理产生随机缺失数据的轴承数据集,缺失率以此为5%、10%、15%和20%,然后为每个缺失属性选择并生成训练样本。缺失轴承数据集通过改进后基于缺失数据的BP网络训练,并估值之后,得到恢复后的完整数据集,再将恢复后的数据集用模糊C均值进行聚类分析,最终得到隶属度矩阵U(c×n),其中,c=4代表四个聚类中心,表示轴承数据集的四种类别,分别为正常、内圈故障、外圈故障、滚动体故障,n表示轴承数据样本的数目,隶属度矩阵中的值uik表示每个轴承数据样本隶属于每个故障类别的程度。由隶属度矩阵可以得到每个轴承样本隶属于四个故障类别的隶属度值,通过比较四个隶属度值的大小判定轴承样本的故障类别,即属于隶属度值最大的类别。例如,一个数据样本隶属于正常状态的程度为0.1、隶属于内圈故障的程度为0.8、隶属于外圈故障的程度为0.02、隶属于滚动体故障的程度为0.08,则可以判定该样本属于内圈故障的隶属度值最大,所以为内圈故障。
为了说明本发明提出的改进BP神经网络估值的缺失数据轴承故障诊断(IBPFCM)方法的有效性,分别同完备数据策略(WDS)、局部距离策略(PDS)、优化完整策略(OCS)和最近原型策略(NPS)针对滚动轴承健康度评价的聚类结果进行对比。分别对模糊聚类评价指标RI、FR、JR和MR进行对比,其中,RI、FR和JR的值越大说明聚类结果越好,MR的值越小说明聚类结果越好。表1为五种方法10次试验所得RI指标的平均值,表2为五种方法10次试验所得FR指标的平均值,表3为五种方法10次试验所得JR指标的平均值,表4为五种方法10次试验所得MR指标的平均值,由表1-4可以看出,本发明所提出的算法IBPFCM与其他四种算法相比能够更好的对轴承的健康度进行评价,尤其当缺失率为5%、10%和15%的时候都能够得到最优的结果。
表1
表2
表3
表4
Claims (2)
1.一种改进BP神经网络估值的缺失数据轴承故障诊断方法,其特征在于,步骤如下:
1)轴承数据预处理:将采集到的滚动轴承的原始数据进行特征提取,选取其中的9个特征,并将确定的滚动轴承数据进行人工随机缺失处理,得到缺失样本;
2)确定和优化训练样本:采用局部距离公式(1)计算每个缺失样本与其他所有样本的相似度,将得到的相似度从大到小排列,为每个缺失属性选取相似度最大的样本作为相应的预训练样本集,再针对每个样本缺失属性的位置对其预训练样本集中数据的相应属性做缺失处理,将处理后的数据集作为训练样本集,训练样本集作为网络的输入,同时每个输入值也作为网络的期望输出Y;局部距离公式如下:
其中,对于缺失数据集 与均是中的数据样本,xia和xib分别是和的第i个属性,s表示样本属性的数目,N表示数据集中样本的总数;
3)初始化网络:根据选取的缺失轴承数据训练样本集来确定网络输入层节点数n=9,隐含层节点数l=14,输出层节点数m=9;初始化输入层,隐含层和输出层神经元之间的连接权值wij,wjk,初始化隐含层阈值a和输出层阈值b,给定学习速率和神经元激励函数,确定最大训练次数M,误差精度ε1;
4)改进后基于缺失数据的BP网络的训练:运用每个缺失属性的训练样本集对改进的BP网络进行训练,得到针对每个缺失属性训练的神经网络,获得相应的权值wij,wjk和阈值a,b:
4a)隐含层输出计算:根据训练样本输入向量输入层和隐含层之间的连接权值wij以及隐含层阈值a,计算隐含层输出H;
且
式中,n表示输入层节点数,l表示隐含层节点数,是训练样本集中的一个数据样本,xi表示数据样本的第i个属性,为输入层节点数目恢复系数,f为隐含层激励函数,激励函数为:
4b)输出层计算输出:根据隐含层输出H,连接权值wjk和阈值b,计算输出层的预测输出O;
式中,l表示隐含层节点数,m表示输出层节点数;
4c)误差计算:根据网络预测输出O和期望输出Y,计算网络预测误差e;
且
式中,m表示输出层节点数,Yk表示期望输出数据样本的第k个属性,Ok表示预测输出数据的第k个属性,为输出层节点数目恢复系数。
4d)权值更新:根据网络预测误差e更新连接权值wij和wjk;
wjk=wjk+ηHjek,j=1,2,...,l;k=1,2,...m; (13)
式中,n表示输入层节点数,l表示隐含层节点数,m表示输出层节点数,x(i)表示数据样本的第i个属性,η=0.1为学习速率。
4e)阈值更新:根据网络预测误差e更新网络节点阈值a,b;
bk=bk+ek,k=1,2,...m (15)
4f)算法终止条件判定:当e<ε1或者训练次数大于最大训练次数M时,得到相应的权值和阈值,转到步骤4g);否则,则返回步骤4a);
4g)利用步骤4f)得到的针对相应缺失属性训练的网络的权值和阈值,对隐含层和输出层节点之间的连接权值wij和wjk,隐含层阈值a和输出层阈值b进行赋值;
5)对缺失属性进行估值:利用训练好的改进BP网络对每个缺失属性进行估值,进而将缺失数据恢复完整,最后得到恢复完整后的滚动轴承数据集:
5a)将得到的wij和隐含层阈值a带入公式(3)计算隐含层输出H;
5b)将得到的隐含层输出H、权值wjk和输出层阈值b带入公式(7),从得到的输出层输出值中获取相应缺失属性的估计值,将整个缺失数据集填补成完整数据集;
6)对数据集进行聚类分析:利用模糊C均值算法对恢复完整后的滚动轴承数据集进行聚类,最终得到轴承数据的故障分类结果。
2.根据权利要求1所述的一种改进BP神经网络估值的缺失数据轴承故障诊断方法,其特征在于:所述的步骤6)的具体步骤为:
6a)初始化参数:设定聚类中心数目c=4,即缺失轴承数据集故障分类数目,迭代最大次数为G;确定模糊加权系数m和迭代终止阈值ε,m的值一般取2,ε的值一般取0.001至0.01之间的数;初始化隶属度矩阵U(0);
6b)计算聚类中心矩阵V:当迭代到第l(l=1,2,...)次时,根据U(l-1),利用式(16)计算聚类中心矩阵V(l);
式中,vi表示聚类中心矩阵V中的第i个中心,uik表示隶属度矩阵中第k个数据样本隶属于第i类的程度,xk表示轴承数据集中的第k个样本;
6c)计算隶属度矩阵:根据V(l),利用式(17)计算隶属度矩阵U(l);
k=1,2,...,n (17)
式中,m为模糊加权系数,n为轴承数据集的样本数目;
6d)迭代终止阈值判定:对于给定的阈值ε,如果max|U(l+1)-U(l)|≤ε,或者迭代次数l>G,则迭代终止,否则l=l+1,转到步骤6b),最终,得到聚类中心矩阵V和隶属度矩阵U;
6f)通过判断每个数据样本隶属于每个故障类别的程度,来对缺失的轴承数据进行故障诊断。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510172600.0A CN104751229B (zh) | 2015-04-13 | 2015-04-13 | 一种改进bp神经网络估值的缺失数据轴承故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510172600.0A CN104751229B (zh) | 2015-04-13 | 2015-04-13 | 一种改进bp神经网络估值的缺失数据轴承故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104751229A true CN104751229A (zh) | 2015-07-01 |
CN104751229B CN104751229B (zh) | 2017-05-24 |
Family
ID=53590874
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510172600.0A Active CN104751229B (zh) | 2015-04-13 | 2015-04-13 | 一种改进bp神经网络估值的缺失数据轴承故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104751229B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095962A (zh) * | 2015-07-27 | 2015-11-25 | 中国汽车工程研究院股份有限公司 | 一种基于bp人工神经网络的材料动态力学性能预测方法 |
CN106323636A (zh) * | 2016-08-16 | 2017-01-11 | 重庆交通大学 | 栈式稀疏自动编码深度神经网络的机械故障程度特征自适应提取与诊断方法 |
CN106610584A (zh) * | 2015-10-27 | 2017-05-03 | 沈阳工业大学 | 一种基于神经网络与专家系统的再制造工序质量控制方法 |
CN107273429A (zh) * | 2017-05-19 | 2017-10-20 | 哈工大大数据产业有限公司 | 一种基于深度学习的缺失值填充方法及系统 |
CN107480787A (zh) * | 2017-08-10 | 2017-12-15 | 西安建筑科技大学 | 一种基于bp神经网络的露天矿旋回式破碎机故障诊断方法 |
CN107655690A (zh) * | 2017-10-30 | 2018-02-02 | 王才旺 | 基于人工神经网络的电机轴承故障诊断方法 |
CN107729943A (zh) * | 2017-10-23 | 2018-02-23 | 辽宁大学 | 信息反馈极限学习机优化估值的缺失数据模糊聚类算法及其应用 |
CN108805193A (zh) * | 2018-06-01 | 2018-11-13 | 广东电网有限责任公司 | 一种基于混合策略的电力缺失数据填充方法 |
CN108985515A (zh) * | 2018-07-24 | 2018-12-11 | 国网河南省电力公司电力科学研究院 | 一种基于独立循环神经网络的新能源出力预测方法及系统 |
CN109240270A (zh) * | 2018-10-09 | 2019-01-18 | 宁波大学 | 一种基于假设缺失数据迭代估计误差的动态过程监测方法 |
CN109292567A (zh) * | 2018-02-28 | 2019-02-01 | 武汉大学 | 一种基于bp神经网络的电梯故障预测方法 |
CN109973331A (zh) * | 2019-05-05 | 2019-07-05 | 内蒙古工业大学 | 一种基于bp神经网络的风力发电机组风轮叶片故障诊断算法 |
CN110705181A (zh) * | 2019-10-13 | 2020-01-17 | 重庆交通大学 | 基于卷积长短时记忆循环神经网络的滚动轴承剩余寿命预测方法 |
CN111140986A (zh) * | 2019-12-23 | 2020-05-12 | 珠海格力电器股份有限公司 | 空调系统的运行状态检测方法、装置、存储介质及空调 |
CN111259927A (zh) * | 2020-01-08 | 2020-06-09 | 西北工业大学 | 基于神经网络和证据理论的火箭发动机故障诊断方法 |
CN112434929A (zh) * | 2020-11-19 | 2021-03-02 | 广西电网有限责任公司电力科学研究院 | 一种配电变压器抗风险能力评估方法 |
CN113469107A (zh) * | 2021-07-16 | 2021-10-01 | 辽宁大学 | 一种融合空间密度分布的轴承故障诊断方法 |
CN113469106A (zh) * | 2021-07-16 | 2021-10-01 | 辽宁大学 | 一种缺失数据的knn区间型轴承故障诊断方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102063109B (zh) * | 2010-11-29 | 2012-09-05 | 株洲南车时代电气股份有限公司 | 一种基于神经网络的地铁列车故障诊断装置及其方法 |
CN104299035A (zh) * | 2014-09-29 | 2015-01-21 | 国家电网公司 | 基于聚类算法和神经网络的变压器故障诊断方法 |
-
2015
- 2015-04-13 CN CN201510172600.0A patent/CN104751229B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102063109B (zh) * | 2010-11-29 | 2012-09-05 | 株洲南车时代电气股份有限公司 | 一种基于神经网络的地铁列车故障诊断装置及其方法 |
CN104299035A (zh) * | 2014-09-29 | 2015-01-21 | 国家电网公司 | 基于聚类算法和神经网络的变压器故障诊断方法 |
Non-Patent Citations (2)
Title |
---|
张军 等: "滚动轴承振动诊断的BP神经网络方法", 《轻工机械》 * |
高佃波: "基于BP 神经网络和振动测量的轴承故障诊断", 《长沙交通学院学报》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105095962B (zh) * | 2015-07-27 | 2017-07-28 | 中国汽车工程研究院股份有限公司 | 一种基于bp人工神经网络的材料动态力学性能预测方法 |
CN105095962A (zh) * | 2015-07-27 | 2015-11-25 | 中国汽车工程研究院股份有限公司 | 一种基于bp人工神经网络的材料动态力学性能预测方法 |
CN106610584A (zh) * | 2015-10-27 | 2017-05-03 | 沈阳工业大学 | 一种基于神经网络与专家系统的再制造工序质量控制方法 |
CN106323636A (zh) * | 2016-08-16 | 2017-01-11 | 重庆交通大学 | 栈式稀疏自动编码深度神经网络的机械故障程度特征自适应提取与诊断方法 |
CN107273429A (zh) * | 2017-05-19 | 2017-10-20 | 哈工大大数据产业有限公司 | 一种基于深度学习的缺失值填充方法及系统 |
CN107480787A (zh) * | 2017-08-10 | 2017-12-15 | 西安建筑科技大学 | 一种基于bp神经网络的露天矿旋回式破碎机故障诊断方法 |
CN107480787B (zh) * | 2017-08-10 | 2020-12-11 | 西安建筑科技大学 | 一种基于bp神经网络的露天矿旋回式破碎机故障诊断方法 |
CN107729943B (zh) * | 2017-10-23 | 2021-11-30 | 辽宁大学 | 信息反馈极限学习机优化估值的缺失数据模糊聚类算法及其应用 |
CN107729943A (zh) * | 2017-10-23 | 2018-02-23 | 辽宁大学 | 信息反馈极限学习机优化估值的缺失数据模糊聚类算法及其应用 |
CN107655690A (zh) * | 2017-10-30 | 2018-02-02 | 王才旺 | 基于人工神经网络的电机轴承故障诊断方法 |
CN109292567A (zh) * | 2018-02-28 | 2019-02-01 | 武汉大学 | 一种基于bp神经网络的电梯故障预测方法 |
CN108805193B (zh) * | 2018-06-01 | 2022-04-26 | 广东电网有限责任公司 | 一种基于混合策略的电力缺失数据填充方法 |
CN108805193A (zh) * | 2018-06-01 | 2018-11-13 | 广东电网有限责任公司 | 一种基于混合策略的电力缺失数据填充方法 |
CN108985515B (zh) * | 2018-07-24 | 2021-11-26 | 国网河南省电力公司电力科学研究院 | 一种基于独立循环神经网络的新能源出力预测方法及系统 |
CN108985515A (zh) * | 2018-07-24 | 2018-12-11 | 国网河南省电力公司电力科学研究院 | 一种基于独立循环神经网络的新能源出力预测方法及系统 |
CN109240270B (zh) * | 2018-10-09 | 2021-03-09 | 宁波大学 | 一种基于假设缺失数据迭代估计误差的动态过程监测方法 |
CN109240270A (zh) * | 2018-10-09 | 2019-01-18 | 宁波大学 | 一种基于假设缺失数据迭代估计误差的动态过程监测方法 |
CN109973331A (zh) * | 2019-05-05 | 2019-07-05 | 内蒙古工业大学 | 一种基于bp神经网络的风力发电机组风轮叶片故障诊断算法 |
CN110705181A (zh) * | 2019-10-13 | 2020-01-17 | 重庆交通大学 | 基于卷积长短时记忆循环神经网络的滚动轴承剩余寿命预测方法 |
CN110705181B (zh) * | 2019-10-13 | 2022-06-24 | 重庆交通大学 | 基于卷积长短时记忆循环神经网络的滚动轴承剩余寿命预测方法 |
CN111140986A (zh) * | 2019-12-23 | 2020-05-12 | 珠海格力电器股份有限公司 | 空调系统的运行状态检测方法、装置、存储介质及空调 |
CN111259927A (zh) * | 2020-01-08 | 2020-06-09 | 西北工业大学 | 基于神经网络和证据理论的火箭发动机故障诊断方法 |
CN111259927B (zh) * | 2020-01-08 | 2022-08-05 | 西北工业大学 | 基于神经网络和证据理论的火箭发动机故障诊断方法 |
CN112434929A (zh) * | 2020-11-19 | 2021-03-02 | 广西电网有限责任公司电力科学研究院 | 一种配电变压器抗风险能力评估方法 |
CN112434929B (zh) * | 2020-11-19 | 2023-04-18 | 广西电网有限责任公司电力科学研究院 | 一种配电变压器抗风险能力评估方法 |
CN113469106A (zh) * | 2021-07-16 | 2021-10-01 | 辽宁大学 | 一种缺失数据的knn区间型轴承故障诊断方法 |
CN113469107A (zh) * | 2021-07-16 | 2021-10-01 | 辽宁大学 | 一种融合空间密度分布的轴承故障诊断方法 |
CN113469106B (zh) * | 2021-07-16 | 2024-06-04 | 辽宁大学 | 一种缺失数据的knn区间型轴承故障诊断方法 |
CN113469107B (zh) * | 2021-07-16 | 2024-06-04 | 辽宁大学 | 一种融合空间密度分布的轴承故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104751229B (zh) | 2017-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104751229A (zh) | 一种改进bp神经网络估值的缺失数据轴承故障诊断方法 | |
Shao et al. | A novel approach of multisensory fusion to collaborative fault diagnosis in maintenance | |
Guo et al. | Hierarchical adaptive deep convolution neural network and its application to bearing fault diagnosis | |
CN106405640B (zh) | 基于深度信念神经网络的微震信号到时自动拾取方法 | |
Jiang et al. | Two-stage structural damage detection using fuzzy neural networks and data fusion techniques | |
CN111967378B (zh) | 拉剪倾倒式岩溶危岩失稳预警的声发射多前兆方法及装置 | |
CN111967486A (zh) | 一种基于多传感器融合的复杂装备故障诊断方法 | |
CN106599920A (zh) | 一种基于耦合隐半马尔可夫模型的航空轴承故障诊断方法 | |
CN103926526A (zh) | 一种基于改进的rbf神经网络的模拟电路故障诊断方法 | |
CN107016241A (zh) | 基于加噪样本增强深度自编码学习网络的旋转机械寿命阶段识别方法 | |
Di et al. | Ensemble deep transfer learning driven by multisensor signals for the fault diagnosis of bevel-gear cross-operation conditions | |
Liu et al. | Fault diagnosis of rod pumping wells based on support vector machine optimized by improved chicken swarm optimization | |
CN110455512B (zh) | 基于深度自编码器dae的旋转机械多集成故障诊断方法 | |
CN105975744A (zh) | 一种基于d-s证据理论的纺织过程数据融合系统 | |
Petković et al. | Soft methodology selection of wind turbine parameters to large affect wind energy conversion | |
CN114118138A (zh) | 一种基于多标签领域自适应模型的轴承复合故障诊断方法 | |
CN112329520B (zh) | 基于生成对抗学习的货车轴承故障识别方法 | |
CN104537243A (zh) | 一种基于模糊推理结合加权相似度量的心电分类方法 | |
CN103810392B (zh) | 基于支持向量机和rbf神经网络的退化数据缺失插补方法 | |
CN109145951B (zh) | 一种基于贝叶斯网络的隧道衬砌结构服役状况评价方法 | |
Ghiasi et al. | Comparison of seven artificial intelligence methods for damage detection of structures | |
CN116383757A (zh) | 一种基于多尺度特征融合和迁移学习的轴承故障诊断方法 | |
CN109214089A (zh) | 一种桥梁结构健康监测指标的隶属度评估方法 | |
CN114252266A (zh) | 一种基于dbn-svdd模型的滚动轴承性能退化评估方法 | |
Yang et al. | An improved probabilistic neural network with ga optimization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240110 Address after: 905, Building G, Huangjin Times Square, No. 9999 Jingshi Road, Lixia District, Jinan City, Shandong Province, 250000 Patentee after: Zhongchangxing (Shandong) Information Technology Co.,Ltd. Address before: 110000 58 Shenbei New Area Road South, Shenyang, Liaoning. Patentee before: LIAONING University |
|
TR01 | Transfer of patent right |