CN101169623A - 基于核主元分析贡献图的非线性过程故障辨识方法 - Google Patents
基于核主元分析贡献图的非线性过程故障辨识方法 Download PDFInfo
- Publication number
- CN101169623A CN101169623A CN200710158468.3A CN200710158468A CN101169623A CN 101169623 A CN101169623 A CN 101169623A CN 200710158468 A CN200710158468 A CN 200710158468A CN 101169623 A CN101169623 A CN 101169623A
- Authority
- CN
- China
- Prior art keywords
- phi
- alpha
- sigma
- data
- fault
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
- G06F18/21355—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis nonlinear criteria, e.g. embedding a manifold in a Euclidean space
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/243—Classification techniques relating to the number of classes
- G06F18/2433—Single-class perspective, e.g. one-against-all classification; Novelty detection; Outlier detection
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
一种基于核主元分析贡献图的非线性故障辨识方法,包括采集数据、提取特征空间的基、提取非线性主元、故障检测和故障辨识五个步骤,本发明提出了一种新的特征空间基的提取方法,通过基的提取来减少了数据映射到特征空间后的线性冗余,以及在训练样本数量很大的时候减少KPCA的计算量。本发明方法采用贡献图进行故障辨识,对过程的每个变量计算其贡献图和控制限,通过贡献图和控制限的关系判断出各变量失控对故障的发生所承担的责任,此辨识方法克服了输入空间和特征空间不能自由转换所带来的故障辨识的困难。
Description
技术领域
本发明属于故障检测与诊断技术领域,提出了一种基于核主元分析贡献图的非线性过程故障辨识方法。
背景技术
近十年里,多变量统计方法在包括故障检测和故障诊断的过程监控中已经得到了迅速的发展。如主元分析(PCA),偏最小二乘(PLS),独立元分析(ICA)和支持向量基(SVM)等已经在工业过程中得到了广泛的应用。但是这些方法主要解决线性变化问题,基于这些方法的过程监控一般假定过程为线性的,而实际的工业过程都存在不同程度的非线性变化。对于非线性变化较强的过程,应用这些方法进行在线监控将产生较高的误报率。为了解决观测数据的非线性问题,人们就提出了一系列的非线性方法。如一种基于主曲线和神经网络的非线性PCA、一种基于自联想的五层神经网络的非线性PCA以及基于输入训练神经网络的非线性PCA方法等。但是大部分现有的非线性PCA方法都是基于神经网络,必须离线训练,要计算主元的话也必须解决非线性最优化问题。
最近,由Schlkopf等人首次在《Neural Computation》提出一种新的PCA的非线性扩展,即核主元分析(KPCA)。核主元分析是首次通过非线性函数把输入空间映射成特征空间,接着在特征空间中提取主元。它通过引入一个核函数避免执行非线性映射和在特征空间中计算内积。同其他非线性方法相比,它主要的优点是不涉及非线性最优化问题,本质上仅需要线性代数,使用起来同PCA一样简单。但是在KPCA的训练过程中核矩阵的大小是样本数量的平方,当样本数变大,特征值和特征向量的计算就变的很费时间。而且对于现有的KPCA,故障变量的分离和故障的识别仍然是个问题。在PCA方法中,故障变量和监控量之间存在一定的线性关系,可以比较容易的计算变量贡献用于绘制贡献图。如Qin在《JOURNAL OFCHEMOMETRICS》中所提出的在输入空间基于重构的贡献图方法。而对于KPCA,由于输入空间到特征空间的非线性映射函数未知,因此类似PCA中的贡献图无法直接用于KPCA的故障变量辨识。
发明内容
针对上述现有技术存在的不足,本发明提供了一种基于改进的核主元分析的非线性过程故障辨识方法。它在原始KPCA的基础上进行了以下两方面的改进:
1)提出了特征空间基的概念,通过提取基来减少了数据映射到特征空间后的线性冗余,以及在训练样本数量很大的时候减少KPCA的计算量;
2)提出了基于重构的故障辨识策略,首次通过基于重构的贡献图来辨识过程中各变量失控对故障产生所承担的责任。
本发明的故障诊断方法包括如下步骤:
步骤一、采集数据
采集过程相关变量的数据,并对其进行中心化和归一化,即规范化处理。
步骤二、提取特征空间的基
根据特征因子依次判断规范化后的数据与基中数据的线性关系,并去除线性关系较强的数据,其余数据加入基中,最后得到基Ω。为了获得一个能够代表特征空间所有信息的特征空间的基,本发明采用基于几何学思想的方法。假定 是一个特征空间的基,其中Φ(·)表示输入空间映射到特征空间的非线性映射函数。特征空间的任意数据点能够表示如下:
当差值ηi等于0时,则 表明数据点Φ(xi)可由基ΦB线性表示。
通过对方程(2)的最小化求解得到
学习的开始,基中只有一个数据点,即 第i个特征空间的基定义为 i=1,2,…,t1,且满足ni<t1,t1为采集数据的个数。根据方程(3)得到
根据方程(1)和(2)定义特征因子
其中 当Si+1=0时,Φ(xi+1)可由基Φi中数据点线性表示。
本发明引入核技巧k(x,y)=(Φ(x),Φ(y)),从而特征空间的特征提取转化为在输入空间进行计算,方程(4)和方程(5)分别转化为
因此,通过方程(6)和(7),即可在输入空间确定最终的基。学习的开始基中只有一个数据即 当 其中γ0是满足 且预先确定的较小常数,新数据xi+1将被引入建立新基,得到新基Φi+1={Φi,xi+1};否则,这个数据点将被拒绝,没有基扩展,即Φi+1=Φi。按此方法依次对采集的数据进行特征提取,去除与基中数据线性相关性较强的数据,最后得到基Ω,它保留了N个数据。
步骤三、提取非线性主元
KPCA把标准PCA扩展到了非线性数据分布。KPCA的基本思想是通过非线性映射把输入空间的观测数据x映射到特征空间F,接着在特征空间F执行线性PCA。在执行PCA之前,将Ω中的观测数据xk∈Rm,k=1,2,…,N,映射到一个高维特征空间F
xk→Φ(xk) (8)
假定 则映射后数据的协方差矩阵可表示为
对协方差进行特征分解
λv=CFv (10)
其中特征值λ≥0,v∈F≠{0}是λ对应的特征向量。这里,CFv能够表示如下:
其中<x,y>表示x和y之间的内积。因此λv=CFv等价于
λ<Φ(xk),v>=<Φ(xk),CFv>,k=1,2,…,N (12)
且存在系数αi(i=1,2,…,N)使得
结合方程(12)和(13),得到
其中所有k=1,2,…,N。通常,尽管映射函数Φ(·)存在,但计算上并不容易处理。然而,它不需要明确被计算,仅仅需要知道特征空间中两个向量的内积。通过引入核函数k(x,y)=<Φ(x),Φ(y)>避免了执行非线性映射以及在特征空间计算两者的内积。
现在,定义一个N×N维的核矩阵K
[K]ij=Kij=(Φ(xi),Φ(xj))=k(xi,xj) i=1,2,…,N,j=1,2,…,N (15)那么方程(14)的左边能够表示为
由于k=1,2,…,N,方程(16)转化为λKα。方程(14)的右边能够表示为
由于k=1,2,…,N,方程(17)转化为(1/N)K2α。结合方程(16)和(17),得到
λNKα=K2α (18)
其中系数α=[α1,…,αN]T。方程(18)的求解转化为下式特征向量的求解
Nλα=Kα (19)
在应用KPCA之前,应先执行高维空间中的中心化。核矩阵K中心化如下
因此方程(19)转化为
现在,在特征空间执行PCA等价于求解方程(21)。求解后得到特征值λ1≥λ2≥…≥λN对应的特征向量α1,α2,…αN。通过仅保留前p个特征向量来降低问题的维数。按特征空间F中相应向量归一化的要求对α1,α2,…αp进行归一化,即,
<vk,vk>=1,k=1,2,…,p (22)
利用 方程(22)改为
式中αi k,αj k——分别表示第k个特征向量αk的第i行和第j行的值。通过把Φ(x)投影到特征空间F中的特征向量vk上,其中k=1,2,…,p,任意数据x的非线性主元t计算如下
式中k-k=1,2,…,p;
x-要进行非线性主元提取的数据;
xi-训练数据中的第i个数据,i=1,2,…,N;
其中 k=[k(x,x1),…,k(x,xN)]。
步骤四、故障检测
利用Hotelling的T2统计和平方预测误差(SPE)统计进行故障检测,通过控制图即可判断出是否发生故障,当T2统计和SPE统计超出各自的控制限时,认为有故障发生,反之,整个过程正常。
Hotelling的T2统计量和平方预测误差(SPE)统计量的计算如下
T2=tTΛ-1t=[t1,t2,…,tp]Λ-1[t1,t2…,tp]T (25)
其中tk(k=1,2,…,p)由方程(24)决定,Λ-1是训练数据前p个最大正特征值的对角阵的逆。
T2的100α%控制限是利用F分布获得的
其中F(p,n-1,α)是自由度为p和n-p的F分布的上限100α%的临界点。
SPE的100(1-α)%控制限是
其中
且cα是上限100(1-α)%相应的正常偏差。m表示特征空间有效的维数。
步骤五、故障辨识
本发明在特征空间提出了一种KPCA技术的故障识别方法,即贡献图方法。一旦检测到故障发生,采用贡献图进行故障辨识,对过程的每个变量计算其贡献图和控制限,通过贡献图和控制限的关系即可判断出各变量失控对故障的发生所承担的责任。
按照每个已知故障的方向重构故障。首次提出基于重构的非线性贡献图。当Ξi被假定为输入空间的故障方向,Φ(Ξi)是特征空间的故障方向,他们并不一定是真实的故障,可能是也可能不是,因而,按照Φ(Ξi)方向定义重构的测量向量如下
其中
式中 为训练数据协方差矩阵的特征向量,αj为训练数据的中心化核矩阵的特征向量。
因此
且重构的残差向量变为
和
因此SPE指标就被引入了控制限δ2,得到
方程(42)取代了方程(36),得到
令 因为 因此,
其中
考虑到自由故障测量变量服从正态分布且基于重构的分布在二次型中定义,则基于重构的SPE贡献图的控制限定义如下
其中
式中S是x的协方差矩阵。
利用矩阵轨迹特性,可以把方程(48)和(49)重新排列为
方程(50)和(51)括号里的项都是标量且能够根据迹的运算求出。由以上得出
因此,δi 2表示为:
通过各变量的贡献图SPE_t和对应的控制限δi 2,就可判断出各变量对故障所承担的责任。当 时,则第i个变量对故障的发生要承担一定的责任。根据差值的大小判断责任的大小。
具体操作步骤如下:
一、离线训练
(1)采集数据,并进行规范化处理;
(2)对规范化后的数据进行特征提取,得到基Ω;
(3)根据方程(15)对基Ω中的数据计算核矩阵K;
(4)根据方程(20)在特征空间执行中心化,得到中心化核矩阵
(5)根据方程(21)求解特征向量,并根据方程(23)对前p个最大正特征值所对应的特征向量进行归一化;
二、在线监测
(6)获得新数据,并对它进行规范化得到xt∈Rm,由[kt]i=[kt(xt,xi)]计算核向量kt∈R1×N,
其中xi是被规范化的标准操作数据。xi∈Rm,i=1,2,…,N;
(7)计算中心化核向量如下:
其中K由离线建模的第(3)步确定,1t=(1/I)[1,…,1]∈R1×N;
(8)计算非线性主元
(9)计算新数据的检测统计量
(10)确定统计量的控制限;
(11)计算SPE贡献图;
(12)确定SPE贡献图的控制限。
以美国ROCKWELL公司的可编程控制器(PLC)实现基础控制,监控程序用RSView32提供的VBA应用软件编制。RSLinx与监控程序之间通过DDE方式进行双向通讯。把故障辩识结果输出到计算机的系统管理画面,同时把故障辩识得到的变量保存到实时数据库中,为操作者或相关技术工人进行监督操作提供参考指导作用。
本发明——基于改进的核主元分析的非线性过程故障辩识方法的优越性表现为:
1)提出一种新的特征空间基的提取方法;
2)基于重构的故障辩识方法克服了输入空间和特征空间不能自由转换所带来的故障辩识的困难。
附图说明
图1田纳西-伊斯曼过程示意图;
图中:1,2,3,4,5,6,7,8,9,10,11,12,13分别表示为:流1,流2,流3,流4,流5,流6,流7,流8,流9,流10,流11,流12,流13;FC:流量控制;FI:流量指示器;PI:压力指示器;PHL:压力控制;SC:同步回旋加速器;TC:温度控制;TI:温度指示器;LI:液位指示器;LC:液位控制。XC:成分控制;XA,XB,XD,XE分别是成分A分析,成分B分析,成分D分析,成分E分析。
图2田纳西-伊斯曼过程中故障10基于原始KPCA的监测结果;
图3田纳西-伊斯曼过程中故障10基于改进的KPCA的监测结果;
图4田纳西-伊斯曼过程中故障10基于735次采样的改进的KPCA的贡献图;
图5田纳西-伊斯曼过程中故障10基于785次采样的改进的KPCA的贡献图;
图6本发明方法实现的流程图。
具体实例
田纳西-伊斯曼过程
本发明方法应用到了田纳西-伊斯曼过程仿真数据中,并同原始KPCA的检测结果进行对比。田纳西-伊斯曼过程是一个复杂非线性的过程,它是由Eastman化学品公司所创建,其目的是为评价过程控制和监测方法提供一个真实的工业过程。控制结构如图1所示。过程包括五个主要单元:反应器,冷凝器,压缩机,分离器,汽提塔;而且,它包含了八种成分:A,B,C,D,E,F,G和H。四种反应物A,C,D以及E和惰性成分B一起被加进反应器里,形成产品G和H,还有副产品F。田纳西-伊斯曼过程包括22个连续过程测量,12个控制变量,和19个成分测量。如表1所示。除了反应器的搅拌器的搅动速度(因为没有对它进行控制),共52个观测变量用于本研究的监测。过程包括21个预设定的故障,如表2所示。本实例针对故障10进行分析说明。在故障10的情况中,流4的C进料温度是随机变化的。当故障发生时,汽提塔的温度也发生变化,因而导致了汽提塔压力的变化。为了弥补汽提塔温度和压力的变化,汽提塔蒸汽阀通过一个控制回路来控制,从而汽提塔蒸汽流率也改变了。
表1.田纳西-伊斯曼过程中的监测变量
No. | 过程测量 | No. | 成分测量 | No. | 控制变量 |
123456789101112 | A进料(流1)D进料(流2)E进料(流3)总进料(流4)再循环流量(流8)反应器进料速度(流6)反应器压力反应器等级反应器温度排放速度(流9)产品分离器温度产品分离器液位 | 232425262728293031323334 | 成分A(流6)成分B(流6)成分C(流6)成分D(流6)成分E(流6)成分F(流6)成分A(流9)成分B(流9)成分C(流9)成分D(流9)成分E(流9)成分F(流9) | 424344454647484950515253 | D进料量(流2)E进料量(流3)A进料量(流1)总进料量(流4)压缩机再循环阀排放阀(流9)分离器灌液流量(流10)汽提器液体产品流量(流11)汽提器水流阀反应器冷却水流量冷凝器冷却水流量搅拌速度 |
13141516171819202122 | 产品分离器压力产品分离器塔底低流量(流10)汽提器等级汽提器压力汽提器塔底流量(流11)汽提器温度汽提器流量压缩机功率反应器冷却水出口温度分离器冷却水出口温度 | 35363738394041 | 成分G(流9)成分H(流9)成分D(流11)成分E(流11)成分F(流11)成分G(流11)成分H(流11) |
表2.田纳西-伊斯曼过程的过程故障描述
No. | 描述 | 类型 | No. | 描述 | 类型 |
1234567891011 | A/C进料比率,B成分不变(流4)B成分,A/C进料比率不变(流4)D的进料温度(流2)反应器冷却水的入口温度冷凝器冷却水的入口温度A进料损失(流1)C存在压力损失-可用性降低(流4)A,B,C进料成分(流4)D的进料温度(流2)C的进料温度(流4)反应器冷却水的入口温度 | 阶跃阶跃阶跃阶跃阶跃阶跃阶跃随机变量随机变量随机变量随机变量 | 12131415161718192021 | 冷凝器冷却水的入口温度反映动态反应器冷却水阀门冷凝器冷却水阀门未知未知未知未知未知流4的阀门固定在稳态位置 | 随机变量慢偏移粘住粘住恒定位置 |
步骤一、采集数据
对于训练数据和实时工况数据采用了三分钟的采样间隔采集数据。故障10的训练数据由480个观测数构成,实时工况数据由960个观测数构成。开始时都是没有故障的,训练数据中,故障都是在第20次采样时引入,实时工况数据中,故障都是在第160次采样引入的。训练数据和实时工况数据中的数据都包含了52个观测变量。随机选取了故障10的训练数据和实时工况数据中控制变量的5组数据分别如表3和表4所示。采集数据后再对其进行规范化处理。
表3.故障10训练数据中控制变量的五组数据
表4.故障10实时工况数据中控制变量的五组数据
步骤二、提取特征空间的基
对规范化后的训练数据进行特征空间基的提取,本例核函数采用3次多项式核函数k(x,y)=(x·y+1)3,γ0=0.05。学习的开始基中只有一个数据点Φ1={x1},引入数据x2,根据方程(6)和(7)计算S2,当 时,数据x2被加进基中,得到新基Φ2={x1,x2},否则,数据x2被拒绝,Φ2=Φ1。接着引入x3,根据方程(6)和(7)计算S3,当 时,数据x3被加入基中,得到新基Φ3={Φ2,x3},否则,数据x3被拒绝,Φ3=Φ2。按此方法依次计算,直到计算完480组数据。最终得到基Ω。
步骤三、提取非线性主元
选取径向基核函数 对基Ω中的数据根据方程(15)进行核计算,得到Gram核矩阵K。接着根据公式(20)对Gram核矩阵K进行中心化处理,得到中心化的核矩阵通过计算出了它的特征值,选取了最大的12个特征值以达到降维的目的。根据方程(21)计算出了选取的12个特征值所对应的特征向量,并根据方程(23)对特征向量进行归一化得到归一化后的特征向量αk,k=1,2,…,12。
参数都确定之后,就可以对规范化后的新数据进行非线性主元提取,根据公式[kt]i=[kt(xt,xi]计算新数据的核向量,其中xt是规范化后的新数据,xi是训练数据,i=1,2,…,316。然后根据方程(55)对其进行了中心化处理,得到中心化核向量最后根据方程(56)计算得到新数据的非线性主元。
步骤四、故障检测
首先通过SPE统计量对田纳西-伊斯曼过程进行检测,根据方程(58)计算SPE统计,接着根据方程(28)确定它的控制限,得到改进的KPCA对田纳西-伊斯曼的监测结果如图3所示,从图中看出SPE统计量在大约170s检测到故障。本发明还利用原始KPCA对故障进行检测,原始KPCA也是在大约170s检测到故障,延时时间是10s。如图2所示。由此可以得出,改进的KPCA能达到原始KPCA同样的效果,但改进的KPCA检测结果基于了更少的样本,减少了计算量。
步骤五、故障辨识
由步骤四的结果得知检测到故障发生,因此采用SPE贡献图确定哪些变量引起了故障的发生。针对每个变量分别计算它们的贡献图和控制限,本发明随机选取了第735次采样的数据对各变量进行计算,首先,对第一个变量A进料(流1)进行计算,看看它对故障所承担的责任,假定故障方向Ξ1=(1,0,0,…,0)T∈R52,根据方程(43)和(54)分别计算得到A进料(流1)的SPE_t和控制限δ1 2;接着对第二个变量D进料(流2)进行计算,看看它对故障所承担的责任,假定故障方向Ξ2=(0,1,0,…,0)T∈R52,根据方程(43)和(54)分别计算得到D进料(流2)的SPE_t和控制限δ2 2;按此方法依次对田纳西-伊斯曼的52个变量进行计算,得到贡献图如图4所示。当各变量的SPE_t>δ2时,那么此变量要对故障的发生承担一定的责任,根据SPE_t-δ2的大小判断责任的轻重。从图中可以看出第1-6,8-13,15-17,19,21,22,24,28,30,32,46-48,52个变量失控。接着又随机选取了第785次采样的数据对个变量进行计算,用上面同样的方法,得到结果如图5所示。从图中知道第1-17,19,21,22,24-26,28,30-33,38,43,45-48,52个变量失控。因此,通过此方法就可以很明显的从图中看出各变量对故障的发生所要承担的责任,就不用盲目的对整个过程的所有变量进行检查。
通过上面的实例,表明了本发明——基于改进的核主元分析的非线性过程故障辨识方法的有效性以及同原始KPCA相比的优越性。它不仅实现了故障的检测,还辨识出各变量失控对故障所承担的责任,减少了设备维修的盲目性。
Claims (5)
1.一种基于核主元分析贡献图的非线性故障辨识方法,其特征在于包括以下步骤:
步骤一、采集数据
采集过程相关变量的数据,并对其进行中心化和归一化,即规范化处理;
步骤二、提取特征空间的基
通过选取一个特征空间的基来提取特征,根据特征因子依次判断规范化后的数据与基中数据的线性关系,并去除线性关系较强的数据,其余数据加入基中,最后得到基Ω;
步骤三、提取非线性主元
通过非线性映射把输入空间的观测数据x映射到特征空间F,接着在特征空间F执行线性主元分析,非线性主元计算公式如下:
式中tk-非线性主元t的第k个元素;
k=1,2,...,p,其中p是非线性主元要提取的个数;
vk-训练数据协方差矩阵的第k特征向量;
αi k-中心化核矩阵的第k个特征向量的第i个元素;
x-要进行非线性主元提取的数据;
xi-训练数据中的第i个数据,i=1,2,…,N;
其中 k=[k(x,x1),...,k(x,xN)];
步骤四、故障检测
利用Hotelling的T2统计和平方预测误差SPE统计进行故障检测,通过控制图即可判断出是否发生故障,当T2统计和SPE统计超出各自的控制限时,认为有故障发生,反之,整个过程正常;
步骤五、故障辨识
当步骤四检测到故障发生,则采用贡献图进行故障辨识,对过程的每个变量计算其贡献图和控制限,通过贡献图和控制限的关系即可判断出各变量失控对故障的发生所承担的责任;当变量的贡献图小于其控制限时,则此变量对故障的发生要承担一定的责任;根据差值的大小判断责任的大小。
2.根据权利要求1所述的一种基于核主元分析贡献图的非线性故障辨识方法,其特征在于步骤二所述特征空间基的求解过程如下:
学习的开始,基中只有一个数据点,即 第i个特征空间的基定义为 i=1,2,...,t1,且满足ni<t1,t1为采集数据的个数,根据方程(3)得到系数向量
根据方程(1)和(2)定义特征因子
其中 当Si+1=0时,Φ(xi+1)可由基Φi中数据点线性表示,
引入核技巧k(x,y)=(Φ(x),Φ(y)),从而特征空间基的提取转化为在输入空间进行计算,方程(4)和方程(5)分别转化为
因此,通过方程(6)和(7),即可在输入空间确定最终的基,学习的开始基中只有一个数据即 当 其中γ0是满足 且预先确定的较小常数,新数据xi+1将被引入建立新基,得到新基Φi+1={Φi,xi+1};否则,这个数据点将被拒绝,没有基扩展,即Φi+1=Φi,按此方法依次对采集的数据进行特征提取,去除与基中数据线性相关性较强的数据,最后得到基Ω,保留了N个数据。
3.根据权利要求1所述的一种基于核主元分析贡献图的非线性故障辨识方法,其特征在于步骤四所述T2统计量和SPE统计量及其各自的控制限计算如下:
T2=tTΛ-1t=[t1,t2,…,tp]Λ-1[t1,t2,…,tp]T (25)
式中t-非线性主元,且tk(k=1,2,…,p)由方程(24)决定;
Λ-1-训练数据前p个最大正特征值的对角阵的逆;
T2统计量的控制限为
式中F(p,n-1,α)-自由度为p和n-p的F分布的上限100α%的临界点;
SPE统计量的控制限为
式中 i=1,2,3,m为特征空间有效的维数;
cα-上限100(1-α)%相应的正常偏差。
4.根据权利要求1所述的一种基于核主元分析贡献图的非线性故障辨识方法,其特征在于步骤五所述贡献图及其控制限的求解过程如下:
按照每个已知故障的方向重构故障,当Ξi被假定为输入空间的故障方向,Φ(Ξi)是特征空间的故障方向,贡献图SPE_t的计算公式为
其中 为训练数据的协方差矩阵的特征向量,αj为中心化核矩阵的特征向量;方程(43)的各项计算如下:
贡献图的控制限δi 2计算公式为:
式中S-x的协方差矩阵。
5.根据权利要求1所述的一种基于核主元分析贡献图的非线性故障辨识方法,其特征在于该方法的具体操作步骤如下:
离线训练:
(1)采集数据,并进行规范化处理;
(2)对规范化后的数据进行特征提取,得到基Ω;
(3)根据方程(15)对基Ω中的数据计算核矩阵K;
(5)根据方程(21)求解特征向量,并根据方程(23)对前p个最大正特征值所对应的特征向量进行归一化;
在线监测:
(6)获得新数据,并对它进行规范化得到xt∈Rm,由[kt]i=[kt(xt,xi)]计算核向量kt∈R1×N,其中xi是被规范化的标准操作数据;xi∈Rm,i=1,2,…,N;
(7)计算中心化核向量如下:
其中K由离线建模的第(3)步确定,1t=(1/I)[1,…,1]∈R1×N;
(8)计算非线性主元
(9)计算新数据的检测统计量
(10)确定统计量的控制限;
(11)计算SPE贡献图;
(12)确定SPE贡献图的控制限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710158468.3A CN101169623B (zh) | 2007-11-22 | 2007-11-22 | 基于核主元分析贡献图的非线性过程故障辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200710158468.3A CN101169623B (zh) | 2007-11-22 | 2007-11-22 | 基于核主元分析贡献图的非线性过程故障辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101169623A true CN101169623A (zh) | 2008-04-30 |
CN101169623B CN101169623B (zh) | 2010-07-07 |
Family
ID=39390272
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200710158468.3A Expired - Fee Related CN101169623B (zh) | 2007-11-22 | 2007-11-22 | 基于核主元分析贡献图的非线性过程故障辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101169623B (zh) |
Cited By (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101776531A (zh) * | 2010-02-10 | 2010-07-14 | 东北大学 | 一种球磨机负荷参数软测量方法 |
CN101470426B (zh) * | 2007-12-27 | 2011-02-16 | 北京北方微电子基地设备工艺研究中心有限责任公司 | 一种故障检测的方法和系统 |
CN102449660A (zh) * | 2009-04-01 | 2012-05-09 | I-切塔纳私人有限公司 | 用于数据检测的系统和方法 |
CN102608303A (zh) * | 2012-03-20 | 2012-07-25 | 天津大学 | 一种橡胶硬度的在线测量方法 |
CN103197663A (zh) * | 2013-03-07 | 2013-07-10 | 北京信息科技大学 | 一种故障预测方法及系统 |
CN103207567A (zh) * | 2013-03-08 | 2013-07-17 | 华北电力大学 | 一种低误报率的改进主元分析过程监测方法及其监测系统 |
CN103336507A (zh) * | 2013-06-24 | 2013-10-02 | 浙江大学 | 基于多模态协同时段自动划分的统计建模与在线监测方法 |
CN103366092A (zh) * | 2013-07-16 | 2013-10-23 | 上海大学 | 基于状态转移的工程风险监测系统及方法 |
CN103576594A (zh) * | 2013-11-11 | 2014-02-12 | 浙江工业大学 | 一种基于张量全局-局部保持投影的间歇过程在线监控方法 |
CN104182623A (zh) * | 2014-08-12 | 2014-12-03 | 南京工程学院 | 一种基于当量变化率计算的热工过程数据检测方法 |
CN104199441A (zh) * | 2014-08-22 | 2014-12-10 | 清华大学 | 基于稀疏贡献图的高炉多工况故障分离方法及系统 |
CN104499001A (zh) * | 2015-01-20 | 2015-04-08 | 重庆科技学院 | 基于特征子空间优化相对矩阵的铝电解槽况诊断方法 |
CN104635724A (zh) * | 2014-12-25 | 2015-05-20 | 重庆科技学院 | 基于动态核独立分量分析的天然气净化过程异常检测方法 |
CN104699077A (zh) * | 2015-02-12 | 2015-06-10 | 浙江大学 | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 |
CN104914854A (zh) * | 2015-05-29 | 2015-09-16 | 东北大学 | 一种基于kpca的工业过程故障诊断方法 |
CN105116323A (zh) * | 2015-08-14 | 2015-12-02 | 江苏科技大学 | 一种基于rbf的电机故障检测方法 |
CN105137354A (zh) * | 2015-08-27 | 2015-12-09 | 江苏科技大学 | 一种基于神经网络电机故障检测方法 |
CN105204496A (zh) * | 2015-09-25 | 2015-12-30 | 清华大学 | 动车组空气制动控制系统传感器故障诊断的方法与系统 |
CN105759217A (zh) * | 2016-02-26 | 2016-07-13 | 江苏快乐电源(涟水)有限公司 | 一种基于可测数据的铅酸蓄电池组在线故障诊断方法 |
CN105911977A (zh) * | 2016-04-22 | 2016-08-31 | 浙江大学 | 嵌套迭代费舍尔判别分析与相对变化结合的故障诊断方法 |
CN107544477A (zh) * | 2017-10-23 | 2018-01-05 | 中国石油大学(华东) | 基于核主元分析的非线性工业过程故障检测方法 |
CN107831743A (zh) * | 2017-10-31 | 2018-03-23 | 东北大学 | 基于可信软标签线性判别分析的工业过程故障诊断方法 |
CN108388221A (zh) * | 2018-03-21 | 2018-08-10 | 江南大学 | 流程工业过程的一种多层模式监控方法 |
CN108664002A (zh) * | 2018-04-27 | 2018-10-16 | 中国石油大学(华东) | 一种面向质量的非线性动态过程监控方法 |
CN109165160A (zh) * | 2018-08-28 | 2019-01-08 | 北京理工大学 | 基于核主成分分析算法的软件缺陷预测模型设计方法 |
CN109491347A (zh) * | 2018-12-18 | 2019-03-19 | 江南大学 | 基于ppls模型的批次运行中操作轨迹的调整方法 |
CN109683594A (zh) * | 2019-01-11 | 2019-04-26 | 河南工学院 | 一种异常变量精确识别与定位方法 |
CN109932904A (zh) * | 2019-03-06 | 2019-06-25 | 东方证券股份有限公司 | 基于特征选择与主元控制的异常状态监测与控制系统 |
CN111382809A (zh) * | 2020-06-01 | 2020-07-07 | 国网江西省电力有限公司电力科学研究院 | 一种基于电机输出功率的隔离开关机械故障诊断方法 |
CN113191219A (zh) * | 2021-04-15 | 2021-07-30 | 华能威宁风力发电有限公司 | 一种风机轴承故障自适应辨识方法 |
CN116577718A (zh) * | 2023-04-13 | 2023-08-11 | 中国电力科学研究院有限公司 | 一种对电容式电压互感器进行在线监测的方法及系统 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112002114A (zh) * | 2020-07-22 | 2020-11-27 | 温州大学 | 一种基于5G-ZigBee通信的机电设备无线数据采集系统和方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5544077A (en) * | 1994-01-19 | 1996-08-06 | International Business Machines Corporation | High availability data processing system and method using finite state machines |
CN1655082A (zh) * | 2005-01-27 | 2005-08-17 | 上海交通大学 | 基于核主元分析的非线性故障诊断的方法 |
CN100480926C (zh) * | 2006-11-23 | 2009-04-22 | 浙江大学 | 一种基于小波分析的工业生产过程故障诊断系统及方法 |
-
2007
- 2007-11-22 CN CN200710158468.3A patent/CN101169623B/zh not_active Expired - Fee Related
Cited By (55)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101470426B (zh) * | 2007-12-27 | 2011-02-16 | 北京北方微电子基地设备工艺研究中心有限责任公司 | 一种故障检测的方法和系统 |
CN102449660B (zh) * | 2009-04-01 | 2015-05-06 | I-切塔纳私人有限公司 | 用于数据检测的系统和方法 |
CN102449660A (zh) * | 2009-04-01 | 2012-05-09 | I-切塔纳私人有限公司 | 用于数据检测的系统和方法 |
US8744124B2 (en) | 2009-04-01 | 2014-06-03 | Curtin University Of Technology | Systems and methods for detecting anomalies from data |
CN101776531A (zh) * | 2010-02-10 | 2010-07-14 | 东北大学 | 一种球磨机负荷参数软测量方法 |
CN101776531B (zh) * | 2010-02-10 | 2013-07-10 | 东北大学 | 一种球磨机负荷参数软测量方法 |
CN102608303A (zh) * | 2012-03-20 | 2012-07-25 | 天津大学 | 一种橡胶硬度的在线测量方法 |
CN102608303B (zh) * | 2012-03-20 | 2014-10-22 | 天津大学 | 一种橡胶硬度的在线测量方法 |
CN103197663A (zh) * | 2013-03-07 | 2013-07-10 | 北京信息科技大学 | 一种故障预测方法及系统 |
CN103197663B (zh) * | 2013-03-07 | 2015-05-20 | 北京信息科技大学 | 一种故障预测方法及系统 |
CN103207567A (zh) * | 2013-03-08 | 2013-07-17 | 华北电力大学 | 一种低误报率的改进主元分析过程监测方法及其监测系统 |
CN103207567B (zh) * | 2013-03-08 | 2015-06-24 | 华北电力大学 | 一种低误报率的改进主元分析过程监测方法及其监测系统 |
CN103336507B (zh) * | 2013-06-24 | 2015-08-19 | 浙江大学 | 基于多模态协同时段自动划分的统计建模与在线监测方法 |
CN103336507A (zh) * | 2013-06-24 | 2013-10-02 | 浙江大学 | 基于多模态协同时段自动划分的统计建模与在线监测方法 |
CN103366092A (zh) * | 2013-07-16 | 2013-10-23 | 上海大学 | 基于状态转移的工程风险监测系统及方法 |
CN103366092B (zh) * | 2013-07-16 | 2016-12-28 | 上海大学 | 基于状态转移的工程风险监测系统及方法 |
CN103576594B (zh) * | 2013-11-11 | 2015-12-02 | 浙江工业大学 | 一种基于张量全局-局部保持投影的间歇过程在线监控方法 |
CN103576594A (zh) * | 2013-11-11 | 2014-02-12 | 浙江工业大学 | 一种基于张量全局-局部保持投影的间歇过程在线监控方法 |
CN104182623A (zh) * | 2014-08-12 | 2014-12-03 | 南京工程学院 | 一种基于当量变化率计算的热工过程数据检测方法 |
CN104199441A (zh) * | 2014-08-22 | 2014-12-10 | 清华大学 | 基于稀疏贡献图的高炉多工况故障分离方法及系统 |
CN104199441B (zh) * | 2014-08-22 | 2017-03-01 | 清华大学 | 基于稀疏贡献图的高炉多工况故障分离方法及系统 |
CN104635724A (zh) * | 2014-12-25 | 2015-05-20 | 重庆科技学院 | 基于动态核独立分量分析的天然气净化过程异常检测方法 |
CN104635724B (zh) * | 2014-12-25 | 2017-02-22 | 重庆科技学院 | 基于动态核独立分量分析的天然气净化过程异常检测方法 |
CN104499001A (zh) * | 2015-01-20 | 2015-04-08 | 重庆科技学院 | 基于特征子空间优化相对矩阵的铝电解槽况诊断方法 |
CN104499001B (zh) * | 2015-01-20 | 2017-03-29 | 重庆科技学院 | 基于特征子空间优化相对矩阵的铝电解槽况诊断方法 |
CN104699077A (zh) * | 2015-02-12 | 2015-06-10 | 浙江大学 | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 |
CN104914854A (zh) * | 2015-05-29 | 2015-09-16 | 东北大学 | 一种基于kpca的工业过程故障诊断方法 |
CN104914854B (zh) * | 2015-05-29 | 2017-05-10 | 东北大学 | 一种基于kpca的工业过程故障诊断方法 |
CN105116323A (zh) * | 2015-08-14 | 2015-12-02 | 江苏科技大学 | 一种基于rbf的电机故障检测方法 |
CN105116323B (zh) * | 2015-08-14 | 2017-10-17 | 江苏科技大学 | 一种基于rbf的电机故障检测方法 |
CN105137354A (zh) * | 2015-08-27 | 2015-12-09 | 江苏科技大学 | 一种基于神经网络电机故障检测方法 |
CN105137354B (zh) * | 2015-08-27 | 2018-04-06 | 江苏科技大学 | 一种基于神经网络电机故障检测方法 |
CN105204496A (zh) * | 2015-09-25 | 2015-12-30 | 清华大学 | 动车组空气制动控制系统传感器故障诊断的方法与系统 |
CN105204496B (zh) * | 2015-09-25 | 2018-01-12 | 清华大学 | 动车组空气制动控制系统传感器故障诊断的方法与系统 |
CN105759217A (zh) * | 2016-02-26 | 2016-07-13 | 江苏快乐电源(涟水)有限公司 | 一种基于可测数据的铅酸蓄电池组在线故障诊断方法 |
CN105759217B (zh) * | 2016-02-26 | 2020-01-07 | 江苏快乐电源(涟水)有限公司 | 一种基于可测数据的铅酸蓄电池组在线故障诊断方法 |
CN105911977A (zh) * | 2016-04-22 | 2016-08-31 | 浙江大学 | 嵌套迭代费舍尔判别分析与相对变化结合的故障诊断方法 |
CN105911977B (zh) * | 2016-04-22 | 2018-10-26 | 浙江大学 | 嵌套迭代费舍尔判别分析与相对变化结合的故障诊断方法 |
CN107544477B (zh) * | 2017-10-23 | 2019-05-31 | 中国石油大学(华东) | 基于核主元分析的非线性工业过程故障检测方法 |
CN107544477A (zh) * | 2017-10-23 | 2018-01-05 | 中国石油大学(华东) | 基于核主元分析的非线性工业过程故障检测方法 |
CN107831743A (zh) * | 2017-10-31 | 2018-03-23 | 东北大学 | 基于可信软标签线性判别分析的工业过程故障诊断方法 |
CN107831743B (zh) * | 2017-10-31 | 2019-12-10 | 东北大学 | 基于可信软标签线性判别分析的工业过程故障诊断方法 |
CN108388221B (zh) * | 2018-03-21 | 2019-06-21 | 江南大学 | 流程工业过程的一种多层模式监控方法 |
CN108388221A (zh) * | 2018-03-21 | 2018-08-10 | 江南大学 | 流程工业过程的一种多层模式监控方法 |
CN108664002A (zh) * | 2018-04-27 | 2018-10-16 | 中国石油大学(华东) | 一种面向质量的非线性动态过程监控方法 |
CN109165160A (zh) * | 2018-08-28 | 2019-01-08 | 北京理工大学 | 基于核主成分分析算法的软件缺陷预测模型设计方法 |
CN109491347A (zh) * | 2018-12-18 | 2019-03-19 | 江南大学 | 基于ppls模型的批次运行中操作轨迹的调整方法 |
CN109491347B (zh) * | 2018-12-18 | 2020-04-10 | 江南大学 | 应用于青霉素发酵过程中的批次运行中操作轨迹的调整方法 |
CN109683594A (zh) * | 2019-01-11 | 2019-04-26 | 河南工学院 | 一种异常变量精确识别与定位方法 |
CN109932904A (zh) * | 2019-03-06 | 2019-06-25 | 东方证券股份有限公司 | 基于特征选择与主元控制的异常状态监测与控制系统 |
CN111382809A (zh) * | 2020-06-01 | 2020-07-07 | 国网江西省电力有限公司电力科学研究院 | 一种基于电机输出功率的隔离开关机械故障诊断方法 |
CN111382809B (zh) * | 2020-06-01 | 2020-11-03 | 国网江西省电力有限公司电力科学研究院 | 一种基于电机输出功率的隔离开关机械故障诊断方法 |
CN113191219A (zh) * | 2021-04-15 | 2021-07-30 | 华能威宁风力发电有限公司 | 一种风机轴承故障自适应辨识方法 |
CN113191219B (zh) * | 2021-04-15 | 2022-11-01 | 华能威宁风力发电有限公司 | 一种风机轴承故障自适应辨识方法 |
CN116577718A (zh) * | 2023-04-13 | 2023-08-11 | 中国电力科学研究院有限公司 | 一种对电容式电压互感器进行在线监测的方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101169623B (zh) | 2010-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101169623B (zh) | 基于核主元分析贡献图的非线性过程故障辨识方法 | |
CN100565403C (zh) | 一种非线性过程故障诊断方法 | |
Zhang et al. | A comparison and evaluation of key performance indicator-based multivariate statistics process monitoring approaches | |
CN101446831B (zh) | 一种分散的过程监测方法 | |
CN108062565B (zh) | 基于化工te过程的双主元-动态核主元分析故障诊断方法 | |
CN104699077B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障变量隔离方法 | |
CN105955219A (zh) | 基于互信息的分布式动态过程故障检测方法 | |
CN101403923A (zh) | 基于非高斯成分提取和支持向量描述的过程监控方法 | |
CN109143995B (zh) | 质量相关慢特征分解的闭环系统精细运行状态监测方法 | |
CN104536439B (zh) | 一种基于嵌套迭代费舍尔判别分析的故障诊断方法 | |
CN109507972B (zh) | 基于分层式非高斯监测算法的工业生产过程故障监测方法 | |
CN101308385A (zh) | 基于二维动态核主元分析的非线性过程故障检测方法 | |
CN104714537A (zh) | 一种基于联合相对变化分析和自回归模型的故障预测方法 | |
CN107153409A (zh) | 一种基于缺失变量建模思路的非高斯过程监测方法 | |
CN108469805A (zh) | 一种基于动态性最优选择的分散式动态过程监测方法 | |
CN108375965A (zh) | 一种基于多变量块交叉相关性剔除的非高斯过程监测方法 | |
CN110009126B (zh) | 基于pls模型与pca贡献度融合的在线报警分析方法 | |
CN108830006B (zh) | 基于线性评价因子的线性-非线性工业过程故障检测方法 | |
CN108181893B (zh) | 一种基于pca-kdr的故障检测方法 | |
CN114757269A (zh) | 一种基于局部子空间-邻域保持嵌入的复杂过程精细化故障检测方法 | |
Ge | Improved two-level monitoring system for plant-wide processes | |
CN111639304A (zh) | 基于Xgboost回归模型的CSTR故障定位方法 | |
CN114611067A (zh) | 一种基于典型变量相异性分析的化工过程缓变故障检测方法 | |
CN112904810B (zh) | 基于有效特征选择的流程工业非线性过程监测方法 | |
Fei et al. | Online process monitoring for complex systems with dynamic weighted principal component analysis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100707 Termination date: 20101122 |