[发明专利]一种基于特征空间的双核多向偏最小二乘的发酵过程质量预测方法有效

专利信息
申请号: 201410075308.2 申请日: 2014-03-03
公开(公告)号: CN103810396B 公开(公告)日: 2016-10-19
发明(设计)人: 王普;王锡昌;高学金;齐咏生 申请(专利权)人: 北京工业大学
主分类号: G06F19/00 分类号: G06F19/00
代理公司: 北京思海天达知识产权代理有限公司 11203 代理人: 刘萍
地址: 100124 *** 国省代码: 北京;11
权利要求书: 查看更多 说明书: 查看更多
摘要: 一种基于特征空间的双核多向偏最小二乘的发酵过程质量预测方法,既属于生物发酵领域,又属于信息科学领域,具体涉及在发酵过程中用于预测在线测量中难以及时取得的质量变量。具体步骤包括:获取训练数据,将数据展开并进行标准化,对标准化的测量数据和质量数据进行核空间中的特征提取,以提取后的特征空间为基进行高维投影,建立偏最小二乘模型,采集实时数据预测质量数据在高维空间中的投影,将高维空间中的质量数据处理为真实的预测值。本发明达到了降低整体模型非线性度,提高模型的预测精度的能力,对实现发酵过程的质量预测、控制和优化控制有重要意义。
搜索关键词: 一种 基于 特征 空间 多向 最小 发酵 过程 质量 预测 方法
【主权项】:
一种基于特征空间的双核多向偏最小二乘的发酵过程质量预测方法,其特征在于包含以下步骤:步骤(1):获取训练数据:采集发酵过程中可在线测量的测量变量和难以在线测量的质量变量形成训练数据,具体的一次数据采样表示形式如下:X=(x1,x2,...,xnx)    (1)Y=(y1,y2,...,yny)    (2)其中,x1,x2,...,xnx表示发酵过程中在线可实时测量的nx个测量变量,即测量点;同理,y1,y2,...,yny表示发酵过程中在线不可实时测量的ny个质量变量,经过多次采样形成数据Xp,nx,r和Yp,ny,r,其中p表示批次代数,用以区分间歇过程不同的批次,即批次数;r表示采样数;经过测量的数据组成三维矩阵数据:X3D和Y3D,矩阵规格为批次数×测量点数×采样数;步骤(2):将X3D和Y3D按照AT方法进行数据展开和标准化,记录取得的方差和均值,处理后的数据为X2D和Y2D;步骤(3):分别对X2D和Y2D进行核空间中的特征提取,依次得到特征基SX和SY;步骤(4):计算测量数据X2D和质量数据Y2D到核空间中的投影:选取高斯核函数,利用核技巧将X2D、Y2D以SX和SY为基分别投影到核空间中,得到KX=K(X2D,SX)和KY=K(Y2D,SY);步骤(5):利用得到的数据Kx和KY建立偏最小二乘预测模型:利用偏最小二乘算法建立数据KX到KY的映射,并获得PLS映射系数B,具体表现形式为:KY=KXB    (3)步骤(6):在线测量和处理数据,具体步骤为:步骤(6.1):在线采集新的测量数据Xnew,其数据格式与公式(1)中的X相同;步骤(6.2):将采集到的新数据Xnew按照建模时刻的AT处理方法进行处理和标准化;步骤(6.3):利用核技巧将Xnew以SX为基投影到核空间中,得到KXnew=K(Xnew,SX);步骤(7):将KXnew输入到已建立好的PLS模型中,利用公式(3)得到KYnew,即KYnew=KXnewB;步骤(8):将KYnew进行逆映射重构,并逆标准化,得到预测值Ynew,具体步骤为:步骤(8.1):利用SY构成新的矩阵数据,其中SY的数据格式为:其中,Vmn表示SY中的向量矩阵,一共有m个行向量,且每一个行向量的元素个数为n;其中m表示经过特征提取的特征基SY中行向量的个数,n=ny;经过计算得到的KYnew的数据格式为:KYnew=(Ky1,Ky2,...,Kym)    (5)令:L=KYnew    (6)对于L中的每一个元素Li,i为自然数且i∈[1,m],Li=Kyi;又令:<mrow><mi>Lmi</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>L</mi><mn>1</mn></msub><mo>-</mo><msub><mi>L</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><msub><mi>L</mi><mn>2</mn></msub><mo>-</mo><msub><mi>L</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>L</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>L</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>步骤(8.2):计算公式Ynew=‑0.5(VmiTVmi)‑1VmiT(Lmi‑VmisI)    (10)其中VmiT为Vmi的转置,I为所有元素为1的n维列向量;步骤(8.3):按照AT方法的数据处理方法逆向计算出Ynew在标准化之前的实际预测数据Ypre;步骤(9):若生产过程未结束,则重复步骤(6)至(8)。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京工业大学,未经北京工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201410075308.2/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top