优先级多目标稳定化约束模型预测控制_何德峰(1)

合集下载

多目标最优化模型

多目标最优化模型

多目标最优化模型多目标最优化是一种将多个目标函数优化问题组合在一起的方法,旨在找到一个让所有目标函数达到最优的解。

这种方法广泛应用于工程、经济学和决策科学等领域,因为在现实世界中,很少有问题只涉及一个目标。

通过解决多目标最优化问题,我们可以在平衡各种需求和限制条件的基础上做出更好的决策。

在多目标最优化问题中,我们需要同时考虑多个冲突的目标函数。

这些目标函数可以是相互独立的,也可以存在相互依赖关系。

例如,对于一个制造公司来说,我们可能希望同时最小化生产成本和最大化产量,这两个目标是相互矛盾的。

当我们试图减少成本时,产量可能会受到影响,而当我们试图提高产量时,成本可能会增加。

在解决多目标最优化问题时,我们需要定义一个衡量目标函数的目标向量。

这个向量通常包含所有目标函数的值,通过改变决策变量的值,我们可以在目标向量中找到不同的点。

我们的目标是找到一个解,使得目标向量达到最优,即找到一个无法通过改变决策变量的值而得到更好结果的点。

多目标最优化问题的解可以有多个,这些解通过一种称为帕累托前沿的概念呈现。

帕累托前沿是指在不改变任何目标函数值的前提下,无法找到另一个解使得一些目标函数值变得更好的解。

换句话说,帕累托前沿是指在一个多目标最优化问题中,无法一次达到所有目标函数的最优值,因为它们往往是相互冲突的。

解决多目标最优化问题的方法有很多,包括传统的数学编程方法和启发式算法。

在数学编程方法中,我们可以使用多目标规划模型来定义和求解问题。

这种方法的优点是准确性和可解释性高,但在面对大规模和复杂问题时效率较低。

另一种方法是使用启发式算法,如遗传算法、模拟退火算法和粒子群优化算法等。

这些算法通过模拟生物进化和物理过程,逐步解空间并逐渐改进解的质量。

启发式算法的优点是能够在较短的时间内找到满足要求的解,但无法保证最优解。

除了解决问题的方法外,还有一些问题需要考虑。

首先,我们需要定义目标函数,这是一个非常关键和困难的任务。

有界扰动多变量Hammerstein系统输入到状态稳定模型预测控制

有界扰动多变量Hammerstein系统输入到状态稳定模型预测控制

有界扰动多变量Hammerstein系统输入到状态稳定模型预测控制何德峰;余世明【摘要】考虑具有状态和控制约束的有界未知扰动多变量Hammerstein系统,提出一种具有输入到状态稳定和有限L2增益性能的鲁棒非线性模型预测控制策略.基于多变量线性子系统H∞控制律,滚动预测非线性代数方程的解算误差,继而在线优化计算满足系统约束条件的预测控制量.利用输入到状态稳定性概念和L2增益思想,建立闭环系统关于该扰动信号具有鲁棒稳定性和L2增益的充分条件,使闭环系统不仅满足系统约束,而且对不确定扰动输入和解算误差具有鲁棒性.最后以工业聚丙烯多牌号切换过程控制为例,仿真验证本文算法的有效性.【期刊名称】《控制理论与应用》【年(卷),期】2019(036)004【总页数】8页(P605-612)【关键词】Hammerstein模型;模型预测控制;约束控制;输入到状态稳定性;L2增益【作者】何德峰;余世明【作者单位】浙江工业大学信息工程学院,浙江杭州 310023;浙江工业大学信息工程学院,浙江杭州 310023【正文语种】中文1 引言Hammerstein模型是一类由静态非线性输入函数和动态线性模型串联组成的非线性系统.该模型结构简单,但能充分刻画pH中和、空气分离、聚烯烃牌号切换等工业过程的非线性特性[1–4].因此,针对Hammerstein系统的非线性模型预测控制(nonlinear model predictive control,NMPC)算法得到了较多研究[3–16],总体上分为整体策略和两步策略.整体策略直接使用非线性函数和线性模型描述模型预测控制(model predictive control,MPC)最优控制问题,并直接优化控制变量,具有统一处理系统输入与状态约束的优点,但为建立名义闭环系统的稳定性,需要计算“三要素”等稳定性约束条件[5–6],增加了NMPC在线求解的计算量.两步策略[7–14]则根据模型的串级结构,先设计线性子系统的控制输入(即Hammerstein系统的中间量MPC),再求解非线性代数方程与解饱和运算计算实际控制量,具有计算量小且实现方便的优点.进一步,文献[16]考虑解算误差和实际控制输入加权对系统性能的影响,提出了一种改进的两步策略,并结合解算误差建立名义闭环系统的稳定性充分条件.实际系统普遍存在不确定扰动,而MPC是基于有限时域目标函数的控制策略,当系统同时存在约束和扰动时,其可行性和稳定性可能丢失[17],恶化闭环系统的控制性能.输入到状态稳定性(input-to-state stability,ISS)理论[18]为一类持续扰动不确定非线性系统的稳定性分析提供有效工具,近年已用于NMPC的鲁棒性研究.例如,文献[19]引入区域ISS概念,分析约束不确定系统NMPC的吸引域;文献[20]使用双模策略建立NMPC的ISS充分条件;而文献[21]结合H∞控制方法建立NMPC的ISS和有限L2增益的充分条件;等.尽管上述结果原理上也适用于Hammerstein系统鲁棒NMPC整体策略设计,但由于该系统两步法NMPC不采用终端罚函数和终端约束,因此不能用于不确定多变量Hammerstein系统鲁棒NMPC两步策略研究. 在文献[16]的成果基础上,本文应用ISS理论研究输入时变的多变量约束不确定Hammerstein系统鲁棒NMPC策略.采用两步法设计鲁棒NMPC策略,首先建立线性子系统的H∞控制策略,再通过对非线性代数方程解算误差做多步预测,进而在线滚动优化关于线性子系统H∞控制律的跟踪性能,计算满足系统实际控制和状态约束的鲁棒预测控制律.进一步,利用ISS理论和L2增益概念,建立保证闭环系统鲁棒渐近稳定和有限L2增益的解算误差的上界.因此,本文策略不仅可以处理系统各种约束而且对解算误差具有鲁棒性能.最后,以聚丙烯多牌号切换过程控制[16,24]为例,验证本文结果的有效性和实用性.2 问题描述考虑离散时间不确定多变量Hammerstein系统状态空间模型:其中:x∈Rn,u∈Rm和w∈Rp分别是系统的状态向量、控制输入和不确定扰动输入;v∈Rr和z∈Rl分别是系统的中间变量和辅助输出向量;A,B1,B2和H为已知的恰当维数矩阵;g是连续可微的非线性函数向量,描述系统输入和中间变量之间的静态非线性环节,满足g(0)=0.假设原点为系统的平衡点,且本文仅限于状态反馈控制器设计.定义系统的状态和控制输入约束如下:其中Cx⊂Rn和Cu⊂Rm分别是包含原点为内点的紧凸集.为保证存在满足约束(2)的控制输入量,考虑系统中间变量约束[11]其中Cv⊂Rr是包含原点为内点的紧凸集.进一步,考虑不确定扰动输入w,满足其中扰动集W是关于有界持续扰动信号的集合.本文目标是针对有界持续扰动输入(4),设计多变量约束Hammerstein系统的鲁棒预测控制器,使对应闭环系统相对于该扰动输入满足输入到状态稳定和有限L2增益定义1 给定有限常数γ>0和鲁棒不变紧集S⊆Rn,其中S内含原点.如果存在一个有限项ξ(x0),对任意初始条件x0∈Rn,满足ξ(0)=0和不等式条件其中wt∈W(t=0,1,···,N),则称不确定系统(1)和(4)在集S上具有有限L2增益[18]. 定义2 考虑不确定非线性系统xt+1=f(xt,wt),如果存在K∞类函数γ和KL类函数β,使得对任意输入wt∈W和初始条件x0∈Rn,在任意时刻t系统的解xt=x(t;x0,w)存在,且满足如下不等式:其中w={w0,···,wt−1},则称该系统具有输入到状态稳定性(ISS)[18].引理1 不确定非线性系统xt+1=f(xt,wt)是输入到状态稳定的,系指存在一个连续的正定函数V(x)满足其中:∀x∈Rn,w∈Rp,函数α1,α2和α3为K∞类函数,δ为K类函数,并称V(x)为系统的ISS–Lyapunov函数[18].3 鲁棒NMPC策略设不确定线性子系统的一个鲁棒控制律为,定义时刻t非线性代数方程解算误差令xt+i|t,ut+i|t和vt+i|t别为系统(1)在时刻t对时刻t+i的状态、输入和解算误差的预测变量,则其中:∀i=0,1,···,N−1,N>0为预测时域.将等式组(9)简写为由于扰动w不可测量,对约束不确定系统(1)–(4),采用名义解算误差模型定义N步滚动优化目标函数.其中R>0为控制输入加权矩阵,则在时刻t定义有限时域最优控制问题应用非线性数值规划算法或随机搜索优化算法求解优化问题(12),得最优控制序列.但MPC采用滚动优化与前景控制策略,即在时刻t施加控制量,在时刻t+1重复求解优化问题(12),得到.注1 为了保证反馈控制律是不确定线性子系统的鲁棒控制律,令为该线性系统的一个H∞控制律[21].其中对称正定矩阵P满足Riccati不等式组且矩阵A+B1K+B2K2稳定,其中:显然,Hammerstein系统MPC律与不确定线性子系统xt+1=Axt+B1vt+B2wt的鲁棒控制律密切相关.为此,令H∞控制律(13)为Hammerstein系统线性子系统的鲁棒控制律,并据此定义Hammerstein系统MPC律为则对应闭环系统为注2 最优控制问题(12)采用名义解算误差模型定义,此时得到的预测控制律及其闭环系统(15)在扰动作用下不一定满足系统约束(2)–(3).对此,可采用微分对策(Differential game)原理或紧缩(tightening)状态约束集[22],修改定义最优控制问题(12),保证该优化问题存在时间迭代优化可行性,即预测控制律在任意时刻总是存在.为描述清晰,本文假设最优控制问题(12)在任意时刻都是优化可行的,进而分析闭环系统(15)关于不确定扰动(4)的输入到状态稳定性和有限L2增益性能.定理1 考虑闭环系统(15),如果存在一个鲁棒不变集S⊆Rn,对任意状态x∈S,解算误差∆v满足不等式则该系统在集S上关于不确定扰动(4)具有输入到状态稳定性和有限L2增益性能. 证考虑线性子系统xt+1=Axt+B1vt+B2wt的鲁棒控制律(13),定义如下集合:由于约束集Cx和Cv是包含原点的紧凸集,且vrb是x的连续函数,故集Sx总是存在且非空.考虑Riccati不等式组(14)的一个正定矩阵解P,定义正定函数V(x)=xTPx 及其一个水平集S={x∈Rn:V(x)6r}⊆Sx,则存在r>0使集S存在且非空.定义线性子系统的一个Hamiltonian函数其中x+=Ax+B1v+B2w.令.由H∞控制理论[23]可知,线性子系统xt+1=Axt+B1vt+B2wt存在鞍点(v∗,w∗)=(Kx,K2x)[22],其中,增益K和K2见注1所示.将v∗和w∗代入函数(17),整理得其中矩阵令,则和.将K和K2代入Π,则由文献[21]和不等式(14)可得,Π=HTH−P+ATPΘA<0.令ε为一充分小的正数,满足Π6−εI,则代入式(18)可得将不等式(16)代入式(19),得Z(x,v∗,w∗)6−εxTx对任意x∈S成立.将Hamiltonian函数在鞍点(v∗,w∗)做一阶泰勒阶数展开,得其中O(∥s∥)表示∥s∥的高阶无穷小量,则结合式(14)可得,对任意状态x∈S,函数(17)满足结合不等式(21),沿闭环系统(15)的轨迹对V(x)作差分运算,得令δ(∥w∥)=∥w∥2和α3(∥x∥)=xT(εI+HTH)x,则δ和α2为K∞类函数,且进一步,函数V(x)是x的正定函数,存在K∞类函数α1和α2满足α1(∥x∥)6V(x)6α2(∥x∥),则根据引理1可得,闭环系统(15)在集S上对扰动输入(4)具有ISS.为证明闭环系统(15)在集S上对扰动输入(4)具有有限L2增益,考虑任意时间区间[0,N],其中N>0.对任意初始状态x0∈S,在上述区间对不等式(22)进行累加,得即从而有根据定义1可知,闭环系统(15)在集S上对扰动输入(4)具有有限L2增益. 证毕.注3 定理1表明,可以通过判断每个时刻的解算误差条件(16),建立有界扰动多变量闭环Hammerstein系统的鲁棒稳定性,从而不再使用整体法策略的终端约束条件[5–6].尽管条件(16)与文献[16]定理1中条件在形式上具有相似性,但文献[16]中的条件是基于名义线性子系统的二次型最优控制律设计得到,对于不确定扰动作用,通常不能保证闭环线性子系统的鲁棒渐近稳定性[21].尽管如此,结合条件(16)和文献[16]定理1中条件,本文对多变量约束Hammerstein系统在无扰动和有扰动作用下MPC闭环稳定性分析给出了一种形式统一的充分条件.推论1 考虑闭环系统(15),如果存在一个鲁棒不变集S⊆Rn,对任意状态x∈S,解算误差∆v满足不等式其中λmax(M)表示矩阵M的最大特征值,则该系统在集S上关于不确定扰动(4)具有输入到状态稳定性和有限L2增益性能.证条件(27)是不等式(16)成立的充分条件,则结合上述定理结论得到该推论. 证毕. 下面给出约束Hammerstein系统NMPC算法的实施步骤.步骤1 初始化设计参数(H;R;N);步骤2 离线求解Riccati不等式组(14),并计算鲁棒控制律(13);步骤3 读入时刻t的状态测量值xt,应用数值优化算法或进化算法求解优化问题(12).当时刻t解算误差满足条件(16)时,得当前时刻的最优控制序列,并将首个控制量作用于系统(1);步骤4 令t=t+1,转步骤3.注4 在Hammerstein系统现有的两步法NMPC策略中[7–14],MPC控制律是在第1步针对线性子系统的虚拟输入即中间变量设计,需要通过解饱和代数运算得到实际控制输入,此时无法考虑非线性代数方程解算误差对系统约束和性能的影响;而本文MPC控制律是在第2步直接对实际控制输入计算,显式地考虑了该解算误差对系统满足约束和控制性能的要求.其次,稳定性条件(16)或(27)与预测时域N无关,从而在保证闭环系统稳定的条件下,调整预测时域以减小优化问题(12)的在线计算量.最后,由于本文算法直接对控制输入设计,可以处理更一般的系统约束,因此可以在优化过程中直接对实际控制输入做加权处理,以利于系统安全生产.4 聚丙烯多牌号切换控制考虑环管反应器聚丙烯牌号切换过程,对比现有Hammerstein系统两步法鲁棒MPC策略[13],验证本文算法的有效性.为书写简便,将两步法鲁棒MPC策略[13]简记为con–RMPC,而用new–RMPC表示本文策略.假设环管反应器为连续搅拌反应器,则聚丙烯质量连续时间状态空间模型为[16,24]其中:状态向量x=[x1x2]T=[lgMIcEtc]T,中间变量v=[g1(u)g2(u)]T,控制输入u=[T CH2/Cm Cm2/Cm]T,非线性函数g1(u)=lgMIi和g2(u)=Eti分别为在式(28)–(30)中,常数τ表示丙烯在反应器内的平均停留时间(τ=2h),MIc和Etc表示聚合物的累积熔融指数和累积乙烯含量,MIi和Eti表示聚合物的瞬时熔融指数和瞬时乙烯含量;T,CH2,Cm和Cm2分别为反应温度、氢气浓度、丙烯浓度和乙烯浓度;模型参数ki和rj(i=1,···,6;j=1,2)通常随生产牌号变化,可用工业数据辨识计算[3,24].在工业生产中,累积熔融指数和累积乙烯含量通过离线化验或利用软测量技术获得[24],通常存在不确定测量误差.设聚丙烯累积质量参数的测量误差为±10%,并令其测量误差为系统(28)的不确定扰动输入w=[w1w2]T,其中w1=ω1lgMIc和w2=ω2Etc及−0.16ω1,ω260.1,即实际累积熔融指数和累积乙烯含量动态变化过程状态方程为其中v1和v2由式(28)给定.进一步,令控制周期Ts=0.5h,则聚丙烯牌号切换过程质量不确定离散时间状态空间模型为选择聚丙烯均聚牌号A,B和无规共聚牌号C及切换顺序A→B→C,并假设生产装置在第10小时由A→B切换,在第30小时由B→C切换.采用双层控制结构作为切换控制组态[3],即牌号切换控制器作为上层控制器,聚丙烯装置回路控制器作为底层控制器,其中牌号切换控制器的控制量将作为回路控制器的设定值,实现聚丙烯多牌号切换生产过程的控制.注意,底层控制通常采用PID控制器直接控制聚丙烯生产过程的回路,在本文策略中并不考虑聚丙烯生产过程回路的底层控制器设计,故在该仿真中假设聚丙烯生产过程的底层控制器具有理想化的控制效果.此外,均聚牌号A和B 生产工艺不含乙烯输入量,即Cm2/Cm对应的控制输入始终为零.进一步,3种牌号规格及切换过程约束可见表1,其中上标c,i和∆分别表示累积量、瞬时量和增量约束.此外,不确定扰动系数ω1和ω2分别为[−0.1,0.1]的随机数.在仿真中,现有两步法鲁棒NMPC(con–RMPC)策略[13]的预测时域N=4,二次型性能指标中状态加权矩阵Q=I2,中间量加权矩阵R=5I2,采用Min-max策略定义不确定线性子系统的MPC控制器,并用MATLAB函数fminimax求解.考虑相同预测时域取值,本文两步法鲁棒NMPC(new–RMPC)策略的参数为H=diag{1,0.7}和R=diag{1,1010,1010},则求解Riccati不等式组(14),得一个正定对称矩阵解和H∞控制律增益分别为P=4I和K=−0.75I,并采用MATLAB函数fmincon优化问题(12),则图1–3给出了分别采用两种鲁棒MPC策略的聚丙烯牌号切换过程控制结果,其中:图1为牌号切换过程累积质量变化曲线,图2为牌号切换过程瞬时质量变化曲线,图3为牌号切换过程控制输入曲线.在图1–3中:setpoint表示牌号切换过程控制的设定值;new–RMPC表示本文所提出鲁棒预测控制算法;con–RMPC表示常规鲁棒预测控制算法;limits表示对应变量的上下限约束.图1 牌号切换过程累积质量曲线Fig.1 Cumulative quality pro files of the grade transition process图2 牌号切换过程瞬时质量曲线Fig.2 Instantaneous quality pro files of the grade transition process图3 牌号切换过程控制输入曲线Fig.3 Control input pro files of the grade transition process对比2种Hammerstein系统两步法鲁棒MPC策略,由图2上子图点线和表1可知,采用con–RMPC策略实施聚丙烯牌号由A→B切换控制时,瞬时熔融指数在切换时刻超出约束范围;进一步,由图1–3点线可知,当con–RMPC策略控制聚丙烯牌号由B→C切换时,熔融指数和乙烯含量存在较大的稳态偏差,这将严重影响目标牌号产品C的质量,而造成该稳态误差的主要原因是con–RMPC策略计算得到的控制分量乙烯与丙烯比值过小(如图3下子图点线所示).与此同时,在本文两步法鲁棒MPC策略控制下,聚丙烯的整个切换过程都是严格满足系统的约束条件,且具有较小的稳态偏差.如果以聚丙烯物性规格指标的±5%作为切换目标牌号的合格检验标准,则由图1–3可以看出,本文两步法鲁棒MPC牌号控制尽管存在累积质量的不确定测量误差,但在均聚牌号和无规共聚牌号两次切换过程中,聚合物的MI和Et离开起始生产牌号的合格品区后都先进入废料生产区,再都能最终进入目标牌号的合格品区,同时两次切换过程的(操作)控制变量都满足切换过程的约束.表1 3种牌号规格及切换过程约束[16]Table 1 Three grade specifications and constraints of transition process牌号规格MI/(g/10min)Et/%T/KCH2/%Cm2/Cm/%A 2.7—343.150.050—B 39—343.150.330—C 10 2.4 343.15 0.236 1.99[2.4,42.0]c—[341.15,345.15][0.02,0.35]—A→B[2.0,45.0]i ——[−0.1,0.10]∆ —[8.0,42.0]c [0.0,2.5]c [341.15,345.15][0.20,0.35] [0.0,2.5]B→C [6.0,45.0]i [0.0,2.6]i—[−0.1,0.10]∆ [−0.1,0.10]∆表2给出了本文策略与现有两步法鲁棒MPC策略作用下的牌号切换过程过渡时间对比结果.尽管在现有两步法鲁棒MPC策略作用下,均聚牌号切换过程(A→B)的过渡时间与本文策略作用下的过渡时间相近,但在现有两步法鲁棒MPC策略控制下的无规共聚牌号切换过程(B→C)将无法保证进入合格品生产过程(如图1点线所示);相反,在本文策略作用下,均聚牌号和无规共聚牌号切换过程都能够快速进入目标牌号的合格品生产过程.表2 两种策略作用下牌号切换过渡时间Table 2 Transient times of grade transition by the two schemes两种策略A→B B→C本文 9.0 h 7.0 h文献[13] 9.5 h∞但需要指出的是,在均聚牌号切换过程中,现有两步法鲁棒MPC策略计算的控制输入比本文方法更加平滑,如表3给出的两种策略作用下各个控制量的偏差平方和,这将有利于底层控制器的安全平稳运行.表3 两种策略作用下控制量偏差平方和Table 3 Square sum of control deviations by the two schemes两种策略T/KCH2/Cm/%Cm2/Cm/%本文24.0500 0.0433 0.1912文献[13]4.3×10−74.7×10−41.6×10−105 结语以聚丙烯树脂牌号切换过程为研究背景,针对不确定多变量约束Hammerstein系统,提出了具有扰动输入到闭环状态稳定的NMPC策略.该策略以线性子系统H∞控制为基础,集成了Hammerstein系统两步法和整体法NMPC设计优点,可以处理系统各种实际约束以及在优化过程中考虑系统实际输入加权以平缓系统运行.同时该策略显式地考虑非线性代数方程解算误差对实际控制输入满足系统约束的影响,建立了保证闭环系统鲁棒稳定和有限L2增益的解算误差的界.最后通过对聚丙烯牌号切换过程的仿真控制,验证了本文策略的有效性和实用性.参考文献:【相关文献】[1]BHANDARI N,ROLLINS D.Continuous-time Hammerstein nonlinear modeling applied to distillation.AIChE Journal,2004,50(2):530–533.[2]LAWRYNCZUK M.Nonlinear predictive control of dynamic systems represented by winer-hammerstein models.Nonlinear Dynamics,2016,86(2):1193–1214.[3]HE Defeng,ZOU Tao,YU Li.Receding horizon control of grade transition of large-scale polypropylene production plants.CIESC Journal,2010,61(2):405–412.(何德峰,邹涛,俞立.大型聚丙烯生产装置牌号切换滚动时域控制.化工学报,2010,61(2):405–412.) [4]FENG Shaohui,QIAN Feng.The application of model prediction control based on Hammerstein model in Chemical Process.Automation in Petro-ChemicalIndustry,2004(3):25–27.(冯少辉,钱锋.基于Hammerstein模型的预测控制在化工过程中的应用.石油化工自动化,2004(3):25–27.)[5]KHANI F,HAERI M.Robust model predictive control of nonlinear processes represented by wiener or Hammerstein models.Chemical Engineering Science,2015,129:223–231. [6]BLOEMEN H H J,BOOM T J J V,VERBRUGGEN H B.Modelbased predictive control for Hammerstein-wiener systems.International Journal of Control,2001,74(5):482–495.[7]ZHUXF,SEBORGDE.NonlinearpredictivecontrolbasedonHammerstein models.Control Theory&Applications,1994,11(5):564–575.[8]DING B C,XI Y G,LI S Y.Stability analysis on predictive control of discrete-time systems with input nonlinearity.Acta Automatica Sinica,2003,29(6):827–834.[9]XU Xiangyuan,MAO Zongyuan.The analysis and reseaerch of predictive control based on Hammerstein model.Control Theory&Applications,2000,17(4):529–532.(徐湘元,毛宗源.基于Hammerstein模型的预测控制的分析与研究.控制理论与应用,2000,17(4):529–532.)[10]DING BC,PING XB.Dynamic output feedback model predictive control for nonlinear systems represented by Hammerstein-Wiener model.Journal of ProcessControl,2012,22(9):1773–1784.[11]DING Baocang,LI Shaoyuan.Design and analysis of Hammerstein nonlinear controlsystems with constraints.Control and Decision,2003,18(1):24–28.(丁宝苍,李少远.具有约束的Hammerstein非线性控制系统的设计与分析.控制与决策,2003,18(1):24–28.)[12]DING B C,XI Y G.A two-step model predictive control of Hammersteinmodel.International Journal of Robust and Nonlinear Control,2006,16(7):353–367.[13]NIU Yongxiao,DING Baocang,SUN Hexu.Robust stability of twostep predictive control for systems with input nonlinearitis.Control and Decision,2006,21(4):457–462.(牛永肖,丁宝苍,孙鹤旭.输入非线性系统的两步法预测控制的鲁棒稳定性.控制与决策,2006,21(4):457–462.)[14]ZHANG H T,LI H X,CHEN G R.Dual-model predictive control algorithm for constrained Hammerstein systems.International Journal of Control,2008,81(10):1609–1625.[15]LI Yan,MAO Zhizhong,WANG Yan,et al.Predictive control of Hammerstein-Wiener nonlinearity based on polytopic terminal region.Acta Automatica Sinica,2011,37(5):629–638.(李妍,毛志忠,王琰,等.基于多面体终端域的Hammerstein-Wiener非线性预测控制.自动化学报,2011,37(5):629–638.)[16]HE Defeng,YU Li.Nonlinear predictive control of constrained Hammerstein systems and its research on simulation of polypropylene grade transition.Acta Automatica Sinica,2009,35(12):1558–1563.(何德峰,俞立.约束Hammerstein系统非线性预测控制及在聚丙烯牌号切换中的仿真研究.自动化学报,2009,35(12):1558–1563.)[17]XI Yugeng,LI Dewei.Fundamental philosophy and status of qualitative synthesis of model predictive control.Acta Automatica Sinica,2008,34(10):1225–1234.(席裕庚,李德伟.预测控制定性综合理论的基本思路和研究现状.自动化学报,2008,34(10):1225–1234.)[18]KHALIL HK.Nonlinear Systems[M].3rd Edition.New York:Prentice Hall,2002.[19]MAGNI L,RAIMONDO DM,SCATTOLINI R.Regional input-tostate stability for nonlinear model predictive control.Regional Inputto-State Stability for Nonlinear Model Predictive Control,2006,51(9):1548–1553.[20]LAZAR M,PENA D M,HEEMELS WPMH,et al.On input-to-state stability of min-max nonlinear model predictive control.Systems&Control Letter,2008,57(1):39–48.[21]HE Defeng,JI Haibo,ZHENG Tao.Nonlinear H∞robust predictive control with bounded persistent disturbances.Acta Automatica Sinica,2008,34(2):215–219.(何德峰,季海波,郑涛.持续有界扰动下的非线性H∞鲁棒预测控制.自动化学报,2008,34(2):215–219.)[22]HE Defeng,DING Baocang,YU Shuyou.Review of fundamental properties and topics of model predictive control for nonlinear systems.ControlTheory&Applications,2013,30(3):274–288.(何德峰,丁宝苍,于树友.非线性系统模型预测控制若干基本特点与主题回顾.控制理论与应用,2013,30(3):274–288.)[23]LIN W,BYRNES CI.H∞Control of discrete-time nonlinear systems.IEEE Transactions on Automatic Control,1996,41(4):494–510.[24]RICHARDS J R,CONGALIDIS J P.Measurement and control of polymerization puters and Chemical Engineering,,2006,30(10/12):1447–1463.。

具有输入约束和扰动补偿的四旋翼无人机姿态稳定模型预测控制

具有输入约束和扰动补偿的四旋翼无人机姿态稳定模型预测控制

具有输入约束和扰动补偿的四旋翼无人机姿态稳定模型预测控

郑潇;何德峰;王巧灵;马家腾
【期刊名称】《计算机测量与控制》
【年(卷),期】2024(32)5
【摘要】为了解决四旋翼无人机飞行过程中的姿态控制问题,考虑受到执行器输入约束和外界未知扰动影响,设计了一种具有输入约束和扰动补偿的四旋翼无人机姿态稳定模型预测控制方法;设计简化的四旋翼无人机姿态动力学模型,降低控制器设计的复杂程度,设计带有输入约束的控制器,模拟饱和输入现象,实现饱和输入下的四旋翼姿态稳定控制;设计风扰观测器,实现对外部未知扰动的估计,有效跟踪外部持续扰动,并由此设计扰动补偿律;围绕代价函数设计带有扰动补偿律的最优控制律,作用于四旋翼姿态系统,实现四旋翼无人机姿态的稳定控制;最后进行数值仿真,设置风扰观测器参数λ_(i)为0.25,预测时域N_(p)为10,控制时域N_(c)为9,仿真测试文章方法与不带观测器的非线性预测控制方法(NMPC),验证文章控制方法的有效性和优越性。

【总页数】8页(P94-101)
【作者】郑潇;何德峰;王巧灵;马家腾
【作者单位】浙江工业大学信息工程学院
【正文语种】中文
【中图分类】TP273
【相关文献】
1.一种基于内模补偿的四旋翼无人机姿态系统的渐近跟踪控制律的实现
2.输入饱和与姿态受限的四旋翼无人机反步姿态控制
3.具有自适应增益调节的四旋翼无人机姿态期望补偿RISE控制
4.具有时变扰动的四旋翼无人机有限时间预定性能控制
5.基于扰动观测器的四旋翼无人机自适应姿态控制方法
因版权原因,仅展示原文概要,查看原文内容请购买。

目标规划模型

目标规划模型

目标规划模型目标规划是一种多目标决策方法,旨在寻找一个可行的目标向量,这个向量最好满足一组优先级排序的目标。

目标规划模型可以用来解决多目标决策问题。

目标规划模型通常包括以下几个要素:决策者的目标向量、决策变量、约束条件和目标函数。

决策者的目标向量是指决策者对决策问题中各个目标的优先级排序。

在目标规划模型中,通常将目标向量表示为一个具有多个元素的向量,每个元素表示各个目标的权重。

决策变量是可以被决策者调整的变量,在目标规划模型中,在决策变量的取值范围内寻找一个可行的解。

决策变量的具体取值将影响各个目标的实现程度。

约束条件是对决策变量的限制条件。

这些限制条件可能是由于资源有限,或由于业务规则等原因导致的。

约束条件是确保决策方案可行和符合实际情况的必要条件。

目标函数是目标规划模型的核心部分。

目标函数是一个由决策变量和目标向量构成的函数,表示决策方案对各个目标的实现程度。

目标函数的含义是在满足约束条件的前提下,最大化或最小化目标向量中的各个元素。

目标规划模型的解决方法通常有两种:基于罚函数的解法和基于切比雪夫距离的解法。

基于罚函数的解法通过引入罚函数,将目标规划问题转化为单目标规划问题,然后使用传统的单目标规划方法求解。

基于切比雪夫距离的解法则通过计算决策方案与目标向量之间的切比雪夫距离,将目标规划问题转化为一个单目标规划问题。

目标规划模型的求解过程通常包括以下几个步骤:确定决策变量、建立目标函数、建立约束条件、确定目标权重、求解目标规划模型。

目标规划模型具有以下几个优点:可以考虑多个目标,能够灵活地适应不同的决策需求;可以根据决策者的需求制定不同的目标权重,不受固定的优先级限制;可以通过引入不同的解决方法,得到不同的结果,提供更多的选择。

总之,目标规划模型是一种多目标决策方法,可以用于解决多目标决策问题。

它通过优化决策方案和目标向量之间的关系,寻找一个满足决策者需求的最优解。

目标规划模型具有灵活性和鲁棒性等优点,是现代决策科学中的重要工具之一。

多目标决策理论与方法

多目标决策理论与方法

多目标决策理论与方法
多目标决策理论与方法是一种能够同时考虑多个目标和多个决策变量的决策模型和决策方法。

它的主要目标是在决策过程中寻找一种最优解,使得在给定的约束条件下,多个目标都能够得到最优的平衡。

多目标决策理论与方法主要包括以下几个方面:
1. 多目标决策模型:多目标决策模型是描述多个目标之间的关系和权衡的数学模型。

它可以通过建立目标函数和约束条件,并确定决策变量的取值来寻找最优解。

2. 多目标优化方法:多目标优化方法是为了求解多目标决策模型而设计的算法和技术。

它包括传统的多目标线性规划、多目标非线性规划、多目标动态规划等方法,以及基于进化算法的多目标优化方法如多目标遗传算法、多目标粒子群优化算法等。

3. 多目标决策分析:多目标决策分析是为了帮助决策者对多个目标进行评估和权衡的方法。

它可以通过建立决策树、层次分析法、模糊综合评价等方法来帮助决策者进行决策分析和决策评估。

4. 多目标决策支持系统:多目标决策支持系统是为了帮助决策者进行多目标决策而设计的计算机应用系统。

它通过结合多目标优化方法和决策分析方法,提供
了一种交互式的决策支持环境,帮助决策者进行决策方案的生成、评估和选择。

总之,多目标决策理论与方法为决策者提供了一种全面、系统和有效的决策工具,能够考虑多个目标和多个变量之间的复杂关系,帮助决策者做出更好的决策。

第十章 多目标优化方法简介

第十章 多目标优化方法简介
q * 2
hv ( X ) 0(v 1, 2, , p)
求解上述问题得到的设计方案既考虑了目标函 数的重要性,又最接近完全最优解,因此,它是原
多目标优化问题的一个更加理想、更加切合实际的
相对最优解。
(3)功效系数法
每个分目标函数 f k ( X ) 都可以用一个对应的功效系
数 k (0 k 1)来表示该项设计指标的好坏,规定:
分层序列法及宽容分层序列法
分层序列法的基本思想是将多目标优化问 题式中的J个目标函数分清主次,按其重要程度 逐一排除,然后依次对各个目标函数求最优解。 不过后一目标应在前一目标最优解的集合域内寻 优。
现在假设f1(x)最重要,f2 (x)其次,f3 (x)再其次,…。 首先对第一个目标函数f1(x)求解,得最优值
D ( k ) x f i min x f i x f i max x
xD
min f k x (k )
统一目标法
统一目标法又称综合目标法。它是将原多目标 优化问题,通过一定方法转化为统一目标函数或综 合目标函数作为该多目标优化问题的评价函数,然 后用前述的单目标函数优化方法求解。
(1)加权组合法
加权组合法又称为线性加权法或加权因子法。
即在将各个分目标函数组合为总的“统一目标函数”
的过程中,引入加权因子,以平衡各指标及各分目
标间的相对重要性以及他们在量纲和量级上的差异,
因此,原目标函数可写为:
min f X wk f k X
k 1 q
s.t.
gu ( X ) 0(u 1,2, , m)
多目标优化问题概述
实际的工程设计和产品设计问题通常有多个设
计目标,或者说有多个评判设计方案优劣的标准。

基于扰动预测的网联车鲁棒协同巡航预测控制

基于扰动预测的网联车鲁棒协同巡航预测控制

, %% ,
浙江工业大学学报
第 !" 卷
率和驾驶安 全 性 有 着 重 要 意 义 & )#=A* 相 较 于 自 适 应 巡航控制!+NDXMET;OULEP;O8IMU89"+77#"7+77 不 仅可以通过传感器 获 取 前 车 速 度 和 位 置 信 息"还 可 以通过 3"3 通信获取更多附加信息"如前车当前的 加速度信息&有效利用这些信息可以很好地提高车 辆队列的整体 性 能& 在 网 联 车 协 同 控 制 中"如 何 有 效利用 3"3 的附加信 息 提 高 控 制 性 能 已 经 受 到 了 学者们的广泛关注 & )%=#$* 模型预测控制!W8N;9XU;= NEOMET;O8IMU89"W_7#作 为 一 种 先 进 的 控 制 技 术 已 被广泛应用于智能 网 联 车 控 制 中"与 传 统 控 制 算 法 相比"它能够处 理 饱 和 约 束 和 多 目 标 控 制 问 题& 然 而传统的 W_7 对 模 型 较 为 依 赖"当 出 现 模 型 外 的 扰 动 时 "控 制 器 可 能 违 反 约 束 "无 法 保 证 原 本 的 控 制 性能 & )##=#"* 为了解 决 这 一 问 题"提 出 了 可 以 显 式 处 理扰动影 响 的 鲁 棒 W_7 算 法& 目 前 鲁 棒 W_7 主 要分为 CEI=CDY W_7 和 ML<; W_7& 前 者 考 虑 扰 动作用下性能指标最差的 情 况 来 求 解 优 化 变 量 & )#A* 然而考虑到该算法在线求解的优化问题较为复杂" 在实际应 用 中 实 时 性 可 能 无 法 得 到 保 证& 而 ML<; W_7 在线计算时"无须显式处理扰动 的 影 响&ML<; W_7 在离线阶段已经计算好扰动的鲁棒 不变 集"之 后通过扰动的鲁棒不变集对名义系统的约束进行紧 缩处理"以确保实际系统在扰动下仍满足约束 & )#%=#!* 在车辆控制过程中 扰 动 是 不 可 避 免 的"因 此 针 对 使 用鲁棒 W_7 处理扰 动 已 经 有 了 诸 多 研 究 成 果",EI 等 针 )#B* 对 +77 系 统 前 车 未 来 速 度 未 知 的 情 况"首 先通过常数预测的 方 法 得 到 前 车 的 速 度 预 测 值"然 后根据已知的物理 约 束 条 件"通 过 离 线 计 算 得 到 一 个更紧凑 的 扰 动 不 变 集"达 到 降 低 保 守 性 的 目 的& ,L8等 通 )#>* 过 设 计 比 例 多 重 积 分 观 测 器 来 估 计 扰 动"采用前馈和干 扰 补 偿 的 方 式 组 成 综 合 控 制 律 来 控制 _,-!_U;N;O;PP8U=9;DN;U=H8998Z;U#通 信 拓扑 下 的 车 队 "保 证 车 辆 在 扰 动 下 达 到 控 制 目 标 &

基于收缩块的高速列车多目标滚动时域速度规划

基于收缩块的高速列车多目标滚动时域速度规划

0引言
在能源价格日益上涨与环保压力日益增大的今 天,铁路运输以其运输能力强、安全舒适、快速准时、 节能环保和全天候运输等优势,已成为运输方式未 来的发展趋势⑴。至2018年底,中国铁路营业里程 已达到13.1万公里,高速铁路营业里程已达到2. 9 万公里,占全世界高速铁路总里程的66%以上,位 居世界第一⑵。近5年来我国铁路旅客发送量平 均增长率为9. 89% ,2018年我国铁路运输系统完成 旅客发送量33.75亿人⑵。但逐年增加的发送量 给铁路运输带来了各种不可忽视的问题,如能源消 耗的增加和运行延误的增多。2018年国家铁路能 源消耗折算为标准煤1624. 21万吨⑵,2018年中国 高速铁路电力消耗240.9亿kWh⑶。巨大的能源 消耗,不仅给运营企业带来巨大的压力,也向环境问 题日益严峻的社会发出了挑战。在列车能耗中,大
的保证。因此,自动列车驾驶(automatic train opera­ tion, ATO)系统成为了如今HST运行控制系统研究 的热点⑺。ATO系统能够自动调整HST的车速达 到最佳车速,负责HST的启动、牵引、制动和停车等 工况下的自动控制。ATO系统具有提高HST运行 准点性、能效性和舒适性的特点,能有效避免人为失 误导致的事故发生⑺。一般来说,ATO系统分为最 优轨迹规划和速度跟踪控制上下两层。通过优化上 层的轨迹规划算法,得到更优的参考轨迹,是一种主 要的降低列车能耗的方式⑻。HST在ATO的操控 下,其运行性能与上层规划轨迹紧密相关,规划一条 具有高准点性、能效性、舒适性和安全性的速度轨 迹,对提高列车的运行性能至关重要。
高技术通讯2021年第31卷 第6期:615-627
doi :10. 3772/j. issn. 1002-0470. 2021.06. 006
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

1832





第 28 卷
用目标函数在线切换方法, 设计多目标预测控制器; 文献 [15-16] 结合理想点法提出了一类稳定化多目标 MPC 策略. 这些方法可有效地处理多目标 MPC 问题, 但当系统运行模式和目标函数及其优先级发生变化 时, 这些方法都需要重新设计相应优化参数. 显然, 这 需要高级人员参与操作完成, 控制系统无法自行完成 控制器组态更新. 近年来, 针对多目标函数广泛存在的优先级、 不 可公度和冲突问题, 文献 [17] 在模块多变量控制[18] 的 基础上提出了字典序多目标 MPC 策略, 并成功应用 于水箱液位等过程控制
0 引

法[1,3,6] , 即通过权系数将多目标控制问题转换为单目 标问题, 并以权系数的大小表示目标函数的相关重要 性. 但是, 目前并没有完整的权系数调整规则, 需要在 实际应用中试凑确定. 虽然函数加权法设计简单, 但 权系数无法显式处理各目标函数的优先级、 不可公度 和冲突等问题. 为此, 近年来国内外相关学者提出了 一些新的多目标 MPC 策略. 例如: 文献 [7] 将多目标 MPC 问题转换为离线计算目标函数加权系数问题, 进 而将多目标控制问题转换为单目标问题; 进一步, 文 献 [8] 给出了多目标线性 MPC 闭环系统的稳定性结 论和多参数规划计算方法; 文献 [9-12] 将各目标函数 的优先级表示为混合逻辑整数规划问题, 通过整数规 划算法在线计算多目标预测控制量; 文献 [13-14] 采
模型预测控制 (MPC) 具有统一处理系统约束、 性能指标优化和多变量控制问题的能力, 在过程、 机 械和航天等领域都有成功的应用[1-3] . 在预测控制器 设计中, 目标函数是必不可少的要素之一, 通常用来 表征控制系统的性能, 如超调量、 调节时间和稳态误 差等. 然而, 随着控制系统及其功能的复杂化, 目标函 数除常见的设定值目标外, 还有如生产效率、 节能降 耗和控制量约束等经济性目标. 这些目标函数具有不 同的重要性即优先级和一定的冲突性, 但缺乏统一的 度量标准即不可公度
其中: ������������ 为当前 ������ 时刻的初始状态值; ������������∣������ 和 ������������∣����� + ������ 时刻的状态变量和控制 变量的预测值; ������������ ∈ ������������������ 为当前 ������ 时刻优化序列, 即
Stabilizing constrained model predictive control with prioritized multiobjectives
HE De-feng, SONG Xiu-lan, HUANG Hua
(College of Information Engineering, Zhejiang University of Technology, Hangzhou 310023, China. Correspondent: HE De-feng, E-mail: hdfzj@) Abstract: For the different prioritization problem of objective functions, a multi-objectives model predictive control scheme with guaranteed stability is proposed for constrained multi-variables linear invariant-time systems. Firstly, some results on lexicographic optimal solutions on the multi-objectives predictive control problem are given based on the multi-objective optimization theory. Then the multi-objectives predictive control problem is reformulated as the lexicographic one that is able to deal with the different prioritization of objective functions. Then, ingredients of terminal constraints, terminal penalty functions and local state feedback laws are used to prove the asymptotical stability of the multi-objectives predictive control closed-loop system. Finally, a simulation example is exploited to illustrate the effectiveness of the results obtained. Key words: model predictive control;prioritization;multi-objective control;stability;lexicographic programming
首先, 基于多目标优化理论给出多目标预测控制问题的字典序最优解结果, 并在此基础上考虑目标函数的优先级, 重 新将多目标预测控制问题定义为字典序多目标预测控制问题; 然后, 采用终端约束、 终端罚函数和局部状态反馈律 等三要素, 证明多目标预测控制闭环系统是渐近稳定的; 最后, 通过一个仿真实例验证了所提出方法的有效性. 关键词: 模型预测控制;优先级;多目标控制;稳定性;字典序规划 中图分类号: TP273 文献标志码: A
(5)
极小化每个目标函数, 若多目标优化问题 (4) 可行, 即 式 (5) 存在, 则根据滚动优化原理定义预测控制器为
������mpc = ������∗ 0∣������ , ������ = 0, 1, ⋅ ⋅ ⋅ . ������
(6)
下面结合多目标优化理论[22] , 对多目标预测控 制问题 (4) 定义如下基本概念. 定义 1 化序列 ������������ ∈ ������ 定义 2 可行解 ������∗ ������ , 等式成立:
������������ = {������0∣������ , ������1∣������ , ⋅ ⋅ ⋅ , ������������ −1∣������ }.
. 进一步, 文献 [20] 总结了
现有性能指标、 约束及其不可行问题, 给出了更一般 性的字典序多目标 MPC 策略; 基于此策略, 文献 [21] 给出了城市污水处理系统的多目标预测控制器设 计. 然而, 现有字典序多目标 MPC 侧重于应用研究, 缺乏对稳定性和鲁棒性等理论性质的分析, 目前尚未 见到有关字典序多目标 MPC 稳定性和鲁棒性等理论 成果的报道. 稳定性是闭环控制系统的基本理论性质 之一, 也是控制器能有效运行的基本条件, 因此有必 要研究字典序多目标 MPC 闭环系统的稳定性问题. 本文考虑多变量离散时间线性定常系统, 针对优 化目标函数的不同优先级问题, 提出一种稳定化多目 标约束预测控制策略. 采用多目标优化理论, 定义多 目标预测控制问题的字典序最优解[22] . 在此基础上, 首先考虑各目标函数的优先级, 重新将多目标预测控 制问题定义为字典序多目标预测控制问题. 进一步, 采用 MPC 稳定性理论中的终端约束、 终端罚函数和 局部状态反馈 “三要素” 法[2,23] , 给出了多目标预测控 制闭环系统渐近稳定的充分条件. 最后, 通过一个实 例仿真验证了本文结果的有效性.
[19]
标 ������������1 优先级高于目标 ������������2 ; 整数 ������ ⩾ 1 为预测时域;
������������ (������) 为终端罚函数; 函数 ������������ (������, ⋅), ������������ (⋅, ������) 和 ������������ (������) 分
������������
(4a)
s.t. ������������+1∣������ = ������������������∣������ + ������������������∣������ , ������0∣������ = ������������ , ������������+1∣������ ∈ ������, ������������∣������ ∈ ������, ������ = 1, 2, ⋅ ⋅ ⋅ , ������ − 1. (4b)
则 ������������ 称为该优化问题的一个可行解. 给定多目标优化控制问题 (4) 的一个 如果对于任意其他可行解 ������������ , 使得如下不 (7)
则 ������∗ ������
称为该优化问题的最优解.
显然, 求解多目标优化控制问题 (4) 的最优解是 很困难的, 而且当子目标函数之间有冲突时, 最优解 可能不存在. 因此, 对于控制问题 (4), 通常没有绝对 或唯一的最优解. 定义 3 给定优化控制问题 (4) 的一个可行解 (8)
������������ (������ )∣������∗ < ������������ (������ )∣������������ , ������ = 0, 1, ⋅ ⋅ ⋅ , ������, ������
对于优化控制问题 (4), 如果多目标优
������������
及其驱动的闭环系统满足约束 (4b),
������∗ ������ ,
1 问题描述与预备
考虑如下多变量离散时间线性定常系统:
������������+1 = ������������������ + ������������������ , ������ = 0, 1, ⋅ ⋅ ⋅ .
相关文档
最新文档