基于A RM A X 模型自适应预测函数控制.pdf
《基于A RM A X 模型自适应预测函数控制.pdf》由会员分享,可在线阅读,更多相关《基于A RM A X 模型自适应预测函数控制.pdf(6页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、文章编号:100220411(2000)0520431206基于ARMAX模型自适应预测函数控制张泉灵王树青(浙江大学工业控制技术国家重点实验室杭州310027)摘要:本文提出了基于ARMAX模型的自适应预测函数控制,该算法的特点是占用内存少,计算速度快,并具有较强的鲁棒性.ARMAX模型参数是通过带遗忘因子的递推最小二乘算法在线辨识得到.仿真结果表明,该控制算法比P I D控制具有更好的控制品质.关键词:ARMAX模型,模型预测控制,预测函数控制,自适应控制中图分类号:TP13文献标识码:B1引言R ichalet和Culter等人在70年代后期提出的基于模型的预测控制经过20多年的发展,在
2、工业过程控制中得到了广泛的成功应用.到目前为止,已有以下几类典型的模型预测控制算法:R ichalet等人提出的基于有限脉冲响应(F IR)模型的模型预测启发控制(M PHC)1,并在PVC生产、蒸馏过程及电厂里得到了较好的应用;Culter2等人提出的基于有限阶跃响应(FSR)模型的动态矩阵控制(DM C);Garcia和M orari3提出的内模控制;Clarke4,5等人提出的基于受控自回归积分滑动平均模型(CAR I MA)的广义预测控制(GPC),GPC带有自校正机制,可在线修正模型参数;在80年代后期,R ichalet6等人提出了第三代模型预测控制算法,即预测函数控制(PFC),
3、并在工业机器人的快速高精度跟踪控制中得到了成功的应用.与其它模型预测控制不同,PFC认为控制输入的结构是确保控制系统性能的关键.由于实际生产过程的数学模型往往难以确知或者由于工作情况改变或环境变化而造成被控对象特性的改变.对于这类对象,常规控制往往难以对付,而自适应控制却能较好地解决这个问题.本文提出了基于ARMAX模型自适应预测函数控制,它结合了自适应控制可在线修改模型参数和预测函数控制快速跟踪能力的优点,仿真结果表明,该控制算法比P I D控制具有较快的跟踪能力及较强的鲁棒性.2基于ARMAX模型自适应预测函数控制PFC仍属于模型预测控制的范畴,因此它仍具有模型预测控制的三大特征:预测模型
4、;滚动优化;误差校正.预测模型用来预测过程的未来输出,PFC对模型没有特殊的要求,它可以是任意结构.本文采用的是ARMAX模型.2.1ARMAX模型被控对象用ARMAX模型来描述,如式(1)所示:第29卷第5期2000年10月 信 息 与 控 制Information and ControlVol.29,No.5Oct.,2000收稿日期:1999-09-22基金项目:浙江省重点科研项目资助A(z-1)y(k)=B(z-1)u(k)+C(z-1)(k)(1)其中:A(z-1)=1+a1z-1+a2z-2+anaz-naB(z-1)=b1z-1+b2z-2+bnbz-nbC(z-1)=1+c1z
5、-1+c2z-2+cncz-nc其中(k)是均值为零,方差为 2的白噪声.2.2预测函数控制2.2.1输出预测预测函数控制的模型预测输出ym(k)由两部分组成,一部分为模型自由响应yl(k),它仅仅依赖于过去时刻的控制量及输出量,与当前时刻及将来的控制量无关.另一部分为模型的受迫输出yf(k),它是当前时刻起加入控制作用后新增加的模型响应.与其它模型预测控制不同,预测函数控制认为控制输入的结构是确保控制性能的关键.在预测函数控制中,控制作用认为是若干个已知基函数ubn的线性组合,即:u=(k+i)=Nn=1nubn(i),i=0,1,H-1(2)其中N是基函数个数,ubn(i)为基函数在t=i
6、Ts时的值,Ts为采样周期,H为优化时域的长度,n为线性组合系数.其中基函数的选择取决于设定值的性质,通常采用多项式基函数的形式,当设定值在被控区间里其变化率小于或等于某一值时,控制输入的结构可以取一个基函数,当设定值在被控区间里其变化率大于某一值时,控制输入的结构可以取两个基函数.由上分析有:ym(k)=yl(k)+yf(k)(3)由于(k)是均值为零,方差为 2的白噪声,且不可测量.在推导模型的预测输出时可以设(k)=0,根据式(1)模型最合适的预测输出为:ym(k)=nai=1(-ai)ym(k-i)+nbj=1bju(k-j)(4)为了在k时刻预测得到k+i时刻(i=1,2,H)的模型
7、输出值,可根据式(3),(4)推导得到:y1(k+ik)=naj=1(-aj)y0m(k+i-j)+nbj=1bju0(k+i-j)i=1,2,H(5)其中:y0m(k+i-j)=ym(k+i-j),ij+1y1(k+i-j),ij+1u0(k+i-j)=u(k+i-j),ij0,ijyf(k+ik)=ij=1cju(k+i-j)i=1,2,H(6)其中:c1=b1cj=m in(j-1,na)k=1(-ak)cj-k+bj,2jnbcj=m in(j-1,na)k=1(-ak)cj-k,nbjH(7)234信息与控制29卷当取一个基函数时,基函数可以取阶跃函数,即:u(k+i)=u(k)=1
8、,i=0,1,H-1;当取两个基函数时,基函数可以取阶跃函数和斜坡函数,即:u(k+i)=1+i32,i=0,1,H-1.2.2.2滚动优化预测函数控制的参考轨迹可以采用各种形式.对于一个稳定的系统,通常采用一阶指数形式,如式(8)所示:yr(k+i)=c(k+i)-i(c(k)-y(k)(8)式中yr是参考轨迹;c是设定值;,=e(-Ts?Tr),Tr是95的参考轨迹响应时间;y(k)为k时刻的过程输出.预测函数控制的优化算法中,未来控制作用是通过迭代,最优化,约束限制等一系列步骤来实现的.最常用的方法是在优化点上参考轨迹和过程预测输出的误差平方和最小化.其最小化的标准式如式(9)所示:m
9、inJp=Hi=1(yr(k+i)-y(k+i)2(9)y(k+i)=ym(k+ik)+e(k+i)(10)式中y(k+i)为预测过程输出;ym(k+ik)是k+i时刻的模型输出;e(k+i)是未来误差.2.2.3误差校正实际情况下,由于模型失配,二次输入及噪声的影响,模型输出与过程输出之间常存在着误差.为此,需要对未来优化时域中的误差进行预测,预测的方法有多种,其中可以取未来的误差为:e(k+i)=y(k)-ym(k)(11)其中ym(k)为k时刻的模型输出.2.3模型参数的在线辨识参数辨识方法必须假定一种模型结构,通过极小化模型与过程之间的误差准则函数来确定模型的参数.本文考虑式(1)结构
10、的模型,通过极小化广义误差和函数来确定模型的参数,即最小二乘类参数辨识法.随着数据的增加,最小二乘法将出现所谓的“数据饱和”现象.这是由于增益矩阵K(k)随着k的增加将逐渐趋近于零,以致递推算法慢慢失去修正能力之故.考虑到需要在线辨识模型参数,为了克服“数据饱和”现象,我们采用带遗忘因子的递推最小二乘算法.针对式(1),且C(z-1)=1的情况下具体算法如下:设h(k)=-y(k-1),-y(k-2),-y(k-na),u(k-1),u(k-2),u(k-nb)(12)=a1,a2,ana,b1,b2,bnb(13)初始条件:P(0)=a2I,a是充分大的宵数,I是(na+nb)(na+nb)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于A RM 模型自适应预测函数控制 基于 模型 自适应 预测 函数 控制
限制150内