基于BP神经网络的故障诊断方法模板

1、：ZHEJIANGUNIVERSITYOFSCIENCE&TECHNOLOGY智能控制基础研究生课程设计报告题目基于BP神经网络的故障诊断方法学院机械与汽车工程学院专业班级车辆工程学号学1601852020学生姓名李跃轩指导教师武晓莉完成日期2016年12月10日目录1设计概述2.1.1 研究对象介绍2.1.2 设计内容及目标2.2设计原理、方法及步骤3.2.1 基于BP算法的神经网络模型32.2 神经网络信息融合故障诊断步骤43结果及分析6.3.1 数据仿真.6.3.2 结果分析8.4设计小结9.参考文献1Q附录程序111设计概述1.1 研究对象介绍信息融合是多源信息综合处理的一项新

2、技术，是将来自某一目标（或状态）的多源信息加以智能化合成，产生比单一信息源更精确、更完全的估计和判决。信息融合所处理的多传感器信息具有更为复杂的形式，可以在不同的信息层次上出现。多传感器信息融合的优点突出地表现在信息的冗余性、容错性、互补性、实时性和低成本性。神经网络是由大量互联的处理单元连接而成，它是基于现代神经生物学以及认知科学在信息处理领域应用的研究成果。它具有大规模并行模拟处理、连续时间动力学和网络全局作用等特点，有很强的自适应学习和非线性拟合能力，从而可以替代复杂耗时的传统算法，使信号处理过程更接近人类思维活动。柴油机故障具有相似性，故障与征兆的关系不明确，具有较强的模糊性，故障特征

3、相互交织，柴油机故障诊断是一个复杂的问题。综合柴油机故障的特点以及神经网络的优势，采用基于BP神经网络的多传感器信息融合技术对柴油机机械故障进行诊断。1.2 设计内容及目标设计内容：针对传统故障诊断方法存在的诊断准确性不高的问题，提出了BP神经网络信息融合的方法，实现对柴油机的机械故障诊断。由多个传感器采集信号，分别经过快速傅里叶变换后获得故障频域特征值，再经BP神经网络对柴油机进行故障局部诊断，能够对相应传感器的不同故障类型做出一个准确地分类，最终完成对汽轮机机械故障的准确诊断。实验结果表明，该方法克服了单个传感器的局限性和不确定性，是一种有效的故障诊断方法。采用方法：通过BP神经网络进行局

4、部诊断，最终判定故障及故障类型。基于BP神经网络多传感器信息融合，故障诊断方法是特征层状态属性融合，并利用MATLAB仿真。2设计原理、方法及步骤基于神经网络多传感器信息融合故障诊断方法是特征层状态属性融合，也就是特征层联合识别方法，多传感器检测系统为识别提供了比单传感器更多的有关目标(状态)的特征信息，增大了特征空间维数。本文运用神经网络多传感器信息融合方法对机械设备运行状态进行诊断识别，是基于这样一种思想：设备运行状态与其各种征兆参数(温度、压力、电压、电流、振动信号等)之间存在着因果关系，而这种关系之复杂是难用公式表达的，由于神经网络所具有的信息分布式存储方式、大规模自适应并行处理、高度

5、的容错能力等是其可用于模式识别的基础，特别是其学习能力、容错能力和高度的非线性映射能力对机械设备运行状态的不确定性模式识别具有独到之处。2.1基于BP算法的神经网络模型本文采用的是3层BP神经网络模型，由输入层、隐层和输出层构成，图1所示为一个典型的三层BP神经网络模型图1一个简单的BP神经网络模型网络的前馈意义在于每一层节点的输入仅来自前面一层节点的输出。对于输入信号，先前向传播到隐层节点，经过激活函数后，再把隐层节点的输出信息传播到输出节点，最后得到输出结果。(1)输入层节点i(i=1,2,.,n),其输出Oi等于输入Xi,将变量值传送到第二层。(2)隐层节点j(j=1,2,.,p),其输

6、入Ij,输出Oj分别为：nIj=£叫Oj+Hj(2-1)i1Oj=f(Ij)=1/1+exp(-Ij)(2-2)式中，切ji为隐层节点j与输入层节点i之间的权值，句为隐层节点j的偏置，f为sigmoid函数，其表达式为：f(x)=1/1+exp(-x)(2-3)(3)输出层节点k(k=1,2,.,m),其输入1,输出yk分别为:pIk=£防0wk(2-4)j4yk=f(Ik)=1/1+exp(Ik)(2-5)式中，6勾为输入层节点k与隐层节点j之间的连接权值，讥为输出层节点k的偏置。对于给定的训练样本(xp1,xp2,.,xpn),p为样本数(p=1,2,.,P),网络输出

7、与训练目标之间的均方误差可表示为：1 p/、E=£Ep(2-6)ppm_12/cr、Ep=-Z(tpijpi)(2-7)2 11其中p为样本数，tpi为第p个样本的第l个输出单元的目标输出结果，yp为第p个样本的第l个输出单元的网络运算结果。BP网络训练的过程包括网络内部的前向计算和误差的反向传播，其目的就是通过调整网络内部连接权值使网络输出误差最小。对于多层前馈网络中输入层与隐层之间、隐层与输出层之间连接权值利用BP算法调整。2.2神经网络信息融合故障诊断步骤神经网络知识表示是一种知识的隐式表示，知识表现为网络的拓扑结构和连接权值，采用神经网络技术的专家系统，由于神经网络是一种信息

8、存储和处理统一的网络系统，因此，在采用神经网络技术的专家系统中，知识的存储与问题求解过程中的推理过程均在系统的神经网络模块中进行，是推理机和知识库的统一。首先从已有的设备特征信号提取特征数据，经过数据预处理（归一化处理）后作为神经网络输入，从已知的故障结果提取数据作为神经网络输出，构建BP神经网络，利用已有的特征数据和已知的故障结果数据形成的训练样本集对构建的BP神经网络进行训练和网络自学习，使BP神经网络的权值、阀值与已知的故障结果之间存在对应关系达到期望的故障结果输出。当BP神经网络训练完毕后，就可以利用训练成功的BP神经网络进行故障诊断。故障诊断的过程如下：1）将故障样本输入给输入层各节

9、点，同时它也是该层神经元的输出。2）由式（2-2）求出隐层神经元的输出，并将其作为输出层的输入。3）从式（2-5）求得输出层神经元的输出。4）由阈值函数判定输出层神经元的最终输出结果。柴油机的故障诊断首先从待诊断的故障信号中提取数据并进行数据预处理，而后将待诊断故障数据输入训练成功的神经网络。利用神经网络信息融合进行故障诊断步骤如图2所示:压力传感器待测原件信号采集特征提取特征关联神经网络分类器诊断结果，温度传感器图2信息融合进行故障诊断步骤3结果及分析3.1 数据仿真以某种柴油机中的4个实际故障样本为例，每个故障样本有5个故障特征值，因此选取网络的输入节点为5。将样本输入到神经网络模型中，每

10、个输出节点代表一种故障类型，则这4种故障类型和一个正常状态共需5个输出节点与之相对应，因此网络的输出节点为5。表3-1训练数据。表3-1训练数据样本序号故障特征值X温度压力电压电流振动rp10.00000.00000.00130.00000.000020.00740.03350.00150.00320.010630.00430.02230.00000.00470.005340.00560.03150.00120.00270.008650.00630.02050.00200.00330.006660.00490.03000.00180.00260.007870.55200.31070.25810

11、.30940.231680.54520.27930.26110.29880.203690.55020.24580.27170.31150.2347100.54060.28620.26450.30090.2145110.56040.30280.28350.31250.2225120.55250.31250.27780.30050.2164130.24620.15080.09470.09640.0999140.25350.10610.09680.09710.0810150.26650.08940.09370.09440.0908160.25560.14650.09880.10000.1000170

12、.26280.14320.308990.09650.0899180.25620.09990.09000.09990.0998190.66150.52510.51950.47101.000200.67380.44130.52250.47320.9667210.66650.47490.52550.47690.9758220.66890.46520.52000.47560.9899230.67000.47480.52350.48251.0000240.66230.48320.52680.46690.9887251.00001.00000.98121.00000.8206260.97970.97771

13、.00000.99600.8206270.98460.97270.98470.98570.7600280.98000.98250.98350.98870.8000290.99650.99340.98620.98990.7800301.00000.99120.99380.99610.8014在表3-1训练数据中样本序号16是正常状态；样本序号712是故障1状态；样本序号1318是故障2状态；样本序号1924是故障3状态；样本序号2530是故障4状态。表3-2是测试数据。表3-2测试数据样本序号故障特征值X温度压力电压电流振动振p10.54930.26260.26590.30880.222120.

14、00310.02350.00050.00300.004530.99200.98990.99790.99370.797940.67040.49720.52350.47410.979150.25720.10060.09580.09810.089060.98000.98250.98350.98870.8000根据故障诊断的特点建立BP神经网络模型，网络的输入层节点、隐层节点和输出层节点个数分别为5、60和5,设置系统误差为1e-3,学习率为1,最大迭代次数为1000次，其中输入层到隐层之间的激活函数为双曲正切函数，隐层到输出层之间的激活函数设置为logsig的S型传递函数，在输出端设置的阈值为0.8

15、5,即网络输出值大于0.85的置1,小于等于0.85的置0来处理，利用MATLAB神经网络工具箱建立BP神经网络模型3.2 结果分析在本文中将正常状态设置为10000,故障1设置为01000,故障2设置为00100,故障3设置为00010,故障4设置为00001,输出阈值设置为0.85,即在输出值大于0.85即为1,否则为0o根据输出数据的分析判断属于哪种故障，表3-3为仿真输出。表3-3仿真输出样本序号网络输出y正常故障1故障2故障3故障410.00061.00080.0021-0.00960.001121.01390.0007-0.01150.0021-0.00243-0.0014-0.0

16、0190.0010-0.00190.992340.0410-0.0090-0.01590.9811-0.01445-0.0014-0.00201.0039-0.00270.00156-0.00020.00040.00250.00090.9970将表3-2中“样本序号1”的5个故障特征值输给网络的输入层节点，则网络输出层节点与其对应的输出为表3-3中“样本序号1”所在行的5个输出值,其中只有y2=1.0008>0.85,其他5个输出均远小于0.85,所以，网络故障诊断的结果为故障1;表3-2中“样本序号2”的5个故障特征值输给网络的输入层节点，则网络输出层节点与其对应的输出为表3-3中“样

17、本序号2”所在行的5个输出值，其中只有y=1.0139>0.85,其他5个输出均远小于0.85,所以，网络故障诊断的结果为正常状态；表3-2中“样本序号3”的5个故障特征值输给网络的输入层节点，则网络输出层节点与其对应的输出为表3-3中“样本序号3”所在行的5个输出值，其中只有y4=0.9923>0.85,其他5个输出均远小于0.85,所以，网络故障诊断的结果为故障4,其他类型的故障诊断以此类推。4设计小结因为学习的专业是车辆工程，因此选择了对柴油机故障诊断这一研究对象。此篇设计经过查阅文献资料，对BP神经网络的故障诊断及Matlab仿真有了一定的认识。柴油机故障诊断是一个非常复杂

18、的过程，其故障与征兆关系不是很明确，存在一种非线性映射关系。BP神经网络具有良好的学习能力，利用BP网络信息融合方法，能够在足够多的样本数目前提下网络保证良好的容错性和鲁棒性是比较好的。所以在故障诊断的过程中，神经网络信息融合发挥其联想记忆和分布并行处理能力，不仅能够诊断出已有的故障还能对故障进行一个预测，从而满足柴油机故障断的要求。人工神经网络信息融合为智能诊断提供了良好的方法，为自适应学习和决策高度智能化控制系统提供了强有力的基础，并具有广泛的应用潜力和发展前景。参考文献1朱大奇，于盛林.基于D-S证据理论的数据融合算法及其在电路故障诊断中的应用J.电子学报,2002,30(2):221-

19、223.2王万良.人工智能及其应用(第3版)(附光盘)(BZ)M.高等教育出版社，2016.3朱大奇，于盛林.电子电路故障诊断的神经网络数据融合算法J.东南大学学报(自然科学版),2001,31(2):87-90.4何友，王国宏.多传感器信息融合及应用M.北京：电子工业出版社,2000.5张绪锦，谭剑波，韩江洪.基于BP神经网络的故障诊断方法J.系统工程理论与实践,2002,22(6):61-66.6 Y.G.Lei,Z.J.He,Y.Y.Zi,Q.Hu,FaultdiagnosisofrotatingmachinerybasedonmultipleANFIScombinationwithGa

20、sJ.Mech.Syst.SignalProcess.2007(21)2280W294.7 V.T.Tran,F.AlThobiani,A.Ball,AnapproachtofaultdiagnosisofreciprocatingcompressorvalvesusingTeager-KaiserenergyoperatoranddeepbeliefnetworksJ.ExpertSyst.Appl.2014(41)4113-4122.附录程序clcclear%喻入训练数据（包括正常数据和故障数据）%”练样本中，每一列分别代表：“温度”，“压力”，“电压”，“电流”和“振动信号”；每一行分别

21、代表一个样本p11=0.00000.00000.00130.00000.0000'；p12=0.00740.03350.00150.00320.0106'p13=0.00430.02230.00000.00470.0053'p14=0.00560.03150.00120.00270.0086'p15=0.00630.02050.00200.00330.0066'p16=0.00490.03000.00180.00260.0078'%正常数据p21=0.55200.31070.25810.30940.2316'p22=0.54520.279

22、30.26110.29880.2036'p23=0.55020.24580.27170.31150.2347'p24=0.54060.28620.26450.30090.2145'p25=0.56040.30280.28350.31250.2225'p26=0.55250.31250.27780.30050.2164'%故障1p31=0.24620.15080.09470.09640.0999'p32=0.25350.10610.09680.09710.0810'p33=0.26650.08940.09370.09940.0908'

23、;p34=0.25560.14650.09880.10000.1000'p35=0.26280.14320.08990.09650.0899'p36=0.25620.09990.09000.09990.0998'%故障2p41=0.66150.52510.51950.47101.0000'p42=0.67380.44130.52250.47320.9667'p43=0.66650.47490.52550.47690.9758'p44=0.66890.46520.52000.47560.9899'p45=0.67000.47480.5235

24、0.48251.0000'p46=0.66230.48320.52680.46690.9887'%故障3p51=1.0001.00000.98121.00000.8206'p52=0.97970.97771.00000.99600.7759'p53=0.98460.97270.98470.98570.7600'p54=0.98000.98250.98350.98870.8000'p55=0.99650.99340.98620.98990.7800'p56=1.00000.99120.99380.99610.8014'%故障4p=p

25、11p12p13p14p15p16p21p22p23p24p25p26p31p32p33p34p35p36p41p42p43p44p45p46p51p52p53p54p55p56;将该数据集设置%将所有的数据：包括正常数据和异常数据放在一个数据集中，为输入神经网络的特征值%喻被检测对象类别进行编码（设置目标值，即标签，用二进制类型表示）%normalt11=10000't12=10000't13=10000't14=10000't15=10000't16=10000't21=01000't22=01000't23=01000

26、9;t24=01000't25=01000't26=01000'%fault1t31=00100't32=00100't33=00100't34=00100't35=00100't36=00100'%fault2t41=00010't42=00010't43=00010't44=00010't45=00010't46=00010'%fault3t51=00001't52=00001't53=00001't54=00001't55=00001&#

27、39;t56=00001'%fault4t=t11t12t13t14t15t16t21t22t23t24t25t26t31t32t33t34t35t36t41t42t43t44t45t46t51t52t53t54t55t56;p,s1=mapminmax(p);%俞入样本归一化归一化的范围是(-1,1),s1记录归一化的参数net=newff(minmax(p),60,5,'tansig','purelin','trainlm');%新建BP网络，p为样本输入：p=martric(5*30)分别表示5个传感器，30个样本%隐层神经元个数为60,根据输出目标t确定输出层神经元个数为5,%从输入层到隐层的激励函数为双曲正切，隐层到输出层的激励函数为线性函数，训练方法利用LM(Levenberg-Marquardt)算法进行网络参数sita=W,b的更新%吸置网络训练参数net.trainParam.show=10;%设置数据显示刷新频率，学习次刷新一次图象net.trainParam.epochs=1000;%最大训练次数net.trainParam.goal=1e-5;%设置训练误差net=init(net);%网络初始化nettr=train(net,p,t);%训练网络%testingp1=0.54930.