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

广义预测控制理论

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）｝分别表示系统的输入、输出和白噪声序列，‘刊。广义预测控制使用如下的二次目标函数

N

2 N

M

J =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 q

nc

q‘是后移算子；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^q4nb

A也(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 J

y*(k j k) = %(k j) (' g i q」)：u(k j -1), N^ j < N2 (2.2.12)

i =0