【技术实现步骤摘要】
本专利技术涉及一种数据驱动的过程监测方法，尤其涉及一种基于多子空间正交典型相关分析的过程监测方法。
技术介绍
1、近年来，由于工业大规模化，生产过程中采集到了大量的数据。在这种背景下，数据驱动的过程监测愈加得到重视并被广泛研究。然而，由于工业过程中闭环控制的存在，并非每个故障都会导致产品质量的恶化。根据故障是否影响产品质量，故障可分为三类：质量相关、质量半相关和质量无关。当质量相关的故障发生时，需要立即诊断故障以采取措施使其恢复正常。从成本效益的角度来看，当与质量无关的故障发生时，不需要花费大量的精力和资源来立即采取行动。为了判断故障是否影响产品质量，人们对质量相关过程监测方法进行了广泛的研究，例如，偏最小二乘(partial least squares，pls)、主成分回归分析(principal component regression，pcr)、典型相关分析(canonical correlationanalysis，cca) 算法等等。此外，由于闭环控制系统的存在，生产过程往往是动态的，有些传感器采集到的数据具有较强的时间序列相关性，而有些传感器采集到的数据具有较弱的时间序列相关性，甚至没有时间序列相关性。如果在过程监测模型中忽略时间尺度上的序列相关性，会导致监测结果不理想，增加误检和漏检的风险。因此，考虑到时间序列的相关性，可以根据过程知识、专家知识、统计数据分析和工程经验将过程变量划分为时间序列相关子空间和时间序列不相关子空间，同时设置时滞对时间序列相关子空间进行矩阵增广以考虑时序相关性。另外，一些过程变量与质量变量
技术实现思路
1、本专利技术所要解决的主要技术问题是：在考虑采样数据时间序列相关性的基础上，提取采样数据的质量相关信息和质量不相关信息，并在此基础上实施过程监测。具体来讲，本专利技术首先将原始信息空间划分为四个不同的子空间；其次，根据不同子空间中的数据特征，利用正交cca、pca方法在不同的子空间中提取特征分别建立对应的监测模型，进而构建各子空间内监测统计量；最后，采用贝叶斯融合策略构建综合质量相关监测统计量和综合质量无关监测统计量以检测故障是否发生。
2、本专利技术解决上述技术问题所采用的技术方案为：一种基于多子空间正交典型相关分析的过程监测方法，包括以下所示步骤：
3、步骤(1)：采集生产过程正常运行状态下的样本，组成训练数据集x∈rn×m，y∈rn×p，其中n为训练样本数，m为过程测量变量数，p为质量测量变量数，r为实数集，rn×m表示n×m 维的实数矩阵，rn×p表示n×p维的实数矩阵；
4、步骤(2)：对矩阵x中各列向量实施z-score标准化处理，得到标准化后的矩阵：其中xi∈rm×1(i＝1,2,…,n)，yi∈rp×1(i＝ 1,2,…,n)；
5、步骤(3)：采用lasso回归度量过程变量与质量变量的相关性，并采用自相关函数(acf) 度量过程变量的自相关性，从而将过程数据集划分为四个子空间：质量相关-序列相关子空间：a＝[a1,a2,…,an]t∈rn×a，质量无关-序列相关子空间：b＝[b1,b2,…,bn]t∈rn×b，质量无关-序列不相关子空间：c＝[c1,c2,…,cn]t∈rn×c，质量相关-序列不相关子空间：d＝ [d1,d2,…,dn]t∈rn×d。其中a+b+c+d＝m。具体的实施过程如下所示：
6、步骤(3.1)：采用lasso算法构建回归模型，作为输入，yj(j＝1,2,…,p)作为输出。 lasso方法估计参数如下：
7、
8、∑i|βi|≤t (1)
9、其中，k＝1,2,…,n；i＝1,2,…,m。t是调整参数。对于所有的t，当时，α可忽略。
10、步骤(3.2)：由式(1)可得β＝[β1,β2,…,βm]t，则满足下式的变量可以被挑选为质量相关的过程变量
11、
12、其中d<m且ε取为一个极小的正数。
13、步骤(3.3)：对于每个输出变量yj(j＝1,2,…,p)，通过lasso方法获得对应的质量相关过程变量集，并取并集作为整体质量相关变量集将剩余变量作为质量无关变量集
14、步骤(3.4)：自相关函数(acf)度量了xt(i)与xt+k(i)之间的相关性，其中k为滞后。对于 x(i)，其滞后k阶的自相关式为：
15、acf(k)＝ck/c0 (3)
16、
17、
18、置信区间：
19、
20、对于每个过程变量，若其acf在置信区间内，则认为该变量没有时间序列相关性；否则，认为该变量具有时间序列相关性。因而，根据acf将训练数据集划分为序列相关子空间和序列不相关子空间
21、步骤(3.5)：根据下式，将过程数据集划分为四个子空间：
22、a＝xq∩xd (7)
23、b＝xu∩xd (8)
24、c＝xu∩xn (9)
25、d＝xd∩xn (10)
26、其中，a＝[a1,a2,…,an]t∈rn×a为质量相关-序列相关子空间，b＝[b1,b2,…,bn]t∈rn×b为质量无关-序列相关子空间，c＝[c1,c2,…,cn]t∈rn×c为质量无关-序列不相关子空间，d＝ [d1,d2,…,dn]t∈rn×d为质量相关-序列不相关子空间，且a+b+c+d＝m。
27、步骤(4)：对质量相关-序列相关子空间a∈rn×a和质量无关-序列相关子空间b∈rn×b，进行矩阵增广，将时滞设置为2，从而得到新的数据集和
28、步骤(5)：在数据扩展后的质量相关-序列相关子空间中以及在质量相关-序列不相关子空间中，分别构建正交cca模型。以数据扩展后的质量相关-序列相关子空间为例，具体的实施过程如下所示：
29、步骤(5.1)：在数据扩展后的质量相关-序列相关子空间中，a'∈r(n-2)×3a和y'∈r(n-2)×3p，构造矩阵并对k进行svd分解：
30、k＝r∑vt (11)
31、其中则典型相关变量系数矩阵
32、步骤(5.2)：构造矩阵h＝jλl-1，并对h进行svd分解：
33、ht＝u[ξ,0][vy,vn]t (12)
34、从而得到了质量相关投影矩阵vy和质量无关投影矩阵vn；
35、步骤(5.3)：然后，依据如下公式
36、
37、
本文档来自技高网...【技术保护点】
1.一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，所述步骤(3)中将过程数据集划分为四个子空间具体实施过程如下所示：
3.根据权利要求1所述的一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，所述步骤(5)中建立相应的正交CCA模型的具体实施过程如下所示：
4.根据权利要求1所述的一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，所述步骤(6)中建立相应的PCA模型的具体实施过程如下所示：
5.根据权利要求1所述的一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，所述步骤(7)中根据贝叶斯融合策略构建两个概率型指标BICq和BICu的具体实施过程如下所示：
【技术特征摘要】
1.一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，所述步骤(3)中将过程数据集划分为四个子空间具体实施过程如下所示：
3.根据权利要求1所述的一种基于多子空间正交典型相关分析的过程监测方法，其特征在于，所述步骤(5)中建立相应的正交cca模...
【专利技术属性】
技术研发人员：宋冰，王炯，侍洪波，陶阳，谭帅，张天清，
申请(专利权)人：华东理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。