主权项 |
一种由数据差异驱动的间歇过程自学习动态优化方法,其特征在于包括下述步骤:(1)针对操作完整的间歇过程,按批次收集待优化变量和最终质量或产率指标;(2)对所述步骤(1)中采集的数据,以批次为变量进行主元分析并在主元模式图中剔除奇异点,使得所有数据点在一个可信度之内;(3)将剔除奇异点后的操作数据在时间轴上间隔划分成N段;将每一个时间间隔所包含的各批次数据表达为一个连续变量C<sub>i</sub>,i=1,2,…N,这些变量被称为时段变量;所述时段变量的值由待优化变量在一个特定时间区间的各个批次数据所组成;由多批次时段变量构成的数据矩阵称为时段变量矩阵,记为L<sub>p×N</sub>,p是批次数;(4)将步骤(3)中所对应的每个批次质量或收率指标,称为指标变量Y<sub>p×1</sub>;所述指标变量的值是由p个批次最终质量或收率形成的连续变量;(5)根据步骤(3)和步骤(4)中形成的时段变量矩阵L<sub>p×N</sub>和指标变量Y<sub>p×1</sub>,分别计算协方差矩阵S<sub>LL</sub>和联合协方差矩阵S<sub>LY</sub>;(6)对协方差矩阵S<sub>LL</sub>和联合协方差矩阵S<sub>LY</sub>做主元分析并得到PLS系数向量F<sub>i</sub>,i=1,2,…N;(7)对步骤(6)中的PLS系数变量元素,按符号大小进行分类,定义作用符号如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>s</mi><mi>i</mi><mi>g</mi><mi>n</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = "}"><mtable><mtr><mtd><mrow><mo>+</mo><mn>1</mn></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>F</mi><mi>i</mi></msub><mo>≥</mo><mi>e</mi></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mo>|</mo><msub><mi>F</mi><mi>i</mi></msub><mo>|</mo><mo><</mo><mi>e</mi></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>F</mi><mi>i</mi></msub><mo>≤</mo><mo>-</mo><mi>e</mi></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mi>N</mi></mrow>]]></math><img file="FDA0000858207840000011.GIF" wi="805" he="238" /></maths>其中e是对噪声的阈值限定,sign(i)是第i个时段对应的PLS系数符号;(8)计算出各个时段变量的均值和标准差;按以下摄动量计算公式建立已采集批次数据各时段变量的初始优化策略:J<sub>i</sub>=M<sub>i</sub>+sign<sub>i</sub>×3σ<sub>i</sub>此处的J<sub>i</sub>,M<sub>i</sub>和σ<sub>i</sub>分别是第i个时段变量的优化目标值,均值和标准差;(9)将步骤(8)中所得各时段的优化目标值,按时段顺序i=1,2,…N组合成一条针对整个批次过程的初始优化变量曲线;(10)收集新批次的时段变量C<sub>i</sub>(k+1),i=1,2,…N和指标变量数据Y(k+1),利用以下递归公式,计算更新协方差矩阵S<sub>LL</sub>(k+1)和联合协方差矩阵S<sub>LY</sub>(k+1):S<sub>LL</sub>(k+1)=λS<sub>LL</sub>(k)+C(k+1)<sup>T</sup>C(k+1)S<sub>LY</sub>(k+1)=λS<sub>LY</sub>(k)+C(k+1)<sup>T</sup>Y(k+1)其中C(k+1)=[C<sub>1</sub>(k+1),C<sub>2</sub>(k+1)…C<sub>N</sub>(k+1)],0<λ<1是对现有协方差矩阵的遗忘因子;当λ=1时,代表没有数据从旧的协方差矩阵中去除;(11)对协方差矩阵S<sub>LL</sub>(k+1)和联合协方差矩阵S<sub>LY</sub>(k+1)做主元分析并得到PLS系数向量F<sub>i</sub>(k+1);(12)利用如下递归公式,计算新批次数据下各时段变量的均值和标准差:M<sub>i</sub>(k+1)=M<sub>i</sub>(k)+[C<sub>i</sub>(k+1)‑M<sub>i</sub>(k)]/(k+1)σ<sub>i</sub>(k+1)=σ<sub>i</sub>(k)+[C<sub>i</sub>(k+1)‑M<sub>i</sub>(k)]×[C<sub>i</sub>(k+1)‑M<sub>i</sub>(k+1)];(13)按步骤(8)中的摄动量公式计算新批次数据下各时段变量的优化目标值;(14)将步骤(13)中所得各时段的优化目标值,按时段顺序i=1,2,…N结合构成一个新的优化变量曲线;(15)判断是否有新数据更新,有新数据转步骤(10),继续做自学习更新运算,否则结束学习过程。 |