CN110579967A - 基于同时降维和字典学习的过程监控方法 - Google Patents
基于同时降维和字典学习的过程监控方法 Download PDFInfo
- Publication number
- CN110579967A CN110579967A CN201910902691.7A CN201910902691A CN110579967A CN 110579967 A CN110579967 A CN 110579967A CN 201910902691 A CN201910902691 A CN 201910902691A CN 110579967 A CN110579967 A CN 110579967A
- Authority
- CN
- China
- Prior art keywords
- dictionary
- spe
- data
- fault
- dictionary learning
- 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
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于同时降维和字典学习的过程监控方法,包括离线字典学习和在线故障监测两个步骤,在离线字典学习阶段中提出一种同时降维和字典学习方法进行字典学习,在线故障监测阶段包括故障检测、模式识别和故障隔离三个功能,对测试数据进行故障检测,若为故障数据则进行故障诊断以对故障发生的位置进行判断,若非故障数据,进行模式识别。本发明在离线学习阶段通过同时降维和字典学习方法学习到的投影和字典能同时解决数据高维和多模态特点的问题,在线故障监测阶段通过投影矩阵对测试数据构造SPE统计量,能减小在线监测过程的计算复杂度,提高在线监测的实时性。本发明能够学习并保留原始数据更多的空间信息,提升字典的表示和辨识能力。
Description
技术领域
本发明涉及过程监控领域,更具体的说,涉及一种基于同时降维和字典学习的过程监控方法。
背景技术
随着工业过程设备的自动化程度和集成化程度的增高,其结构也越来越复杂,随之影响设备稳定运行的因素也在不断增加,这使得工业过程发生故障的可能性越来越大,且故障的复杂度也越来越高,故工业过程监控引起了学术界和工业界的极大的兴趣。由于传感器成本的降低及其技术的快速发展和先进计算机技术的应用,使得现代工业系统应用了大量的传感器去获得更加丰富的过程信息,获得更准确的数据,而利用大量传感器获得过程数据往往维度很高,虽然这些高维数据包含了大量的信息,但是由于各个变量之间的存在很强的相关性,使得这些高维数据包含很多冗余信息。同时工业过程设备往往工作在多个正常操作模态下,即工况条件受到设定值、操作条件、季节变化和工业设备的老化等因素的影响,这增加了数据处理的复杂性。数据的高维度和多模态特性是目前工业过程监控的两大问题。
现有技术通过降维和字典学习两个独立的过程来解决上述问题,即首先对传感器的数据进行降维,然后用降维后的数据进行字典学习。然而这种预训练的降维过程会使得字典学习的训练数据结构是固定的,而这种固定的低维数据结构会导致字典学习过程无法有效提取原始数据最好的特征,降低了字典的表示和辨识能力,并且这种预训练的降维过程也无法保证投影矩阵能提取出原始数据的空间结构。
发明内容
有鉴于此,本发明提出了一种基于同时降维和字典学习的过程监控方法,能够学习并保留原始数据更多的空间信息,提升字典的表示和辨识能力,解决了现有技术中字典学习过程无法有效提取原始数据最好的特征,字典的表示和辨识能力较低,降维过程无法保证投影矩阵能提取出原始数据的空间结构的技术问题。
本发明提供了一种基于同时降维和字典学习的过程监控方法,包括:
利用历史数据进行离线字典学习,得到全局结构化字典、降维投影矩阵和控制阈值;
利用全局结构化字典、降维投影矩阵和控制阈值进行在线故障监测;
所述离线字典学习包括同时降维和字典学习和计算控制阈值,其中,同时降维和字典学习包括以下步骤:
1)固定P和D,利用正交匹配跟踪方法计算以下公式的解Xi,从而更新稀疏编码矩阵X:
其中||·||F表示矩阵的Frobenius范数,P∈Rp×m表示学习到的降维投影矩阵,其中p(p<m)表示降维后样本的维数,Y=[Y1,...,Yi,...,Yc]表征样本,且i=1,...,c,c为模态的类数,表示第i类模态的样本,且D=[D1,...,Di,...,Dc]∈Rp×k是被学习到的全局结构化字典,表示与第i类相关的特定子字典,且X=[X1,...,Xi,...,Xc]∈Rk×n表示样本Y的稀疏编码矩阵,表示样本Yi在全局结构化字典D下的稀疏编码矩阵,且j=1,...,c,表示样本Yi在子字典Di下的稀疏编码矩阵,是一个选择向量,且选择向量Ei能够通过Di=DEi从全局字典D选出Di,同样,通过求出与第i类局部字典相关第i类稀疏编码
2)固定P和X,利用以下公式计算全局结构化字典D的原子dl,从而更新全局结构化字典D:
其中E/i=[E1,...,Ei-1,Ei+1,...,Ec], l=1,...,ki,
3)固定X和D,利用以下公式计算降维投影矩阵P,从而更新降维投影矩阵P:
Pt=Pt-1+γ(U(1:d,:)-Pt-1)
其中,Pt=U(1:d,:),t为迭代次数,U(1:d,:)为奇异值矩阵∑对应的最小d个向量,所述U和∑利用求得,其中γ为小常数控制P的平稳,
可选的,所述控制阈值为SPE统计量的控制线,按照以下公式计算控制线SPEc:
其中,α为置信水平,{SPE1 2,..,SPEn 2}为{x1,...,xn}的SPE统计量,γ是带宽系数,((·)表示为核方程。
可选的,所述在线故障监测包括以下步骤:
1)利用离线字典学习过程得到的降维投影矩阵P按照以下公式计算数据y的SPE统计量SPEnew,
SPEnew=||(I-PTP)xnew||
其中,Xnew为新的数据样本;
2)将SPEnew和控制线SPEc进行比较,若SPEnew超过SPEc,则该数据为故障数据,则进行故障诊断;否则,该数据为正常数据,则进行模式识别。
可选的,在所述故障诊断中发生故障的维度K利用下式进行诊断:
SPEnew>SPEc
其中,ξi∈Rm是一个方向向量,它表示单位矩阵的第i列,yi=y-ξifi为y第i维重构样本,其中为第i维的故障幅值,表示的伪逆,P为离线字典学习过程得到的降维投影矩阵。
可选的,在所述模式识别中利用以下公式计算数据y的模态类别I:
SPEnew<SPEc
其中,x=[x1;...;xi;...;xc]为稀疏编码矩阵,利用正交匹配跟踪方法求得公式的解,xi表示与子字典Di相关的系数子矢量,P为离线字典学习过程得到的降维投影矩阵,表示与第i类相关的特定子字典。
与现有技术相比,本发明之技术方案具有以下优点:本发明的基于同时降维和字典学习的过程监控方法包括离线字典学习和在线故障监测两个步骤,在离线字典学习阶段中提出一种同时降维和字典学习方法进行字典学习,在线故障监测阶段包括故障检测、模式识别和故障隔离三个功能,首先对测试数据进行故障检测,若为故障数据则进行故障诊断,否则对其进行模式识别。本发明能够学习并保留原始数据更多的空间信息,提升字典的表示和辨识能力,解决了现有技术中字典学习过程无法有效提取原始数据最好的特征,字典的表示和辨识能力较低,降维过程无法保证投影矩阵能提取出原始数据的空间结构的技术问题。本发明在离线学习阶段通过同时降维和字典学习方法学习到的投影和字典能同时解决数据高维和多模态特点的问题,在线故障监测阶段通过投影矩阵对测试数据构造SPE统计量,能减小在线监测过程的计算复杂度,提高在线监测的实时性。
附图说明
图1为本发明基于同时降维和字典学习的过程监控方法的流程图。
具体实施方式
以下结合附图对本发明的优选实施例进行详细描述,但本发明并不仅仅限于这些实施例。本发明涵盖任何在本发明的精神和范围上做的替代、修改、等效方法以及方案。
为了使公众对本发明有彻底的了解,在以下本发明优选实施例中详细说明了具体的细节,而对本领域技术人员来说没有这些细节的描述也可以完全理解本发明。
在下列段落中参照附图以举例方式更具体地描述本发明。需说明的是,附图均采用较为简化的形式且均使用非精准的比例,仅用以方便、明晰地辅助说明本发明实施例的目的。
图1示意了本发明基于同时降维和字典学习的过程监控方法的流程图。所述基于同时降维和字典学习的过程监控方法,包括以下步骤:
本发明提供了一种基于同时降维和字典学习的过程监控方法,包括:
利用历史数据进行离线字典学习,得到全局结构化字典、降维投影矩阵和控制阈值;
利用全局结构化字典、降维投影矩阵和控制阈值进行在线故障监测;
所述离线字典学习包括同时降维和字典学习和计算控制阈值,其中,同时降维和字典学习包括以下步骤:
1)固定P和D,利用正交匹配跟踪方法计算以下公式的解Xi,从而更新稀疏编码矩阵X:
其中||·||F表示矩阵的Frobenius范数,P∈Rp×m表示学习到的降维投影矩阵,其中p(p<m)表示降维后样本的维数,Y=[Y1,...,Yi,...,Yc]表征样本,且i=1,...,c,c为模态的类数,表示第i类模态的样本,且D=[D1,...,Di,...,Dc]∈Rp×k是被学习到的全局结构化字典,表示与第i类相关的特定子字典,且X=[X1,...,Xi,...,Xc]∈Rk×n表示样本Y的稀疏编码矩阵,表示样本Yi在全局结构化字典D下的稀疏编码矩阵,且j=1,...,c,表示样本Yi在子字典Di下的稀疏编码矩阵,是一个选择向量,且选择向量Ei能够通过Di=DEi从全局字典D选出Di,同样,通过求出与第i类局部字典相关第i类稀疏编码
2)固定P和X,利用以下公式计算全局结构化字典D的原子dl,从而更新全局结构化字典D:
其中E/i=[E1,...,Ei-1,Ei+1,...,Ec], l=1,...,ki,
3)固定X和D,利用以下公式计算降维投影矩阵P,从而更新降维投影矩阵P:
Pt=Pt-1+γ(U(1:d,:)-Pt-1)
其中,Pt=U(1:d,:),t为迭代次数,U(1:d,:)为奇异值矩阵∑对应的最小d个向量,所述U和∑利用求得,其中γ为小常数控制P的平稳,
所述控制阈值为SPE统计量的控制线,按照以下公式计算控制线SPEc:
其中,α为置信水平,{SPE1 2,...,SPEn 2}为{x1,...,xn}的SPE统计量,γ是带宽系数,((·)表示为核方程。
所述在线故障监测包括以下步骤:
1)利用离线字典学习过程得到的降维投影矩阵P按照以下公式计算数据y的SPE统计量SPEnew,
SPEnew=||(I-PTP)xnew||
其中,Xnew为新的数据样本;
2)将SPEnew和控制线SPEc进行比较,若SPEnew超过SPEc,则该数据为故障数据,则进行故障诊断;否则,该数据为正常数据,则进行模式识别。
在所述故障诊断中发生故障的维度K利用下式进行诊断:
SPEnew>SPEc
其中,ξi∈Rm是一个方向向量,它表示单位矩阵的第i列,yi=y-ξifi为y第i维重构样本,其中为第i维的故障幅值,表示的伪逆,P为离线字典学习过程得到的降维投影矩阵。
在所述模式识别中利用以下公式计算数据y的模态类别I:
SPEnew<SPEc
其中,x=[x1;...;xi;...;xc]为稀疏编码矩阵,利用正交匹配跟踪防范求得公式的解,xi表示与子字典Di相关的系数子矢量,P为离线字典学习过程得到的降维投影矩阵,表示与第i类相关的特定子字典。
上述过程的推导过程如下:
一、离线字典学习
考虑到工业过程中,我们从m个传感器上获得n个样本,样本表示成数据矩阵Y=[y1,...,yn]∈Rm×n,我们根据训练样本对应的标签信息将Y分成c个不同的类别,则表示为Y=[Y1,...,Yi,...,Yc],且i=1,...,c,表示第i种模态的样本,且
定义SDRDL的目标方程如下:
s.t.PPT=I. (1)
其中||·||F表示矩阵的Frobenius范数,P∈Rp×m表示学习到的降维投影矩阵,其中p(p<m)表示降维后样本的维数。D=[D1,...,Di,...,Dc]∈Rp×k是被学习到的全局结构化字典,其中表示与第i类相关的特定子字典,且X=[X1,...,Xi,...,Xc]∈Rk×n表示样本Y的稀疏编码矩阵。其中表示样本Yi在全局字典D下的稀疏编码矩阵,且j=1,...,c。表示样本Yi在子字典Di下的稀疏编码矩阵。为了让方程看起来更加简洁,通过引入选择向量且选择向量Ei能够通过Di=DEi从全局字典D选出Di,同样,通过求出与第i类局部字典相关第i类稀疏编码
由于目标方程(1)对于多个优化变量P,D,X而言是非凸,但是将其中的两个固定求解另外一个变量时,则方程(1)就变成了凸函数。这种交替优化的方法在字典学习中被广泛使用,我们的模型使用这种方法对P,D,X进行交替迭代求解,下面将详细介绍SDRDL方法的优化过程。
(1)更新稀疏编码X
首先固定P、D,方程(1)可以被进一步简化:
我们将方程(2)逐步分解成对Xi优化,此时其他所有的Xj(j≠i)被固定,这方程(2)可以进一步简化为:
为了能让式子能够更简洁的表示,可以将上面的式子进一步简化,如下:
显然,方程(4)是一个典型的LASSO问题,我们通过正交匹配跟踪(OMP)方法求解Xi。
(2)更新字典D
固定字典P、X,方程(1)可被进一步简化:
显而易见,上述式子中包含了字典D的c个子字典,所以我们可以上述问题分成c个单独的问题分别进行求解,当更新子字典Di(i=1,...,c)时,其他子字典Dj(j≠i)被固定。则方程(5)进一步进行简化为:
由于方程(6)只包含F-norm,所以优化该式子属于凸优化过程,则我们可以通过求解解析解找到该式的最优解。
为了便于更加方便的求解Di的解析解,首先使用一个简单的式子表示:
我们让其中则可以将Di分成一列列的原子进行单独更新,则更新时,则其他列被固定,即都是固定的,则方程(7)可转换为:
这里的E/i=[E1,...,Ei-1,Ei+1,...,Ec] 这里令则方程进一步转化为;
由于(9)只有唯一的一个变量dl,为了求得dl的解析解,我们可以让这求得结果如下:
因为,如果我们可以得到dk的解析解为:
由于dk Tdk=1;则
否则,则原子dl保持不更新。通过上述过程更新字典中的每一个原子,最后将整个字典进行更新。
(3)更新投影矩阵P
为了求解P,我们保持D和X固定,则方程(1)中的目标方程可以被简化为:
为了让式子变得更简单,这里令 这(12)可以转化成:
由于等式(13)的优化是一个正交约束优化过程,即PPT=I,则由所以,方程(13)可以转化成:
s.t.PPT=I (14)
由于不便于我们求解,则将其拆开,即:
则(14)进一步简化为:
为了解决上述最小化问题,我们使用第t-1次迭代求得的投影矩阵Pt-1去近似Pt,这里使用奇异值分解(SVD)技术,即:
然后我们将P更新为奇异值矩阵∑对应的最小d个向量,即Pt=U(1:d,:),Pt为第t次迭代获得的投影矩阵,但是这样更新P的变化可能很大,使得方程(1)中整个系统的优化过程变得不稳定。为了避免这种情况,我们采用如下方式逐渐的更新P,即:
Pt=Pt-1+γ(U(1:d,:)-Pt-1) (17)
这里的γ是一个小的正常数去控制迭代中P的平稳变化。
这里反复的计算上述过程直到达到最大迭代次数或者收敛,则停止迭代,输出P,D。在下面对具体优化过程进行详细阐述。
实际工作中同时降维和字典学习过程按照以下方式进行迭代:
1:输入:训练样本Y=[Y1,...,Yc]∈Rm×n,累计百分比方差(cPV),参数η,γ,最大迭代次数N;
2:初始化:初始化P(0)作为Y的PCA投影,初始化Di(0)作为Yi的KSVD字典,D(0)=[D1 (0),...,Dc (0)]
3:设置J=1.
当J≤N时,固定字典D(J)和投影矩阵P(J),利用OMP方法根据公式(4)按类计算稀疏编码×(J+1);固定稀疏编码X(J+1)和投影矩阵P(J),根据公式(11)按原子计算字典D(J+1);固定稀疏编码X(J+1)和字典D(J+1),根据公式(17)计算投影矩阵P(J+1);
4:J=J+1,转到步骤3
5:当J>N时,输出投影矩阵P,字典D和稀疏编码X。
通过上述可以求得投影矩阵P,然后通过投影矩阵将原始空间划分成主元空间和残差空间。因此:
是主元空间,Y=PX∈Rd×n是X的投影子空间,是残差空间。接下来,我们将计算原域所对应的平方预测误差(SPE)统计量,根据SDRDL模型,SPE统计量构建如下所示,对于给定一个新的数据样本xnew,它可以被分成两个部分,即:
则SPE被用来测量异常数据的变化,其被计算为:
SPEnew=||(I-PTP)xnew|| (20)
为了确定工业过程是否在正常的条件下运行,必须定义合适的控制线,这里SPE统计量的控制线由核密度估计(KDE)所确定,首先,所有训练样本{x1,...,xn}的SPE统计量{SPE1 2,...,SPEn 2}被估计为:
其中K(·)表示为核方程,在本文中,我们选用最常见的高斯方程作为我们进行KDE的核函数,γ是带宽系数,则在置信水平为α下,SPE统计量的控制线SPEc被计算为:
二、在线故障监测
在离线字典学习阶段,本专利可以得到投影矩阵、字典和控制线。接下来,本专利训练得到的投影矩阵和字典对新数据进行监控和分类。这个过程主要包括以下步骤:
第一步,本专利利用投影矩阵求解测试数据y的SPE统计量,如式(20)所示。
第二步,将求解到的SPEnew与控制线SPEc进行比较,如果超过SPEc,则该数据为故障数据。否则,该数据为正常数据。
第三步,若该数据为故障数据,则进行故障诊断。假设ξi∈Rm是一个方向向量,它表示单位矩阵的第i列,则定义yi=y-ξifi为y第i维重构样本,其中fi为第i维的故障幅值。类比式子(20),这里我们构造基于RBC的SPE统计量如下:
其中而重构的任务就是找到fi,使得SPE(yi)最小化,这个最小化是通过SPE(yi)对fi求导后等于0实现的.即
令方程(24)等于0,并求出fi,则其中表示的伪逆。对于单维故障,如果哪一维减去故障后,其重构量变为最小,则这一维即为故障,并对故障进行隔离。即发生故障维度k的数学定义为:
SPEnew>SPEc (25)
若该数据为正常数据,则对该数据进行模式识别。我们先通过Py将测试样本投影到低维空间中,然后在结构化字典上对低维样本Py进行编码。即通过求解下面这个等式进行稀疏编码:
μ是一个正标量,稀疏编码x可能写出x=[x1;...;xi;...;xc],其中xi表示与子字典Di相关的系数子矢量。式(26)可以通过OMP求出其稀疏编码,因此我们可以通过下式求得正常数据的类别l,即:
SPEnew<SPEc (27)
传感器的高维数据包含了大量的信息,但是由于各个变量之间的存在很强的相关性,使得这些高维数据包含很多冗余信息。在故障监控中,一般通过降维消除这些冗余信息,其中使用最广泛的是主成分分析(PCA)、偏最小二乘(PLS),局部嵌入投影(LPP)等多元统计监控方法。这些方法是将高度相关的高维过程数据投影到低维子空间,通过数据变化最大的几个主成分作为数据的主要特征。其中PCA虽然保留了数据的全局结构,但是却破坏了数据间的局部关系,而LPP主要保留过程数据的邻接结构,但是没有明确考虑数据集的全局结构,会导致方差信息的丢失。传统的多元统计监控方法假设工业过程都是在同一种操作工况条件下运行,然而由于设定值,操作条件,季节变化和工业设备的老化等因素的影响,实际工业过程往往工作在多个正常操作模态下。现有的多模态过程监控方法都是针对每一种模态分别建模,而不能用统一结构表达,从而影响了监控效果。近年来,字典学习(DL)在计算机视觉,图像处理等领域已经取得了很大的成功。字典学习方法的基本思想是位于高维环境空间的数据点存在于多个低维子空间,并且相同类的数据点聚集在一起构成低维子空间。因此数据的内在数据结构对字典学习方法的分类效果很重要。字典学习在故障监控方面越来越受欢迎,虽然其无法同时解决大量工业过程中数据存在高维度和多模态特性的问题,但是字典学习良好的表示能力能够很好解决过程数据多模态特性的问题。
本发明基于同时降维和字典学习的过程监控方法联合投影和字典交替迭代求解,该过程能学习到保留原始数据更多的空间信息的投影和结构化的字典。在离线学习阶段通过同时降维和字典学习方法学习到的投影和字典能同时解决数据高维和多模态特点的问题。在在线故障监测阶段,通过投影矩阵对测试数据构造SPE统计量,能减小在线监测过程的计算复杂度,提高在线监测的实时性。
虽然以上将实施例分开说明和阐述,但涉及部分共通之技术,在本领域普通技术人员看来,可以在实施例之间进行替换和整合,涉及其中一个实施例未明确记载的内容,则可参考有记载的另一个实施例。
以上所述的实施方式,并不构成对该技术方案保护范围的限定。任何在上述实施方式的精神和原则之内所作的修改、等同替换和改进等,均应包含在该技术方案的保护范围之内。
Claims (5)
1.一种基于同时降维和字典学习的过程监控方法,其特征在于,包括:
利用历史数据进行离线字典学习,得到全局结构化字典、降维投影矩阵和控制阈值;
利用全局结构化字典、降维投影矩阵和控制阈值进行在线故障监测;
所述离线字典学习包括同时降维和字典学习和计算控制阈值,其中,同时降维和字典学习包括以下步骤:
1)固定P和D,利用正交匹配跟踪方法计算以下公式的解Xi,从而更新稀疏编码矩阵X:
其中||·||F表示矩阵的Frobenius范数,P∈Rp×m表示学习到的降维投影矩阵,其中p(p<m)表示降维后样本的维数,Y=[Y1,...,Yi,...,Yc]表征样本,且i=1,...,c,c为模态的类数,表示第i类模态的样本,且D=[D1,...,Di,...,Dc]∈Rp×k是被学习到的全局结构化字典,表示与第i类相关的特定子字典,且X=[X1,...,Xi,...,Xc]∈Rk×n表示样本Y的稀疏编码矩阵,表示样本Yi在全局结构化字典D下的稀疏编码矩阵,且j=1,...,c,表示样本Yi在子字典Di下的稀疏编码矩阵,是一个选择向量,且选择向量Ei能够通过Di=DEi从全局字典D选出Di,同样,通过求出与第i类局部字典相关第i类稀疏编码
2)固定P和X,利用以下公式计算全局结构化字典D的原子dl,从而更新全局结构化字典D:
其中Eli=[E1,...,Ei-1,Ei+1,...,Ec], l=1,...,ki,
3)固定X和D,利用以下公式计算降维投影矩阵P,从而更新降维投影矩阵P:
Pt=Pt-1+γ(U(1:d,:)-Pt-1)
其中,Pt=U(1:d,:),t为迭代次数,U(1:d,:)为奇异值矩阵∑对应的最小d个向量,所述U和∑利用求得,其中γ为小常数控制P的平稳,
2.根据权利要求1所述的基于同时降维和字典学习的过程监控方法,其特征在于:所述控制阈值为SPE统计量的控制线,按照以下公式计算控制线SPEc:
其中,α为置信水平,{SPE1 2,...,SPEn 2}为{x1,...,xn}的SPE统计量,γ是带宽系数,K(·)表示为核方程。
3.根据权利要求1所述的基于同时降维和字典学习的过程监控方法,其特征在于:所述在线故障监测包括以下步骤:
1)利用离线字典学习过程得到的降维投影矩阵P按照以下公式计算数据y的SPE统计量SPEnew,
SPEnew=||(I-PTP)xnew||
其中,Xnew为新的数据样本;
2)将SPEnew和控制线SPEc进行比较,若SPEnew超过SPEc,则该数据为故障数据,则进行故障诊断;否则,该数据为正常数据,则进行模式识别。
4.根据权利要求3所述的基于同时降维和字典学习的过程监控方法,其特征在于:在所述故障诊断中发生故障的维度K利用下式进行诊断:
SPEnew>SPEc
其中,ξi∈Rm是一个方向向量,它表示单位矩阵的第i列,yi=y-ξifi为y第i维重构样本,其中为第i维的故障幅值,表示的伪逆,P为离线字典学习过程得到的降维投影矩阵。
5.根据权利要求3所述的基于同时降维和字典学习的过程监控方法,其特征在于:在所述模式识别中利用以下公式计算数据y的模态类别I:
SPEnew<SPEc
其中,x=[x1;...;xi;...;xc]为稀疏编码矩阵,利用正交匹配跟踪方法求得公式的解,xi表示与子字典Di相关的系数子矢量,P为离线字典学习过程得到的降维投影矩阵,表示与第i类相关的特定子字典。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910902691.7A CN110579967B (zh) | 2019-09-23 | 2019-09-23 | 基于同时降维和字典学习的过程监控方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910902691.7A CN110579967B (zh) | 2019-09-23 | 2019-09-23 | 基于同时降维和字典学习的过程监控方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110579967A true CN110579967A (zh) | 2019-12-17 |
CN110579967B CN110579967B (zh) | 2020-06-30 |
Family
ID=68813446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910902691.7A Expired - Fee Related CN110579967B (zh) | 2019-09-23 | 2019-09-23 | 基于同时降维和字典学习的过程监控方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110579967B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111340110A (zh) * | 2020-02-25 | 2020-06-26 | 湖南师范大学 | 一种基于工业过程运行状态趋势分析的故障预警方法 |
CN111829425A (zh) * | 2020-08-06 | 2020-10-27 | 厦门航空有限公司 | 民机前缘位置传感器的健康监测方法及系统 |
CN112541554A (zh) * | 2020-12-18 | 2021-03-23 | 华中科技大学 | 基于时间约束的核稀疏表示的多模态过程监控方法及系统 |
CN113129122A (zh) * | 2021-04-30 | 2021-07-16 | 国家电网有限公司 | 财务风险预警审计方法、装置、电子设备、及存储介质 |
WO2021185044A1 (zh) * | 2020-03-19 | 2021-09-23 | 中南大学 | 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质 |
CN113609602A (zh) * | 2021-06-21 | 2021-11-05 | 国网浙江省电力有限公司湖州供电公司 | 一种基于字典学习的卸载式多驱输送机故障监测方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930790A (zh) * | 2016-04-19 | 2016-09-07 | 电子科技大学 | 基于核稀疏编码的人体行为识别方法 |
CN106980848A (zh) * | 2017-05-11 | 2017-07-25 | 杭州电子科技大学 | 基于曲波变换和稀疏学习的人脸表情识别方法 |
US20180121759A1 (en) * | 2016-10-28 | 2018-05-03 | International Business Machines Corporation | Simultaneous feature extraction and dictionary learning using deep learning architectures for characterization of images of heterogeneous tissue samples |
CN108509833A (zh) * | 2017-12-21 | 2018-09-07 | 苏州大学 | 一种基于结构化分析字典的人脸识别方法、装置及设备 |
CN108985177A (zh) * | 2018-06-21 | 2018-12-11 | 南京师范大学 | 一种结合稀疏约束的快速低秩字典学习的人脸图像分类方法 |
CN109766748A (zh) * | 2018-11-27 | 2019-05-17 | 昆明理工大学 | 一种基于投影变换与字典学习的行人再识别的方法 |
US20190180414A1 (en) * | 2017-12-08 | 2019-06-13 | Tata Consultancy Services Limited | Systems and methods for reconstructing super-resolution images under total aliasing based upon translation values |
CN109921799A (zh) * | 2019-02-20 | 2019-06-21 | 重庆邮电大学 | 一种基于聚能量字典学习的张量压缩方法 |
CN110265039A (zh) * | 2019-06-03 | 2019-09-20 | 南京邮电大学 | 一种基于字典学习和低秩矩阵分解的说话人识别方法 |
-
2019
- 2019-09-23 CN CN201910902691.7A patent/CN110579967B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105930790A (zh) * | 2016-04-19 | 2016-09-07 | 电子科技大学 | 基于核稀疏编码的人体行为识别方法 |
US20180121759A1 (en) * | 2016-10-28 | 2018-05-03 | International Business Machines Corporation | Simultaneous feature extraction and dictionary learning using deep learning architectures for characterization of images of heterogeneous tissue samples |
CN106980848A (zh) * | 2017-05-11 | 2017-07-25 | 杭州电子科技大学 | 基于曲波变换和稀疏学习的人脸表情识别方法 |
US20190180414A1 (en) * | 2017-12-08 | 2019-06-13 | Tata Consultancy Services Limited | Systems and methods for reconstructing super-resolution images under total aliasing based upon translation values |
CN108509833A (zh) * | 2017-12-21 | 2018-09-07 | 苏州大学 | 一种基于结构化分析字典的人脸识别方法、装置及设备 |
CN108985177A (zh) * | 2018-06-21 | 2018-12-11 | 南京师范大学 | 一种结合稀疏约束的快速低秩字典学习的人脸图像分类方法 |
CN109766748A (zh) * | 2018-11-27 | 2019-05-17 | 昆明理工大学 | 一种基于投影变换与字典学习的行人再识别的方法 |
CN109921799A (zh) * | 2019-02-20 | 2019-06-21 | 重庆邮电大学 | 一种基于聚能量字典学习的张量压缩方法 |
CN110265039A (zh) * | 2019-06-03 | 2019-09-20 | 南京邮电大学 | 一种基于字典学习和低秩矩阵分解的说话人识别方法 |
Non-Patent Citations (3)
Title |
---|
HIROYUKI KASAI等: "Riemannian joint dimensionality reduction and dictionary learning on symmetric positive definite manifolds", 《2018 26TH EUROPEAN SIGNAL PROCESSING CONFERENCE (EUSIPCO)》 * |
XIANWEI等: "Reconstructible Nonlinear Dimensionality Reduction via Joint Dictionary Learning", 《IEEE TRANSACTIONS ON NEURAL, NETWORKS AND LEARNTNG SYSTEMS》 * |
郑思龙等: "基于字典学习的非线性降维方法", 《自动化学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111340110A (zh) * | 2020-02-25 | 2020-06-26 | 湖南师范大学 | 一种基于工业过程运行状态趋势分析的故障预警方法 |
CN111340110B (zh) * | 2020-02-25 | 2023-04-18 | 湖南师范大学 | 一种基于工业过程运行状态趋势分析的故障预警方法 |
WO2021185044A1 (zh) * | 2020-03-19 | 2021-09-23 | 中南大学 | 基于迁移学习的重金属废水处理过程异常工况智能化监测方法、装置及存储介质 |
CN111829425A (zh) * | 2020-08-06 | 2020-10-27 | 厦门航空有限公司 | 民机前缘位置传感器的健康监测方法及系统 |
CN112541554A (zh) * | 2020-12-18 | 2021-03-23 | 华中科技大学 | 基于时间约束的核稀疏表示的多模态过程监控方法及系统 |
CN112541554B (zh) * | 2020-12-18 | 2024-03-22 | 华中科技大学 | 基于时间约束的核稀疏表示的多模态过程监控方法及系统 |
CN113129122A (zh) * | 2021-04-30 | 2021-07-16 | 国家电网有限公司 | 财务风险预警审计方法、装置、电子设备、及存储介质 |
CN113609602A (zh) * | 2021-06-21 | 2021-11-05 | 国网浙江省电力有限公司湖州供电公司 | 一种基于字典学习的卸载式多驱输送机故障监测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110579967B (zh) | 2020-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110579967B (zh) | 基于同时降维和字典学习的过程监控方法 | |
Kadlec et al. | Data-driven soft sensors in the process industry | |
Zhang et al. | Stationarity test and Bayesian monitoring strategy for fault detection in nonlinear multimode processes | |
Fezai et al. | Online reduced kernel GLRT technique for improved fault detection in photovoltaic systems | |
Hollmén | User profiling and classification for fraud detection in mobile communications networks | |
CN111340110B (zh) | 一种基于工业过程运行状态趋势分析的故障预警方法 | |
Ayodeji et al. | Causal augmented ConvNet: A temporal memory dilated convolution model for long-sequence time series prediction | |
Wang et al. | Time-weighted kernel-sparse-representation-based real-time nonlinear multimode process monitoring | |
Chen et al. | Full decoupling high-order dynamic mode decomposition for advanced static and dynamic synergetic fault detection and isolation | |
Ferwerda et al. | KRLS: A Stata package for kernel-based regularized least squares | |
Kong et al. | Semi-supervised classification of wafer map based on ladder network | |
CN111639304A (zh) | 基于Xgboost回归模型的CSTR故障定位方法 | |
CN109324595B (zh) | 一种基于增量pca的工业监测数据分类方法 | |
CN113420815A (zh) | 半监督rsdae的非线性pls间歇过程监测方法 | |
Othman et al. | Abnormal patterns detection in control charts using classification techniques | |
Idé et al. | Change Detection Using Directional Statistics. | |
CN108009586B (zh) | 封顶概念分解方法及图像聚类方法 | |
CN115470838A (zh) | 一种基于鲁棒动态变分贝叶斯字典学习的过程监测方法 | |
Han et al. | Grounding learning of modifier dynamics: An application to color naming | |
CN114036947A (zh) | 一种半监督学习的小样本文本分类方法和系统 | |
Natsumeda et al. | Feature Selection with Partial Autoencoding for Zero-Sample Fault Diagnosis | |
Wang et al. | Knockoff-Guided Feature Selection via A Single Pre-trained Reinforced Agent | |
Wang et al. | Historical Information-Aided Monitoring of Few-Sample Modes in Industrial Processes With Orthogonal Transferred Projection | |
Donat et al. | Data visualization, data reduction and classifier fusion for intelligent fault detection and diagnosis in gas turbine engines | |
Guo et al. | A Fault Diagnosis Method using Convolutional Neural Network with Visualized Samples |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200630 Termination date: 20210923 |
|
CF01 | Termination of patent right due to non-payment of annual fee |