分布曲线对象的无偏模型预测范文

本站小编为你精心准备了分布曲线对象的无偏模型预测参考范文，愿这些范文能点燃您思维的火花，激发您的写作灵感。欢迎深入阅读并收藏。

《化工学报》2015年第S2期

摘要：

随着系统总体性能要求的提高，分布曲线对象的控制问题成为研究热点。分布曲线对象是通过宏观的操作手段控制微观的输出曲线，不是任意形状的目标曲线都是能无偏跟踪的。本文提出了一种分布曲线对象的无偏模型预测控制算法。该方法首先采用B样条模型进行分布曲线对象的描述，其次通过复合梯形方法对滚动优化命题进行离散化，最后采用二次规划方法进行无偏目标曲线计算，以实现无偏控制。仿真结果表明了算法的有效性。

关键词：

分布曲线对象；模型预测控制；无偏控制；优化；数学模拟

1.引言

在传统的工业控制中，被控变量通常是n维空间中的的一个点。随着控制系统性能要求的提高，面临着曲线控制的需求，被控变量从n维空间的一个点过渡到n维空间中的一条曲线。例如对于聚合产品而言，高分子链和低分子链代表了不同的特性，在传统工业过程控制中常采用平均分子量或者熔融指数等宏观指标作为被控目标[1-3]，但是这些宏观指标无法反映出完整的分子结构，导致聚合产品性能如抗冲击性能、热力学性能等达不到预期的要求。因此为了满足性能要求，提出对于反映完整分子结构信息的分子量分布曲线的控制是十分必要的。分布曲线的控制不仅仅用于聚合反应的分子量分布，还广泛用于其它化工领域，如造纸过程纸张两维质量分布控制、燃烧过程温度场分布控制等等。

分布曲线的控制在近10多年来得到了广泛的关注，WangH[4]利用神经网络建立随机分布模型，对静态系统和动态系统的建模和控制都进行了研究，随后为这种新型随机系统设计鲁棒控制[5]，自适应控制和非线性控制算法[6]，解决了随机系统的建模与控制问题，并且将所推导的控制算法成功运用在造纸的生产过程中。张春雨[7]利用B样条神经网络和线性递归神经网络相结合的方法对管式聚合反应的温度分布进行了建模，在此基础上引入了实时模型误差的策略，提出了扩展的积分平方误差控制指标，解决了管式聚合反应的温度分布控制问题。吴海燕[8]提出当分子量分布曲线采用离散正交多项式来描述时，权值和分子量分布的矩向量存在线性关系，并利用最优算法实现了分子量分布的形状控制。申珊华[9]基于灰箱模型对分子量分布矩值设计预测控制器，利用分子量分布矩值可测的功能完成反馈校正的功能。

分布曲线对象是通过宏观的操作手段控制微观的输出曲线，不是任意形状的目标曲线都是能无偏跟踪的。在现有研究成果中，无偏控制问题仍然没有很好解决。模型预测控制由于其模型预测、滚动优化、反馈校正的鲜明特点，使其成为目前公认处理复杂过程约束控制的有效算法，在石油、化工等过程领域中得到成功应用[10-13]。因此，本文在模型预测控制框架下进行分布曲线对象的无偏控制研究，首先采用B样条模型进行分布曲线对象的描述，其次通过复合梯形方法对滚动优化命题进行离散化，最后通过二次规划方法进行无偏目标曲线计算，以实现无偏控制。

2.分布曲线对象模型

本文采用如下模型进行分布曲线描述

3.面向分布曲线的模型预测控制

3.1预测模型令k时刻已经得到状态的估计值x(k|k)，因此，在M个连续的控制增量作用下，系统在未来P个时刻的权值预测如下所示：

3.2滚动优化在采样时刻k，通过最小化性能指标函数（7）确定最优控制序列，使得未来预测输出曲线尽可能接近期望的目标曲线，同时对控制作用的大小加以约束，避免控制作用变化过于剧烈。

3.3反馈校正由于模型失配、不可测干扰等因素存在，系统的实际运行可能偏离理想的优化结果，反馈校正策略必不可少。上述三步组成一次完整的校正，消除由于模型误差以及不可测扰动所带来的影响，使系统的实际运行尽可能逼近理想的优化结果。

4.无偏目标曲线计算

以宏观的控制变量控制微观的分布曲线，由于对象特性本身的原因，不是任意形状的目标曲线都是可达的，特别是在执行机构有约束的情况下。比如对于分子量分布而言，进料的流量、速率都是限定在一个范围之内，低分子比例与高分子比例不能任意变化。因此，有必要在控制分布曲线时根据对象特性对目标曲线进行分析，对不可达的目标曲线做出调整，放弃部分性能要求，以实现无偏控制。对于分布曲线对象，由公式(3)通过推导可以得到分布曲线的稳态模型。

5.仿真结果

本文针对如下的分布曲线对象设计控制器：通常在预测控制中，优化时域P，一般取近似于过程的上升时间，本文选择P=6。控制时域M的选择应当兼顾快速性和稳定性，本文选择M=3。输出偏差加权阵Q的大小反映了每个时刻对预测输出曲线逼近给定曲线的重视程度，控制增量加权阵R是防止控制增量的剧烈变化，为了简单起见这里Q和R均选取单位阵。

5.1仿真一在本次仿真中，对象无模型失配，且目标曲线为：()1.16()0.31()1.08()taryyyy123γ=D+D+D(18)通过本文提出的无偏目标曲线计算方法，该目标曲线是可达的，稳态误差较小。仿真结果如下图所示，经过8个控制周期，输出曲线已经跟踪上目标曲线。

5.2仿真二在本次仿真中，对象无模型失配，且目标曲线为通过本文提出的无偏目标曲线计算方法，该目标曲线是不可达的，稳态误差较大。选取qss(y)常值函数，计算出稳态误差最小的可达目标曲线，如下图虚线由仿真结果可知，输出分布曲线经过八个控制周期达到了新的目标曲线，为了区分不同区间的重要性，可以修改稳态误差的权值函数qss(y)。

5.3仿真三在本次仿真中，对象存在模型失配，目标曲线可达，在模型存在误差的情况下需要使用反馈校正逼近理想的优化效果。在本例中模型失配是由于模型的基函数数目少于真实对象的数目，导致动态特性与静态特性都存在失配。

6.结论

本文首先采用B样条基函数的描述方式，将曲线控制的动态过程用状态空间来描述，其次利用模型设计预测控制器，使得输出曲线尽可能逼近目标曲线，最后为了分析目标曲线的可达性，推导出模型的稳态表达式，利用二次规划进行求解，解决了分布曲线控制的无偏性问题，下一步工作是针对变形过程开展多分辨率控制。

参考文献

[1]费正顺，胡斌，叶鲁彬等．带路径约束的聚烯烃牌号切换操作优化方法[J]．化工学报，2010，61(4)：893-900．FEIZS,HUB,YELB,etal.optimalgradetransitiOnofpolymerizationprocesswithpathconstraints[J].CIESCJournal，2010，61(4)：893-900．

[2]张伟杰，梁军，费正顺等．具有不确定性的聚烯烃牌号切换过程动态优化[J]．化工学报，2011，62(10)：2797-2804.DOI:10.3969/j.issn.0438-1157.2011.10.018.ZHANGWJ,LIANGJ,FEIZS,etal.Optimalgradetransitioninpolymerizationprocessesunderuncertainty[J].CIESCJournal,2011,62(10)：2797-2804．DOI:10.3969/j.issn.0438-1157.2011.10.018.

[3]王平．淤浆法聚乙烯生产中牌号切换策略研究[D]：[硕士学位论文]．杭州：浙江大学，2006．WANGP.Optimalstrategyofgradetransitioninindustrialslurryhigh-densitypolyethyleneproductionprocess[D].Hangzhou:ZhejiangUniversity,2006.

[4]WangH．Boundeddynamicstochasticdistributionsmodellingandcontrol[M]．London：Springer—Ver—lagLtd．2000．14-20．

[5]WangH．Robustcontroloftheoutputprobabilitydensityfunctionsfordynamicstochasticsystems[C]//Proc．ofthe1998ConferenceonDecisionandControl．Tampa，Madison：Omini，1998．134～150．

[6]WANGH．Controloftheoutputprobabilitydensityfunctionsforaclassofnonlinearstochasticsystems[C]//Proc．oftheIFACWorkshoponAlgorithmsandArchitecturesforRealtimeControl．Cancun，Mexico：Elesvier，1998．152～180．

[7]张春雨．管式聚合反应温度分布建模及控制方法的研究[D]：[硕士学位论文]．北京：北京化工大学，2005．ZHANGCY.Thestudyoftemperaturedistributionmodelandcontrolintubularpolymeiuzationreaction[D].Beijing:BeijingUniversityofChemicalTechnology,2005.

[8]吴海燕．利用神经网络实现聚合反应分子量分布的建模与控制[D]：[硕士学位论文]．北京：北京化工大学，2004．WUHY.Modelingandcontrolingmolecularweightdistributionforpolymerizationprocessesvianeuralnetworks[D].Beijing:BeijingUniversityofChemicalTechnology,2004.

[9]申珊华，曹柳林，王晶.基于分布函数矩的聚合物分子量分布预测控制[J].化工学报，2013,64(12)：4379-4384，DOI:10.3969／j.issn.0438-1157.2013.12.016.SHENSH,CAOLL,WANGJ.PredictivecontrolofmolecularweightdistributioninpolymerizationreactionbasedonmomentofMWD[J].CIESCJournal,2013,64(12)：4379-4384，DOI:10.3969／j.issn.0438-1157.2013.12.016.

[10]钱积新，赵均，徐祖华.预测控制[M].北京,化学工业出版社,2007.QIANJX,ZHAOJ,XUZH.ModelPredictiveControl[M].Beijing:ChemicalIndustryPress,2007.

[11]韩恺，赵均，ZHUYucai,等.一种扰动自适应的鲁棒预测控制算法[J].化工学报，2009，60(7):1730-1738.HANK,ZHAOJ,ZHUYC,etal.ArobustMPCtechniquewithadaptivedisturbancemodel[J].CIESCJournal,2009,60(7):1730-1738.

[12]徐祖华，赵均，钱积新.基于Min-Max的预测控制鲁棒参数设计[J].化工学报，2004，55(4):613-617.XUZH,ZHAOJ,QIANJX.RobustdesignmethodofpredictivecontrllerparameterbasedonMin-Maxruler[J].CIESCJournal,2004,55(4):613-617.

[13]徐祖华，ZHUYucai，赵均,等.新一代的自适应模型预测控制器[J].化工学报，2008，59(5):1207-1215.XUZH,ZHUYC,ZHAOJ,etal.Next-generationadaptiveMPCcontroller[J].CIESCJournal,2008，59(5):1207-1215.

[14]程正兴，李水根，数值逼近与常微分方程数值解[M].西安交通大学出版社,1998.CHENZX,LISG.NumericalapproximationandNumericalSolutionofOrdinaryDifferentialEquations[M].Xi'anJiaotongUniversityPress,1998.

作者：康岳群 徐祖华 赵均 邵之江 单位：浙江大学控制科学与工程学院 工业控制国家重点实验室