[发明专利]一种基于GA-ICA的新型分散式非高斯过程监测方法有效

专利信息
申请号: 201910229753.2 申请日: 2019-03-12
公开(公告)号: CN111695229B 公开(公告)日: 2023-10-17
发明(设计)人: 唐俊苗;童楚东;史旭华 申请(专利权)人: 宁波大学
主分类号: G06F30/20 分类号: G06F30/20;G06F18/2134;G06N3/126
代理公司: 暂无信息 代理人: 暂无信息
地址: 315211 浙江省*** 国省代码: 浙江;33
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 基于 ga ica 新型 分散 式非高斯 过程 监测 方法
【权利要求书】:

1.一种基于GA-ICA的新型分散式非高斯过程监测方法,其特征在于,包括以下步骤:

离线建模阶段包括步骤(1)至步骤(7):

步骤(1):在生产过程对象正常运行状态下,采集n个样本数据x1,x2,…,xn组成训练数据矩阵X=[x1,x2,…,xn]∈Rm×n,并计算均值向量μ=(x1+x2+…+xn)/n以及标准差向量δ∈Rm×1

其中,⊙表示向量(xi-μ)与(xi-μ)中的对应元素相乘、xi∈Rm×1表示第i个样本数据、m为测量变量的总个数、i=1,2,…,n、R为实数集,Rm×n表示m×n维的实数矩阵;

步骤(2):根据公式对样本数据x1,x2,…,xn分别实施标准化处理,从而得到矩阵其中,表示向量(xi-μ)与标准差向量δ中的对应元素相除,为标准化后的数据向量;

步骤(3):根据生产过程对象的组成结构,确定各个生产单元所涉及到的测量变量,以此将m个测量变量分成B个变量子块,并根据这B个变量子块从矩阵中选取相应的行向量分别组成B个子块矩阵

步骤(4):设置分离向量的个数为D,并利用十进制的遗传算法求解矩阵的D个分离向量w1,w2,…,wD、D个混合向量p1,p2,…,pD、以及D个独立成分向量t1,t2,…,tD

步骤(5):根据D个分离向量w1,w2,…,wD,依次对B个变量子块矩阵实施独立成分的提取,具体的实施过程如下所示:

步骤(5.1):初始化d=1;

步骤(5.2):根据步骤(3)中的B个变量子块,从第d个分离向量wd中选取相应的元素分别组成B个变量子块对应的分离向量wd(1),wd(2),…,wd(B);

步骤(5.2):针对b∈{1,2,…,B}的所有取值,分别判断是否满足条件:若否,则依次根据公式与分别计算第b个变量子块对应的独立成分向量td(b)与载荷向量ad(b)后,跳转至步骤(5.3);若是,则将第b个变量子块对应的分离向量w1(b),w2(b),…,wd(b)组成分离矩阵Wb,并将已计算出的第b个变量子块对应的所有独立成分向量与所有载荷向量分别组成独立成分矩阵Tb与载荷矩阵Pb后,再跳转至步骤(5.4);

步骤(5.3):根据公式更新第b个变量子块矩阵后,跳转至步骤(5.5);

步骤(5.4):保留第b个变量子块对应的独立成分矩阵Tb、分离矩阵Wb、以及载荷矩阵Pb

步骤(5.5):判断是否满足条件:d<D?若是,则置d=d+1后,返回步骤(5.2);若否,则得到所有B个变量子块对应的分离矩阵W1,W2,…,WB及载荷矩阵P1,P2,…,PB

步骤(6):先根据公式Λb=TbTbT/(n-1)计算协方差矩阵Λ1,Λ1,…,ΛB后,再根据公式Qb=diag(TbTΛb-1Tb)计算监测统计指标向量Q1,Q2,…,QB,其中diag()表示将矩阵对角线的元素转变成向量的操作,上标号T表示矩阵或向量的转置;

步骤(7):分别对监测统计指标向量Q1,Q2,…,QB中的元素按大小进行降序排列,并将第n/100个最大值分别做为各变量子块的监测统计指标上限CL1,CL2,…,CLB

在线监测阶段执行如下所示步骤:

步骤(8):收集新采样时刻的样本数据x∈Rm×1,并根据公式对x实施标准化处理得到向量

步骤(9):根据步骤(3)中的B个变量子块,可对应地将向量分成B个子向量

步骤(10):先根据公式分别计算各子向量对应的独立成分向量y1,y2,…,yB,再根据公式计算监测指标

步骤(11):针对下标号b的所有取值,判断是否满足条件:若是,则当前采样时刻生产过程对象是正常运行的,返回步骤(8)继续监测下一时刻的样本数据;若否,则当前采样时刻生产过程对象进入异常运行状态。

2.根据权利要求1所述的一种基于GA-ICA的新型分散式非高斯过程监测方法,其特征在于,所述步骤(4)的详细实施过程具体为:

步骤(4.1):设置分离向量的个数D与遗传算法的最大迭代次数U,初始化d=1与u=0,并在区间[-1,1]上随机生成一个N×m维的实数矩阵W,其中N为种群个数;

步骤(4.2):根据公式求解矩阵F∈RN×n,并计算矩阵F中n个列向量的均值向量f=(F1+F2+…+Fn)/n,其中Fi表示矩阵F中的第i列向量;

步骤(4.3):找出均值向量f中元素的最大值fmax,并记录该最大值在均值向量f中所在的位置υ;

步骤(4.4):判断是否满足条件u<U?若是,则置u=u+1后,跳转至步骤(4.5);若否,将实数矩阵W中的第υ行向量W(υ)记录为第d个分离向量wd,并跳转至步骤(4.13);

步骤(4.5):依次根据公式γk=fk/(f1+f2+…+fN)与公式计算选择概率γk与累计概率其中k=1,2,…,N、fk表示均值向量f中的第k个元素;

步骤(4.6):在区间(0,1)上随机生成一个数s后,从中找出满足条件的最小值并将实数矩阵W中的第k行向量保留;

步骤(4.7):重复实施步骤(4.6)N次,从而将保留的N个行向量组成矩阵V∈RN×m,并初始化j=1;

步骤(4.8):在区间(0,1)上随机生成一个数c后,判断是否满足条件:其中表示交叉概率;若是,则按照如下所示公式②对矩阵V中第j行向量Vj与第j+1行向量Vj+1实施交叉操作:

上式中,θj为区间(0,1)上随机生成的一个数;若否,则矩阵V中第j行向量Vj与第j+1行向量Vj+1保持不变;

步骤(4.9):判断是否满足条件:j<N-1?若是,则置j=j+2后,返回步骤(4.8);若否,则得到实施交叉操作后的矩阵V,并初始化k=1;

步骤(4.10):在区间(0,1)上随机生成一个数g后,判断是否满足条件:其中表示变异概率,建议取值若是,则按照如下所示公式③对矩阵V中第k行向量Vk的任意一个元素V(k,ε)实施变异操作:

上式中,ε表示任意选择的元素V(k,ε)在向量Vk中的位置,rk为区间(0,1)上随机生成的一个数,函数h(u)=rk(1-u/U)2;若否,则保持矩阵V中第k行向量V(k)不变;

步骤(4.11):判断是否满足条件:k<N?若是,则置k=k+1后,返回步骤(4.10);若否,则得到实施变异操作后的矩阵V;

步骤(4.12):依据公式VN=W(υ)更新矩阵V中最后一行的行向量VN,其中W(υ)表示实数矩阵W中的第υ行向量,再设置实数矩阵W=V后,返回步骤(4.2);

步骤(4.13):判断是否满足条件:d<D?若是,则依次根据公式和分别计算第d个独立成分向量td、混合向量pd、和更新矩阵再置d=d+1与u=0,并在区间[-1,1]上随机生成一个N×mb维的实数矩阵W后,返回步骤(4.2);若否,则得到所有D个分离向量w1,w2,…,wD、D个混合向量p1,p2,…,pD、以及D个独立成分向量t1,t2,…,tD

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于宁波大学,未经宁波大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/201910229753.2/1.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top