广义预测控制原理-江苏科技大学

广义预测控制理论1引言预测控制思想主要是在70年代形成的，进人80年代后，随着模型算法控制（MAC）的问世，相继出现了动态矩阵控制（DMC）、扩展时域预测自适应控制（EPSAC）等结构各异的预测控制算法，这些算法分别基于有限脉冲响应和有限阶跃响应模型，算法简单，容易实现，1984年，Clarke及其合作者在上述算法的基础上，提出了广义预测控制（GPC）思想及基本方法，GPC基于参数模型，引入了不相等的预测水平和控制水平，使系统设计更灵活。

由于广义预测控制具有预测模型、滚动优化和反馈校正三个基本特征，因而具有优良的控制性能和鲁棒性，被认为是具有代表性的预测控制算法之一并被广泛应用于过程工业中。

近年来，广义预测控制吸引了众多学者对其进行研究国际上，各大控制会议和杂志对它也非常关注，近10年来的美国控制会议（ACC）、IEEE决策与控制会议（CDC）和国际自动控制联合会（IFAC）世界大会几乎每年都有关于预侧控制的专题分组及以预侧控制为主题的工作讨论会，1995年在韩国又召开了关于预测控制的国际讨论会，在广义预测控制方面也发表了不少综述文献和著作。

2广义预测控制2.1广义预测控制的基本算法GPC采用如下CARIMA模型来描述系统A(z ' )y(t) = B(z ')u (t -1) C (z (t) / ：其中A（z」），B（z'）,C（z'）分别是阶数位的n a,n b,n c的z J的多项式，A（z‘）和C（z‘）是首一多项式，｛u（t）｝,｛y（t）｝,「（t）｝分别表示系统的输入、输出和白噪声序列，‘刊。

广义预测控制使用如下的二次目标函数N2 NMJ =E{ ' [y(t k) 一，(t k)]2亠二；[：u(t k -1)]2} (2.1.2) k 出1 k =1其中，N1,N2分别为最小、最大预测长度，N M为控制长度，满足关系仁N「N2,N M乞N2且当k>N2时看，假定u（t k-T，■为控制加权序列,(2.1.1)■'(t k)是经柔化后的参考值，在GPC 中，不要求对象输出直接跟踪设定值'■, 只要求y(t)沿着参考轨迹到达设定值■。

极小化目标函数J ,并根据滚动优化的原则，得控制律为u(t) =u(t -1) g ( • - f)其中g 为一行向量「为一参考序列向量,f 为由已知输入和输出组成的数 据向量。

GPC 控制方法的具体推导应合理选择 GPC 中的N I ，N 2,N M 以及’，可使GPC 取得较好的控制性能，其它的一些预测控制的方法可以认为是 GPC 的特殊情况， 例如：当N 1 =1，N 2 = N M = N, ■ =0时，相当于 Richalet 提出的IDCOM ，当 弘二也小皿"时相当于性能指标中不加权的 GMV 控制算法。

2.2广义预测控制的预测模型在GPC 中，采用最小方差控制中所用的受控自回归积分滑动平均模型 (CARIMA )来描述受到随机干扰的对象：A(q 」)y(k) =B(q[q 』u(k) C(q ) (k)其中：A(q J ) =1 yq 」B(q 」)二b 。

bq‘ 他4" C(qJ) =c ° qq/G c qncq‘是后移算子；y(k)q 4 = y(k -1) ； =1-q‘为差分算子；(k)是一个独立的随机噪声序列，为研究方便，如若假设 d =1，则模型可简化为：则k j 时刻系统模型为:A(q 4)y(k j) =B(q 4)u(k j -1) C (qj)因为y(k - j)中含有未知信息，因此引入Diophantine 方程获得系统在k j 时刻的输出预测值。

Diopha ntine 方程：(2.1.3)(2.2.1)(2.2.2)(2.2.3)仁A(q 二)E j(q»： q」F j(q」) (2.2.4)其中：E j(q」)二e j,。

• e j,i q‘•乳如丄亠F i(q1^ f j,。

f j,i q‘f j^q"E j和F j由A(q」)和预测长度j唯一确定，由(131.3)、(131.4)可化简得到如下方程：y(k j^BE j.：u(k j -1) F j y(k) E j (^ j)，从而得到GPC预测模型为：y M(k j)=G j.：u(k j -1) F j y(k)(2.2.5)其中G j 二BE j 二(j j) =g j,0 g j,1q J - g j,nb^q4nbA也(2.2.6)因此，对于未来k - j时刻的输出估计只使用k时刻之前的输出以及我们根据最优性能指标确定的输入来确定即可。

式(2.2.2)可简化为：A(q°)y(k) =B(q 冷丸(k-1) C(q J (k)(2.2.7)其中A(q4)=A(q[(1-q')=1 Nq_n a = n a 1, =1,a n^ -a n a,a^ a^a i4,^ i 乞n a，则k 时刻对k - j 时刻的误差可记为：~(k+j|k)=y(k+j)-y(k+jk),j 釘(2.2.8)使预测误差的方差：J-E{~(k+jk)}(2.2.9)最小的j步最优预测y*(k jk)由下列差分方程给出:C(q[y*(k jk)二F j(q')y(k) G j(q「：u(k j T)(2.2.10)此时最优预测误差表示为：~ (k • j k) -E j(q ) (k • j)(2.2.1 1假设C(q」)=1此时式(2210)可简化为:j Jy*(k j k) = %(k j) (' g i q」)：u(k j -1), N^ j < N2 (2.2.12)i =0j/其中，％(k j) = F j (q 」)y(k) [G j (q 」)-(、gq 」)u(k j-1)]j」式(2.2.12)中的y/k - j )由过去的控制输入和输出决定，而 (g i q-L \ =u(k j -1)i =0由现在和未来的控制输入决定。

式(2.2.12)用向量和矩阵的形式表示是：y = y 1 G = u(2.2.13)其中，y * 二[y *(k N 1 k),y *(k N 1 1k), y *(k N 2k)]Ty 二[%(k N 1), Y 1(k N 1 1),%(k N 2)]T•心二[ ：u(k), . ：u(k 1), . ：u(k N u -1)]TpN^g N严…g 。

0…0 IgN tg N^Jg 0G =I ：: • --gq卫叫」9叫4g N^Nu _(N 2亠制趣u其中N 1,N 2分别称为最小和最大预测步长，N u 为控制步程(N u ：：： N 2)，G 中参数 均为被控对象开环阶跃响应系数。

根据最优预报可知： y(k j)二 y *(k jk) E(q 」)(k j) (2.2.14) E{y(k j)} =y *(k jk)(2.2.15)在GPC 中，为了将输出值y(k)按一定响应速度平滑地过渡到由参数轨迹确定的 期望值w ，参考轨迹通常可取为如下的一阶滞后(一阶平滑)模型：y r (k)二 y(k) (2.2.16) Y r (k j) -y(k j _1)(1 )w j =1,2，…,N 2其中[0,1]为输出柔化系数(2.2.17)2.3滚动优化GPC采用的是对输出误差和控制增量加权的二次型性能指标:J = E{( - y)T(y「- y) ■ . ：u T. ：u}(2.3.1)其中,y r 二[y r(k NJ,y r(k N1 1), y r(k N2)]T(2.3.2)y =[y(k N1), y(k N1 1), y(k 2)]丁(2.3.3)■为控制增量加权系数，将式(1.321)对.心求极值，可得最优控制律为:T _1 T=u = (G G 川1) G (y r -y1)(2.3.4)则当前的控制作用为:u(k) =u(k-1) [1,0/ ,0](G T G J)」G T(y r-y1) (2.3.5) 2.4在线辨识与反馈校正GPC控制算法只使用一个控制模型，通过在线反馈校正来保证其准确的预测。

将对象(227)改写为：y(k) =-A(q')y(k) B(qJ・：u(k -1) (k) (241) 其中，A(q」)二A(q」)-1 则有，y(k)二T(kp (k)(2.4.2)其中，(k) =[-y(k-1) -y(k -n a)」u(k-1) • ：u(k - 门厂1)]丁V - [a a n_ b o b n b]T可用带遗忘因子的递推最小二乘法(RLS)来估计模型参数二值:^k)二？(k -1) K(k)[y(k) - ：(k)彳(k -1)](2.4.3)K(k) =P(k -1) ：(k)L： (k)P(k) ：(k)」]' (2.4.4) 1P(k) [I —K(k)「(k)]P(k —1) (2.4.5) 其中，0 ：：：「1是遗忘因子，一般可取0.95」叮；K(k)是权因子，P(k)为正定协方差阵。

控制启动时需要设置参数向量和协方差阵P的初值，通常设定方法为：”P(0)=G2|, a为充分大的实数」(2.4.6) 邂0)=名，5为充分小的实向量总结GPC算法计算机实现步骤如下：1) 根据被控对象及其控制要求给出n a,门匕小仆也‘肌，■和〉；2) 设置初值P(0)和彳(0)，输入初始数据；3) 读取y(k)，用带遗忘因子的RLS递推估计；4) 用辨识得到的参数代替式(227)中的A和B，递推求解Diophantine方程得到F j和G j ；5) 构造向量和矩阵G，并计算(G T G"I)‘ ；6) 计算(G T G」I)°G T的第一行；7) 按式(2.3.5)计算并施加控制u(k)；8) 返回3)继续循环。