显式预测控制
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
显式模型预测控制综述
亚历山德罗·阿莱西奥和阿尔贝托·本波拉德
摘要:
显式模型预测控制解决了模型预测控制的主要缺点,即需要解决在线的数学程序来估计控制动作。在一些文献中,这个估计阻止了MPC的应用,因为在采样时间内解决最优问题的计算技术比较困难或根本不可行,或者是计算机代码实现数值求解器时导致软件认证问题,尤其是存在安全性的关键问题。
显式MPC在给定的感兴趣的操作条件范围内允许其离线计算优化问题。显式MPC通过利用多参数的编程技术离线计算出的最优控制动作可以作为状态的显式功能和参考矢量,使上线操作简化为简单的函数求值。这样一个函数在大多数情况下是分段仿射的,所以MPC控制器可以映射成线性增益的查找表。
本文调查显式MPC在科学文献上的主要贡献。首先回顾基本概念和MPC 方案的问题,复习解决显式MPC问题的主要方法包括新颖而简单的次优实用方法以降低显式形式的复杂性。本文结尾评论了一些未来的研究方向。
关键词:模型预测控制,显式解,多参数编程,分段仿射控制器,混合系统,最小—最大控制。
1模型预测控制
模型预测控制的控制动作是通过解决在每个采样时刻的有限时域开环最优控制问题而得到的。每个优化产量的最优控制序列,只有第一步应用到过程中:在下个时间步长中,采用最新的可用状态信息作为最优控制新的初始状态通过偏移一个时间范围重复计算。因此,MPC也被成为后退或滚动时域控制。
该解决方案依赖于过程的动态模型,遵守所有的输入和输出(状态)的约束,
并优化其性能指标。这通常表示为一个二次或线性准则,因此,对于线性预测模型产生的优化问题可以转换为一个二次规划(QP )或线性规划 (LP ),而对于混合预测模型,所得到的优化问题可以转换为将在下一章节的讲述的混合整数二次或线性规划(MIQP/ MILP )。因此,MPC 和常规控制的主要区别在于,后者的控制函数在离线情况下预先计算的。MPC 在工业应用上的成功是由于其在处理有许多操纵的控制流程,控制变量及系统对它们的约束的能力。
过程控制通常是由差分方程的系统建模的
(1)((),())x t f x t u t += (1)
其中,()n
x t ∈
是状态矢量, ()m
u t ∈
是输入矢量。假设一个简单的(0,0)0f =,
控制和状态序列要满足的约束条件
()x t ∈℘,()u t μ∈ (2)
m
μ⊆
,n χ⊆在它们的内部封闭集包含原点。假设控制目标引导状态的起源,
MPC 解决了如下的约束监管问题。假设在当前时刻t 的状态量x(t)的完整测量是可使用的,那么下面的有限视距优化监管问题可以解决。
1
0(()):min ((),())()N N z
k x t l x k u k F xN -=P +∑ (3a)
..s t 1(,),0,...,1k k k x f x u k N +==- (3b)
0()x x t = (3c)
,0,...,1k U u k N μ∈=- (3d)
,1,...,1k x k N χ∈=- (3e)
,N N x χ∈ (3f)
(),,...,1k k u u x k N N κ==- (3g) 其中,z ∈
是优化变量的向量,01['...']'N z u u -=,u mN =(一般情况下,z 包括
命令输入和额外的优化变量),闭集N χχ⊆的选择,终端成本F,和终端增益κ,确保MPC 方案的闭环稳定性。在每个时间步长t ,k x 表示在t+k 时刻的预测状态向量,通过将输入序列01,...,k u u -应用到模型(1)中获得k x ,起始状态0()x x t =。预测范围N>0,u N 是输入范围(1U N N ≤≤),“≤”表示特定组件的不平等。N 是有穷的,如果f,l 和F 是连续的,μ是紧凑的,在式(3a)里有最小值。在每个时间步长t ,问题(())N D x t 的求解方法是通过求解程序发现的
min (,())
..(,())0,z
q
h z x t s t g z x t g ≤∈
(4)
从(3)求出最优控制序列*(())z x t ,只有第一个输入被应用在系统(1)
0()*(())u t u x t = (5)
在时刻t+1基于新的状态x(t+1),再次计算优化问题(3)。
MPC 基本设置(3)可以在不同的情况下专门设置。主要取决于预测模型,性能指数和使用的终端条件。 1.1线性模型和二次成本
一个有限域的最优问题(3)的二次阶段的成本是通过设置在(3a )中的公式(6)制定的。
(,),()k k k k k k N N N l x u x Qx u Ru F x x Px '''=+= (6)
其中0,0,Q Q R R ''=≥=>和0P P '=≥为适当维数权重矩阵。式(3b )是一个确定性的线性离散时间预测模型
(,)k k k k f x u Ax Bu =+ (7)
在
式
(3g)
中()x Kx κ=,
μ
,
χ
是多面集,例如
{}}{m i n
m a x
m i n m a x
:,:m
n
u u u u x x
x x
μχ=∈
≤≤=∈≤≤,N χ也是多面集。如果
代入1
10()k k j k k j j x A x t A Bu ---==+∑,式(4)就变为二次规划问题:
11
(,())()()()22
h z x t z Hz x t C z x t Yx t ''''=
++ (8a) (,())()g z x t Gz W Sx t =-- (8b)
其中,0H H '=>,C,Y,G,W,S 为适当维数的矩阵[18].注意Y 是不需要计算
*(())z x t ,它只影响(8a)的最优值。
1.2线性模型和线性成本
范数∞-和1-用来测量性能
其中,,nR m
nQ n
nP n
R Q P ⨯⨯⨯
∈
∈∈
,使用1.1节中所描述的相同设置。遇到范数∞-
等价于
,1,...,,0,...,1,1,...,,1,...,1,1,...,u i k k R u x i k k Q x i N N P
R u i n k N Q x i n k N P x i n εεε≥±==-≥±==-≥±= (10) 其中,在(10)中的下标i 表示第i 行,问题(3)可以被映射为线性规划(LP )[12]
1
(,())[1...1 1...10...0]u u
N mN N
h z x t z -=∂ (11a)
(,())()g z x t Gz W Sx t =-- (11b)
其中,1u N N ∂=-+是最后一次输入的数就是在预测时域重复的次数,
10110[,...,,...,...]u N u u u x x N N z u u εεεε--'''是优化向量,G,W,S 是从模型矩阵A,B(10)的权重Q,R,P 得到的。约束集,,N μχχ和增益K 。1 - 规范的情况下同样可以通
过引入松弛变量,,u i x i x i
ik k ik k iN N R u Q x P x εεε≥±≥±≥±来处理。