故障诊断学与可靠性工程作业(吐血分享).docx
《故障诊断学与可靠性工程作业(吐血分享).docx》由会员分享,可在线阅读,更多相关《故障诊断学与可靠性工程作业(吐血分享).docx(26页珍藏版)》请在三一文库上搜索。
1、机械故障诊断学及可靠性工程MECHANICAL FAULT DIAGNOSTICS ANDRELIABILITY ENGINEERING作业题目基于冗余第二代小波包变换和邻域粗 糙集以及支持向量机的机械故障诊断作者姓名黎原作者学号S12080203030学科专业机械设计及理论 指导教师董敏教授2013年6月机械故障诊断学及可靠性工程基于冗余第二代小波包变换和邻域粗糙集以及支持向量机的机械故障诊断硕士研究生:黎原硕士生学号:S12080202030导师:董敏教授学科专业:机械设计及理论所在单位:机械工程学院Mechanical fault diagnostics and Reliability
2、EngineeringMECHANICAL FAULT DIAGNOSIS BASED ONREDUNDANT SECOND GENERATION WAVELETPACKET TRANSFORM, NEIGHBORHOODROUGH SET AND SUPPORT VECTOR MACHINEby Li YiianSupervisor: Professor Dong MinYanshan UniversityMav,2013J 7摘要本文研究的冗余第二代小波包变换(RSGWPT),邻域粗糙集和支持向量机 (SVM) (NRS)关于故障检测,属性简约和模式分类中的应用。在此基础上,提 出了基于R
3、SGWPT, NRS和SVM的机械故障诊断的一种新方法,它利用RSGWPT从 小波包系数的统计特征中提取故障特征参数构成特征向量的,然后通过NRS法使 性能简约从而取得的主要的功能特性,最后这些主要特点被输入到SVM实现故障 模式分类。所提出的方法用于变速箱和汽油发动机气门列车的故障诊断,实验结 果表明,该方法可以提取故障特征,具有更好的分类能力,同时在保证分类准确 率的情况下,减少了很多的冗余功能,从而提高分类效率,达到更好的分类性能。关键词:冗余第二代小波包变换;邻域粗糙集;支持向量机;属性约简;故障诊 断I目录摘要I正文11引言12 RSGWPT和故障特征提取33NRS和故障特征选择54
4、SVM和故障特点分类85实验装置和数据采集105.1 变速箱试验105.2 汽油发动机气门列车实验106结果与讨论116.1 案例1:变速箱故障诊断116.2 案例2:汽油发动机气门的故障诊断136.3 讨论147结论17致谢18参考文献19H正文正文1引言随着越来越多的高品质的生产需求,有效的机械故障诊断已获得越来越多的 关注一般来说,在故障诊断系统中的两个重要步骤:第一个是故障的特 征提取的信号处理的,它决定了设备中是否发生了故障,第二步骤包括在先前步 骤中获得的特性的基础上的模式分类。目前故障诊断方法的有效性进行评估时, 精度和速度都是非常重要的因素,所以寻找一个准确且快速的方法进行故障
5、诊断 是一个重要的问题4, 5。总之,机械故障诊断的研究的目的是要找到一种方法, 通过使用从收集到的振动信号或声音发射信号中提取出来的故障特征来准确评 估故障类型。初期的故障特征往往是薄弱的,并且隐藏在背景信号中,所以很难用传统的 信号处理方法来检测它们。第二代小波变换 信GWT),由WimSweldens提出,是 一种新的小波构成方法,该方法采用了提升方案6, 7。它可以被看作是经典的 离散小波变换的替代实现。SGWT主要特征的是它提供了一个完全的空间域变换 的解释,而不是传统的基于频域的结构。SGWT的时间-频率分辨率随分解的水 平不同而变化。它在高频子带提供了良好的时间分辨率和较差的频率
6、分辨率,在 低频子带,它提供良好的频率和较差的时间分辨率。为了在故障特性经常存在的 高频子带中获得更高的分辨率,第二代小波包变换(SGWPT)已经建成,并因此 在每个层次的细节系数进一步分解从而得到的近似和细节分量4, 8-10 o不幸的 是,SGWPT没有时间不变性。延迟信号的分解结果,将原始信号的时移版本。这 可能会导致有用信息的特征提取与故障诊断故障的损失。冗余提升方案具有时间 不变性并且克服了提升方案的缺点,摆脱了小碎步和算子和更新算预测的零填 充,这使得各级的近似和细节信号与原始信号的长度是相同的11T3。基于冗 余提升方案小波包变换不仅能提供更详细的局部时频信号的描述,还能抑制频率
7、 混叠成分的分析结果,因为它没有分割和合并步骤的分解和重建阶段14。因此, 统计功能从变换系数RSGWPT提取有更大的能力来检测故障信号。一般来说,震动信号是从 监控机械设备中获得的。A / D转换后,采样的 振动数据分解使用RSGWPTo从每一个所得到的子带的小波包系数,可以计算出 统计特性来描述每一个尺度的信号的特性。这些统计功能,可以直接在故障诊断 应用,但在实践中,从实验数据中提取的某些功能通常是不完善和冗余,甚至对 方不兼容,因此给机械故障诊断带来了很多问题,如高计算复杂度,低的识别速 度和识别效果弱。因此,调查故障诊断方法使用较少的属性值,而不错过任何故 障信息是非常有意义的。粗糙
8、集理论是解决不确定的,不精确的和模糊的信息问题的一种有用的数学 工具,这方面的知识是关于域名COM分区和粒度特性。粗糙集理论的特点是创造 近似描述对象进行数据分析,优化和认可,它不需要域或先验知识。因此,使用 这种方法可以评估的各种属性的重要性,并没有额外的知识,除了保留一些关键 的属性所需提供的数据5。目前,粗糙集方法已经应用在许多领域,如机器故 障诊断,股市预测,决策支持系统,医疗诊断,数据过滤和软件工程15T8。然而,经典的粗糙集模型可以被用来处理类别特征是离散值,因此机械故障 诊断系统中基于粗糙集的特征选择,因此有必要引入一个离散化算法分区值域实图1.提出的故障诊断程序的流程图。时值的
9、变量分成几个时间间隔,然后把它们看作是类别特征。近年来许多的数值 属性离散化方法已经提出,包括等距离方法,频率相同的方法,最大燧法等16。 显然,离散化的数值属性,可能会导致信息丢失,因为不考虑数值离散化值的隶 属度19, 20。为了解决这个问题,引入一个邻域粗糙集模型,可以处理分类和 数值属性,一些属性约简算法基于此模型已建成20-22。由于数值属性离散化 的情况下的故障信息,邻域粗糙集模型的属性约简,可以完美地保存原始数据。对于分类,文献中有很多关于特定应用的最优分类,例如齿轮故障检测,复 式压缩机,滚子轴承等等23-26。支持向量机(SVM)是一种新型的基于统计学 习理论的机器学习方法,
10、该方法最初是由Vapnik介绍,其中有许多吸引人的特 点,包括其理想的分类能力有少量的非现场故障数据样本和有前途的实证表现线性和高维模式识别27-29L根据上述原则,本文提出了一个新的机械设备故障诊断方法,它利用了构建 在冗余提升方案和SGWPT基础之上的RSGWPT来提取故障特征参数,然后运用NRS 方法使故障特征向量的属性约简从而获得主要特性参数,最后这些主要特点输入 到SVM实现故障模式分类。故障诊断方法的流程图建议如图1所示。本文的其余部分安排如下:在第2节中,统计参数的小波包变换系数作为故 障特征由RSGWPT所提取。在第3节,通过NRS法,特征向量被减少,一些关键 功能也被选择了进
11、行分类。在第4节中,我们使用支持向量机来完成基于选定的 功能基础上的故障模式分类。在第5及6中,所提出的方法应用于汽油发动机的 齿轮箱和气门传动诊断不同的状态,结果已验证建议的诊断方法的有效性。最后 的结论将在第7节中进行了简要叙述。2 RSGWPT和故障特征提取为了完成RSGWPT的构建,本节介绍SGWPT和冗余提升方案。首先SGWPT的分解和重建阶段说明如下。SGWT的分解阶段包括三个步骤: 分裂,预测和更新。在分解阶段,根据下面的公式计算每个子带的系数:X/+i,2 =X/l0-P(X/je)Xi+i,i =X/e + U(X/+i,2)< =X'2o-P(X|,2%)=X
12、|,2% + U(X/+,2»i)SGWT的重建阶段的是分解阶段的一个相反的过程,其中包括逆更新步骤, 逆预测步骤和合并步骤。在重建阶段,将要被重构子带系数被保留,然后其他子 带系数被设定为零。最后,重构的结果,由公式如下得到:=X/+ 2仆1_1 U(X(+.21Al)X1,2b =Xj+i+ P(X1,2%)Xl2i(2i) = Xl2ieX|,ie =Xi + ij -U(X + i,2)X.o =X/ + 1 2+P(X1,ie)X/,1 ) = X/eX:(2i + 1) =X,io图2 . RSGXVPT的分解阶段其次,在冗余提升方案,小碎步将被丢弃。假设P和U】分别代表
13、1级冗余 提升方案的预测和更新算子,则户和U1的系数通过用零11填充P和U初始算 子P;和U,得到。pli=,p?.oo 虑.就凤-12(-12,-12f-lq = ng, u?,,,叫 _2。,021-12f-12'-1运用冗余提升计划和SGWPT, RSGWPT很容易被构造。RSGWPT在1级的预测 步骤和更新步骤通过使用P和U1进行预测,表示如下:'X/十,2 =X/,1-P'-1(X,1)X|+1,1 =1+户1(为 + 1,2)< 为十 i2+i =X/2i-P/+1(X/2i)X1+12+1_i =方力 + U'+1 (X/+胃+1)RSGWP
14、T的重建阶段可以从其分解阶段得到,并且可以由下述公式表达:J 1,2' =(V2)(X/+1 2i+i_1-U,+1(X/+1 2i+i)+X/+1 2i+i +P/*1(Xf+1 2r+i_1-U,*1(Xf+1 2r+i)X/,1=(l/2)(X/+1,1-U/+,(X/+1,2)+Xf+1,2+P/+,(X/+1,1-U,+ 1(X/+1,2)RSGWPT的分解和重构的阶段分别由图中2和图3表示。由于没有RSGWPT分解阶段的分割操作,所有各级的近似和细节系数与输入 信号具有相同的长度。因此,RSGWPT的分解结果具有时间不变性,并且保持完 美的原始信号的信息。此外,它还具有抑制
15、频率混叠性能14,这是非常有用的 故障特征提取与故障诊断。使用RSGWPT对采样的振动数据分解,从每一个所得 到的子带的小波包系数,九个统计特性包括峰值,均值,标准差,均方根,形状 系数,偏度,峰度,峰值因数和脉冲指数指的是在每一个尺度来描述的信号的分 布特性。如果RSGWPT的分解水平是1,即,子带的小波包的数目为* 和n统计特 性从各子带中提取的,那么所形成的故障的特征向量有乘以n功能的参数。根 据这些统计特性,可以得到的特征向量,并且,故障的原因和类别的特征向量的 映射关系,也可构建。3 NRS和故障特征选择基于粗糙集的特征选择的动机是选择一个最小的属性子集,它作为整个属性 集具有相同的
16、特征。古典粗糙集模型只适合处理的分类功能。虽然作为数字数据, 它广泛存在于机械故障诊断的实际应用,但有必要转换数值属性分类属性。研究 人员通常采用离散化算法值域实值变量分成几个间隔分区,然后把它们看作是类 别特征。这种转换可能不仅会导致信息丢失,而且也使得获得的结果在很大程度 上取决于离散效应。为了解决这个问题,邻里粗糙集模型被提出20, 21。下面介绍了它的基本 概念和原理。在信息系统U,A中,U是一个非空有限集合,被称为总体的样本XLX2,., XN), A是一组属性Al, A2, AM(也称为功能),它表明了样品特性。具体来说,U, A也称为决策表,如果A二CD,其中C为条件属性和D为决
17、策属性 集。给定一个xi GU, BC功能空间B中的xi的邻域被定义为:5(七)=化H引其中6是邻域的大小,而/是一个距离函数。特别的,如果B是分类的功能 集,则6=0。一般归结距离被定义:一 nZe%-%)”_k=l一其中xik和xjk分别是xi和xj在功能为k时的值。P是一个不小于1的实 数。根据所选择的不同的P值得到三种距离函数:(1)当P=1时,为曼哈顿距离 /I; (2)当P=2时,为欧氏距离/2; (3)当P=8时,为切比雪夫距离,3。SB(xi)是以样本xi为中心的邻域信息点,它受两个重要因素的影响。一个是所使用的距离函数,该函数决定邻域的形状,另一个是阈值8,用于控制邻 域的大
18、小。如果属性集A中同时包含数值属性B1和分类属性B2,那么由Bl, B2 和B1PIB2得出的样本x的邻域被定义为:鼠(x)= k /(占用) 5,% £ u4(K)二句金式为)=。,玉名叫(x)= k(x,xj K 5c 殳(x,x) = 0从上面的定义中,我们可以看到,样本领域内在类别特征的数值 的功能很有多相同的值,其距离比阈值3低。总体中的领域关系式N 可以写成相关矩阵(旬%也的形式,如Eq中所示。邻域关系中绘制在一 起的对象在距离上具有相似性或不可分辨性。=A g*0, otherwise如果在系统中有一个属性在总值范围内生成一个邻域关系,那么 这个信息系统被称为邻域信息系
19、统,记为NIS=U,A,N。具体来说, 如果在系统中有两种属性集,即人=(20口,其中C为条件属性集, D是决策属性集,那么至少存在一个条件属性,它包括总体中的邻域 关系,邻域信息系统也被称为邻域决策系统,它由NDT= U, CUD,N表示。假设xi,X2,xm是命令1-N的对象子集,8 b(Xx) 是由属性子集BcC产生的样本我的邻域信息。相对于属性B的定值 D的下限和上限的约数被定义为:=3 n 用。=5瓦 X, 其中NpX = %怜8(Z)=*,七 £(/,NrX = Maa)cX r。,七 £定值的下限约数也被称为定值正区域,用POSB(D)表示。对 象的子集邻域地
20、区均属于决策类。与此相反,位于邻域子集的边界区 域(BN(D)= NbD-NbD)的样本来自于一个以上的决策类中。因此邻 域模型将样本分为正区域和决策类边界地区。正区域中的样品可分为 没有不确定性的决策类和不能被确定的边界区域分类,后者很容易被 错误分类。因此,边界区域越大,条件属性的特征越不明显。为了表 征功能的显著性,DB的依赖程度被定义为:其中I I是一个集合的基数。Yb(D)反映B近似D的能力,它可 以用来测量特征子集的重要性。Sig/,5,O)=(0九(0运用的上述属性评价方法,本文采用了基于邻域粗糙集模型的属 性约简算法计算正无穷的数值。从一个空集开始,该算法计算每一轮 的所有其他
21、属性中该属性的比重,并具有将比重最大的属性加入到还 原集的功能。直到该属性的比重在其余所有属性中的比重接近于零, 该算法才停止。该算法被描述如下:步骤1:给定一个邻域决策系统NDT=<U,A,N>作为输入,其中U是 一个非空有限的样本集Xi,X2,xu; A=CUD,其中C为条件 属性集,D是决策属性集;步骤2:对于Va£A,计算Na的领域关系;步骤3:让。一red,其中ied是包含所选属性的区域;步骤4:属性的子集的比重可以由SIG功能定义。对于每一个af A-red,计算其 Sig(ai,red,D)= Yredua(D)- yred(D);步骤5:属性选择的目的是要
22、搜索属性的一个子集,使得它的分类问 题在所选择的特征空间中具有最大的一致性,属性评价的措施之一是 选择一个满足下式Sig(ak,red,D尸max(Sig(ai,red,D)的加值。步骤6:如果Sig(ak,red,D)>0, redUakred,那么程序转到步骤4; 否则回到ied的还原效果。图3. RSGWPT的重建阶段4 SVM和故障特点分类SVM是一种统计机器学习技术,已被广泛应用在模式识别领域。让 (%, y) 1 = 12,N表示一个实验样本集s,每一个样本%属于集合y”-1,1。 支持向量机的目标是要找到一个超平面来划分S,例如,与同一个标签上的所有 点的超平面的同一侧上,
23、而这两个类A和B之间与和超平面的距离达到最大化。 两个数据集的最优超平面的一个例子如图4所示。如该图4所示,圆圈和方块分别代表两类相同的样本点0 H是一个分离平面。 H1和H2是分别平行于H的平面,并且通过这些级别中离H最近的采样点。H1和 H2之间的距离被定义为余量。具有泛化误差最小的最佳分离平面是一个不仅能 正确将所有样本点分为这两级,而且但也在H1和H2之间的留有最大分差。应用 核函数,支持向量机可以用于非线性分类任务中。其基本思想是使用非线性变换 到一个高维特征空间的变换输入向量,然后在特征空间中进行线性分离。为了构 建一个非线性支持向量分类模型,内积(xi,x)的被核函数K(xi,x
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 故障诊断 可靠性工程 作业 吐血 分享
链接地址:https://www.31doc.com/p-12080368.html