系统最优化及控制.ppt课件.ppt

1、最优控制理论选用教材：xxx 编著 最优控制理论 科学出版社教学参考书：符曦编著 系统最优化及控制 机械工业出版社 解学书 最优控制理论与应用 清华大学出版社第一章第一章 绪绪 论论第二章第二章 数数 学学 准准 备备第三章第三章 用变分法求解最优控制问题用变分法求解最优控制问题第四章第四章 极小值原理及其应用极小值原理及其应用第五章第五章 线性二次型问题的最优控制线性二次型问题的最优控制第六章第六章 动态规划法动态规划法第一章第一章 绪绪 论论 1-1最优控制发展简史 最优控制是系统设计的一种方法。它所研究的中心问题是如何选择控制信号才能保证控制系统的性能在某种意义下最优。一：最优控制的发展

2、第二次世界大战以后发展起来的自动调节原理，对设计与分析单输入单输出的线性定常系统是有效的；然而近代航空及空间技术的发展对控制精度提出了很高的耍求，并且被控制的对象是多输入多输出的，参数是时变的。面临这些新的情况建立在传递函数基础上的自动调节原理就日益显出它的局限性来。这种局限性首先表现在对于时变系统，传递函数根本无法定义，对多输入多输出系统从传递函数概念得出的工程结论往往难于应用。由于工程技术的需要，以状态空间概念为基础的最优控制理论渐渐发展起来。最优控制理论是现代控制理论的核心，20世纪50年代发展起来的，已形成系统的理论。最优控制理论所要解决的问题是：按照控制对象的动态特性，选择一个容许控

3、制，使得被控对象按照技术要求运转，同时使性能指标达到最优值。二：研究最优控制的方法从数学方面看，最优控制问题就是求解一类带有约束条件的泛函极值问题，因此这是一个变分学的问题：然而变分理论只是解决容许控制属于开集的一类最优控制问题，而在工程实践中还常遇到容许控制属于闭集的一类最优控制问题，这就要求人们研究新方法。在研究最优控制的方法中，有两种方法最富成效：一种是苏联学者庞特里雅金提出的“极大值原理”；另一种是美国学者贝尔曼提出的“动态规划”。极大值原理是庞特里雅金等人在1956至1958年间逐步创立的，先是推测出极大值原理的结论，随后又提供了一种证明方法。动态规划是贝尔曼在1953年至1958年

4、间逐步创立的，他依据最优性原理发展了变分学中的哈密顿-雅可比理论，构成了动态规划。由于电子计算机技术的发展，使得设计计算和实时控制有了实际可用的计算工具，为实际应用些更完善的数学方法提供了工程实现的物质条件，高速度、大容量计算机的应用，一方面使控制理论的工程实现有了可能，另一方面又提出了许多需要解决的理论课题，因此这门学科目前是正在发展的，极其活跃的科学领域之一。最优控制理论在一些大型的或复杂的控制系统设计中，已经取得了富有成效的实际应用。目前很多大学在自动控制理论课程中已经开始适当增加这方面的内容，而对于自动控制方面的研究生则普遍作为必修课程。求解最优控制问题，可以采用解析法或数值计算法1-

5、2 最优控制问题的实例 例11月球上的软着陆问题 飞船靠其发动机产生一与月球重力方向相反的推力u(t)，以使飞船在月球表面实现软着陆，要寻求发动机推力的最优控制规律，以便使燃料的消耗为最少。设飞船质量为m(t)，高度为h(t)，垂直速度为v(t)，发动机推力为u(t)，月球表面的重力加速度为常数g。设不带燃料的飞船质量为M，初始燃料的总质量为F初始高度为h0，初始的垂直速度为v0，那么飞船的运动方程式可以表示为：)()()()()()()(tkutmtmtugtvtvth初始条件 FMmvvhh)0()0()0(00终端条件 0)(0)(fftvth性能指标是使燃料消耗为最小，即 约束条件)(

6、0tu)(ftmJ 达到最大值 我们的任务是寻求发动机推力的最优控制规律u(t),它应满足约束条件，使飞船由初始状态转移到终端状态，并且使性能指标为极值(极大值)。例12拦截问题在某一惯性坐标系内，设拦截器质心的位置矢量和速度矢量为：LxLx 目标质心的位置矢量和速度矢量为：F(t)为拦截器的推力MxMx MLMLxxvxxx 则拦截器与目标的相对运动方程为：ctFmtmtFtavvx)()()()(其中a(t)是除控制加速度外的固有相对加速度，是已知的。初始条件为：000000)()()(mtmvtvxtx终端条件为：任意)(0)(fftvtxefmtm)(从工程实际考虑，约束条件为)(ma

7、x)(0tFtF如果我们既要求拦截过程的时间尽量短，又要求燃料消耗尽量少，则可取性能指标：fttdttFcJ0)(1为最小 综上所述，所谓最优防天拦截问题，即选择满足约束条件的控制F(t),驱使系统从初始状态出发的解，在某个时刻满足终端条件，且使性能指标为极值（极小值）。1-3最优控制问题的提法 在叙述最优控制问题的提法之前，先讨论一些基本概念。1：受控系统的数学模型一个集中参数的受控系统总可以用一组一阶微分方程来描述，即状态方程，其一般形式为：),(),()(ttutXftXTnxxxX,21是n维状态向量 Tpuuuu,21为p维控制向量),(),(ttutXf为n维函数向量),()(),

8、(),()(),(),()(),(),()(),(),()(),(),()(),(),(),(),(),(),(),(),(),()(2121212122121121ttutututxtxtxfttutututxtxtxfttutututxtxtxfttutXfttutXfttutXfttutXftXpnnpnpnn2：目标集如果把状态视为n维欧氏空间中的一个点，在最优控制问题中，起始状态（初态）通常是已知的，即)0()(0XtX而所达到的状态（末态）可以是状态空间中的一个点，或事先规定的范围内，对末态的要求可以用末态约束条件来表示：0),(0),(11ffffttxgttxg满足末态约束的状

9、态集合称为目标集，记为M，即：0),(,0),(,)();(21ffffnffttxgttxgRtxtxM至于末态时刻，可以事先规定，也可以是未知的。有时初态也没有完全给定，这时，初态集合可以类似地用初态约束来表示。3：容许控制在实际控制问题中，大多数控制量受客观条件的限制，只能在一定范围内取值，这种限制通常可以用如下不等式约束来表示：piuutui2,1)(0max或上述由控制约束所规定的点集称为控制域U，凡在t0-tf上有定义，且在控制域U内取值的每一个控制函数u(t)均称为容许控制。4：性能指标通常情况下，最优控制问题的性能指标形如：fttffdtttutxFttxJ0),(),(),(

10、（）（其中第一项是接近目标集程度，即末态控制精度的度量，称为末值型性能指标。第二项称为积分型性能指标，它能反映控制过程偏差在某种意义下的平均或控制过程的快速性，同时能反映燃料或能量的消耗。5：最优控制的提法已知受控系统的状态方程及给定的初态),(),()(ttutXftX)0()(0XtX规定的目标集为M，求一容许控制u(t)U,t t0,tf,使系统从给定的初态出发，在tf t0时刻转移到目标集M，并使性能指标 fttffdtttutxFttxJ0),(),(),(（）（为最小。这就是最优控制问题。如果问题有解，记为u*(t),t t0,tf,则u*(t)叫做最优控制（极值控制），相应的轨线

11、X*(t)称为最优轨线（极值轨线），而性能指标J*=J（u*()）则称为最优性能指标。1-4最优控制的应用类型设计最优控制系统时，很重要的一个问题是选择性能指标，性能指标按其数学形式可分为如下三类：1）积分型性能指标 fttdtttutXFJ0),(),(这样的最优控制问题为拉格朗日问题。2）终值型性能指标),(ffttXJ这种性能指标只是对于系统在动态过程结束时的终端状态提出了要求，而对于整个动态过程中系统的状态和控制的演变未作要求。这样的最优控制问题为迈耶尔问题。3）复合型性能指标 fttffdtttutXFttXJ0),(),(),(这样的最优控制问题为波尔扎问题。通过适当变换，拉格朗日

12、问题和迈耶尔问题可以相互转换。按控制系统的用途不同，所选择的性能指标不同，常见的有：1：最小时间控制fttfdtttJ0102：最小燃料消耗控制粗略地说，控制量u(t)与燃料消耗量成正比，最小燃料消耗问题的性能指标为：fttdttuJ0|)(|3：最小能量控制设标量控制函数u2(t)与所消耗的功率成正比，则最小能量控制问题的性能指标为：fttdttuJ0)(24：线性调节器给定一个线性系统，其平衡状态X(0)=0，设计的目的是保持系统处于平衡状态，即这个系统应能从任何初始状态返回平衡状态。这种系统称为线性调节器。线性调节器的性能指标为：fttniidttxJ012)(加权后的性能指标为：ftt

13、niiidttxqJ012)(对u(t)有约束的性能指标为：fttTTdttRututQXtXJ0)()()()(21式中Q和R都是正定加权矩阵。一般形式，有限时间线性调节器性能指标：fttTTffTdttRututQXtXtPXtXJ0)()()()(21)()(21无限时间线性调节器性能指标：0)()()()(21tTTdttRututQXtXJP0，Q0，R0，均为对称加权矩阵。5：线性跟踪器若要求状态X(t)跟踪或尽可能接近目标轨迹Xd(t),则这种系统称为状态跟踪器，其相应的性能指标为：fttTdTddttRututXtXQtXtXJ0)()()()()()(21Q0，R0，均为对称

14、加权矩阵。若要求系统输出y(t)跟踪或尽可能接近目标轨迹yd(t),则这种系统称为输出跟踪器，其相应的性能指标为：fttTdTddttRututytyQtytyJ0)()()()()()(21Q0，R0，均为对称加权矩阵。除了上述几种应用类型外，根据具体工程实际的需要，还可以选取其他不同形式的性能指标，在选取性能指标时需注意：1）应能反映对系统的主要技术条件要求2）便于对最优控制进行求解3）所导出的最优控制易于工程实现第二章第二章 数数 学学 准准 备备 2-1函数极值问题一：多变量函数极值问题设二元函数f（x1,x2），在点（x1*,x2*）处有极值f（x1*,x2*）的必要条件为：0),(

15、),(*2*11211xxxxxfxf0),(),(*2*12212xxxxxfxff（x1*,x2*）取极小值的充分条件为：0)(2)(22*21*21*222111xfxxfxfxxxxxx0)(21*2,122212111xxffffxxxxxxxxxx或*22212111xxxxxxxxxxffffF正定),(),(*2*12121211xxxxxfxxf),(),(*2*12121221xxxxxxfxxf),(),(*2*12221222xxxxxfxxf其中上述结论可以推广到自变量多于两个的情形 设n 个变量的多元函数f（x1，x2，xn），若f（x）在x*处有极小值，其必要条件

16、为：0,*2*1221nxxxxfxfxfxF充分条件为：2222122222212212212212*)()()()()()()()()(nnnnnxxxxfxxxfxxxfxxxfxxfxxxfxxxfxxxfxxfF为正定矩阵。二：有约束条件的函数极值问题 设二元函数f（x1,x2），x1和x2必须满足下列方程：g（x1,x2）0 为求函数f（x1,x2）的极值，并找出其极值点（x1*,x2*），作一辅助函数拉格朗日函数：),(),(),(212121xxgxxfxxL式中为辅助变量，称为拉格朗日乘子。函数f（x1,x2）求极值问题，转变为无约束条件函数求极值问题（拉格朗日乘子法），其存

17、在极值的必要条件为 021LxLxLxL或0111xgxfxL0222xgxfxL0),(21xxgL同样，用拉格朗日乘子法可以求有约束条件的n元函数的极值。设n元函数为f（x1，x2，xn），有m个约束方程 0),(21nixxxgi1，2，m（nm）),(),(),(211212121nimiinmnxxxgxxxfxxxL作拉格朗日函数：函数L有极值的必要条件为：01111xgxfxLimii02122xgxfxLimii0),(2111nxxxgL0),(2122nxxxgL0),(21nmmxxxgL01nimiinnxgxfxL2-2泛函极值问题一.无条件约束的泛函极值问题设函数x

18、（t）在 t0,tf 区间上连续可导 定义下列形式的积分dtttxtxFJftt),(),(0J的值取决于函数x（t），称为泛函 1：始端时刻t0和终端时刻tf都给定时的泛函极值 设 dtttxtxFJftt),(),(0函数x*(t)使J为极小 令：)()()(*ttxtx式中是一个很小的参数，（t）是一个连续可导的任意函数 dttttxttxFxJftt),()(),()()(0*其取极小值的必要条件为：00)(xJ上式为J(x)取极小值的必要条件J(x)为极大、极小，通常可根据系统的物理性质来判断。00)(22xJJ(x)取极小值的充分条件0)()()()()(000000ffffftt

19、ttttttttxFdtxFdtddtxFtdxFdtxFdtxFtxFtxJJ(x)取极值的必要条件为：0 xFdtdxF0)(0fttxFt欧拉方程横截条件00)(xJ由必要条件不同函数F的欧拉方程为：),(ttxF0 xF)(txF 022xxF),(ttxF 0222txFxxF)(),(txtxF0222xFtxFxxF xtxtxttxtxF),(),(),(),(0tx当t0和tf给定时，根据x(t0),x(tf)是固定的或自由的各种组合，可导出边界条件（1）固定始端和固定终端0)(0ftttx(t0)=x0,x(tf)=xf 故边界条件为：x(t0)=x0,x(tf)=xf X

20、(t)X1(t)X2(t)X3(t)t0tft0)(0fttxFt由横截条件(2)自由始端和自由终端 00txF0ftxFX(t)t0tft（3）自由始端和固定终端00txFx(tf)=xf X(t)t0tft（4）固定始端和自由终端x(t0)=x0 0ftxFX(t)t0tft极小值的充分条件：00)(22xJ0222222220dtxFxxFxFftt02222220dtxFxxFxxFxFttf故J(x)取极小值的充分条件：222222xFxxFxxFxF为正定 例1 设性能指标为：dttxxJ)(2122边界条件为：x(1)=1,x(2)=2,求J为极值时的x*(t)解 22),(tx

21、xtxxF由欧拉方程 0 xFdtdxF0)21(2txdtd21*)(CtCtx根据边界条件，x(1)=1,x(2)=2 32)(*ttx22000222222txFxxFxxFxF正半定，J(x)为极小值 2：未给定终端时刻的泛函极值问题 若始端时刻t0给定，始端状态x(t0)固定或沿规定的边界曲线移动；而终端时刻tf自由，终端状态x(tf)自由或沿规定的曲线移动，这类最优控制问题称之为未给定终端时刻的泛函极值问题。设系统性能指标:dtttxtxFJftt),(),(0式中t0是已知的，tf未给定，x(t0)给定或未给定)()()(*ttxtx)()()(*ttxtx)(*ffftttdt

22、tttxttxFJffttt),()(),()()(*0J取极值的必要条件为：。00J0)(),(),()(*0fffftttttxtxFdtxFxFtf上式第二项分部积分*0*0*0)()(fffttttttdtxFdtdtxFtdtxF于是有：0)(),(),()()(*0*0fffftttttttxtxFxFtdtxFdtdxFtff得J(x)取极值得必要条件为 0 xFdtdxF欧拉方程 横截条件。0)(),(),()(*0fffftttttxtxFxFtfdttttxttxFJffttt),()(),()()(*0由横截条件可推出各种情况下的边界条件：1）给定始端和自由终端 X(t0

23、)t0tftX*(t)X(tf)tf*此时，x(t0)=x0,(t0)=0,(tf)和(tf)自由 可得边界条件与横截条件为：x(t0)=x0 0*fttxF0),(),(*fffttxtxF由于最优轨线x*(t)的tf即是最优时刻tf*,上式可写为：00)(xtx0,fttxxF。0ftxF0)(),(),()(*0fffftttttxtxFxFtf2）给定始端x(t0)=x0和终端有约束x(tf)=C(tf)X(t0)t0tftX*(t)X(t)tf*C(t)()()(*ttxtx代入)()(fftCtx)()()()(*ffffffffttCttxttttx上式对求偏导，并令0)()()

24、()(*ffffttxtCt可得边界条件与横截条件为：00)(xtx)()(fftCtx。0),()()(ftttxxFxFtxtc0)(),(),()(*0fffftttttxtxFxFtf（3）终端x(tf)固定，始端有约束x(t0)=(t0)X(tf)t0tftX*(t)X(t)tf*(t)边界条件与横截条件为：ffxtx)()()(00ttx0),()()()(0ttxxFtxtxF从以上讨论可以看出，不论边界情况如何，泛函极值都必须满足欧拉方程，只是在求解欧拉方程时，对于不同边界情况，应采用不同的边界条件与横截条件。00)(xtx ffxtx)(ffxtx)(00ttxF00)(xt

25、x 0fttxF00ttxF 0fttxF00)(xtx 0fttxF0,ftttxxF00)(xtx )()(fftCtx0)()(fttFxFtxtc)()(00ttx ffxtx)(0)()(0ttFtxtxFtf固定x(t0)固定x(tf)固定x(t0)自由x(tf)固定tf固定x(t0)固定x(tf)自由x(t0)自由x(tf)自由tf自由x(t0)固定x(tf)自由x(t0)固定x(tf)约束x(t0)约束x(tf)固定无条件约束的泛函极值问题中的边界条件和横截条件列表例2 求使性能指标 dtxJftt021)1(2为极小时的最优轨线x*(t)。设x(0)=1，x(tf)=C(tf

26、)，C(tf)=2-t,tf未给定。解 显然，所给出的性能指标就是x(t)的弧长，也就是说，要求从x(0)到直线C(t)的弧长未最短。tx(t)c(t)2x(0)x*(t)x(t)021)1(),(2xtxxF欧拉方程为：0 xFdtdxF0)1(0212xxdtd。cxx21)1(22221accxax battx)(这是一个x(t0)固定，x(tf)约束情况下的极值问题。由边界条件 x(t0)=x（0）b=1,x(t)=at+1 横截条件 0)()(fttFxFtxtc0)1()1(121212fttxxxx解得 1)(ftx 1a。1)(*ttxtx(t)c(t)2x(0)x*(t)x(

27、t)0由边界条件)()(fftCtxfftt2121*ft2202*2121)1(dtxJtx(t)c(t)2x(0)x*(t)x(t)03：向量函数泛函极值问题 在上面所讨论的公式中，都假定x是1维变量，但是，所有公式都可推广到n维变量的情况 设性能指标 dttxxFJftt),(0式中)()()()(21txtxtxtxn则欧拉方程为 0 xFdtdxF式中 nxFxFxFxF21。nxFxFxFxF21对于始端时刻t0和终端时刻tf都给定时，横截条件 0)(0fttTxFt式中)()()()(21ttttn对于未给定终端时刻tf时的横截条件为：（1）给定始端和终端有约束：0)()(ftt

28、TFxFtxtc（2）给定终端和始端有约束。0)()(0 ttTFtxtxF二有约束条件的泛函极值问题 在实际问题中，对应泛函极值的最优轨线x*(t)通常不能任意选取，而受着各种约束。求泛函在等式约束下的极值，称为条件泛函极值问题。1.代数方程约束设 dtttxtxFJftt),(),(0约束方程 0),(txGnRx mRG nm 构造增广泛函dttxGtxxFJTttaf),(),(0mR令纯量函数),(),(),(txGtxxFtxxLT 00dtxxJfttTLTxLTxLa分部积分000ffttTxLttTTxLdtdTxLxdtGx由于x,相互独立，为使上式成立，应同时满足下述欧拉

29、方程，约束方程和横截条件：欧拉方程：0 xLdtdxL 0 xFdtdTxGxF约束方程：0),(txG横截条件：00fttTxLx利用横截条件，根据始端状态x(t0)和终端状态x(tf)的不同情况，可以导出具体的边界条件和横截条件，其讨论过程和结论与无约束条件的泛函极值问题相同。2：微分方程约束设 dttxxFJftt),(0约束条件：0),(txxGmRG 设纯量函数),(),(),(txxGtxxFtxxLT欧拉方程 0 xLdtdxL 0TxGdtdxFdtdTxGxF约束条件 0),(txxG横截条件 00fttTxLx00fttTxGTxFx3：积分方程约束设 dttxxFJftt

30、),(0约束方程 cdttxxGftt),(0mRG c为一常数 设),()(txxGtZ则 0)(),(tZtxxG0)(0tZctZf)(令 zxxzxx)(ZGFLT欧拉方程 0 xLdtdxL0 xLdtdxL0约束方程),()(txxGtZ0)(0tZctZf)(横截条件 00fttTxLx可见，对于有约束条件的泛函极值问题，可采用拉格朗日乘子法将其转化为无约束条件的泛函极值问题进行求解。在不同边界条件情况下，欧拉方程不变，只是边界条件及横截条件不同。第三章第三章 用变分法求解最优控制问题用变分法求解最优控制问题设系统状态方程:),(tuxfx 性能指标:fttffdttuxFttx

31、J0),(),(式中 pnRuRx和F为纯量函数 最优控制问题就是寻求最优控制)(*tu及最优状态轨迹)(*tx使性能指标J取极值.一.初始时刻 及始端状态 给定,给定,终端自由 0t)(0txft构造增广泛函 fttTffadtxtuxftuxFttxJ0),(),(),(令哈密尔顿函数:),(),(),(tuxftuxFtuxHT则 fttTffadtxtuxHttxJ0),(),(0)()()()(0dtxxHuuHxxHxxJTTTTTttttTaff注意到:fffttTttTttTxdtxdtx0000)(0tx0)()()()(0dtxHuuHxxHxxaJTTTttttTff为使

32、上式成立,应同时满足下列方程:欧拉方程(伴随方程)xH状态方程 Hx 控制方程 0uH横截条件 fftfttTxtxtxxx)()(0)(00对于两端固定的情况下横截条件 ffxtxxtx)(,)(00例 1设系统状态方程为)()(tutxtx）（)(tx的边界条件为0)(,1)0(ftxx求最优控制)(tu使下列性能指标 ftdtuxJ022)(21为最小 解:作哈密尔顿函数)()(2122uxuxH欧拉方程 x控制方程 0,0uuH状态方程 uxxHx,xHxxx)0()(221)0()12()12(2212222tttteexeex)0()12()12(221)0()(2212222tt

33、tteexee消除u由边界条件 0)(,1)0(ftxxfffftttteeee2222)12()12()0(得最优控制)12()12()12()12(2212222222tttttttteeeeeeeeuffffxx1111二.初始时刻 及始端状态 给定,给定,终端约束.0t)(0txft设终端约束方程为 0)(),(ffftxMttxMqRM 构造增广泛函:fttTfTfadtxtuxftuxFtxMvtxJ0),(),()()(fttTfTfdtxtuxHtxMvtx0),()()(式中 qRv0)()()()(0dtxHuuHxxHxvxMxJTTTttttTTaffJ取极值的必要条件

34、是 正则方程 欧拉方程状态方程xHHx控制方程 0uH边界条件和横截条件,)(00 xtx0)(ftxMfttTfvxMxt)()(三.初始时刻 及始端状态 给定,自由,终端约束 0t)(0txft设终端约束为 0),(ffttxM构造增广泛函 fttTffTffadtxtuxHttxMvttxJ0),(),(),(0aJ得J取极值的必要条件为:正则方程 xHHx控制方程 0uH边界条件和横截条件,)(00 xtx0),(ffttxMfttTfvxMxt)()(0fttTtMvtH用变分法求解最优解的必要条件 性能指标 fttffdttuxFttxJ0),(),(),(),(),(tuxftu

35、xFtuxHT系统方程),(tuxfx 约束条件,)(00 xtx0),(ffttxM正则方程 Hx xH控制方程 0uH条件边界条件和横截条件tf给定终端固定终端自由终端约束tf自由 终端固定终端自由终端约束,)(00 xtxfffxtxtxM)(0)(00)(xtx)()(fftxt00)(xtx00)(xtx00)(xtx00)(xtx0)(ftxMftTfvxMxt)()(ffxtx)(ffttH)()()(fftxtffttH)(0),(ffttxMfttTfvxMxt)()(fttTftMvttH)(例2 已知系统状态方程为 1)0(),(xtux 求最优控制)(*tu使性能指标

36、10222)(dtuxeJt为最小 解 本题为 ft给定，终端自由的情况 uuxeHt)(222正则方程:Hx xH控制方程 0uH得 ux txe22tteuuue22)2(202tteuuxe22)2(22消除u02xxx ttecectx)21(2)21(1)(ttececxu)21(2)21(1)21()21(边界条件与横截条件 00)(xtx0)1(,1)0()()(xtxtff求得 221)12()12(12ec22222)12()12()12(eec最后得最优控制)12()12(1)()21(22)21(22*tteeeetu0591.07957.14142.04142.2tte

37、e例3设系统的状态方程为 21xx uxx220)0(1x0)0(2x性能指标 2022221212)2(215)2(21dtuxxJ终端约束条件 15)2(5)2(21 xx试求使 minJ的最优控制 解 本题为 给定2ft终端受约束的 最优解问题)(2122212uxxuHuxu22212)(2122212)2(215)2(21xx015)2(5)2(21xxM正则方程 uxxxxHx2221,11111)(,0,ctxHxH1221222)(,cectxHt控制方程 122)(00cectuuuHtutxtx)()(22123221)(cecectxtt4123121)(ctcecect

38、xtt)(2122212uxxuHuxu22212)(2122212)2(215)2(21xx015)2(5)2(21xxM边界条件和横截条件 05.0,0)0(432cccx05.0321ccc)()()()(111fTffftvtxMtxt1115)2()2(cvx)()()()(222fTffftvtxMtxt1222252)2()2(cecvx代入),2(1x)2(2x123221)(cecectxtt4123121)(ctcecectxtt22212)2(215)2(21xx015)2(5)2(21xxM255.1255.0315)2(5)2(1543705.005.03222143

39、222121432221321432vcececvccececxxccececcccccc解得 73.01c13.02c73.013.0)(12*ttecectu例 4设系统状态方程为 ux 边界条件 0)(,1)0(ftxx试确定最优控制)(*tu使 ftfdtutJ0221为极小 解 这是 ft自由,终端固定的最优解问题 uuH221正则方程 uxHx1,0cxH控制方程 uuuH00应用边界条件 11)0(2cx21ctcxffttH)(1)()()(212ffftuttu2)(1)()(2122fffttt2,21uc12)(21tctcx22,0)(ffttx2)(*tu22*ft例

40、 5设控制对象方程为 ux 0)0(xx终端时刻 ft自由,终端固定 0)(ctxf求)(*tx和)(*tu 使得 ftdtxxJ022)(为极小.解 本题 ft自由,终端固定 uuxuxxH2222uxHxxxH2uuuH2020 xu 0 xxxuuxux ttecectx21)(uecectxtt21)(由边界条件和横截条件 0)0(xxftfttH)(0)()()()(22fffftuttutx)()()(222ffftxtutx0)()()()(fffftxtxtxtx021cc021021ccxcc故 01xc 02c或 01c02xc u2于是最优轨线和最优控制为:当 00cx

41、textx0*)(textxtu0*)()(当 00cx textx0*)(textxtu0*)()(0)(ctxf由可求出终端时刻tf*例6 磁场控制的直流电动机如图所示 MconstIaRfLfUf数学模型 fuxx100010 xkyf0边界条件 0)0()0()0(021xxx00)(1tx性能指标 102tfdtuJ1t给定 试求在t1时间内由x(0)转移到x(t1),并使控制能量具有极小值时的控制输入(励磁电压)uf*,最优性能指标J*和最优轨线x(t)*,解:这是tf 给定,x(tf)固定的最优控制问题 fffuxuxxuH221222112正则方程 11110cxH212122

42、ctcxH控制方程 0fuH212212102ctcuuff代入状态方程得322122141)()(ctctcdttutxf4322312141121)()(ctctctcdttxtx代入给定边界条件310124tc212102tc03c04c则最优控制为 21*2121ctcuf210310*612tttuf最优性能指标 1031202*12tftdtuJ最优轨线 ttttx21023102*660221033101*32ttttx第四章 极小值原理及其应用 用古典变分法解最优控制问题时,假定u(t)不受限制,从而得到最优控制应满足 0uH实际上在工程问题中,控制变量总有一定的限制.设控制变

43、量被限制在某一闭集内 u即u(t)满足 0),(),(ttutxG满足限制条件的u(t)称为容许控制，由于u不能是任意的,0uH的条件已不存在 4-1.连续时间系统的极小值原理设系统状态方程为:),(),()(ttutxftx初始条件,)0()(0pnRuRxxtx为有界闭集,不等式约束为,0),(),(ttutxGG为m维连续可微的向量函数,pm 系统从x0转移到终端状态x(tf)，tf未给定,终端状态x(tf)满足等式约束 0),(ffttxMM为q 维连续可微向量函数,nq 性能指标:fttffdtttutxFttxJ0),(),(),(最优控制问题就是要寻找最优容许控制u(t)使J为极

44、小 令 0)()()(0ttut0)(),(),()()(),(),()(0221tZttutxGtZtztztztZmT且于是,系统方程为:0)(0)()(),(),(00002ttzxtxtxGZtxfx终端时刻tf 未给定,终端约束 0),(ffttxM要求确定最优控制 u使性能指标 fttffdttttxFttxJ0),(),(),(为极小引入拉格朗日乘子向量及,写出增广性能指标泛函dtZtxGxtxftttxFttxMvttxJTTttffTffaf),(),(),(),(),(),(20令哈密而顿函数为),(),(),(txftxFtxHT拉格朗日纯量函数),(),(),(2Ztx

45、GxtxHtzxxTT则 fttffTffadttzxxttxMvttxJ0),(),(),(对J取一阶变分得*0*)()()()()(*fffttTTfttTTfttTadtzzdtddtdxxdtdxtxxMvxttMvtJ令 0aJ可得增广性能指标泛函取极值的必要条件为 欧拉方程 0 xdtdx00zdtddtd横截条件:000)(0)(*ffffttttttTttTTzxvxMxtMvtxx把的表达式代入欧拉方程:00)(zdtddtdxGxHT横截条件:*)()()(*ffttTfttTfvxMxttMvttH由欧拉方程和横截条件知,最优轨线 00*zz以上为使性能指标J取极值的必要

46、条件,为使性能指标为极小,还必须满足维尔斯特拉斯函数沿最优轨线非负的条件,即:0)()()()(),(),(*zzzxxxtzxxtzxxET或:),(),(:0),(),(),(),(*tuxHtuxHtxHtxHxtzxxxtzxxET即上式表明,沿最优轨线函数H相对最优控制u*(t)取绝对极小值,这是极小值原理的一个重要结论.00-*TTuGuHGH)(0)(0由上式表明,在有不等式约束的情况下,沿最优轨线 0uH不再成立定理:(极小值原理)设系统的状态方程为),(),()(ttutxftx控制u(t)是有第一类间断点的分段连续函数,属于p维空间中的有界闭集,满足不等式约束:0),(),

47、(ttutxG在终端时刻tf 未知的情况下,为使状态自初态 00)(xtx转移到满足边界条件 0),(ffttxM的终态,并使性能指标 fttffdtttutxFttxJ0),(),(),(达极小值.设哈密而顿函数为),(),(tuxftuxFHT则最优控制u*(t),最优轨线x*(t)和最优伴随向量*(t)必须满足下列条件:(1).沿最优轨线满足正则方程:TxGxHHx)(式中是与时间t无关的拉格朗日乘子向量,其维数与G相同,若G中不包含x,则:xH(2)横截条件及边界条件:0),()(0)(),()()(00ffttTttTfttxMxtxvtMttuxHvxMxtff(3)在最优轨线x*

48、(t)上与最优控制u*(t)相对应的H函数取绝对极小值,即),(),(*tuxHtuxH并且沿最优轨线,下式成立 TuGuH)(上述条件与不等式约束下的最优控制的必要条件相比较,横截条件及端点边界条件没有改变,仅 0uH这一条件不成立,而代之以与最优控制相对应的函数为绝对极小,其次是正则方程略有改变,仅当G中不包含x时,方程才不改变.当 t0和x(t0)给定,根据tf给定或自由,x(tf)给定,自由或受约束等不同情况下所导出的最优解必要条件列表如下:tf给定 性能指标 终端状态 正则方程 极值条件 边界条件与横截条件 固定 自由 约束 fttfdttuxFtxJ0,)(xHtuGtuxftux

49、FHxGxHHxTT则若0),(),(),()(),(),(min*tuxHHtuxHHHHuffxtxxtx)()(00)()()(00fftxtxtxftTffvxMxttxMxtx)()(0)()(00tf给定 性能指标 终端状态 正则方程 极值条件 边界条件与横截条件 固定 自由 约束 fttdttuxFJ0,xHtuGtuxftuxFHxGxHHxTT则若0),(),(),()(),(),(min*tuxHHtuxHHHHu未知)()()(00ffftxtxxtx0)()(00ftxtxftTffvxMttxMxtx)()(0)()(00tf给定 性能指标 终端状态 正则方程 极值条

50、件 边界条件与横截条件 固定 自由 约束)(ftxJffxtxxtx)()(00)()()(00fftxtxtxftTffvxMxttxMxtx)()(0)()(00 xHtuGtuxftuxFHxGxHHxTT则若0),(),(),()(),(),(min*tuxHHtuxHHHHutf自由性能指标 终端状态 正则方程 极值条件 边界条件与横截条件 固定 自由 约束 fttfdttuxFtxJ0,)(),(),()(tuxftuxFHxGxHHxTT),(),(min*tuxHHtuxHHHHuffffttHxtxxtx)()()(00ffffttHtxtxtx)()()()(00fTftT