基于量子遗传算法优化支持向量机的轴承故障诊断方法
【专利摘要】本发明的基于量子遗传算法优化支持向量机的轴承故障诊断方法，包含以下步骤：（1）采集轴承的振动信号；（2）计算无量纲指标；（3）基于云模型量子遗传算法优化支持向量机的模型参数C和σ；（4）对支持向量机模型进行训练；（5）利用支持向量机模型进行故障诊断；（6）输出轴承故障诊断结果。本发明的轴承故障诊断方法具有诊断准确率高的优点，为解决轴承故障诊断问题提供一种新的方法。
【专利说明】
基于量子遗传算法优化支持向量机的轴承故障诊断方法
技术领域
[0001] 本发明涉及轴承故障诊断方法，特别是涉及一种基于量子遗传算法优化支持向量 机的轴承故障诊断方法。
【背景技术】
[0002] 随着科学技术的进步，机械设备内部结构越来越复杂，功能越来越强大，智能化水 平越来越高。轴承是机械的核心部件。如果机械设备的轴承在运转过程中出现故障，可能导 致整个生产过程中断，不但影响企业的经济效益，而且可能危害工作人员的人身安全，甚至 带来灾难性的后果。为了避免机械轴承发生故障造成损失，就必须对机械轴承的运行状况 进行监测，能够在发生故障之前进行预测出故障或者对早期故障的快速诊断，以便采取相 应的措施避免造成重大损失。
[0003] 目前轴承故障诊断主要方法有基于神经网络、人工免疫、专家系统，模糊逻辑、遗 传算法、证据理论、支持向量机等诊断方法。现有的轴承故障诊断方法存在算法的训练时间 长且诊断的正确率有待提高的问题。
[0004] 针对现有技术中的存在问题，结合支持向量机和量子遗传算法的优点，提供一种 训练时间短、诊断正确率高的基于量子遗传算法优化支持向量机的轴承故障诊断方法，显 得尤为重要。
【发明内容】
[0005] 本发明的目的在于避免现有技术中的不足之处而提供一种训练时间短、正确率高 的基于量子遗传算法优化支持向量机的轴承故障诊断方法。
[0006] 本发明的目的通过以下技术方案实现：
[0007] 提供一种基于量子遗传算法优化支持向量机的轴承故障诊断方法，包括有以下步 骤：
[0008] 步骤一、采集轴承的振动信号
[0009] 利用测振系统采集轴承的振动信号；
[0010] 步骤二、计算无量纲指标
[0011] 对所采集的振动信号计算出5个无量纲指标:波形指标、脉冲指标、裕度指标、峰值 指标及峭度指标，得到训练样本数据集和测试样本数据集；
[0012] 步骤三、优化支持向量机的模型参数(：和〇
[0013] 利用云模型量子遗传算法(CQGA)对支持向量机的模型参数(：和〇进行寻优，得到最 优的参数(：和〇的取值；
[0014] 步骤四、对支持向量机进行训练；
[0015] 将最优的参数(：和〇的取值代入支持向量机模型中，利用训练样本数据对支持向量 机模型进行训练；
[0016] 步骤五、利用支持向量机模型进行故障诊断；
[0017] 将测试样本数据输入到训练好的支持向量机模型中进行故障诊断；
[0018] 步骤六、输出轴承故障诊断结果。
[0019] 优选的，步骤二中，计算5个无量纲指标所用公式如下：
[0025] 式中，|l|表不平均幅值，Xrms表不均方根，Xmax表不最大方根，Xr表不方根幅值，X表 示振动幅值，P(X)表示振动幅值的概率密度函数。
[0026] 优选的，步骤五中，将采集获得振动信号通过步骤二计算获得的无量纲指标作为 测试样本集。
[0027] 更优选的，支持向量机的核函数为径向基函数，径向基函数表示为
[0029] 式中，〇是核参数，即核函数的宽度。
[0030] 优选的，步骤三中采用云模型量子遗传算法对支持向量机参数(：和〇进行寻优，其 具体的寻优方法如下：
[0031] Stepl按照一维二进制编码方法随机产生初始种群，产生初始种群： 尸⑴= {"??_，/〇，n是种群的规模，//,(/ = 1，2,…㈨为种群第t代的一个个体，每个个体 采用量子比特幅编码的染色体结构表示如下：
[0033]其中，1为每个基因的量子比特数，第t代种群的第i个基因的基因云记为C=(EXl， 已取，1^)，其中&^取，1^为第1个基因的基因云的三个特征值；
[0034] Step2对种群P( t)中的所有个体进行测量，得到所有个体的值；
[0035] Step3利用支持向量机计算种群P(t)所有个体的适应度值；
[0036] Step4记录最优个体及其适应度值；
[0037] Step5判断是否满足结束条件，若满足结束条件，则输出最优的参数(：和〇的值，否 则继续计算；
[0038] Step6通过调整种群基因云的参数EXl，Em，Hei来优化子代种群产生的策略，利用 进化方程9 = Eru(pi-xi)+Exi(p-xi)+Hei(pj-xi)自动进行调整量子门的旋转角，按照量子的 旋转门对种群个体进行更新，获得子代种群P(t);其中：E Xl，Em，Hei为第i个基因的基因云 的三个特征值;P为全局极值;？1为个体所在种群极值;P」为邻域种群极值；
[0039] Step7t = t+1，转至Step2，直到最大迭代次数。
[0040]其中，Step2得到所有个体的值包含(：和〇的取值。
[0041]优选的，步骤四中将步骤二获得的无量纲指标作为训练样本集作为输入值训练支 持向量机。
[0042]本发明的有益效果：
[0043] 本发明的基于量子遗传算法优化支持向量机的轴承故障诊断方法对通过测振系 统采集轴承的振动信号，其中振动信号包括轴承的正常状态振动信号及故障状态振动信 号，故障状态的振动信号又包括轴承磨损、轴承外裂、轴承外裂的振动信号，并计算出脉冲 指标、裕度指标、峰值指标、波形指标、峭度指标，得到训练样本数据集和测试样本数据集； 利用云模型量子遗传算法(CQGA)优化支持向量机的模型参数(：和 〇，获得最优的参数(：和〇的 取值;将最优的参数(：和〇的取值代入支持向量机模型中，利用训练样本数据对支持向量机 模型进行训练;将测试样本集输入到训练的支持向量机轴承故障诊断模型中进行测试，得 到诊断结果。
[0044] 无量纲指标对于振动信号中的扰动不敏感，能很好地反映机械轴承的运行状态， 而且几乎不受振动信号绝对水平和工况、载荷等变化的影响。因此利用无量纲指标进行故 障诊断的性能较为稳定。
[0045]支持向量机在小样本分类中具有明显的优势，采用基于云模型量子遗传算法对支 持向量机模型的惩罚参数C和核函数参数〇进行优化，避免模型参数选择的盲目性，提高模 型的诊断精度。
【附图说明】
[0046] 利用附图对本发明做进一步说明，但附图中的内容不构成对本发明的任何限制。
[0047] 图1是本发明的基于量子遗传算法优化支持向量机的轴承故障诊断方法的流程示 意图。
[0048] 图2是本发明的基于量子遗传算法优化支持向量机的模型参数(：和〇的流程示意 图。
【具体实施方式】
[0049] 结合以下实施例对本发明作进一步说明。
[0050] 实施例1
[0051] 参见图1，本实施例的基于量子遗传算法优化支持向量机的轴承故障诊断方法包 括有以下步骤：
[0052]步骤一、采集振动信号
[0053]利用测振系统采集轴承的振动信号；
[0054] 步骤二、计算无量纲指标
[0055]利用所采集的振动信号计算出5个无量纲指标:波形指标、脉冲指标、裕度指标、峰 值指标及峭度指标，得到训练样本数据集和测试样本数据集；
[0061] 式中，|l|表不平均幅值，Xrms表不均方根，Xmax表不最大方根，Xr表不方根幅值，X表 示振动幅值，P(X)表示振动幅值的概率密度函数。
[0062]步骤三、获得最优模型参数(：和〇
[0063]本实施例的支持向量机故障诊断模型中的参数C和〇取值对支持向量机模型泛化 性能起着关键作用。为了得到最优的支持向量机的模型参数(：和〇的取值，利用云模型量子 遗传算法对支持向量机的模型参数(：和 〇进行寻优。
[0064]具体的寻优方法如下：
[0065] Stepl按照一维二进制编码方法随机产生初始种群，产生初始种群： /^)=丨尤'~,/,：}，11是种群的规模，乂(_/ = 1，2^)为种群第七代的一个个体，每个个体 采用量子比特幅编码的染色体结构表示如下：
[0067] 其中，1为每个基因的量子比特数，第t代种群的第i个基因的基因云记为C=(EXl， 已取，1^)，其中&^取，1^为第1个基因的基因云的三个特征值；
[0068] Step2对种群P(t)中的所有个体进行测量，得到所有个体的值，得到所有个体的值 包含C和〇的取值。
[0069] Step3利用支持向量机计算种群P(t)所有个体的适应度值；
[0070] Step4记录最优个体及其适应度值；
[0071] Step5判断是否满足结束条件，若满足结束条件，则输出最优的参数(：和〇的值，否 则继续计算；
[0072] Step6通过调整种群基因云的参数EXl，Em，Hei来优化子代种群产生的策略，利用 进化方程9 = Eru(pi-xi)+Exi(p-xi)+Hei(pj-xi)自动进行调整量子门的旋转角，按照量子的 旋转门对种群个体进行更新，获得子代种群P(t);其中：EXl，Em，Hei为第i个基因的基因云 的三个特征值;P为全局极值; ？1为个体所在种群极值;P」为邻域种群极值；
[0073] Step7t = t+1，转至Step2，直到最大迭代次数。
[0074]步骤四、将步骤三获得最优的参数(：和〇的取值代入支持向量机模型中，将步骤二 获得的无量纲指标作为训练样本集作为输入值，对支持向量机进行训练；
[0075]步骤五、将步骤三获得的测试样本集输入到训练好的支持向量机中进行故障诊 断；
[0076]步骤六、输出轴承故障诊断结果。
[0077]支持向量机是建立在VC维理论和结构风险最小化原理基础上，在解决小样本、非 线性和高维模式识别等问题中有着很多独特的优势.支持向量机的基本原理是首先将实际 问题通过一个非线性映射巾：Rd-H，把样本空间映射到一个高维的特征空间中（Hilbert空 间），然后在这个高维的特征空间中构造最优分类超平面。设有线性不可分的训练样本集 {(xi，yi)，（X2，y2)，…，（x n，yn)}，其中，xiGX = Rd，X为输入空间，输入空间中的每一个xi由d 维特征向量组成，yiGY= {-1，1}，i = l，???，n。通过非线性映射巾：Rd-H得到的高维特征空 间的训练样本集表示为
[0078] {(伞（xi)，yi)，（伞（X2)，y2)，???，（伞（xi)，yi)}
[0079]在高维特征空间构造最优分类超平面的约束条件为：
[0080] yi[ ? ? (J) (xi)+b]-l+Ci^:〇
[0081] 其中L为松弛变量，i = l，2，...，n
[0082]求解最优分类超平面问题可转化为优化问题：
[0085] 其中，《为权系数向量，C(C>0)为规则化(Regularization)参数，b为阈值，表示 向量的内积，li为松弛变量。
[0086] 通过引入拉格朗日函数求解，利用Karush-Kuhn-Tucher(KKT)条件，构造最优分类 超平面可以转化为对偶二次规划问题，表示如下：
[0089]求解上述对偶二次规划问题得到分类的决策函数为：
[0091 ]其中，为确定最优分类超平面的参数。
[0092]根据泛函的有关理论和Mercer条件，可以用核函数K(Xi，Xj)代替高维特征空间中 的内积函数[4>(Xl) ?巾（&)]，得到最终的分类决策函数为：
[0094]支持向量机只能解决二值类分类问题，而故障诊断是多类分类问题，这里采用一 对多方法解决。一对多方法的主要思想:对于K类问题构造出K个两类分类器，在构造第i个 分类器时，将第i类的训练样本作为正类，其他训练样本作为负类。对给定一个测试样本进 行分类时，将其分别输入到K个分类决策函数.
来计算得到的 K个函数值，该测试样本的类别被判定为K个函数值中最大的那个值所属的类别。
[0095] 实施例2
[0096]本实施例的主要技术方案与实施例1基本相同，在本实施例中未作解释的特征，采 用实施例1中的解释，在此不再进行赘述。本实施例中，
[0097]将采集获得振动信号通过步骤二计算获得的无量纲指标作为步骤五的测试样本 集。
[0098]最后应当说明的是，以上实施例仅用于说明本发明的技术方案说明而非对权利要 求保护范围的限制。本领域的普通技术人员参照较佳实施例应当理解，并可以对本发明的 技术方案进行修改或者等同替换，但属于本发明技术方案的实质相同和保护范围。
【主权项】
1. 基于量子遗传算法优化支持向量机的轴承故障诊断方法，其特征在于:包括有以下 步骤： 步骤一、采集轴承的振动信号 利用测振系统采集轴承的振动信号； 步骤二、计算无量纲指标 对所采集的振动信号计算出5个无量纲指标:波形指标、脉冲指标、裕度指标、峰值指标 及峭度指标， 得到训练样本数据集和测试样本数据集； 步骤三、优化支持向量机的模型参数(：和〇 利用云模型量子遗传算法对支持向量机的模型参数(：和〇进行寻优，得到最优的参数C 和〇的取值； 步骤四、对支持向量机进行训练； 将最优的参数(：和〇的取值代入支持向量机模型中，利用训练样本数据对支持向量机模 型进行训练； 步骤五、利用支持向量机模型进行故障诊断； 将测试样本数据输入到训练好的支持向量机模型中进行故障诊断； 步骤六、输出轴承故障诊断结果。2. 根据权利要求1所述的基于量子遗传算法优化支持向量机的轴承故障诊断方法，其 特征在于:步骤二中，计算5个无量纲指标所用公式如下：式中，|X|表不平均幅值，Xrms表不均方根，Xmax表不最大方根，Xr表不方根幅值，X表不振 动幅值，P(X)表示振动幅值的概率密度函数。3. 根据权利要求1所述的基于量子遗传算法优化支持向量机的轴承故障诊断方法，其 特征在于:步骤五中，将采集获得振动信号通过步骤二计算获得的无量纲指标作为测试样 本集。4. 根据权利要求1所述的基于量子遗传算法优化支持向量机的轴承故障诊断方法，其 特征在于：所述支持向量机的核函数为径向基核函数，其 中，〇是核参数，即核函数的宽度。5. 根据权利要求1所述的基于量子遗传算法优化支持向量机的轴承故障诊断方法，其 特征在于:步骤三中采用云模型量子遗传算法对支持向量机参数(：和 〇进行寻优，其具体的 寻优方法如下： Stepi按照一维二进制编码方法随机产生初始种群，产生初始种群：作）=!p:.， n是种群的规模，40=1,2:,???,?)为种群第t代的一个个体，每个个体采用量子比特幅编码 的染色体结构表示如下：其中，1为每个基因的量子比特数，第t代种群的第i个基因的基因云记为 0=江^3取，1^)，其中&^取，1^为第1个基因的基因云的三个特征值； Step2对种群P(t)中的所有个体进行测量，得到所有个体的值； Step3利用支持向量机计算种群P(t)所有个体的适应度值； Step4记录最优个体及其适应度值； Step5判断是否满足结束条件，若满足结束条件，则输出最优的参数(：和〇的值，否则继 续计算； Step6通过调整种群基因云的参数EXl，Em，Hei来优化子代种群产生的策略，利用进化 方程9 = Em(pi-xi)+Exi(p-xi)+Hei(pj-xi)自动进行调整量子门的旋转角，按照量子的旋转 门对种群个体进行更新，获得子代种群P(t);其中：E Xl，Em，Hei为第i个基因的基因云的三 个特征值;P为全局极值;个体所在种群极值;P:为邻域种群极值； Step7 t = t+l，转至Step2,直到最大迭代次数。 其中，Step2得到所有个体的值包含C和〇的取值。6. 根据权利要求1所述的基于量子遗传算法优化支持向量机的轴承故障诊断方法，其 特征在于:步骤四中将步骤二获得的无量纲指标作为训练样本集作为输入值训练支持向量 机。
【文档编号】G01M13/04GK106053067SQ201610352279
【公开日】2016年10月26日
【申请日】2016年5月24日
【发明人】朱兴统, 熊建斌, 许波
【申请人】广东石油化工学院