删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于二阶段自适应多模型的聚合釜温度控制

本站小编 Free考研考试/2020-04-15

王振雷 1 , 毛福兴 1 , 王昕 2
1.华东理工大学 化工过程先进控制和优化技术教育部重点实验室, 上海 200237;
2.上海交通大学 电工与电子技术中心, 上海 200240

收稿日期: 2015-08-25
基金项目: 国家自然科学基金面上基金资助项目(21376077);国家自然科学基金优秀青年基金资助项目(61222303);上海市自然科学基金资助项目(14ZR1410000,14ZR1421800);上海市重点学科建设基金资助项目(B504);流程工业综合自动化国家重点实验室开放课题基金资助项目(PALN201404)
作者简介: 王振雷(1975-),男,教授。E-mail:wangzhen_l@ecust.edu.cn


摘要:针对丙烯腈聚合过程的强时滞和较大参数不确定等特性,该文提出一种基于二阶段自适应多模型的广义预测控制方法。该方法首先根据系统的参数范围,建立多个自适应模型,应用最小二乘算法分别进行参数估计。再利用各自适应模型的参数估计值和预报误差计算模型的权值,将各参数估计值加权求和得到最终参数估计值。将该参数估计值作为参数的真值,利用广义预测控制算法确定各时刻的控制作用。仿真结果显示:该方法能使系统未知参数快速收敛到真值,同时系统的动态性能和对理想温度的跟踪精度较常规多模型自适应控制有明显的提高。
关键词: 多模型 二阶段自适应 广义预测控制 聚合釜 自适应控制
Temperature control of an acrylonitrile polymerization kettle using multiple models with second level adaptation
WANG Zhenlei1, MAO Fuxing1, WANG Xin2
1.Key Laboratory of Advanced Control and Optimization for Chemical Processes of the Ministry of Education, East China University of Science and Technology, Shanghai 200237, China;
2.Center of Electrical and Electronic Technology, Shanghai Jiao Tong University, Shanghai 200240, China


Abstract:A generalized predictive control method was developed from multiple models with second level adaptation for temperature control the acrylonitrile polymerization process which has long time delays and large parameter uncertainties. Several adaptive models are designed for the system parameter ranges with the parameters estimated by a recursive least squares algorithm. Then, the model weights are calculated based on the parameter estimates and the prediction error of each model. Then, the parameter estimates are used as the true values of the parameters to determine the control action via the generalized predictive control algorithm. Simulation results show that this method enables a system with unknown parameters to quickly converge to the true value. The system performance of the system and the tracking accuracy of the ideal temperature are significantly improved compared with conventional multiple model adaptive control.
Key words: multiple modelssecond level adaptationgeneralized predictive controlpolymerization kettleadaptive control
腈纶是一种具有广阔发展前景的人造纤维,素有“人造羊毛”之称,其主要成分为聚丙烯腈,由于聚丙烯腈大分子结构的特点,腈纶纤维具有良好的弹性、强度、耐热等性能[1]。丙烯腈聚合过程是腈纶生产工艺中的关键工序,丙烯腈聚合反应受到温度、聚合物浓度、反应液pH值等的影响,其中温度是至关重要的因素,直接影响产品的质量和产量,甚至工艺的正常安全生产[2]。温度过低,分子运动减弱,反应慢,各单体转化率降低;温度过高,聚合物的分子量下降,使纤维质量下降,同时也加剧了聚合釜的腐蚀。因此,对聚合釜内温度进行测量与控制以达到工艺要求,是丙烯腈聚合过程最重要的控制问题之一。
早期的聚合釜基本上采用比例-积分-微分(PID)控制方案[3],由于聚合釜是一个典型的具有惯性、时滞、非线性、时变等特点复杂控制对象,常规的PID控制难以适应,控制过程响应慢、超调严重、稳定性差。为了改善系统的控制性能,许多更加先进有效的控制算法应用到聚合釜温度控制中,如广义预测控制[4]、模糊PID[5]、基于模糊神经模型的预测控制[6]等,这些方法较常规PID控制在一定程度上提高了控制效果。这些方法主要是降低了经典PID对模型精度的要求,对聚合釜温度系统的强时滞、非线性、慢时变等进行了分析与处理。但这些方法中,没有考虑在实际工作环境下聚合釜系统不可避免地存在的扰动干扰和较大的参数不确定。
二阶段自适应多模型(multiple models with second level adaptation,MMSLA)是在经典多模型自适应控制(multiple models adaptive control,MMAC)[7-8]基础之上进一步拓展的一类多模型自适应控制算法,主要用于处理存在较大参数不确定或时变系统,最早由Zhuo等[9]和Narendra等[10]提出,其核心思想是采用多个自适应模型描述参数未知的被控对象,利用多个自适应模型进行参数估计,确定控制器的参数、控制作用。作为初期的研究,目前对该方法的研究[9, 11-13]仅仅针对参数未知、连续确定性系统。在研究中往往需要假设系统状态可测,且主要处理线性定常或参数缓慢跳变对象,对参数持续变化和扰动噪声缺乏处理,在实际工程中难以应用。
为此,本文提出了一种基于受控自回归积分滑动平均(controlled auto-regressive integrated moving average,CARIMA)模型的多模型二阶段自适应方法。该方法用一类随机噪声描述在实际工作环境下丙烯腈聚合反应过程中存在的不可避免的一些扰动。将系统的机理模型经线性化、离散化以后用CARIMA模型描述。在任意时刻,直接根据对象输入输出及噪声数据,通过上述MMSLA过程在线辨识获得对象参数的估计值,再利用广义预测控制(generalized predictive control,GPC)算法确定控制作用,从而实现MMSLA与GPC的结合,并将其应用到丙烯腈聚合釜温度控制。
1 对象模型聚丙烯腈纤维是一种高分子长链合成聚合物形成的人造纤维。日常生活中的腈纶就是用丙烯腈与其他单体合成的共聚物,再经过纺丝等一系列后续工艺制作成的合成纤维,这些工艺中最关键的就是丙烯腈的聚合过程。丙烯腈聚合反应[14]是以丙烯腈CH2=CH-CN (acrylonitrile,AN)为第一单体,醋酸乙烯酯CH3COOCH=CH2 (vinyl acetate,VA)为第二单体,采用水溶性NaClO3·NaHSO3氧化还原体系为诱导剂,巯基乙醇(β-ME)为链转移剂,在聚合釜内40 ℃~60 ℃的温度条件下进行,生产工艺流程如图 1所示。反应物按生产比例配置混合后通过3条管道连续输送到聚合釜,单体混合物主要为89.05%的丙烯腈和10.95%的醋酸乙烯;NaClO3混合液主要为2.0%的NaClO3,少量CuSO4的水溶液;第三条管道主要为10.65%的NaHSO3,3.72%的NaNO3以及(β-ME)的水溶液,同时还配有含0.7%的NaOH的脱盐水,每条管道都有精确的流量计检测流量大小[15]。在该过程中对温度的控制至关重要,直接影响到产品的质量和产量。
图 1 丙烯腈聚合反应工艺流程示意图丙烯腈聚合反应是一个放热反应过程,在反应中,温度升高,聚合反应的速度加快,同时反应放出的热量也增大,使反应釜系统的温度持续升高,因此,就釜内温度变化来说,该过程具有正反馈性质,必须采取有效的措施,吸收反应所放的能量[15],将聚合釜温度维持在理想温度环境下。
图 1 丙烯腈聚合反应工艺流程示意图
图选项





聚合釜内温度系统模型,按热量传递情况,可分为2个部分,即聚合釜内物料温度及热量变化模型和聚合釜夹套内冷却水温度及热量变化模型[15]。聚合釜内物料温度及热量变化主要包括由反应的焓变、反应物料吸收的热量及聚合釜与夹套间交换的热量,其微分方程如式(1)所示,即[16]
$\begin{align} & V\rho {{c}_{p}}\frac{dT}{dt}=KMV\text{-}\Delta H+ \\ & {{F}_{f}}\rho {{c}_{p}}{{T}_{f}}\text{-}T+UA{{T}_{j}}\text{-}T. \\ \end{align}$ (1)
聚合釜夹套水能量变化主要由冷却水吸收的热量和与聚合釜间能量交换引起,有式(2)所示微分方程,即
$\begin{align} & {{V}_{j}}{{c}_{pj}}{{\rho }_{j}}\frac{d{{T}_{j}}}{dt}= \\ & {{F}_{j}}{{\rho }_{j}}{{c}_{pj}}{{T}_{jf}}\text{-}{{T}_{j}}+UA{{T}_{j}}\text{-}T.~ \\ \end{align}$ (2)
单体浓度变化微分方程如式(3)所示,即
$\begin{align} & \frac{dM}{dt}=\frac{{{F}_{f}}}{V}{{M}_{f}}-\frac{F}{V}M- \\ & 2f{{k}_{d}}IM-{{k}_{d}}MM.~ \\ \end{align}$ (3)
引发剂浓度变化的微分方程式(4)所示,即
$~\frac{dI}{dt}=\frac{{{V}_{f}}}{V}{{I}_{\text{f}}}-\frac{F}{V}I-{{k}_{d}}I.$ (4)
聚合物浓度变化方程式(5)所示,即
$\frac{dP}{dt}~=\frac{{{F}_{f}}}{V}{{I}_{f}}-\frac{F}{V}I-{{k}_{t}}MM.$ (5)
其中:T为釜内物料平均温度;Tf聚合釜进料温度;Tj为夹套冷却水平均温度;Tjf为夹套冷却水入口温度;V为釜内反应液体积;Vj为夹套体积;ΔH反应焓变;ρ反应釜内混合物平均密度;ρj为冷却水密度;cp反应混合物平均比热容;cpj冷却水比热容;UA热交换系数;kd,kp,kt分别为链引发、链增长、链终止速率常数。kd,kp,kt有如式(6)所示关系,即
$\begin{align} & ~{{k}_{d}}={{A}_{d}}\cdot exp-{{E}_{d}}RT, \\ & {{k}_{p}}={{A}_{p}}\cdot exp-{{E}_{p}}RT, \\ & {{k}_{d}}={{A}_{t}}\cdot exp-{{E}_{t}}RT. \\ \end{align}$ (6)
其中:Ad,Ap,At分别为频率因子;Ed,Ep,Et分别为链引发、链增长、链终止速率活化能;R为摩尔气体常数。
通过进一步简化、线性化处理,并加入随机噪声信号模拟聚合釜的进料等扰动,因此,丙烯腈聚合过程的温度动态特性离散化以后可用式(7)所示CARIMA模型描述,即[4]
$\begin{align} & A{{z}^{-1}}yk= \\ & {{z}^{-d}}B{{z}^{-1}}uk+\frac{C{{z}^{-1}}vk}{\Delta }. \\ \end{align}$ (7)
其中:z-1是后移算子;Δ=1-z-1为差分算子;ykuk分别是系统的输出和输入,即聚合釜内反应混合物的平均温度T和聚合釜夹套冷却水的体积流量Fjd是系统的时滞;vk为均值为零的不相关随机序列,表示系统中随机扰动的影响;A,B,Cz-1的多项式,其中AC是首一多项式,BHurwitz多项式。A,B,C的阶次分别为na,nb,nc,已知;且Az-1,Bz-1Cz-1互质,有如式(8)所示形式,即
$\begin{align} & ~A{{z}^{-1}}=1+{{a}_{1}}{{z}^{-1}}+\ldots +{{a}_{{{n}_{a}}}}{{z}^{-{{n}_{a}}}}, \\ & B{{z}^{-1}}={{b}_{0}}+{{b}_{1}}{{z}^{-1}}+\ldots +{{b}_{{{n}_{b}}}}{{z}^{-{{n}_{b}}}}, \\ & C{{z}^{-1}}=1+{{c}_{1}}{{z}^{-1}}+\ldots +{{c}_{{{n}_{c}}}}{{z}^{-{{n}_{c}}}}. \\ \end{align}$ (8)
将式(7)写成如式(9)所示形式,即
$\begin{align} & \Delta A{{z}^{-1}}yk= \\ & B{{z}^{-1}}\Delta uk-d+C{{z}^{-1}}vk. \\ \end{align}$ (9)
并令:
$\bar{A}z=1-zAz=1+\bar{a}z+\ldots +\bar{a}z.~$
等价描述为:
$\begin{align} & y\left( k \right)=-{{{\bar{a}}}_{1}}yk-1-\ldots -{{{\bar{a}}}_{{{n}_{a}}}}yk-{{n}_{{\bar{a}}}}_{_{n}}-1+ \\ & {{b}_{0}}\Delta uk-d+\ldots +{{b}_{{{n}_{b}}}}\Delta uk-{{n}_{b}}-d+ \\ & {{c}_{1}}vk-1+\ldots +{{c}_{{{n}_{c}}}}vk-{{n}_{c}}+vk.\text{ } \\ \end{align}$ (10)
其中:${{{\bar{a}}}_{i}}={{a}_{i}}-{{a}_{i-1}},1\le i\le {{n}_{a}},{{{\bar{a}}}_{{{n}_{a}}+1}}=-{{a}_{{{n}_{a}}}}$
2 多模型集的构建由前述丙烯腈聚合反应过程的描述中可知,由于聚合釜的非线性和复杂性等,使得所建立的数学模型与实际对象有较大的模型摄动。另一方面,由于单体丙烯腈对人体有害,聚合釜是一个密封装置,不容易对反应过程中实时数据进行采集。因此,在对聚合釜温度进行控制时,使用MMSLA,其多模型集的构建分两步进行。
2.1 第一阶段MMSLA的第一步是建立含多个参数可调的自适应模型的多模型集。如前所述,经典MMAC中,即基于切换和基于切换与调节的多模型。切换机制[7]中,多模型集是多个固定模型;切换与调节[12, 17]机制其模型集中既有固定模型又有参数可调的自适应模型,M个固定模型分布在参数不确定空间,自适应模型包含1个可重新赋初值和1个自由运行的自适应模型。在这2种方法中,模型数量与参数向量的维数呈指数增长关系[12]。MMSLA的多模型集由多个自适应模型组成,各个模型的结构完全相同,仅仅是参数的初值不同,参数初值分布在未知参数的确定空间边界,在选择初值时,必须保证对象的真实参数位于各个自适应模型初值所组成的凸壳中[9]。将各自适应模型按式(10)写成如式(11)所示向量形式,即
$yk={{\Phi }^{Tk}}-1\theta +vk.$ (11)
其中:ΦkRn是(k-1)时刻的回归向量,由该时刻的输入输出以及噪声数据组成,即
$\begin{align} & {{\Phi }^{T}}k-1=[-yk-1,\ldots ,-yk-{{n}_{\bar{a}n}}-1, \\ & \Delta uk-d,\ldots ,~\Delta uk-{{n}_{b}}-d, \\ & vk-1,\ldots ,vk-{{n}_{c}}], \\ & \theta =[{{{\bar{a}}}_{1}},\ldots ,{{{\bar{a}}}_{{{n}_{{\bar{a}}}},{{b}_{0}},\ldots ,{{b}_{{{n}_{b}}}},{{c}_{1}},\ldots ,{{c}_{{{n}_{c}}}}}}]. \\ \end{align}$
其中:θRn,是所要估计的对象系统未知参数向量;ykR,是k时刻的系统输出。根据参数估计值${{{\hat{\theta }}}_{p}}$,可得到相应的估计输出$\hat{y}k$,即
${{{\hat{y}}}^{k}}={{\phi }^{T}}k-1\hat{\theta }k-1+vk.$ (12)
设系统参数的估计误差$\tilde{\theta }k\text{-}1={{\theta }_{p}}-\hat{\theta }k-1$,则可得系统的误差e(k),即
$ek={{\phi }^{T}}k-1\tilde{\theta }k-1.$ (13)
同时有,$ek=yk-\hat{y}k$。参数估计应用工程上较为适用的带遗忘因子的递推增广最小二乘算法(RELS),在估计对象模型时,同时估计噪声模型,即
$\begin{align} & \tilde{\theta }k=\hat{\theta }k-1+ \\ & Kkyk-{{\phi }^{T}}k\hat{\theta }k-1, \\ & Kk=Pk-1\phi k{{\phi }^{T}}kPk\text{0}1. \\ & \phi k+{{\mu }^{-1}}, \\ & Pk=\frac{1}{\mu }I-Kk{{\phi }^{T}}kPk-1. \\ \end{align}$ (14)
其中:0<μ<1为遗忘因子;P0=aIa为足够大常数。
在本文中丙烯腈聚合釜温度控制系统的MMLSA多模型集由N=n+1自适应模型组成[12],且每个模型的自适应律均为式(14)所述的RELS算法,各模型的初值θi0分布在对象未知参数的确定参数空间边界,并确保θp0θi0构成的凸壳中。
2.2 第二阶段MMSLA的核心思想是[10, 12]:如果任一时刻下各模型的自适应参数值是对象参数θp的估计,则这些模型的一个凸组合也是对象参数θp的一个估计,即有
${{\hat{\theta }}_{p}}=\sum\limits_{i=1}^{N}{{{a}_{i}}{{{\hat{\theta }}}_{i}}k}.$ (15)
其中:$\sum\limits_{i=1}^{N}{{{a}_{i}}=1}$,且αi≥0同时Han[6]有定理:只要对象参数θp在初始时刻k=0位于各自适应模型的初始条件θi0所构成的凸壳中,并且各模型参数自适应律为式(14),则在任意时刻k>0,θp位于θik所构成的凸壳中,即可保证参数能收敛到真值[18]
MMSLA的第二阶段主要是各模型的权系数估计,由第一阶段得到各个自适应模型的估计参数;再根据各个自适应模型的输出预测误差对其相应系数αi进行估计,再按式(15)构成的模型参数的估计值。
式(15)在k时刻有矩阵形式
$\begin{align} & {{{\hat{\theta }}}_{1}}k{{{\hat{\theta }}}_{2}}k\cdots {{{\hat{\theta }}}_{N}}k\text{a}= \\ & \Theta k\text{a}=\hat{\theta }p. \\ \end{align}$ (16)
其中:α=α1α2…αNT∈RN; Θk∈Rn×N,各列θik是第i个自适应模型对的参数估计。同时在(k-1)时刻有
$\Theta k-1a={{{\hat{\theta }}}_{p}}$ (17)
将式(16)与(17)相减,可得
$\Delta \Theta k\alpha =\Theta k-\Theta k-1\alpha =0.$ (18)
矩阵Θk-Θk-1中每一列为${{{\hat{\theta }}}_{i}}k-{{{\hat{\theta }}}_{i}}k-1$,由参数自适应律式(14),有
${{{\hat{\theta }}}_{i}}k-{{{\hat{\theta }}}_{i}}k-1=-\frac{Pk-1\phi k}{{{\phi }^{T}}kPk-1\phi k+\mu }eik.$ (19)
将矩阵Θk-Θk-1记为Mk,则有
$Mka=0.$ (20)
同时考虑αi的约束$\sum\limits_{i=1}^{N}{{{a}_{i}}=1}$,则有
$Mka=0.$ (21)
其中:$~Wk=_{{{l}^{T}}}^{Mk}\in {{R}^{N\times N}};\text{ }{{l}^{T}}=\left[ 1\text{ }1\ldots 1 \right]\in {{R}^{N}};b=0\text{ }0\ldots 1{{~}^{T}}\in {{R}^{N}}$。则可得αi的自适应律为[12]
${{{\hat{a}}}^{k}}+1={{{\hat{a}}}^{k}}-\frac{{{\mu }_{opt}}{{W}^{T}}k\tilde{b}~k}{1+{{\left\| Wk \right\|}^{2}}}.$ (22)
其中:$\tilde{b}~k=\tilde{b}k-b,\tilde{b}k=Wk\hat{a}k;\hat{b}k$$\hat{a}k$分别是bak时刻的估计值;μopt是自适应增益。
因此,在任时刻下,α用其估计值$\hat{a}k$代替,并有任意时刻对象的参数估计值为
${{{\hat{\theta }}}_{p}}k=\sum\limits_{i=1}^{N}{{{{\hat{a}}}_{i}}k{{\theta }_{i}}k.}$ (23)
即为对象真实参数估计。
3 控制器设计在本文中,丙烯腈聚合釜的控制器根据MMSLA获得对象参数的估计值,以GPC算法确定控制作用。GPC是Clarke等在保持最小方差自校正控制的模型预测、最小方差控制、在线辨识等原理的基础上,汲取了动态矩阵控制、模型算法控制中的多步预测和滚动优化策略,提出来的一种经典的模型预测控制算法,已经成功应用于一些模型精度不高、时滞等复杂工业过程控制中[19]
GPC中的对象如式(7)描述的CARIMA模型,并且在k时刻的优化性能指标有如式(24)所示形式[20],即
$\begin{align} & minJk=E\sum\limits_{i=1}^{N}{yk+j-{{\omega }_{r}}k+{{j}^{2}}+} \\ & \sum\limits_{i=1}^{N}{\lambda j\Delta uk+j-{{1}^{2}}}. \\ \end{align}$ (24)
其中:E(·)表示求数学期望;j为预测步长;ωr(k)为对象的输出参考值;N1,N2分别是优化时域的起始值和终止值,一般N1=1,N2=N即优化时域;NU为控制时域;λj为控制加权系数,为了简便一般为常值λ
为了使上述性能指标最优,需要确定j步后的输出预测值yk+j,并引入如下Diophantine方程
$\begin{align} & C{{z}^{-1}}=\Delta A{{z}^{-1}}{{E}_{j}}{{z}^{-1}}+{{z}^{-d}}{{F}_{j}}{{z}^{-1}}, \\ & {{G}_{j}}{{z}^{-1}}={{E}_{j}}{{z}^{-1}}B{{z}^{-1}}. \\ \end{align}$ (25)
其中:ΔAz-1Bz-1Cz-1系数由前述MMSLA估计得到,且有
$\begin{align} & {{E}_{j}}{{z}^{-1}}=1+\sum\limits_{i=1}^{j=1}{{{e}_{j,i{{z}^{-i}}}}}, \\ & {{F}_{j}}{{z}^{-1}}=1+\sum\limits_{i=1}^{{{n}_{a}}-1}{{{f}_{j,i}}{{z}^{-i}}}, \\ & {{G}_{j}}{{z}^{-1}}=1+\sum\limits_{i=1}^{{{n}_{b}}d-1}{{{g}_{j,i}}{{z}^{-i}}}.~ \\ \end{align}$
最优预测值可由式(26)得到即
$\begin{align} & \hat{y}k+jk= \\ & {{G}_{j}}{{z}^{-1}}\Delta uk+j-1+{{F}_{j}}{{z}^{-1}}yk. \\ \end{align}$ (26)
并有最终控制量u(k)为
$\begin{array}{*{35}{l}} ~uk=uk-1+1\text{ }0\ldots 0\cdot \\ {{G}^{T}}G+\lambda {{I}^{-1}}{{G}^{T}}{{y}_{r}}-f.~ \\\end{array}$ (27)
其中:
$\begin{align} & G=\begin{matrix} {{g}_{1}} & \ldots & 0 \\ \vdots & \ddots & \vdots \\ {{g}_{N}} & \ldots & {{g}_{N-NU+1}}_{N\times NU} \\\end{matrix}, \\ & {{y}_{r}}=yk+1,\ldots ,yk+{{N}^{T}}, \\ & f={{f}_{1}}k+1,\ldots ,{{f}_{N}}k+{{N}^{T}},{{f}_{i}}k+i={{z}^{i-1}}[{{G}_{i}}{{z}^{-1}}-{{z}^{-i-1}}{{g}_{i,i-1}}- \\ & \ldots -{{g}_{i,0}}]\Delta uk+{{F}_{i}}yk. \\ \end{align}$ (28)
由上述内容可得该聚合釜温度系统的机构框图如图 2所示。
图 2 MMSLA控制结构
图选项





4 仿真分析由前面所述丙烯腈聚合过程聚合釜温度系统机理分析可知,在丙烯腈聚合过程中,聚合釜温度系统的数学模型描述如式(29)所示,即
$\begin{align} & ~yk+{{a}_{1}}yk-1+{{a}_{2}}yk-2={{b}_{1}}uk-1+ \\ & {{b}_{2}}uk-2+\frac{vk+{{c}_{1}}vk-1+{{c}_{2}}vk-2}{\Delta }\text{ }. \\ \end{align}$ (29)
其中:由参数a1,a2,b1,b2,c1c2组成的参数向量为θ,各参数值未知;系统的时滞d=1;v(k)表示Gauss白噪声序列,在仿真中由Matlab函数randn生成。在仿真中,分别设计单个自适应模型的GPC控制器,7个模型的经典MMAC的GPC控制器以及7个自适应模型MMSLA的GPC控制器进行仿真。
根据前述机理模型式(1)—式(5)及文[15],得到理想环境下参数,以周期为120 s采样得到离散化后的值为
$\begin{align} & \theta _{p}^{*}= \\ & -0.97\text{ }0.04-0.036\text{ }-0.013\text{ }0.5\text{ }{{0.4}^{T}}. \\ \end{align}$ (30)
考虑在实际工业生产条件下,存在较大扰动,如反应物混合物的密度ρ发生变ρmax=950,ρmin=900; 进料流量 F,流量F存在扰动可能由于流量泵的输入电流不稳定会发生Fmax=26,Fmin=24[15]。因而参数具有不确定性,且各参数的范围为:a1∈-3 -0.5,a2∈0 0.2,b1∈-0.5 0,b2∈-0.5 0,c1∈0 -2以及c2∈0 1。
并由此设计7个参数初值各异,而模型结构完全一致的自适应模型构成该系统的MMSLA多模型集,其参数初值选择具体情况如式(31)所示,即
$\begin{align} & {{\Sigma }_{1}}:{{\theta }_{1}}0=-1.5,0.1,-0.25,-0.25,1.0,{{0.5}^{T}}, \\ & {{\Sigma }_{2}}:{{\theta }_{2}}0=-3,0.001,-0.01,-0.01,0.001,{{0.01}^{T}}, \\ & {{\Sigma }_{3}}:{{\theta }_{3}}0=-0.5,0.2,-0.5,-0.5,2.0,{{0.5}^{T}}, \\ & {{\Sigma }_{4}}:{{\theta }_{4}}0=-3.0,0.001,-0.5,-0.001,2.0,{{0.001}^{T}}, \\ & {{\Sigma }_{5}}:{{\theta }_{5}}0=-0.5,0.2,-0.001,-0.5,1.0,{{0.001}^{T}}, \\ & {{\Sigma }_{6}}:{{\theta }_{6}}0=-3.0,0.001,-0.5,-0.001,0.001,{{1.0}^{T}}, \\ & {{\Sigma }_{7}}:{{\theta }_{7}}0=-0.5,0.2,-0.001,-0.001,2.0,{{1.0}^{T}}. \\ \end{align}$ (31)
系数向量α的初始值为:α0=[0.4,0.1,0.1,0.1,0.1,0.1,0.1]T
同时,取上述参数初值,设计文[12]所述经典的切换与调节MMAC的多模型集,其中模型Σ1为自由运行自适应模型,模型Σ2为可赋初值的自适应模型,模型Σ36为固定模型,在任意时刻,有如式(32)所示性能指标[20],即
$\begin{align} & {{J}_{i}}k=\gamma {{e}_{i}}^{2}k+\eta \sum\limits_{j=1}^{{{L}_{1}}}{{{\rho }^{j}}{{e}_{i}}^{2}k-j}, \\ & i=1,2,\ldots ,m+2. \\ \end{align}$ (32)
其中:γ,η,ρ均为误差权系数,分别取值为0.9、0.1和0.98;L1为数据长度。任意时刻如果Ji<J2,则θi赋给θ2,并将θ2作为实际参数值,通过与MMSLA一样的GPC算法确定控制作用,进行仿真。
此外,单个自适应模型的参数初值可任取,在仿真中取:θ00=[0.001,0.001,0.001,0.001,0.001,0.001]T其他参数,三者的取值完全一致,分别为:控制步长L=400;优化时域N=N2=10;控制时域NU=4;柔化因子a=0.4。
各自适应模型均采用带遗忘因子的ELS进行辨识估计,遗忘因子μ=0.98。
分别在较小扰动环境(噪声方差σ2=0.01)和较大扰动环境(噪声方差σ2=1)下进行仿真。当 σ2=0.01时,通过仿真,得到的输出量的变化轨迹和参数辨识结果分别如图 3图 4所示。在图 3图 4中,y(k)表示系统输出,ωr(k)为幅值为60的参考信号,即为聚合釜内混合反应物的理想平均温度。
图 3 噪声方差为σ2=0.01系统仿真的输出
图选项





图 4 噪声方差为σ2=0.01时参数辨识
图选项





图 34可知,当扰动较小时,无论是使用单个自适应模型,还是多个自适应模型,系统都能跟随参考值。当使用MMSLA时,跟踪效果最好,稳定时输出更平稳;使用经典的MMAC时,效果次之,但明显都较单个自适应模型好。参数辨识的效果更明显,当使用MMSLA辨识时,参数能收敛到真值,且收敛速度快,参数收敛平滑;使用经典MMAC时,由于切换,参数动态过程波动大,且有跳变,但参数收敛速度较单个自适应模型快。
当聚合釜工作环境扰动较大时(σ2=1),由仿真得到的系统输出响应以及参数辨识如图 5图 6所示。由仿真结果可知,当扰动噪声较大时,使用MMSLA和经典MMAC控制时,系统都能较好地跟踪参考值,但MMSLA效果更好,由于存在较大噪声扰动,稳定时输出仍不可避免地存在较小范围的波动;而使用单个自适应模型时,系统不能满足较高精度的跟踪,跟踪误差较大,且频繁波动。
图 5 噪声方差为σ2=1时系统的仿真输出
图选项





图 6 噪声方差为σ2=1时参数辨识
图选项





5 结 论本文综合分析了丙烯腈聚合反应过程中聚合釜温度控制系统的时滞、非线性、参数较大不确定等复杂特性。并针对该温度控制系统,设计了基于二阶段自适应多模型的广义预测控制器。通过建立多个自适应模型在线、实时地对系统参数进行估计,通过协调利用各个模型提供的系统信息进行辨识,提高了模型未知参数收敛到系统参数真实值的速度及系统的动态性能。通过仿真显示,当使用相同数量的模型,相同的进料流量波动等扰动下使用MMSLA时,较使用经典的MMAC时系统未知参数收敛速度显著提高,同时输出温度对丙烯腈聚合的理想温度的跟踪精度明显提高,且都明显优于单个模型的自适应控制。

参考文献
[1] Journal of Central South University(Science and Technology), 41(2):649-654.--> Chan K, Bui T N N, Woo Y Y. Fabrications and electrochemical properties of two-phase activated carbon nanofibers from electrospinning[C]//Proceedings of the IEEE Nanotechnology Materials and Devices Conference. Gyeongju, South Korea:IEEE Nanotechnology Council, 2005:412-413.
[2] Journal of Central South University(Science and Technology), 41(2):649-654.-->丁东湖, 金建祥. 聚氯乙烯温度控制方案的设计及应用[J]. 化工自动化及仪表, 2004, 31(2) : 15–16.DING Donghu, JIN Jianxiang. A study on application and control strategies of PVC temperature system[J]. Control and Instruments in Chemical Industry, 2004, 31(2) : 15–16.(in Chinese)
[3] Journal of Central South University(Science and Technology), 41(2):649-654.-->Nemtoi G, Onu N A. Thermal effect of some polymerization and polycondensation reactions[J]. Ploymeric Ammonium Quaternary Salts. European Ploymer Journal, 2000, 36(12) : 2679–2680.
[4] Journal of Central South University(Science and Technology), 41(2):649-654.--> YE Yanfei, WANG Bailin, SHAO Mingheng, et al. Improved generalized predictive control in polymerization process[C]//Proceedings of the Second International Conference on Mechanic Automation and Control Engineering. Huhehaote:Inner Mongolia University of Technology, 2011:1354-1357.
[5] Journal of Central South University(Science and Technology), 41(2):649-654.--> SUN Zhou, GUO Liji. Modeling and fuzzy controller design for ACR Polymerizing-Kettles[C]//Proceedings of the 8th IEEE International Conference on Control and Automation. Xiamen:IEEE Control System Society, 2010:1688-1693.
[6] Journal of Central South University(Science and Technology), 41(2):649-654.--> Jalili-Kharaajoo M. Predictive control of a continuous stirred tank reactor based on neuro-fuzzy model of the process[C]//Proceedings of the SICE Annual Conference. Fukui, Japan:The Society of Instrument and Control Engineers, 2003:3277-3282.
[7] Journal of Central South University(Science and Technology), 41(2):649-654.-->Morse A S. Supervisory control of families of linear set-point controllers:Part I. Exact matching[J]. IEEE Transactions on Automatic Control, 1996, 41(10) : 1413–1431.DOI:10.1109/9.539424
[8] Journal of Central South University(Science and Technology), 41(2):649-654.-->Narendra K, Balakrishnan S. Adaptive control using multiple models[J]. IEEE Transactions on Automatic Control, 1997, 42(2) : 171–187.DOI:10.1109/9.554398
[9] Journal of Central South University(Science and Technology), 41(2):649-654.-->ZHUO Han, Narendra K S. New concepts in adaptive control using multiple models[J]. IEEE Transactions on Automatic Control, 2012, 42(1) : 171–187.
[10] Journal of Central South University(Science and Technology), 41(2):649-654.--> Narendra K S, ZHUO Han. Discrete-time adaptive control using multiple models[C]//Proceedings of the American Control Conference. San Francisco, USA:American Automatic Control Council, 2011:2921-2926.
[11] Journal of Central South University(Science and Technology), 41(2):649-654.--> Pandey V K, Ka I, Mahanta C. Multiple models and second level adaptation for a class of nonlinear systems with nonlinear parameterization[C]//Proceedings of the 9th International Conference on Industrial and Information Systems. Gwalior, India:ABV-Indian Institute of Information Technology and Management, 2014:1-6.
[12] Journal of Central South University(Science and Technology), 41(2):649-654.-->Narendra K S, ZHUO Han. The changing face of adaptive control:The use of multiple models[J]. Annual Review in Control, 2011, 35(1) : 1–12.DOI:10.1016/j.arcontrol.2011.03.010
[13] Journal of Central South University(Science and Technology), 41(2):649-654.--> ZHUO Han, Narendra K S. Second level adaptation using multiple models[C]//Proceedings of the American Control Conference. San Francisco, USA:American Automatic Control Council, 2011:2350-2355.
[14] Journal of Central South University(Science and Technology), 41(2):649-654.-->唐振波. 丙烯腈水相沉淀聚合生产工艺研究[J]. 应用化工, 2011, 40(6) : 1069–1072.TANG Zhenbo. Research on productive process of precipitation polymerization of acrylonitrilein water phase[J]. Applied Chemical Industry, 2011, 40(6) : 1069–1072.(in Chinese)
[15] Journal of Central South University(Science and Technology), 41(2):649-654.-->刘忠涛, 基于LMI的丙烯腈聚合釜温度控制问题研究[D]. 哈尔滨:哈尔滨工业大学, 2010. LIU Zhongtao. Research on Control of the Acrylonitrile Polymerization Reactor Temperature Based on LMI[D]. Harbin:Harbin Institute of Technology, 2010.(in Chinese).
[16] Journal of Central South University(Science and Technology), 41(2):649-654.-->DING Yongsheng, KUANG Xiaoliang, HAO Rong, et al. An intelligent cooperative decoupling controller for coagulation bath in polyacrylonitrile carbon fiber production[J]. IEEE Transactions on Control Systems Technology, 2013, 21(2) : 467–479.DOI:10.1109/TCST.2012.2186454
[17] Journal of Central South University(Science and Technology), 41(2):649-654.-->Narendra K S, Driolle O A. Adaptive control using multiple models, switching, and tuning[J]. International Journal of Adaptive Control and Signal Processing, 2003, 17(2) : 87–102.DOI:10.1002/(ISSN)1099-1115
[18] Journal of Central South University(Science and Technology), 41(2):649-654.-->Narendra K S, CHENG Xiang. Adaptive control of discrete-time system using multiple models[J]. IEEE Transactions on Automatic Control, 2000, 45(9) : 1669–1686.DOI:10.1109/9.880617
[19] Journal of Central South University(Science and Technology), 41(2):649-654.-->席裕庚. 预测控制[M].北京: 国防工业出版社, 1993.XI Yugeng. Predictive Control[M].Beijing: National Defense Industry Press, 1993.(in Chinese)
[20] Journal of Central South University(Science and Technology), 41(2):649-654.-->李小田, 王昕, 王振雷, 等. 一种基于多模型切换的阶梯式广义预测控制算法[J]. 化工学报, 2012, 63(1) : 193–196.LI Xiaotian, WANG Xin, WANG Zhenlei, et al. A stair-like generalized predictive control algorithm based on multiple models switching[J]. CIESC Journal, 2012, 63(1) : 193–196.(in Chinese)

相关话题/控制 系统

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于最优控制理论的电动汽车机械式自动变速器换档控制
    方圣楠,宋健,宋海军,台玉琢,TRUONGSinhNguyen清华大学汽车安全与节能国家重点实验室,北京100084收稿日期:2015-11-16基金项目:清华大学汽车安全与节能国家重点实验室开放基金资助项目作者简介:方圣楠(1990—),女,博士研究生。通讯作者:宋健,教授,E-mail:daes ...
    本站小编 Free考研考试 2020-04-15
  • 感应同步器测角系统误差测试及补偿
    李海霞,张嵘,韩丰田清华大学精密仪器系,北京100084收稿日期:2015-05-22基金项目:国家自然青年科学基金项目(61203227)作者简介:李海霞(1977—),女,助研。E-mail:li-hx03@mails.tsinghua.edu.cn摘要:工程实用的高精度感应同步器测角系统因安装 ...
    本站小编 Free考研考试 2020-04-15
  • 基于重复滑模观测器的直线电机系统干扰估计
    李方,叶佩青,严乐阳,张辉清华大学机械工程系,摩擦学国家重点实验室,精密超精密制造装备与控制北京市重点实验室,北京100084收稿日期:2015-10-26基金项目:国家科技支撑项目(2015BAI03B00);北京市科技计划项目(Z1411qhdxxb(zrkxb)-56-6-626514015) ...
    本站小编 Free考研考试 2020-04-15
  • 卫星高速数传系统多码率融合LDPC编码器设计
    葛广君,殷柳国清华大学航天航空学院,清华信息科学与技术国家实验室(筹),北京100084收稿日期:2015-12-02基金项目:国家自然科学基金资助项目(61132002,91538203);清华大学自主科研计划(20121088030)作者简介:葛广君(1989-),男,博士研究生通讯作者:殷柳国 ...
    本站小编 Free考研考试 2020-04-15
  • 动态不确定因果图用于复杂系统故障诊断
    赵越1,董春玲2,张勤1,21.清华大学核能与新能源技术研究院,先进核能技术协同创新中心,先进反应堆工程与安全教育部重点实验室,北京100084;2.清华大学计算机科学与技术系,北京100084收稿日期:2015-09-17基金项目:国家自然科学基金资助项目(61273330,61402266)作者 ...
    本站小编 Free考研考试 2020-04-15
  • 请问信号与系统这项专业课的参考书目是哪一本?
    提问问题:请问信号与系统这项专业课的参考书目是哪一本?学院:电子与通信工程学院、人工智能学院提问人:13***09时间:2019-09-1910:37提问内容:请问信号与系统这项专业课的参考书目是哪一本?如何得到历年真题?谢谢回复内容:推荐书目请查看相关学院网站今年我校不再提供历年真题 ...
    本站小编 天津师范大学 2019-11-27
  • 控制科学与工程
    提问问题:控制科学与工程学院:控制科学与工程学院提问人:13***46时间:2018-09-2010:10提问内容:山东的专升本学生(2019年可7月可正常取得本科毕业证)报考可以吗,贵院有没有录取过专升本的考生。回复内容:入学前取得本科毕业证既可以 ...
    本站小编 山东大学 2019-11-26
  • 控制工程(专)
    提问问题:控制工程(专)学院:控制科学与工程学院提问人:18***97时间:2015-09-2516:01提问内容:2016年控制工程(专)招生多少个?推免生有多少?回复内容:我校2016年各专业免试生的人数在10月25日公布。各专业招生人数已经在我校网站公布。 ...
    本站小编 山东大学 2019-11-26
  • 控制理论与控制工程
    提问问题:控制理论与控制工程学院:控制科学与工程学院提问人:11***om时间:2015-09-2512:54提问内容:请问控制理论与控制工程专业的《自动控制原理》教材是哪本书?作者?出版社?回复内容:我校没有指定参考书目,请根据我校的考试大纲复习即可。 ...
    本站小编 山东大学 2019-11-26
  • 控制工程
    提问问题:控制工程学院:控制科学与工程学院提问人:dy***om时间:2015-09-2211:06提问内容:您好,贵院今年的控制工程专业硕士拟招收多少统考生?多少推免生?回复内容:控制工程专业学位一志愿招生约43人,具体接收推免生数以10月份的公示为准 ...
    本站小编 山东大学 2019-11-26