1.一种基于GA‑ICA的新型分散式非高斯过程监测方法,其特征在于,包括以下步骤:离线建模阶段包括步骤(1)至步骤(7):
步骤(1):在生产过程对象正常运行状态下,采集n个样本数据x1,x2,…,xn组成训练数m×n m据矩阵X=[x1,x2,…,xn]∈R ,并计算均值向量μ=(x1+x2+…+xn)/n以及标准差向量δ∈R×1:
m×1
其中,⊙表示向量(xi‑μ)与(xi‑μ)中的对应元素相乘、xi∈R 表示第i个样本数据、m为m×n测量变量的总个数、i=1,2,…,n、R为实数集,R 表示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;
T
步骤(6):先根据公式Λb=TbTb/(n‑1)计算协方差矩阵Λ1,Λ1,…,ΛB后,再根据公式T ‑1Qb=diag(TbΛb Tb)计算监测统计指标向量Q1,Q2,…,QB,其中diag()表示将矩阵对角线的元素转变成向量的操作,上标号T表示矩阵或向量的转置;
步骤(7):分别对监测统计指标向量Q1,Q2,…,QB中的元素按大小进行降序排列,并将第n/100个最大值分别做为各变量子块的监测统计指标上限CL1,CL2,…,CLB;
在线监测阶段执行如下所示步骤:
m×1
步骤(8):收集新采样时刻的样本数据x∈R ,并根据公式 对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为种群个数;
N×n
步骤(4.2):根据公式 求解矩阵F∈R ,并计算矩阵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行向量保留;
N×m
步骤(4.7):重复实施步骤(4.6)N次,从而将保留的N个行向量组成矩阵V∈R ,并初始化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)上随机生成的2
一个数,函数h(u)=rk(1‑u/U) ;若否,则保持矩阵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。