CN108375965A - 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 - Google Patents
一种基于多变量块交叉相关性剔除的非高斯过程监测方法 Download PDFInfo
- Publication number
- CN108375965A CN108375965A CN201810280156.8A CN201810280156A CN108375965A CN 108375965 A CN108375965 A CN 108375965A CN 201810280156 A CN201810280156 A CN 201810280156A CN 108375965 A CN108375965 A CN 108375965A
- Authority
- CN
- China
- Prior art keywords
- matrix
- vector
- formula
- error
- new
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS], computer integrated manufacturing [CIM]
- G05B19/41885—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS], computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32339—Object oriented modeling, design, analysis, implementation, simulation language
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Abstract
本发明公开一种基于多变量块交叉相关性剔除的非高斯过程监测方法,旨在将多变量块之间的交叉相关性考虑进分布式过程建模与过程监测中,从而实施更加可靠而有效的分布式非高斯过程监测。本发明方法首先初步根据各生产单元测量变量的归属,将所有测量变量划分成多个变量子块;其次,利用回归模型将每个变量子块中与其他变量子块之间的交叉相关信息剔除出去;最后,利用剔除交叉相关性后的误差实施基于独立成分分析算法的建模与非高斯过程监测。相比于传统方法,本发明方法的主要创新在于利用回归模型将不同多变量子块之间的交叉相关性考虑进来,以交叉相关性剔除后的误差做为新的监测对象。
Description
技术领域
本发明涉及一种数据驱动的过程监测方法,尤其涉及一种基于多变量块交叉相关性剔除的非高斯过程监测方法。
背景技术
实时监测生产过程运行状态对于保证生产安全与产品质量的稳定性具有重要性意义,生产实际中通常采取的过程监测技术手段就是我们所熟知的故障监测方法,旨在及时地甄别出系统出现的异常状态。近年来,随着工业大规模化与“大数据”化建设的推进,生产过程中可以很容易地采集到大量的样本数据,但是无法为其建立精确的机理模型。在这种应用背景下,传统基于机理模型的故障监测方法就遇到了发展的瓶颈,而数据驱动的过程监测方法逐步得到重视与广泛的研究。在这一研究领域里,多变量统计过程监测得到了最多的研究与关注。这其中,又当以主成分分析(Principal Component Analysis,PCA)与独立成分分析(Independent Component Analysis,ICA)算法为最主流的两种建模与故障监测方法。两者的共同之处在于通过挖掘训练数据的潜在特征,将过程监测任务转换成一个单分类的模式识别问题。
然而,通过现有文献资料与专利技术文件中可以发现,基于ICA的过程监测方法能取得的效果通常优越于基于PCA的过程监测方法。这主要是因为ICA能够挖掘出训练数据的非高斯特征成分,更能揭示事物的本质,同时更能适应于解决非高斯过程的故障监测问题。考虑到现代工业过程的大规模性与复杂特性,生产过程的采样数据一般不符合高斯分布假设,基于ICA的非高斯过程监测方法更值得深入研究。此外,由于现代工业过程对象通常有诸多生产子单元组成,直接对对象整体采用一个故障监测模型取得的监测性能往往差强人意。近年来,分布式的建模与故障监测方法得到了广泛的研究与应用。一般来讲,分布式的过程监测方法能够在建模时简化对过程对象分析的难度,而且多模型的思路通常泛化性能优越于单个模型。基于分布式ICA的非高斯过程监测方法实施的技术难点主要在于如何将所有的测量变量划分成多个变量子块,一种最简单最直接的划分多变量块的方法就是按照测量变量对应的生产单元归属实施变量块划分。
由于现代工业过程各生产子单元之间不是单纯的串联或并联链接关系,相互之间还存在复杂交叉关系,比如产品的回流进入其他生产单元,以及利用不同生产单元间的反馈信号实施控制等等。这种交叉关系体现在采样数据上就是不同变量子块间的数据集存在交叉相关性。因此,直接按照生产单元分成的变量子块的做法,再直接建立多个分布式的过程监测模型不能将不同生产单元之间的交叉相关性考虑进来。因此,传统的分布式的非高斯过程监测方法还有待进一步的改进。可想而知,若是按照生产单元实施分布式故障监测时,能将不同生产单元之间的交叉相关性考虑进来,将会提升相应过程监测方法的性能。可是,若是直接根据生产过程的机理知识分析这种交叉关联关系,需要深入细致的分析过程对象的组成结构情况。这在实际过程中是不可取的,也违背了数据驱动过程监测方法从数据角度出发的优势与理念。
发明内容
本发明所要解决的主要技术问题是:如何将多变量块之间的交叉相关性考虑进分布式过程建模与过程监测中,从而实施更加可靠而有效的分布式非高斯过程监测。具体来讲,本发明方法首先初步根据各生产单元测量变量的归属,将所有测量变量划分成多个变量子块;其次,利用回归模型将每个变量子块中与其他变量子块之间的交叉相关信息剔除出去;最后,利用剔除交叉相关性后的误差实施基于ICA算法的建模与非高斯过程监测。
本发明解决上述技术问题所采用的技术方案为:一种基于多变量块交叉相关性剔除的非高斯过程监测方法,包括以下步骤:
(1)采集生产过程对象中所有生产单元正常运行状态下的样本数据,组成训练数据集X∈Rn×m,并对其进行标准化处理,得到均值为0、标准差为1的新数据矩阵其中,n为训练样本总数,m为过程对象所有测量变量的个数,R为实数集,Rn×m表示n×m维的实数矩阵。
(2)根据各测量变量的生产单元归属,将m个测量变量分成C个不同的变量子块,其中C表示过程对象中生产单元的个数。
值得注意的是,根据过程对象生产单元划分变量子块时,本发明方法不要求对各变量子块做精确的划分。在某些测量变量归属的生产单元可有两个及以上时,可以随意划分该测量变量所属的生产单元。
(3)依据C个不同的变量子块将矩阵对应地分成C个不同的子矩阵X1,X2,…,XC,其中c=1,2,…,C,mc为第c个变量子块中变量个数且满足条件∑mc=m。
(4)将第c个子矩阵Xc做为回归模型的输出,同时将Y=[X1,…,Xc-1,Xc+1,…XC]做为回归模型的输入,利用偏最小二乘算法(Partial Least Square,PLS)建立输入-输出之间的回归模型,具体的实施过程如下所示:
①置k=1与Z=Y后,初始化向量uk为矩阵Xc的第一列;
②根据公式wk=Zuk/||Zuk||、sk=Zwk、qk=Xc Tsk/(sk Tsk)分别计算得到系数向量wk、得分向量sk、和系数向量qk,其中||Zuk||表示计算向量Zuk的长度,上标号T表示矩阵或者向量的转置;
③根据公式unew=Xcqk/qk 2计算向量unew;
④判断是否满足条件||uk-unew||<10-6?若否,则置uk=unew后返回步骤②;若是,则执行⑤;
⑤根据公式pk=ZTsk/(sk Tsk)计算得到第k个投影向量pk∈R(m-1)×1,并保留向量pk、向量wk、和系数向量qk;
⑥判断矩阵Yk=skpk T中的最大元素是否大于0.01?若是,根据公式Z=Z-skpk T更新矩阵Zi后执行步骤⑦;若否,则得到投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2,…,wk]、和向量Q=[q1,q2,…,qk],并执行步骤⑧;
⑦判断k<m-mc?若是,则置k=k+1后,返回步骤②;若否,则得到最终的投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2,…,wk]、和向量Q=[q1,q2,…,qk];
⑧根据公式br=Wr(Pr TWr)-1Qr T计算在保留r个投影向量前提下的回归系数向量br,其中r=1,2,…,k,Pr、Wr、和Qr分别为矩阵P、W、和Q中前1至r列向量组成的矩阵;
⑨根据公式计算输出Xc的预测值后,计算保留不同个数投影向量所对应的模型预测误差其中diag[]表示将矩阵对角线的元素单独组成向量的操作,sum{}表示求取向量中各元素之和;
⑩根据MSE1,MSE2,…,MSEk的数值变化情况找出预测误差不再发生显著变化的收敛点位置,并将其下标号对应的具体数值作为偏最小二乘模型需保留的投影变量个数K,其中K≤k;
根据K确定最终的输入Y与输出Xc之间的回归模型,即:
Xc=YWK(PK TWK)-1QK T+Ec=YΘc+Ec (1)
上式中,回归矩阵Θc=WK(PK TWK)-1QK T,为回归模型误差;
(5)重复步骤(4)直至得到C个PLS模型,并将模型误差E1,E2,…,EC合并成一个误差矩阵E=[E1,E2,…,EC]∈Rn×m。
(6)将估计误差矩阵E作为新的训练数据矩阵,对其中的每一列实施标准化处理,得到均值为0,标准差为1的新数据矩阵
(7)对实施白化处理,将其转换成数据矩阵F∈Rn×M,并保留白化转换矩阵G∈Rm ×M,其中M为白化转换向量个数,具体的实施过程如下所示:
①计算的协方差矩阵
②求解Ψ所有M个大于0.001的特征值η1≥η2≥…≥ηM所对应的特征向量β1,β2…,βM;
③根据公式对进行白化处理,得到白化后的数据矩阵F∈Rn×M,其中,白化转换矩阵G=HΦ-1/2、H=[β1,β2…,βM]、以及矩阵Φ为η1,η2,…,ηM组成的对角矩阵。
(8)利用ICA算法为数据矩阵F∈Rn×M建立非高斯模型,保留分离矩阵W与混合矩阵A,具体的实施过程如下所示:
①初始化i=1;
②设置向量di为M×M维单位矩阵中的第i列;
③按照公式di←E{Fg(di TF)}-E{h(di TF)}di更新向量di,其中函数g(u)=tanh(u)、函数h(u)=[sech(u)]2、u=di TF表示函数自变量;
④对更新后的向量di依次按照下式进行正交标准化处理:
di←di/||di|| (3)
⑤重复步骤③~④直至向量di收敛,并保存向量di;
⑥判断i<M?若是,置i=i+1后,重复步骤②~⑥;若否,则执行步骤⑦;
⑦将得到的所有M个向量d1,d2,…,dM组成矩阵Ω=[d1,d2,…,dM]∈RM×M,并按照如下所示公式计算初始分离矩阵W0∈Rm×M与初始混合矩阵A0∈Rm×M:
W0=HΦ-1/2Ω
A0=HΦ1/2Ω (4)
⑧计算A0中每一列向量的长度,分别记为L1,L2,…,LM,并将L1,L2,…,LM按照数值大小进行降序排列得到l1,l2,…,lM,那么保留的独立成分个数k为满足下列条件的最小值:
⑨将A0中列向量长度最大的μ个列向量组成混合矩阵A∈Rm×μ,同时从W0中取出与A对应的列向量组成新分离矩阵W∈Rm×μ。
(9)根据公式与分别计算独立成分矩阵S与残差矩阵
(10)根据公式I2=diag[SST]与分别计算训练数据对应的监测统计量I2与SPE,并找出I2与SPE中最大值分别记录为与SPElim。
上述步骤(1)~(10)为本发明方法的离线建模过程,如下所示步骤(11)~(17)为在线过程监测过程。
(11)收集新采样时刻的数据样本x∈R1×m,对其实施与步骤(1)中相同的标准化处理得到新数据向量
(12)根据步骤(2)中的C个不同的变量子块,对应地将样本向量对应地分成C个不同的子向量x1,x2,…,xC。
(13)调用步骤(4)中建立的第c个PLS模型,根据如下所示公式计算出多变量块相关性剔除后的误差ec:
ec=xc-yΘc (6)
上式中,y=[x1,…,xc-1,xc+1,,xC]。
(14)重复步骤(13)直至得到C个误差e1,e2,…,eC,并将这些误差合并成一个误差向量e=[e1,e2,…,eC]。
(15)对误差e实施与步骤(6)中相同的标准化处理得到新向量
(16)根据如下所示公式计算监测统计指标I2与SPE:
上式中,上标号T表示矩阵或者向量的转置。
(17)判断是否满足条件:且SPE≤SPElim?若是,则当前监测时刻系统处于正常运行状态,返回步骤(11)继续实施在线故障监测;若否,则当前采样数据来自故障工况。
与传统方法相比,本发明方法的优势在于:
本发明方法的主要创新在于利用回归模型将不同多变量子块之间的交叉相关性考虑进来,以交叉相关性剔除后的误差做为新的监测对象。本发明方法利用的误差不仅能够反映出多变量块之间交叉相关关系是否发生变化,而且还能体现各变量块自身的交叉相关性特征。因此,本发明方法理应具备更优秀的故障监测性能,是一种更为优选的非高斯过程监测方法。
附图说明
图1为本发明方法实施流程图。
图2为TE过程结构流程图。
图3为交叉相关性示意图。
图4为TE过程第20类故障的监测详情对比图。
具体实施方式
下面结合附图与具体的实施案例对本发明方法进行详细的说明。
如图1所示,一种基于多生产单元变量交叉相关解耦策略的故障监测方法。下面结合一个具体的工业过程的例子来说明本发明方法的具体实施过程,以及相对于传统基于ICA的非高斯过程监测方法的优越性。
应用对象是来自于美国田纳西-伊斯曼(TE)化工过程实验,原型是伊斯曼化工生产车间的一个实际工艺流程。如图2所示,TE过程的生产流程比较复杂性,包含5个主要的生产单元:反应器、冷凝器、分离塔、汽提塔、和压缩机。TE过程已作为一个标准实验平台被广泛用于故障监测研究,整个TE过程包括22个测量变量、12个操作变量、和19个成分测量变量。本次选取如表1所示的33个过程变量实施离线建模与在线故障监测,接下来结合该TE过程对本发明具体实施步骤进行详细的阐述。
表1:TE过程监测变量。
序号 | 变量描述 | 序号 | 变量描述 | 序号 | 变量描述 |
1 | 物料A流量 | 12 | 分离器液位 | 23 | D进料阀门位置 |
2 | 物料D流量 | 13 | 分离器压力 | 24 | E进料阀门位置 |
3 | 物料E流量 | 14 | 分离器塔底流量 | 25 | A进料阀门位置 |
4 | 总进料流量 | 15 | 汽提塔等级 | 26 | A和C进料阀门位置 |
5 | 循环流量 | 16 | 汽提塔压力 | 27 | 压缩机循环阀门位置 |
6 | 反应器进料 | 17 | 汽提塔底部流量 | 28 | 排空阀门位置 |
7 | 反应器压力 | 18 | 汽提塔温度 | 29 | 分离器液相阀门位置 |
8 | 反应器等级 | 19 | 汽提塔上部蒸汽 | 30 | 汽提塔液相阀门位置 |
9 | 反应器温度 | 20 | 压缩机功率 | 31 | 汽提塔蒸汽阀门位置 |
10 | 排空速率 | 21 | 反应器冷却水出口温度 | 32 | 反应器冷凝水流量 |
11 | 分离器温度 | 22 | 分离器冷却水出口温度 | 33 | 冷凝器冷却水流量 |
首先,利用TE过程正常工况下的采样数据离线建立故障监测模型,包括以下步骤:
(1):收集生产过程正常运行状态下的数据样本,组成训练数据集X∈R960×33,并对每个变量进行标准化处理,得到均值为0,标准差为1的新数据矩阵
表2:多生产单元变量子块划分结果。
生产单元名称: | 变量子块中的变量编号: |
反应器 | 1,2,3,6,7,8,9,21,23,24,25,32 |
冷凝器 | 33 |
分离塔 | 4,15,16,17,18,19,26,30,31 |
汽提塔 | 10,11,12,13,14,22,28,28 |
压缩机 | 5,20,27 |
(2):根据各测量变量的生产单元归属,将33个测量变量分成5个不同的变量子块,相应的变量子块划分结果如表2所示。
(3):依据5个不同的变量子块将矩阵对应地分成5个不同的子矩阵X1,X2,…,X5。
(4):将第c个子矩阵Xc做为回归模型的输出,同时将Y=[X1,…,Xc-1,Xc+1,…X5]做为回归模型的输入,利用PLS算法建立输入-输出之间的回归模型:Xc=YΘc+Ec。
(5):重复步骤(4)直至得到5个PLS模型,并将模型误差E1,E2,…,E5合并成一个误差矩阵E=[E1,E2,…,E5]。
为验证多变量块相关性剔除特征,特将E5与输入Y=[X1,X2,X3,X4]之间的交叉相关性示意图显示于图3中。从图3中可以发现,压缩机生产单元对应的误差E5与其他生产单元的测量变量间的时序相关性几乎全部为0,即体现出了交叉相关解耦。
(6):将估计误差矩阵E作为新的训练数据矩阵,对其实施标准化处理。
(7):对实施白化处理,将其转换成数据矩阵F∈R960×30,并保留白化转换矩阵G∈R33×30。
(8):利用ICA算法为数据矩阵F∈Rn×M建立非高斯模型,保留分离矩阵W与混合矩阵A。
(9):根据公式与分别计算独立成分矩阵S与残差矩阵
(10):根据公式I2=diag[SST]与分别计算训练数据对应的监测统计量I2与SPE,并找出I2与SPE中最大值分别记录为与SPElim。
利用TE过程第20类故障变量工况的采样数据作为测试数据,实施在线故障监测。值得指出的是,该测试数据集前160个数据为正常工况,故障工况从161个数据采样时刻引入。
(11):收集新采样时刻的数据样本x∈R1×33,对其实施与步骤(1)中相同的标准化处理得到新数据向量后。
(12):根据步骤(2)中的5个不同的变量子块,对应地将样本向量对应地分成5个不同的子向量x1,x2,…,x5;
(13):调用步骤(4)中建立的第c个PLS模型,根据如下所示公式计算出多生产单元变量交叉相关解耦后的误差ec。
(14):重复步骤(10)直至得到5个误差e1,e2,…,e5,并将这些误差合并成一个误差向量e=[e1,e2,…,e5]。
(15):对误差e实施与步骤(6)中相同的标准化处理得到新向量。
(16):计算监测统计指标I2与SPE。
(17):判断是否满足条件:且SPE≤SPElim?若是,则当前监测时刻系统处于正常运行状态,返回步骤(11)继续实施在线故障监测;若否,则当前采样数据来自故障工况。
如图4所示,本发明方法以及传统ICA方法在该测试故障上的监测详情对比图,可以很明显地发现本发明方法的SPE监测统计量的故障漏报率明显低于传统ICA方法。因此,本发明方法是一种更优越的非高斯过程监测方法。
上述实施案例只用来解释说明本发明的具体实施,而不是对本发明进行限制。在本发明的精神和权利要求的保护范围内,对本发明做出的任何修改,都落入本发明的保护范围。
Claims (4)
1.一种基于多变量块交叉相关性剔除的非高斯过程监测方法,其特征在于,包括以下步骤:离线建模阶段的实施过程如下所示:
步骤(1):采集生产过程对象中所有生产单元正常运行状态下的样本数据,组成训练数据集X∈Rn×m,并对其进行标准化处理,得到均值为0、标准差为1的新数据矩阵其中,n为训练样本总数,m为过程对象所有测量变量的个数,R为实数集,Rn×m表示n×m维的实数矩阵;
步骤(2):根据各测量变量的生产单元归属,将m个测量变量分成C个不同的变量子块,其中C表示过程对象中生产单元的个数;
值得注意的是,根据过程对象生产单元划分变量子块时,本发明方法不要求对各变量子块做精确的划分,在某些测量变量归属的生产单元可有两个及以上时,可以随意划分该测量变量所属的生产单元;
步骤(3):依据C个不同的变量子块将矩阵对应地分成C个不同的子矩阵X1,X2,…,XC,其中c=1,2,…,C,mc为第c个变量子块中变量个数且满足条件∑mc=m;
步骤(4):将第c个子矩阵Xc做为回归模型的输出,同时将Y=[X1,…,Xc-1,Xc+1,…XC]做为回归模型的输入,利用偏最小二乘算法(Partial Least Square,PLS)建立输入-输出之间的回归模型:
Xc=YΘc+Ec (1)
上式中,Θc为回归矩阵,为回归模型误差;
步骤(5):重复步骤(4)直至得到C个PLS模型,并将模型误差E1,E2,…,EC合并成一个误差矩阵E=[E1,E2,…,EC]∈Rn×m;
步骤(6):将估计误差矩阵E作为新的训练数据矩阵,对其中的每一列实施标准化处理,得到均值为0,标准差为1的新数据矩阵
步骤(7):对实施白化处理,将其转换成数据矩阵F∈Rn×M,并保留白化转换矩阵G∈Rm ×M,其中M为白化转换向量个数;
步骤(8):利用独立成分分析(Independent Component Analysis,ICA)算法为数据矩阵F∈Rn×M建立非高斯模型,保留分离矩阵W与混合矩阵A;
步骤(9):根据公式与分别计算独立成分矩阵S与残差矩阵
步骤(10):根据公式I2=diag[SST]与分别计算训练数据对应的监测统计量I2与SPE,并找出I2与SPE中最大值分别记录为与SPElim;
在线监测的实施过程如下所示:
步骤(11):收集新采样时刻的数据样本x∈R1×m,对其实施与步骤(1)中相同的标准化处理得到新数据向量
步骤(12):根据步骤(2)中的C个不同的变量子块,对应地将样本向量对应地分成C个不同的子向量x1,x2,…,xC;
步骤(13):调用步骤(4)中建立的第c个PLS模型,根据如下所示公式计算出多变量块相关性剔除后的误差ec:
ec=xc-yΘc (2)
上式中,y=[x1,…,xc-1,xc+1,,xC];
步骤(14):重复步骤(13)直至得到C个误差e1,e2,…,eC,并将这些误差合并成一个误差向量e=[e1,e2,…,eC];
步骤(15):对误差e实施与步骤(6)中相同的标准化处理得到新向量
步骤(16):根据如下所示公式计算监测统计指标I2与SPE:
上式中,上标号T表示矩阵或者向量的转置;
步骤(17):判断是否满足条件:且SPE≤SPElim?若是,则当前监测时刻系统处于正常运行状态,返回步骤(11)继续实施在线故障监测;若否,则当前采样数据来自故障工况。
2.根据权利要求1所述的一种基于多变量块交叉相关性剔除的非高斯过程监测方法,其特征在于,所述步骤(4)中建立输入-输出回归模型的具体实施过程如下所示:
①置k=1与Z=Y后,初始化向量uk为矩阵Xc的第一列;
②根据公式wk=Zuk/||Zuk||、sk=Zwk、qk=Xc Tsk/(sk Tsk)分别计算得到系数向量wk、得分向量sk、和系数向量qk,其中||Zuk||表示计算向量Zuk的长度,上标号T表示矩阵或者向量的转置;
③根据公式unew=Xcqk/qk 2计算向量unew;
④判断是否满足条件||uk-unew||<10-6?若否,则置uk=unew后返回步骤②;若是,则执行⑤;
⑤根据公式pk=ZTsk/(sk Tsk)计算得到第k个投影向量pk∈R(m-1)×1,并保留向量pk、向量wk、和系数向量qk;
⑥判断矩阵Yk=skpk T中的最大元素是否大于0.01?若是,根据公式Z=Z-skpk T更新矩阵Zi后执行步骤⑦;若否,则得到投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2,…,wk]、和向量Q=[q1,q2,…,qk],并执行步骤⑧;
⑦判断k<m-mc?若是,则置k=k+1后,返回步骤②;若否,则得到最终的投影矩阵P=[p1,p2,…,pk]、系数矩阵W=[w1,w2…,wk]、和向量Q=[q1,q2,…,qk];
⑧根据公式br=Wr(Pr TWr)-1Qr T计算在保留r个投影向量前提下的回归系数向量br,其中r=1,2,…,k,Pr、Wr、和Qr分别为矩阵P、W、和Q中前1至r列向量组成的矩阵;
⑨根据公式计算输出Xc的预测值后,计算保留不同个数投影向量所对应的模型预测误差其中diag[ ]表示将矩阵对角线的元素单独组成向量的操作,sum{ }表示求取向量中各元素之和;
⑩根据MSE1,MSE2,…,MSEk的数值变化情况找出预测误差不再发生显著变化的收敛点位置,并将其下标号对应的具体数值作为偏最小二乘模型需保留的投影变量个数K,其中K≤k;
根据K确定最终的输入Y与输出Xc之间的回归模型,即:
Xc=YWK(PK TWK)-1QK T+Ec=YΘc+Ec (1)
上式中,回归矩阵Θc=WK(PK TWK)-1QK T,为回归模型误差。
3.根据权利要求1所述的一种基于多变量块交叉相关性剔除的非高斯过程监测方法,其特征在于,所述步骤(7)中对实施白化处理的具体实施过程如下所示:
①计算的协方差矩阵
②求解Ψ所有M个大于0.001的特征值η1≥η2≥…≥ηM所对应的特征向量β1,β2…,βM;
③根据公式对进行白化处理,得到白化后的数据矩阵F∈Rn×M,其中,白化转换矩阵G=HΦ-1/2、H=[β1,β2…,βM]、以及矩阵Φ为η1,η2,…,ηM组成的对角矩阵。
4.根据权利要求1所述的一种基于多变量块交叉相关性剔除的非高斯过程监测方法,其特征在于,所述步骤(8)中利用ICA算法为数据矩阵F建立非高斯模型的具体实施过程如下所示:
①初始化i=1;
②设置向量di为M×M维单位矩阵中的第i列;
③按照公式di←E{Fg(di TF)}-E{h(di TF)}di更新向量di,其中函数g(u)=tanh(u)、函数h(u)=[sech(u)]2、u=di TF表示函数自变量;
④对更新后的向量di依次按照下式进行正交标准化处理:
di←di/||di|| (3)
⑤重复步骤③~④直至向量di收敛,并保存向量di;
⑥判断i<M?若是,置i=i+1后,重复步骤②~⑥;若否,则执行步骤⑦;
⑦将得到的所有M个向量d1,d2,…,dM组成矩阵Ω=[d1,d2,…,dM]∈RM×M,并按照如下所示公式计算初始分离矩阵W0∈Rm×M与初始混合矩阵A0∈Rm×M:
⑧计算A0中每一列向量的长度,分别记为L1,L2,…,LM,并将L1,L2,…,LM按照数值大小进行降序排列得到l1,l2,…,lM,那么保留的独立成分个数k为满足下列条件的最小值:
⑨将A0中列向量长度最大的μ个列向量组成混合矩阵A∈Rm×μ,同时从W0中取出与A对应的列向量组成新分离矩阵W∈Rm×μ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810280156.8A CN108375965B (zh) | 2018-03-19 | 2018-03-19 | 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810280156.8A CN108375965B (zh) | 2018-03-19 | 2018-03-19 | 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108375965A true CN108375965A (zh) | 2018-08-07 |
CN108375965B CN108375965B (zh) | 2020-06-30 |
Family
ID=63019281
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810280156.8A Active CN108375965B (zh) | 2018-03-19 | 2018-03-19 | 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108375965B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109507972A (zh) * | 2018-12-19 | 2019-03-22 | 中国计量大学 | 基于分层式非高斯监测算法的工业生产过程故障监测方法 |
CN109669412A (zh) * | 2018-12-13 | 2019-04-23 | 宁波大学 | 一种基于新型动态独立成分分析的非高斯过程监测方法 |
CN111694328A (zh) * | 2019-03-12 | 2020-09-22 | 宁波大学 | 一种基于多块独立成分分析算法的分布式过程监测方法 |
CN111695229A (zh) * | 2019-03-12 | 2020-09-22 | 宁波大学 | 一种基于ga-ica的新型分散式非高斯过程监测方法 |
CN111913444A (zh) * | 2019-09-07 | 2020-11-10 | 宁波大学 | 一种基于时序多块建模策略的化工过程监测方法 |
CN113191615A (zh) * | 2021-04-18 | 2021-07-30 | 宁波大学科学技术学院 | 一种基于多块相关成分分析的聚丙烯生产过程异常检测方法 |
CN114415609A (zh) * | 2021-12-22 | 2022-04-29 | 华东理工大学 | 一种基于多子空间划分的动态过程精细化监测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030061185A1 (en) * | 1999-10-14 | 2003-03-27 | Te-Won Lee | System and method of separating signals |
CN101446828A (zh) * | 2008-12-30 | 2009-06-03 | 东北大学 | 一种非线性过程质量预报方法 |
CN103488091A (zh) * | 2013-09-27 | 2014-01-01 | 上海交通大学 | 一种数据驱动的基于动态成分分析的控制过程监控方法 |
CN106092625A (zh) * | 2016-05-30 | 2016-11-09 | 宁波大学 | 基于修正型独立元分析和贝叶斯概率融合的工业过程故障检测方法 |
CN107153409A (zh) * | 2017-06-02 | 2017-09-12 | 宁波大学 | 一种基于缺失变量建模思路的非高斯过程监测方法 |
-
2018
- 2018-03-19 CN CN201810280156.8A patent/CN108375965B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030061185A1 (en) * | 1999-10-14 | 2003-03-27 | Te-Won Lee | System and method of separating signals |
CN101446828A (zh) * | 2008-12-30 | 2009-06-03 | 东北大学 | 一种非线性过程质量预报方法 |
CN103488091A (zh) * | 2013-09-27 | 2014-01-01 | 上海交通大学 | 一种数据驱动的基于动态成分分析的控制过程监控方法 |
CN106092625A (zh) * | 2016-05-30 | 2016-11-09 | 宁波大学 | 基于修正型独立元分析和贝叶斯概率融合的工业过程故障检测方法 |
CN107153409A (zh) * | 2017-06-02 | 2017-09-12 | 宁波大学 | 一种基于缺失变量建模思路的非高斯过程监测方法 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109669412A (zh) * | 2018-12-13 | 2019-04-23 | 宁波大学 | 一种基于新型动态独立成分分析的非高斯过程监测方法 |
CN109669412B (zh) * | 2018-12-13 | 2021-03-26 | 宁波大学 | 一种基于新型动态独立成分分析的非高斯过程监测方法 |
CN109507972A (zh) * | 2018-12-19 | 2019-03-22 | 中国计量大学 | 基于分层式非高斯监测算法的工业生产过程故障监测方法 |
CN111694328A (zh) * | 2019-03-12 | 2020-09-22 | 宁波大学 | 一种基于多块独立成分分析算法的分布式过程监测方法 |
CN111695229A (zh) * | 2019-03-12 | 2020-09-22 | 宁波大学 | 一种基于ga-ica的新型分散式非高斯过程监测方法 |
CN111694328B (zh) * | 2019-03-12 | 2022-03-18 | 宁波大学 | 一种基于多块独立成分分析算法的分布式过程监测方法 |
CN111695229B (zh) * | 2019-03-12 | 2023-10-17 | 宁波大学 | 一种基于ga-ica的新型分散式非高斯过程监测方法 |
CN111913444A (zh) * | 2019-09-07 | 2020-11-10 | 宁波大学 | 一种基于时序多块建模策略的化工过程监测方法 |
CN111913444B (zh) * | 2019-09-07 | 2022-03-18 | 宁波大学 | 一种基于时序多块建模策略的化工过程监测方法 |
CN113191615A (zh) * | 2021-04-18 | 2021-07-30 | 宁波大学科学技术学院 | 一种基于多块相关成分分析的聚丙烯生产过程异常检测方法 |
CN113191615B (zh) * | 2021-04-18 | 2022-05-31 | 宁波大学科学技术学院 | 一种基于多块相关成分分析的聚丙烯生产过程异常检测方法 |
CN114415609A (zh) * | 2021-12-22 | 2022-04-29 | 华东理工大学 | 一种基于多子空间划分的动态过程精细化监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108375965B (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108375965A (zh) | 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 | |
CN101169623B (zh) | 基于核主元分析贡献图的非线性过程故障辨识方法 | |
CN107153409B (zh) | 一种基于缺失变量建模思路的非高斯过程监测方法 | |
CN105955219B (zh) | 基于互信息的分布式动态过程故障检测方法 | |
CN108803520A (zh) | 一种基于变量非线性自相关性剔除的动态过程监测方法 | |
CN108508865B (zh) | 一种基于分散式osc-pls回归模型的故障检测方法 | |
CN106092625B (zh) | 基于修正型独立元分析和贝叶斯概率融合的工业过程故障检测方法 | |
CN108388234A (zh) | 一种基于相关性划分多变量块pca模型的故障监测方法 | |
CN106940808A (zh) | 一种基于改进型主元分析模型的故障检测方法 | |
CN108445867A (zh) | 一种基于分散式icr模型的非高斯过程监测方法 | |
CN105607631B (zh) | 间歇过程弱故障模型控制限建立方法及弱故障监测方法 | |
CN108345284A (zh) | 一种基于两变量块的质量相关故障检测方法 | |
CN113222045B (zh) | 基于加权特征对齐自编码器的半监督故障分类方法 | |
CN107657274A (zh) | 一种基于k‑means的二叉SVM‑tree不平衡数据工业故障分类方法 | |
CN109669415A (zh) | 一种基于结构化典型变量分析的动态过程监测方法 | |
CN108919755A (zh) | 一种基于多块非线性交叉关系模型的分布式故障检测方法 | |
CN106950945A (zh) | 一种基于量纲可变型独立元分析模型的故障检测方法 | |
CN109359662B (zh) | 面向百万千瓦超超临界机组的非平稳分析与因果诊断方法 | |
CN105955214A (zh) | 基于样本时序和近邻相似性信息的间歇过程故障检测方法 | |
CN108960309A (zh) | 一种基于rbf神经网络自相关性剔除的动态过程监测方法 | |
CN108345294B (zh) | 一种基于分散式主元回归模型的故障检测方法 | |
CN108762242A (zh) | 一种基于多块典型相关分析模型的分布式故障检测方法 | |
CN113222046B (zh) | 基于过滤策略的特征对齐自编码器故障分类方法 | |
Wang et al. | Fault monitoring based on adaptive partition non-negative matrix factorization for non-Gaussian processes | |
CN108572639A (zh) | 一种基于主成分自相关性剔除的动态过程监测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20180807 Assignee: Ningbo Science and Technology Innovation Association Assignor: Ningbo University Contract record no.: X2023980033633 Denomination of invention: A Non Gaussian Process Monitoring Method Based on Multivariable Block Cross Correlation Removal Granted publication date: 20200630 License type: Common License Record date: 20230317 |
|
EE01 | Entry into force of recordation of patent licensing contract |