发明名称 基于集成型独立元回归模型的工业过程软测量方法
摘要 本发明涉及一种基于集成型独立元回归模型的工业过程软测量方法,应用于非高斯性工业过程数据。传统的非高斯性软测量回归建模方法需要选择一个的非二次函数,以度量非高斯性大小。然而,不同的工业过程数据或对象会造成实际应用中难以获取足够多的经验知识去指导非二次函数的选择。为此,本发明所涉及的方法通过全面而充分地利用不同的非二次函数来训练得到不同的软测量模型,有效的避免了非二次函数的选择问题。然后,通过加权系数累加得到最终的预测结果,使相应软测量模型的预测精度不再受到非二次函数选择的影响。这大大的提高了软测量模型的预测效果,从而能够对过程中的关键指标或质量指标进行更加精确而可靠的预测。
申请公布号 CN106094786A 申请公布日期 2016.11.09
申请号 CN201610389227.9 申请日期 2016.05.30
申请人 宁波大学 发明人 童楚东;蓝艇
分类号 G05B23/02(2006.01)I 主分类号 G05B23/02(2006.01)I
代理机构 代理人
主权项 一种基于集成型独立元回归模型的工业过程软测量方法,其特征在于,该方法主要包括以下几个步骤:(1)利用集散控制系统收集工业生产过程中容易测量的数据组成软测量模型的输入训练数据矩阵X∈R<sup>n×m</sup>,并对其进行标准化处理使各个过程变量的均值为0,标准差为1,得到新数据矩阵<img file="FSA00001307833600000119.GIF" wi="72" he="56" />其中,n为训练样本数,m为过程测量变量数,R为实数集,R<sup>n×m</sup>表示n×m维的实数矩阵;(2)采用离线分析手段获取与输入训练数据X相对应的产品成分或质量数据组成输出训练数据Y∈R<sup>n×1</sup>,并对其进行标准化处理使各个过程变量的均值为0,标准差为1,得到新数据矩阵<img file="FSA0000130783360000011.GIF" wi="60" he="54" />(3)按照如下所示步骤对数据矩阵<img file="FSA0000130783360000012.GIF" wi="46" he="49" />进行白化处理得到数据矩阵Z∈R<sup>n×M</sup>,其中,M≤m表示矩阵Z中变量个数:①计算<img file="FSA0000130783360000013.GIF" wi="46" he="50" />的协方差矩阵<img file="FSA0000130783360000014.GIF" wi="356" he="64" />其中Φ∈R<sup>m×m</sup>,上标号T表示矩阵转置;②计算矩阵Φ的所有特征值和特征向量,并剔除小于0.0001的特征值及其对应的特征向量,得到特征向量矩阵P=[p<sub>1</sub>,p<sub>2</sub>,…,p<sub>M</sub>]∈R<sup>m×M</sup>以及特征值对角矩阵D=diag(λ<sub>1</sub>,λ<sub>2</sub>,…,λ<sub>M</sub>)∈R<sup>M×M</sup>;③对<img file="FSA0000130783360000015.GIF" wi="46" he="50" />进行白化处理,得到<img file="FSA0000130783360000016.GIF" wi="423" he="54" />(4)设置保留的独立元个数d,利用白化后的输入Z与输出<img file="FSA0000130783360000017.GIF" wi="62" he="52" />选择不同的非二次函数建立起相应的MICR软测量模型,并保存各个模型参数Θ<sub>k</sub>={W<sub>k</sub>,B<sub>k</sub>}以备用,其中,k=1,2,3分别为三种非二次函数的标号,W<sub>k</sub>∈R<sup>d×m</sup>与B<sub>k</sub>∈R<sup>d×1</sup>分别为第k个MICR模型的分离矩阵和回归系数矩阵;(5)利用各个MICR模型对输入训练数据<img file="FSA0000130783360000018.GIF" wi="46" he="50" />进行预测得到相应的预测值<img file="FSA0000130783360000019.GIF" wi="171" he="77" />其中,<maths num="0001"><math><![CDATA[<mrow><msub><mi>Z</mi><mi>k</mi></msub><mo>=</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover><msup><msub><mi>W</mi><mi>k</mi></msub><mi>T</mi></msup><msub><mi>B</mi><mi>k</mi></msub><mo>;</mo></mrow>]]></math><img file="FSA00001307833600000110.GIF" wi="302" he="63" /></maths>(6)采用最小二乘回归方法得到各个MICR模型预测值的权重比<img file="FSA00001307833600000111.GIF" wi="167" he="77" />使平方预测误差<img file="FSA00001307833600000112.GIF" wi="683" he="87" />最小化;(7)收集新的过程容易测量的数据x∈R<sup>m×1</sup>,并对其进行标准化处理得到<img file="FSA00001307833600000113.GIF" wi="190" he="54" />(8)利用各个MICR模型参数分别对<img file="FSA00001307833600000114.GIF" wi="36" he="41" />进行预测得到相应的预测值<img file="FSA00001307833600000115.GIF" wi="162" he="77" />其中,<img file="FSA00001307833600000116.GIF" wi="176" he="63" />为第k个MICR软测量模型预测值;(9)计算对应于当前输入数据<img file="FSA00001307833600000117.GIF" wi="35" he="40" />的预测输出值<img file="FSA00001307833600000118.GIF" wi="477" he="54" />
地址 315211 浙江省宁波市江北区风华路818号