基于混合型优化参数的群智能寻优故障诊断系统技术方案

本发明专利技术公开了一种混合型优化参数的群智能寻优故障诊断系统，用于对田纳西伊斯曼过程进行故障诊断，包括数据预处理模块、主成分分析模块、相关向量机模块以及群智能算法模块。本发明专利技术对田纳西伊斯曼化工过程的重要参数指标进行故障诊断预报，克服已有的化工故障诊断技术仪表测试敏捷性差、难以找到系统最优参数从而影响故障诊断效果的不足，引入群智能算法模块对相关向量机参数寻优，从而得到混合型优化参数的群智能寻优故障诊断系统，实现的田纳西伊斯曼过程故障诊断系统测试敏捷性好、易于找到全局最优解。

Swarm intelligence optimization fault diagnosis system based on mixed optimization parameters

The invention discloses a swarm intelligence optimization fault diagnosis system with mixed optimization parameters, which is used for fault diagnosis of Isman process in Tennessee, including data preprocessing module, principal component analysis module, correlation vector machine module and swarm intelligence algorithm module. The invention makes fault diagnosis and prediction for important parameter indexes of Isman chemical process in Tennessee, overcomes the shortcomings of poor test agility of existing chemical fault diagnosis techniques and instruments, and is difficult to find the optimal parameters of the system, thereby affecting the effect of fault diagnosis, and introduces swarm intelligence algorithm module to optimize the parameters of the correlation vector machine, thereby obtaining the optimal parameters of the correlation vector machine. To the mixed optimization parameter swarm intelligence optimization fault diagnosis system, the Tennessee Eastman process fault diagnosis system has good test agility and is easy to find the global optimal solution.

【技术实现步骤摘要】
基于混合型优化参数的群智能寻优故障诊断系统
本专利技术涉及故障诊断领域、机器学习领域和群智能优化算法领域，尤其涉及一种结合机器学习和群智能优化算法的田纳西伊斯曼过程化工故障诊断系统。
技术介绍
随着工业现代化技术的快速发展，工业自动化系统也不断的升级，越来越多的复杂工业生产过程变得更加可靠，生产变得更加稳定，能够实现更好的控制机器运作。最近几十年当中，微型控制器，传感器，通信总线，集成器等技术的飞速发展。这些元器件变得越来越便宜，性能也越来越稳定可靠。这在很大程度上促进了当代工业自动化系统的发展，促使自动化系统成为工业生产过程中非常重要组成部分。工业生产正向着大型集成化、复杂多样化和精密化的大方向趋势发展。虽然这种趋势能够提高工业生产效率，带来经济、环境和社会效益，但工业生产过程的复杂性也随之增加，受到的干扰因素也逐渐增多，各种各样的故障发生概率也大大增加。一旦出现某个小故障，若不能及时发现故障源，很可能会形成联锁反应，轻则造成机器损坏，经济损失。重则导致人员伤亡，对周围环境造成持续性污染，周围居住的人无法生活。近几年，随着工业过程领域科学技术的革新，工业装置的状态监测、相应设备的监控、控制系统故障检测与诊断等得到快速发展。现如今，工业生产规模不断扩大，生产流程变得越来越复杂。准确、快速有效的故障诊断已成为当前对产品质量和工厂稳定运行至关重要的元素。同时，在工业过程机理结构不断趋向复杂化前提下，生产过程建立准确机理模型增加了难度。因此，促使数据驱动故障检测方法的得到迅速发展。
技术实现思路
为了克服目前已有的故障诊断技术的测试敏捷性差、难以找到系统最优参数从而影响诊断效果的不足，本专利技术的目的在于提供一种系统测试敏捷性好、易于找到系统全局最优参数的群智能寻优故障诊断系统。本专利技术解决其技术问题所采用的技术方案是：一种混合型优化参数的群智能寻优故障诊断系统，用于对田纳西伊斯曼过程进行故障诊断，包括数据预处理模块、主成分分析模块、相关向量机模块以及群智能算法模块，其中：数据预处理模块：田纳西伊斯曼过程的52个变量为数据预处理模块的输入。由于每个变量都有不同的单位，为了防止不同的量纲引起数据量级之间的误差，先对所有数据进行标准化处理，标准化公式如下：其中，mean表示各变量的算术平均值，std表示各变量的标准差，表示输入变量的值，下标i表示第i次检测、j分别表示第j维变量，xij表示标准化后输入变量的值作为输入数据。标准化后的数据为S＝{xi1,xi2,...xi52}。主成分分析模块：通过主成分分析来保证在不降低系统精度的情况下降低系统的复杂度。将标准化后的数据S＝{xi1,xi2,...xi52}进行主成分分析，保留85％的主要成分。相关向量机模块：用于建立诊断系统，提高系统测试敏捷性：在相关向量机中，假设y符合高斯分布，将高斯核函数应用到y(x)，其中yc为y的估计值，κ为函数的宽度参数，可以得到似然估计函数：：其中t＝(t1,t2...tN)T，w＝(w1,w2...wN)T，w表示权重系数，t表示分类类别共N类，P表示概率结果、x为输入数据、y为输出数据，下标i表示第i个数据，上标T表示矩阵的转置。在相关向量机中，权值w可以通过极大似然估计来确定，但是为了防止过拟合，通过高斯先验概率分布来约束参数：其中，α是N+1维超参数，该值的引入导致了系统的稀疏性。根据先验概率和贝叶斯定理，应用拉普拉斯原理计算出近似的后验概率：(1)由于p(w|t,α)∝p(t|w)p(w|α)，对当前固定的α值求出最大可能的权值wMP，采用二阶牛顿法求得wMP。其中，yi＝σ{y(xi；w)}，A＝diag(α0,α1,...αN)。(2)利用拉普拉斯方法，将对数后验概率进行二次逼近，将上式进行两次求导得到：其中，Φ＝[φ(x1),φ(x2)...φ(xN)]T是由多个核函数构成的一个N×(N+1)维度的结构矩阵，并且φ(xi)＝[1,K(xi,x1),K(xi,x2),...,K(xi,xN)]T。B＝diag(β1,β2,..βN)是一个对角矩阵，并且βi＝σ{y(xi)}[1-σ(xi)]。对上式右边取负号再求逆矩阵，即可得到协方差矩阵Σ。(3)利用wMP和Σ，对超参数α进行更新，更新公式如下：其中，γi≡1-αiΣii。最终，可以得到后验协方差为：其中Σ是后验协方差。相关向量机通过引入贝叶斯框架使得系统更加具有稀疏性，在视频追踪、图像检索领域中被广泛应用。同时，该系统的泛化性能较好，测试时间短，适合于在线监测。群智能算法模块：结合差分进化算法和粒子群算法优化相关向量机模块中的函数宽度参数κ，易于找到全局最优参数，实现具体步骤如下：(1)随机初始化种群，设置种群数量N＝100；设置差分进化算法缩放因子F＝0.7，交叉概率CR＝0.5，差分进化算法最大迭代次数50；设置粒子群学习速度c1＝c2＝2，惯性权重ω＝0.67，粒子群算法最大迭代次数100；设置k＝0,i＝0；(2)根据适应度函数来计算候选粒子种群的适应度值，适应度f(x)函数如下所示；其中，表示输出变量的真实值，表示输出变量的预测值，n表示样本个数。(3)根据以下步骤进行差分进化算法变异、交叉、选择操作；(3.1)变异操作：差分进化算法是通过差分的方式进行变异操作的。在当前种群中随机选取两个种群个体进行作差，将他们差向量进行缩放，缩放后的向量与第三个随机选取的种群个体进行求和，生成新个体，更新公式如下所示：Ni(k)＝Lm3(k)+F×(Lm1(k)-Lm2(k)),1≤m1≠m2≠m3≠i≤N(9)其中，F是缩放因子，m1，m2，m3均为1到N之间的整数，k表示迭代次数，Lm1(k)表示第m1个粒子的位置，Lm2(k)表示第m2个粒子的位置，Lm3(k)表示第m3个粒子的位置。通过变异之后，得到第i个变异中间种群Ni(k)。(3.2)交叉操作：根据原始种群和变异得到中间种群，进行个体间的交叉操作，得到交叉中间种群U，更新公式如下所示：其中，CR为交叉概率，处于0到1之间，ni,j(k)表示第i个中间种群第j维变量，li,j(k)表示第i个原始种群第j维变量，ui,j(k)表示交叉中间种群的第i个新种群第j维变量，ni,j(k)表示变异中间种群Ni(k)的第i个新种群第j维变量。(3.3)选择操作：差分进化算法采用了贪心算法，根据适应度值来决定进入下次迭代种群的个体：其中，Ui(k)表示第i个交叉中间种群，Li(k)表示第i个原始种群。(4)如果交叉中间种群的适应度值小于原始种群适应度值，则返回步骤(3)，否则，激活粒子群算法；(5)根据如下公式执行粒子群算法，更新粒子群速度、位置矢量和惯性权重；其中，Vi(k)是粒子在第k次迭代时的速度，Li(k)是粒子在第k次迭代时的位置，Pbest是粒子自身经验的局部最优解；gbest是所有粒子经验的全局最优解，w是惯性权重，c1和c2是学习速率，r1和r2是在0到1之间的随机数，wmax和wmin是惯性权重的最大值和最小值，itermax是最大迭代次数。(6)如果粒子群算法迭代次数小于终止迭代次数，则返回步骤(5)，否则，转到步骤(7)；(7)更新Pbest和gbest作为相关向量机参数的候选解，判断差分进化算法迭本文档来自技高网...

【技术保护点】
1.一种混合型优化参数的群智能寻优故障诊断系统，用于对田纳西伊斯曼过程进行故障诊断，其特征在于：包括数据预处理模块、主成分分析模块、相关向量机模块以及群智能算法模块。

【技术特征摘要】
1.一种混合型优化参数的群智能寻优故障诊断系统，用于对田纳西伊斯曼过程进行故障诊断，其特征在于：包括数据预处理模块、主成分分析模块、相关向量机模块以及群智能算法模块。2.根据权利要求1所述混合型优化参数的群智能寻优故障诊断系统，其特征在于：所述数据预处理模块的输入为田纳西伊斯曼过程的52个变量由于每个变量都有不同的单位，为了防止不同的量纲引起数据量级之间的误差，先对所有数据进行标准化处理，标准化公式如下：其中，mean表示各变量的算术平均值，std表示各变量的标准差，表示输入变量的值，下标i表示第i次检测、j分别表示第j维变量，xij表示标准化后输入变量的值作为输入数据。标准化后的数据为S＝{xi1,xi2,...xi52}。3.根据权利要求1所述混合型优化参数的群智能寻优故障诊断系统，其特征在于：所述主成分分析模块通过主成分分析来保证在不降低系统精度的情况下降低系统的复杂度。将标准化后的数据S＝{xi1,xi2,...xi52}进行主成分分析，保留85％的主要成分。4.根据权利要求1所述混合型优化参数的群智能寻优故障诊断系统，其特征在于：所述相关向量机模块用于建立诊断系统，提高系统测试敏捷性。在相关向量机中，假设y符合高斯分布，将高斯核函数应用到y(x)，其中yc为y的估计值，κ为函数的宽度参数，可以得到似然估计函数：其中t＝(t1,t2...tN)T，w＝(w1,w2...wN)T，w表示权重系数，t表示分类类别共N类，P表示概率结果、x为输入数据、y为输出数据，下标i表示第i个数据，上标T表示矩阵的转置。在相关向量机中，权值w可以通过极大似然估计来确定，但是为了防止过拟合，通过高斯先验概率分布来约束参数：其中，α是N+1维超参数，该值的引入导致了系统的稀疏性。根据先验概率和贝叶斯定理，应用拉普拉斯原理计算出近似的后验概率：(1)由于p(w|t,α)∝p(t|w)p(w|α)，对当前固定的α值求出最大可能的权值wMP，采用二阶牛顿法求得wMP。其中，yi＝σ{y(xi；w)}，A＝diag(α0,α1,...αN)。(2)利用拉普拉斯方法，将对数后验概率进行二次逼近，将上式进行两次求导得到：其中，Φ＝[φ(x1),φ(x2)...φ(xN)]T是由多个核函数构成的一个N×(N+1)维度的结构矩阵，并且φ(xi)＝[1,K(xi,x1),K(xi,x2),...,K(xi,xN)]T。B＝diag(β1,β2,..βN)是一个对角矩阵，并且βi＝σ{y(xi)}[1-σ(xi)]。对上式右边取负号再求逆矩阵，即可得到协方差矩阵Σ。(3)利用wMP和Σ，对超参数α进行更新，更新公式如下：其中，γi≡1-αiΣii。最终，可以得到后验协方差为：其中Σ是后验协方差。5.根据权利要求1所述混合型优化参数的群智能寻优故障诊断系统，其特征在于：所述群智能算法模块结合差分进化算法和粒子群算法优化相关向量机模块中的函数宽度参数κ，找到全局最优...

【专利技术属性】
技术研发人员：刘兴高，何世明，徐志鹏，张泽银，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33