[发明专利]一种自动的步进式有序时段划分方法

权利要求书

1.一种自动的步进式有序时段划分方法，其特征在于，该方法包括以下步骤：步骤1：获取过程分析数据：设一个间歇操作具有J个测量变量和K个采样点，则每一个测量批次可得到一个K×J的矩阵，重复I批次的测量步骤后，得到的数据可以表述为一个三维矩阵其中测量变量为温度、速度、压力、位移等批次运行过程中可被测量的状态参数；

步骤2：数据预处理：将三维矩阵按照采集批次方向展开，即将一个操作批次内的各采样点上的变量按照时间顺序排开得到二维矩阵X(I×KJ)，由K个时间片矩阵X_k(I×J)组成，其中，下标k为时间指标；

设二维矩阵X_k内任意一点的变量为x_k，i，j对该变量进行减均值、除以标准差的标准化处理，其中，下标i代表批次，j代表变量，标准化处理的计算公式如下：

xk,i,j=xk,i,j-x‾k,jsk,j;---(1)]]>

其中：k是时间片指标。是X_k矩阵任一列的均值，s_k，j是X_k矩相应列的标准差，

x‾k,j=1IΣi=1Ixk,i,j,]]>

sk,j=Σi=1I(xk,i,j-x‾k,j)2/(I-1);---(2)]]>

步骤3：时间片PCA建模，该步骤由以下子步骤来实现：

（3.1）对步骤2标准化处理后的每一个时间片矩阵X_k(I×J)执行PCA分解，建立时间片PCA模型，其中PCA分解公式如下：

Xk=TkPkT=Σr=1Jtk,rpk,r;---(3)]]>

其中：t_k，n为正交的主元向量，p_k，n为正交归一化的负载向量，r表示不同的PCA分解方向，上标T表示矩阵的转置；T_k(I×J)代表保留全部主元的得分矩阵，P_k(J×J)代表对应的负载矩阵。

（3.2）选取主元个数，将公式（3）重新表述成如下形式：

Xk=TkPkT+Ek=Σr=1Rtk,rpk,r+Ek---(4)]]>

其中：r表示不同的PCA分解方向；T_k(I×R_k)与P_k(J×R_k)代表负载矩阵分别为保留R_k个主元后的得分矩阵和负载矩阵，E_k为残差矩阵。通过上述变换，多向主元分析法模型将原始数据空间分解为主元空间和残差空间，主元空间内代表主要的系统过程波动信息；这里所保留的主元个数R_k能够反映原过程中90％的过程波动信息。综合考虑全过程，选取出现次数最多的R_k值统一为最终的建模主元个数R，即所有时间片模型保留相同的主元个数。

（3.3）计算残差空间中各时间片k中对应各个批次的SPE指标：

SPE_k，i=e_k，i^Te_k，i      （5）

其中，下标i表示时间片中的批次，e_k，i是对应k时刻第i批次的残差列向量。根据相同时刻上不同批次的SPE值服从带权重系数的χ²分布，从而确定出每个时间点上的控制限Ctr_k，它反应了时间片PCA模型的重构能力。

步骤4：确定基于变量展开模型的SPE指标控制限：从间歇过程初始点开始，依次将下一个时间片与之前的时间片组合在一起并按变量方式展开得到时间块X_v，k(IK×J)，其中下标v代表变量展开方式。对新时间块矩阵进行PCA分析，提取出负载矩阵P_v，k(J×R)。计算其SPE值并根据相同时刻上不同批次的SPE值服从带权重系数的χ²分布，从而确定出每个时间点上的控制限Ctr_k。

步骤5：确定第一时段划分点k^*：比较在相同时间区域内的每个时间点上Ctr_k和Ctr_v，k的大小，如果发现连续三个样本呈现Ctr_v，k＞αCtr_k，那么新加入的时间片对该时间块的PCA监测模型及相应的监测性能都有重大的影响。记加入新时间片前的时刻为k^*。其中，α是依附于Ctr_k的常数，称作缓和因子，它反映的是与时间片模型相比，时间块模型允许监测精度损失的程度。则k^*时刻之前的时间片可认为是一个子时段。

步骤6：过程分析数据更新，确定所有划分时段：根据步骤5中所获得的时刻k^*的指示，移除第一个子时段，把余下的间歇过程数据作为新的输入数据带入到第5步中并重复上述步骤5-6，划分不同时间段，直到没有数据余留。

步骤7：基于时段划分结果的统计建模：根据步骤6时段划分结果，每个时段内的时间片按照变量展开方式组合成子时段代表性建模数据组，X_c(IK_c×J)，其中，下标c是时段指标。每个时段内相似的过程潜在特性可以通过对X_c(IK_c×J)实施PCA分解提取出来：

T_c=X_cP_c

X^c=TcPcT=XcPcPcT---(6)]]>

Ec=Xc-X^c]]>

其中，P_c(J×R_c)是该子时段的PCA负载矩阵，揭示了该子时段内的主要波动方向，R_c代表了该子时段模型保留的PCA主元个数。T_c(IK_c×R_c)代表了该子时段的主元得分矩阵。因此，该时段的主要波动通过T_cP_c^T表征，代表了该子时段内的主要波动信息；剩余为子时段残差矩阵，E_c(IK_c×J)代表了该子时段内的噪声信息。各时间片主元得分，T_k(I×R_c)，可以很容易从子时段得分矩阵T_c(IK_c×R_c)中根据对应的过程时间抽取获得，从而可以据此计算得出各采样时刻上的协方差关系S_k。各时间片残差矩阵E_k(I×J)亦可以从子时段残差矩阵E_c(IK_c×J)中对应获得。

步骤8：计算实时监测统计指标：根据从公式（6）计算的结果中获取的主元得分时间片T_k(I×R_c)和残差矩阵时间片E_k(I×J)可以在每个时刻计算两个监测统计指标：Hotelling-T²统计指标和SPE统计量。

Hotelling-T²统计指标用来测量各采样时刻过程变量偏离正常工况下平均轨迹的距离，并且其在显著性水平α下的控制限计算为：

Ti,k2-(ti,k-t‾k)TSk-1(ti,k-t‾k)~Rc(I-1)I-RcFRc,I-Rc,α---(7)]]>

其中，R_c是该时段PCA模型中所保留的主元个数；t_i，k(R_c×1)是第k时刻第i个批次的主元得分，即对应时间片得分矩阵T_k(I×R_c)的第i行；而是T_k(I×R_c)针对不同批次的均值向量，由于各时间片测量数据已经在数据预处理时中心化为零均值，这里的其实就是零向量。

对于残差子空间，各个时刻不同批次的SPE统计量计算为：

SPEi,k(xi,k-x^i,k)T(xi,k-x^i,k)---(8)]]>

其中，x_i，k代表k时刻第i个批次的过程测量向量，则是由PCA模型重构的对应结果。式（8）中的计算结果在每个时刻都可以构成一个I×1向量[SPE_1，k，SPE_2，k,...,SPE_I，k]^T，并且该向量近似服从加权χ²分布，从而获得每个时刻SPE的监测控制限。

步骤9：基于时段模型的在线过程监测：基于步骤6划分的时段、步骤7建立的时段模型监测系统和步骤8所得的两个监测统计量可以在线监测注塑成型等新运行间歇过程的状态。该步骤由以下子步骤来实现：

（9.1）采集新测量数据及新测量数据预处理：在线监测时，采集到新的过程测量数据x_new(J×1)后，其中，下标new代表新样本，J为测量变量，与步骤1中的测量变量相同；首先需要进行数据预处理。根据步骤2中获得的均值和标准差，根据过程时间的指示调用对应该时刻的均值和标准差对采集到的过程测量数据如公式（1）中所示进行标准化预处理。

（9.2）计算新监测统计量：数据预处理后，根据公式（6）计算的PCA子时段模型，调用对应该新采样时刻所在时段的模型P_c(J×R_c)，按照如下方式计算得到主元得分，估计残差及其对应的Hotelling-T²与SPE两个监测统计指标：

t_new^T=x_new^TP_c

x^new=Pctnew]]>(9)]]>

Tnew2=(tnew-t‾k)TSk-1(tnew-t‾k)]]>

SPEnew=(xnew-x^new)T(xnew-x^new)]]>

其中是新的过程测量数据，是前面根据训练数据获得的T_k的均值向量，S_k是X_k的协方差矩阵。

（9.3）在线判断过程运行状态：实时比较两个监测指标与其各自的统计控制限，若两个监测指标都位于统计控制限之内，表明过程运行正常；若有一个以上监测指标超出正常控制限，表明过程有异常状况发生。

2.根据权利要求1所述自动的步进式有序时段划分方法，其特征在于，所述步骤1中，所述测量变量为以下9个：压力阀开度、流量阀开度、螺杆行程、螺杆速度、注射压力、喷嘴温度、机桶头部温度、机桶中间温度和机桶尾部温度。