《最优控制变分法幻灯片.ppt》由会员分享，可在线阅读，更多相关《最优控制变分法幻灯片.ppt（171页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、最优控制变分法第1页，共171页，编辑于2022年，星期六 一、泛函和变分一、泛函和变分 如绪论中所述，泛函是一个(或一些)函数的函数，它的函数值取决于这个(或这些)函数的选取。更确切地说：如果有一类函数集，对每一个都有一个确定的值与之对应记作 。比如说，旋转体在流体中运动时受到的阻力 取决于它的几何形状。用数学的话来说，就是如果我们把它的几何形状用函数 来描述，那么，是的函数。旋转体可以加工成各式各样的形状，因而函数 同t之间存在着无限多种可能的函数关系，形成一个函数的集合，如图1 1所示。这里便是一个泛函。第2页，共171页，编辑于2022年，星期六第3页，共171页，编辑于2022年，星
2、期六如果要求选择一种形状使旋转体受到的阻力最小，这就意味若要求从这一族函数中寻找出使 达到最小的 来。这就是一个求泛函 的极小值问题，也是一个变分问题。这个问题是牛顿在1686年提出的，它是变分法所解决的最早的一个问题。变分法中有3类基本问题，即拉格郎(Lagrange)问题、马耶耳(Mayer)问题和波尔札(Bo1sa)问题。这3类问题在最优控制问题中都遇到，它们之间的主要区别在于性能泛函的形式不同。（1）拉格朗问题拉格朗问题的泛函表示为式中，是n维矢量；是m维矢量；是独立变量；是 、和 连续函数。绪论中基于性能指标(02)、(010)和(01)的最优控制问题是拉格郎问题的实例。第4页，共1
3、71页，编辑于2022年，星期六(2)马耶耳问题 马耶尔问题的泛函表示为绪论中基于性能指标（09）的最短时间控制问题和基于性能指标（015）的最优推力方向角选择问题就是马耶耳问题的一个特例。(3)波尔扎问题 波尔扎问题的性能泛函是 图1-1所示的最小阻力弹头形状问题就是波尔扎问题的一个实例在高超音速中零攻角下旋转体在流体中受到的压阻力可以精确地表示为 其中，表示旋转体各处半径；是旋转体最大半径；是旋转体的长度；为常数。第5页，共171页，编辑于2022年，星期六从上面的分析可以看出，拉格朗问题的性能泛函是一个积分，马耶尔问题的性能泛函是关于初始状态时间、初始状态和终端时间、终端状态的某个函数，
4、而波尔扎问题的性能泛函则是两者之和。可见，波尔扎问题具有更一般的形式。但是，这3类问题之间常可互相转化。比如，把泛函 改写成 第6页，共171页，编辑于2022年，星期六就把一个马耶耳问题转化成一个拉格郎问题，反之亦然。变分问题是研究泛函的极值路问题；而求某个泛函的极大值等价于求它的负值的极小值。因此，我们主要讨论极小值问题。在讨论泛函的极小值以前，首先回顾普通函数的极小值。给出函数 ，假设它在区间 上连续。如果在该区间上有 ，对于一切 ，有 则 是函数 在区间 上的总体极小值。总体极小值是唯一的，而达到总体极小值 的不一定唯一，比如图1-2所示的函数 在 和 两处都达到总体极小第7页，共17
5、1页，编辑于2022年，星期六第8页，共171页，编辑于2022年，星期六如果 是 的一个局部极小值，则有 且 （1.11）条件(1.11)是局部极小值的必要条件，但不是充分条件。比如 函数 在 处满足条件(1.1 1)，但是 不是极小值，如图13所示。第9页，共171页，编辑于2022年，星期六第10页，共171页，编辑于2022年，星期六如果 且 （1.12）则是一个局部极小值。条件（1.12）是局部极小值的充分条件。给出泛函 ，如果在容许函数（即一切使泛函 存在的函数）中，有函数，对一切容许的 有则 是 的总体极小值。为了定义泛函的局部极小值，需要引出函数空间的“邻域”的概念：考虑区间
6、上两个连续可微的函数 和 ，如果对于每一个，两函数值都互相接近，则函数：和 是在零阶意义下接近的；如果对于每一个 ，不但 值同 值接近，而且它的导数值 值同 值接近；则称函数同函是在一阶意义上接近的。第11页，共171页，编辑于2022年，星期六如果在容许函数的集合函数集合中，对一切同 零阶接近的 ，有 则称 是 泛函的一个强局部极小值；如果对一切同一阶接近的 ，有 则称 是 泛函的一个弱局部极小值弱局部极小值。容易看出，在定义泛函的总体极小值、强局部极小值和弱局部极小值时，函数 是分别同数目逐次减少的函数的集合相比较的。因此，若局部极小值的必要条件也是强局部极小值的必要条件，强局部极小值的必
7、要条件也是总体极小值的必要条件，反之则不尽然。第12页，共171页，编辑于2022年，星期六二、固定端点时间、无约束条件的变分问题二、固定端点时间、无约束条件的变分问题 这一节，我们讨论一类最简单的变分问题，即无约束条件、端点时间固定，只有一个自变量函数的拉格郎问题。通过这个问题来引出欧拉方程和横截条件。求解变分问题，就是要把使泛函达到极值的那个自变量函数找出来，这就需要利用欧拉方程和横截条件。因此，欧拉方程和横截条件是求解变分问题的基础。在推导欧拉方程和横截条件时要使用一个定理，这个定理叫作变分法的基本颈备定理。本节首先介绍基本预备定理，接着推导欧拉方程，然后讨论横截条件，最后讨论泛函取极值
8、的充分条件。第13页，共171页，编辑于2022年，星期六1 1变分法的基本预备定理变分法的基本预备定理 基本预备定理的内容是：如果函数 在区间 上是连续的，而且对于只满足某些一般条件(例如一阶或若干阶可微，在 上的端点处 ，或 且 等)的任意选定的函数，有 (121)则在区间 上 下面我们来证明这个定理。由于函数 ，是任意选定的，因此，可以取 (1.22)其中 的是任一满足条件 第14页，共171页，编辑于2022年，星期六的函数，式中c为某个函数，它在区间 上各点的函数值及 其导数值可以选得任意小。因此，这样选择的 满足定理要求 的一切条件。把方程(122)代入(121)，得 式中被积函数
9、 是非负的。一个非负函数的积分等于 零，只能是被积函数恒等于零，而 是满足某些条件的任意 函数。因此，要使上式成立，必有 进而，如果有 其中 是互相独立，并满足基本预备定理要求的任意选定的函数，那么，可以轮流令除以 外的其它 ()都为零，再使用基本预备定理，可得 ，上述结果可以应用到包含有多个自变量函数的泛函。第15页，共171页，编辑于2022年，星期六2.2.欧拉方程欧拉方程 现在，我们来推导欧拉方程和相应的横截条件。首先讨论固定端点问题，然后讨论未定端点问题。考虑最简单的泛函 (123)的极值。其中 是 二次可微函数；，是变量 和 连函函数，并且有连续二阶偏导数，端点时间 和 固定。首先
10、研究容许函数（或曲线）端点固定的情况，即规定 和 。图14示出了一族容许函数。现在的的问题是要从这一族容许函数(或曲线)中找出使泛函J取极值的函数(或曲线)，即极值函数或极值曲线。第16页，共171页，编辑于2022年，星期六假设 是 极值曲线是任一邻近于它的容许曲线，如图15所示。令 (124)和 (125)是 的连续可微函数，叫做 的变分。是一个小参数。选择不同的 和 就能够把 邻近的容许函数表示出来。因此，式（1.25）表示了 邻近的任意一族容许函数。特别地，当 时，；当 时，。这一族函数都起始于 ，终止于 。第17页，共171页，编辑于2022年，星期六第18页，共171页，编辑于20
11、22年，星期六第19页，共171页，编辑于2022年，星期六将式(125)两边对 求导，可得 (126)将式(125)、(126)代入式(123)，又得 (127)在式(127)中，每选择一个 ，都可作一条 曲 线。选择各式各样的容许的 ，可以作出一族 曲线，如图16所示。因为 所以有 因此，所有 曲线都在 上达到极值。于是，在极值函数上，有 (128)这就是泛函J取极值的必要条件。将式（127）的右边对参数求导，再令结果中的 ，便可求出极值函数 。泛函对参数求导遵守以下法则：假设 是变量 的函数，而又都是独立变量的 可微函数，则 第20页，共171页，编辑于2022年，星期六 将方程(127
12、)两边对求偏导，得到 已知被积函数L具有二次连续偏导数，因此，上式右边的求导和积分可交换顺序，于是得到 利用泛函对参数求导法则，有把它代入上式，则第21页，共171页，编辑于2022年，星期六把它代入上式，则 在上式中令 ，则式(128)便可写成 （1.29）由于式中 和 彼此相关，因此，暂时还不能使用基本预备定理。对上式积分号下第二项使用分积分，可得 (1.210)将式(1210)代入(129)，得到必要条件：第22页，共171页，编辑于2022年，星期六前面假设端点是固定的，即 ，因此，于是，方程(1211)变成 式是连续函数，而 是任意函数，且 ，具备了基本预备定理的全部条件。因此，在极
13、值函 数上，有 第23页，共171页，编辑于2022年，星期六方程(1212)叫作欧拉方程。欧拉方程的积分曲线是一族曲线，这一族曲线叫作极值函数或极值曲线。只有在极值函数上，泛函 才能达到极值。方程(1212)是二阶微分方程，它的解包含两 个积分常数。这两个积分常数要利用边界条件 和来确定。只有满足欧拉方程同时又满足条件的函数，才能在满足给定端点条件下使泛函到达它的极值。因此，求解变分或寻找泛函的极值函数问题归结为求解欧拉方程。在前面我们假设边界条件 和 是预先规定的。现在假设边界条件未规定，和 是任意的。这样的容许函数就更多了，的邻域里包含的函数也更多了。它除包含前面讨论的有公共端点的函数(
14、比如 外，还包含没有公 共端的函数(比如 )，如图17所示。第24页，共171页，编辑于2022年，星期六第25页，共171页，编辑于2022年，星期六不难理解，如果函数 能在未定端点问题中使泛函（123）达到极值，那么，对于同 有公共端点的更狭窄的一类函数来说，必然也能使泛函（1.23）达到极值。也就是说，同样是欧拉方程(212)的解。在这里 和未规定是任意的。因此，必要条件（1.211）包括了下面两个方程：(1.213)和 （1214）上述结果可以解释：既然 和 是任意的，那么，当它们同时为零时，方程（1.211）必须成立，于是得出该式左边第一项，即积分项为零，进而利用变分法的基本预备定理
15、得出式(1.21 3)。另外，当 和 不为零时式(1211)也必须成立，于是得出条件（1.214）。方程（1.213）就是未定端点问题的欧拉方程，方程（1.214）叫做横截条件。这两个方程一起是泛函（1.23）取极值的必要条件。然而，满足必要条件的函数是否确使泛函取得极值，以及其极值究竟是极大值还是极小值还要利用后面将要讨论的充分条件。3对于横截条件的说明对于横截条件的说明 求极值函数 ，不但需要解欧拉方程，还需要正确地使 第26页，共171页，编辑于2022年，星期六第27页，共171页，编辑于2022年，星期六第28页，共171页，编辑于2022年，星期六第29页，共171页，编辑于202
16、2年，星期六（1）固定端点问题 起始条件和终端条件预先规定为 和 。一切容许函数都必须经过两个固定端点 和 ，如图18所示。在这种情况下横截条件(1214)利用 =0来满足。边界条件是 和（2）固定始端、未定终端问题 初始条件固定，终端条件未定，其容许函数如图19所示。这里 0，可任意选择。于是，横截条件(1214)变成 和（3）未定始端、固定终端问题 终端条件规定为 ，初始状态 未规定，其容许函数如图l10所示，相应的横截条件是 和 （4）未定端点问题 端点状态末规定，容许函数如图111所示。假设 、互不相关，条件(1.214)可写成 和又知道 和 任意，于是得到横截条件 和第30页，共17
17、1页，编辑于2022年，星期六例121 求泛函 在边界条件 ，下的极值曲线。解：这个问题的欧换方程是 其通解为 利用边界条件，可求得，。因此，极值曲线是 例122求泛函 的极值曲线。第31页，共171页，编辑于2022年，星期六解：这个问题的欧拉方程是 它的一次积分为引入参数 ，令 ，代入上式可得 式中 。对上式两边微分，得 于是得到 对上式积分，得 或于是得到极值曲线参数方程：，消去 ，得到 这是中心位于纵坐标轴上的一族圆。第32页，共171页，编辑于2022年，星期六 4 4充分条件充分条件 在前面我们根据泛函取极值的必要条件 和基本预备定理导出了欧拉方程和横截条件。如前所述，只有满足必要
18、条件的函数才能使泛函取得极值。但是，究竟是确使泛函达到极值，还是一个驻点?如果是极值，那么，是极大值还是极小值?仅仅依靠必要条件还不能确定。要回答上述问题还必须研究极值的充分条件。我们知道，普通函数求极值时极值的性质是用该函数在一阶导数等于零的点上二阶导数的符号判断的。如果二阶导数为正则极值为极小值；若二阶导数为负，则极值为极大值；如果二阶导数的符号不定则表明是一个驻点（或鞍点）。同样，要确定泛函极值的性质，还需要研究泛函对参数 关于 的二阶导数。根据前面的分析，对于泛函 有 第33页，共171页，编辑于2022年，星期六这也是一个泛函。利用泛函对参数求导法则，可得 令 ，得到第34页，共17
19、1页，编辑于2022年，星期六或写成二次型积分要使 取得极小(极大)值，还必须 这是使 取极小(极大)的另一个必要条件。如果 则J肯定为极小(极大)值。上述条件是使泛函J取极小值。(极 大握i)的自充分条件。容易看出，如果式(1216)中被积函数内的矩阵是正定或 半正定(负定或半负定)的，只要在区间 上被积函数不 恒等于零，它的积分总大于(小于)零。由此可以得出结论，使泛函 第35页，共171页，编辑于2022年，星期六取极小(极大)的充分条件是：在 条件下，二次型矩阵是正定或半正定（负定或半负定）的。小 结 总结以上分析，得到如下结果：给定泛函 其中 是 的二次可微函数，是变量 和 的连续函
20、数，且有连续二阶偏导数，端点时间 固定。那么，使泛函 取极值的函数 必须满足欧拉方程和横截条件 使 极小(极大)值的充分条件是下列二次型矩阵，是正定或半正定（负定或半负定）的。第36页，共171页，编辑于2022年，星期六 例例123123 试求从点到线间弧长最短的曲线。解：这是一个固定始点、未定终点问题。首先建立泛函的数学表达式，定义微分弧长为 ，则 于是，曲线弧长为 这里 ，因此，这个问题的欧拉方程是 它的一次积分为 c为积分常数。把上式展开，可得 ，把上式再积分，得到极值曲线 积分常数a和b利用横截条件来确定。这个问题的横截条件是第37页，共171页，编辑于2022年，星期六第38页，共
21、171页，编辑于2022年，星期六和即和由此可得阿a0，b1。于是我们得到从点 到线 间具有最短弧长的曲线是直线 ，如图112所示。在这个例子中，因此有 ，二次型矩阵半正定。这就从数学上证明了我们求出的极值确是极小值。第39页，共171页，编辑于2022年，星期六例例1.241.24 求泛函在自由端点条件下的极值曲线。解：这个问题的欧拉方程是 欧拉方程的积分曲线是 横裁条件是 ，或 ，解上面两个代数方程，得 ，。因此，极值曲线为 第40页，共171页，编辑于2022年，星期六第41页，共171页，编辑于2022年，星期六布什总统就美国决定退出反导条约发表电布什总统就美国决定退出反导条约发表电视
22、讲话视讲话布什的讲话：早上好，我们刚刚结束了国家安全委员会的一次会议。我们回顾了我与我的朋友俄罗斯总统普京，在过去几个月的多次会议上的会谈情况。美国人民需要放弃1972年达成的反弹道导弹条约。今天，我已正式通告俄罗斯，美国要从这个存在了将近30多年的条约中退出。我们认为，反弹道导弹条约妨碍了我们政府寻求保护民众的新途径所做的努力，阻碍了保护民众免遭恐怖分子和流氓国家的导弹袭击的努力。.第42页，共171页，编辑于2022年，星期六NMD美國國家導彈防御系統介紹美國國家導彈防御系統介紹改迸的預警雷達，它們是NMD系統的眼睛，能預警到4000-4800千米遠的目標。美國除要改進現有部署在阿拉斯加的
23、地地彈預警雷達以及部署在加州与馬薩諸塞州的鋪路爪雷達外，還要在亞洲地區新建一個早期預警雷達。地基雷達是一种X波段、寬頻帶、大孔徑相控陣雷達，將地基攔截彈導引到作戰空域。地基攔截彈是NMD的核心，由助推火箭和攔截器(彈頭)組成，前者將攔截器送到目標鄰近，后者能自動調整方向和高度，在尋找和鎖定目標后与之相撞，將它擊落在太空上。第43页，共171页，编辑于2022年，星期六美国国防部高级官员近日透露，将从今年夏天开始部署用于导弹防御系统的拦截武器，这一计划比白宫原定的“秋季期限”有所提前。美国导弹防御局发言人里克莱纳说：“也许6、7月前后第一枚拦截导弹将进入发射井。”布什政府的导弹防御系统是克林顿提
24、出的国家导弹防御系统的延伸与扩展，从以陆基中段拦截导弹为主发展到陆基、海基和空基多层次导弹防御体系的结合。美国心急火燎地宣布将提前启动该系统，引起了人们的广泛置疑和揣测美急于部署导弹防御系统的背后美急于部署导弹防御系统的背后第44页，共171页，编辑于2022年，星期六中国日报网站消息：据俄罗斯国际文传电讯报道，俄罗斯国防部长伊万诺夫11月29日告诉俄罗斯总统普京，俄罗斯已经成功测试了一个经过改进的反弹道导弹系统。第45页，共171页，编辑于2022年，星期六机动型“白杨-M”的战术技术性能优势显著。导弹采用主动姿态控制系统等先进制导技术，飞行过程中弹体平稳，不会发生翻转。导弹飞行距离超过1万
25、公里时，其命中精度误差小于90米。此外，为了提升突防能力，“白杨-M”具有额外机动能力，可连续变更自己的位置，变化莫测、防不胜防，且能从机动路途中的任意一个点发射，可以说是美国国家导弹防御系统的克星俄试射新型导弹对抗NMD第46页，共171页，编辑于2022年，星期六发射中的以色列“箭”式反导导弹以色列美国成功进行以色列美国成功进行“箭箭”式式反导系统试验反导系统试验中新网7月30日电据法新社报道，以色列和美国星期四(29日)在加利福尼亚州的一个靶场成功地进行了“箭”式反导系统实弹拉截实验。这种反导系统是目前唯一的一种现役反导系统。第47页，共171页，编辑于2022年，星期六发射升空的以色列
26、“箭”式反导导弹以色列时间下午8点25分，一枚美国从伊拉克没收的飞毛腿”导弹从美国海军的一个空战中心发射升空，“箭”式系统的雷达“绿松”发现并确定了这枚导弹的方位，并引导一枚“箭”式导弹拦截了这枚飞毛腿导弹。第48页，共171页，编辑于2022年，星期六中国日报网站消息：遭遇了一连串的败绩后，美国的导弹防御系统试验终于在2月24日重现一丝曙光。当天美军从夏威夷附近海域的舰艇上发射拦截导弹，一举摧毁了“来袭”的模拟敌方导弹。美国国防部导弹防御局发言人理查德莱纳两个多月来第一次宣布了一条“获胜”的消息。他在五角大楼的新闻发布会上兴奋地表示：“我们成功地进行了一次拦截摧毁防御试验。”第49页，共17
27、1页，编辑于2022年，星期六7月6日，日本防卫厅长官石破茂前往首相办公室参加内阁会议。在当天石破茂向内阁会议提交的2004年防卫白皮书中，“专守防卫”和“禁止出口武器”两大军事限制被束之高阁，并声称要耗资100亿美元在7年内建成导弹防御系统日本计划在日本计划在7年内年内建成弹道导弹防御建成弹道导弹防御系统系统第50页，共171页，编辑于2022年，星期六第51页，共171页，编辑于2022年，星期六第52页，共171页，编辑于2022年，星期六第53页，共171页，编辑于2022年，星期六第54页，共171页，编辑于2022年，星期六三、未定终端时间、无约束条件的变分问题 前一节研究泛函 （
28、13一1）的极值时，都假设初始的间 和终端时间 是固定不变的。这问题称为固定端点时间问题。但是，在最优控制中经常还会遇另一类问题，它的终端时间未规定，称为未定终端时间问题。绪论中曾经提到的最短时间问题就是未定终端时间问题的一个实例。假设终端时间 未规定，但不是任意的，它受终端状态约束 ，而 又取决于给定的终端曲线 。规定了终端状态与终端时间之间的关系。如果终端状态分别是 、，则相应的终端时间就只能分别是 、，如图113所示。第55页，共171页，编辑于2022年，星期六第56页，共171页，编辑于2022年，星期六研究未定终端时间问题的任务是：寻找一条连续可微的极值曲线，它由给定的起始点 出发
29、，”在终端时间到达给定的终端曲线 ，即 (132)使性能泛函(131)达到极小(或极大)。仿照固定端点时间问题分析方法，我们用 表示极值函数或曲线，用 （133）表示 与邻近的一族曲线。这里 是 的变分。再用 （134）表示与曲线族 相对应的一切终端时间的集合。表示与极值曲线或最优轨线相对应的终端时间，称为最优终端时间。是 的变分 将式(133)两边对t求导，得到 （135）第57页，共171页，编辑于2022年，星期六将式(133)一(135)代入式(131)，可得 (136)如前所述，泛函J取极值的必要条件是 对于式(136)，可得 （137）上式等号右边第一项相当于固定端点时间问题，利用
30、第二节的结 果，有 第58页，共171页，编辑于2022年，星期六式(137)等号右边第二项，先用积分中值定理，再让它对 求导，然后令 ，则变为 将式(138)、(139)代入式(137)，再令结果等于零，可得到 (1.310)在 终 端处，终端时间和终端状态受终端界线方程(132)约束，因此，变分 同 彼此相关，只允许其中之一可以任意选择，一但确定了其中的一个，则另一个必须按照方程(132)确定。将式(133)、(134)代入式(1，32)，得 （1311）将上式两边对 求偏导，再令 ，得 （1312）其中，于是得到 同 之间的关系，即为 （1313）下面用图解方法对方程(13l 3)稍加解
31、释。在图114中，第59页，共171页，编辑于2022年，星期六第60页，共171页，编辑于2022年，星期六 表示最优轨线，是一条容许的邻近曲线，由图可知：于是得出这个近似等式就是方程(1313)，它同精确等式之间只相差一个高阶无穷小。将方程(1313)代入式(1310)，得 （1.314）在方程(1314)中，同 互不相关，且任意选 择，也是任意的。也就是说，该式左边分别同 、有关的3项是互不相关的，要使该式成立，必然是上述3项 分别等于零。由此得出 第61页，共171页，编辑于2022年，星期六 （1.315）（1.316）（1.317）方程(1315)是未定终端时间问题的欧拉方程。方程
32、(1316)、(1317)是横截条件，式中 应有连续二阶导数，L至少二次连续可微，连续一阶导数。如果始端固定，则 ，初始条件是 ；若 未规定，则始端条件是 。因为可任意选择，则由方程(1317)可得 此外，在求解这类变分问题时还必须在终端时刻满足约束条件(1.32)。第62页，共171页，编辑于2022年，星期六小小 结结 综合以上讨沦，得到如下结果：如果函数 是在给定初始时间从某个容许初始点出发，到达给定终端约束曲线 上的某一点，并且使性能泛函取极值的极值函数，那么，它同最优终端时间 一起，必须满足：1)欧拉方程 2)横截条件第63页，共171页，编辑于2022年，星期六例131 试求在tx
33、平面上由点 到曲线 具有最短弧长的曲线方程(假设 只取正值)。解：由题意，所示问题的性能泛函是 ，未规定参照例123，这个问题的极值曲线方程是联立求解这3个方程，可得 ；于是得到极值曲线和最优终端时间为 和如图l15所示。第64页，共171页，编辑于2022年，星期六第65页，共171页，编辑于2022年，星期六四、无条件变分问题的进一步讨论四、无条件变分问题的进一步讨论 这一节继续研究无约束条件的变分问题，把它从只有一个自变量函数的情况推广到具有 n个自变量函数的泛函。首先讨论欧拉方程和横截条件的矢量形式，然后讨论变分方法。1欧拉方程和横截条件的矢量形式 到这里为止，我们研究的性能泛函只包含
34、一个自变量函数现在讨论性能泛函包含有多个自变量函数的情况，导出它的欧拉方程和横截条件，然后再用状态空间法把它们写成矢量形式。假设性能泛函 (141)是自变量函数 的函数，每一个 都是独立变量的函数，具有连续的二阶导数，把它写成矢量形式，则 (142)式中 是一n维矢量，“T表示矢量或矩阵的转置。(纯量函数 至少二次连续可微。假设起始时间 固定(也可以不规定)，终端时间 受一族给定的终端界线 约束，具有连续一阶导数。第66页，共171页，编辑于2022年，星期六下面，我们采用同前几节类似的方法来推导这类问题的欧拉方程和横截条件。假设最优轨线是 相应的变分是 于是，邻近的曲线族可以表示成 (143
35、)对上式两边求导，得 (144)把它们写成矢量形式 (1.45)(1.46)其中 假设 是最优终端时间，则 (147)格式(1.43)、(1.44)和(1.47)代入式(1.41)，得 (148)第67页，共171页，编辑于2022年，星期六同前面样，在求泛函 对参数 的偏导数时要利用泛函对参数求导法则，在这里利用这个关系，将方程(148)两边对 求导，再令 ，得 (1.49)将上式右边被积函数中每一个方括弧内的第二项使用分部积分，可得 （1.410）第68页，共171页，编辑于2022年，星期六或写成矢量形式 （1.411）在终端时间、终端状态受终端界线族约束，即 （1.412）这里 。因此
36、，与 彼此相关。把式(145)、(147)代入式(1412)，得到 将上式两端对 求偏导，再令 ，整理以后，得 ，(1413)已知泛函J取极值的必要条件是第69页，共171页，编辑于2022年，星期六将式(1，413)代入式(1411)，并令结果等于零，得到泛函(141)取极值的必要条件：（1.414）式中 的每一个分量以及 都浓此互不相关且任意。利用基本预备定理，得到 （1.415）（1.416）（1.417）以及（1.418）第70页，共171页，编辑于2022年，星期六方程(1415)和方程(1416)一(1418)分别是欧拉方程和横截条件的矢量形式。如果终端时间固定，则 ，由方程(14
37、11)可得出终端横截条件 ，（1419）上式的展开形式是 ，如果 彼此 无关且任意，则上式左边 项之和等于零必然是各项分别等于零，于是得出以下 个条件，即第71页，共171页，编辑于2022年，星期六例141 给定性能泛函 式中 。假设 和 固定，在终端时间 满足 。试确定使泛函J取极值的横截条件。解：已知 固定，由方程(1419)，这个问题的终端横截条件是 （1.420）同时，在终端时间还要求满足因此，同 彼此相关，由方程(1420)不能得出的左边两项分别等于零。令 （1.421）代入约束方程（1.421），得将上式两边对 求偏导，再令 ，得到，第72页，共171页，编辑于2022年，星期六
38、上式中 、之一可以任意选择，把它代入式(1420)，得 已知始端 固定，。于是，得到两端边界条件分别为（1423）和 第73页，共171页，编辑于2022年，星期六例例1.421.42 试求从点 到圆的具有最短弧长的曲线方程。解：由题意，这个问题的性能泛函是 欧拉方程是 和欧拉方程的一次积分是 和或 ，其中 ，由此可得，其中，对以上两式分别积分，得 第74页，共171页，编辑于2022年，星期六利用初始条件 ，可得 和 。利用式(1424)，这个问题的终端横截条件是 ，联立求解这两个代数方程，得或 ，于是，我们得到的极值曲线为 或 图136是这个问题的几何说明。不难看出，从点到 的具有最短弧长
39、的曲线方程是第75页，共171页，编辑于2022年，星期六第76页，共171页，编辑于2022年，星期六2 2变分方法变分方法 到目前为止，我们是用微分方法确定泛函极值的，即依据 建立的泛函取极值的必要条件，用的符号确定极值的性质。在这一小节里，我们要介绍确定泛函极值的变分方法；用泛函 的一次变分 等于零建立极值的必要条件，用二次变分 确定极值的性质。考虑性能泛函 （1425）其中 和 固定。假设 是极值曲线或最优轨线，用 和 分别代表 和 ，并定义为 和 的一次变分，那么 （1426）表示xt(6)邻近的任意族容许曲线，且有 （1427）将式(1426)、(1427)代入式(1425)的右边
40、，并在 ，上将 展开成台劳级数，即第77页，共171页，编辑于2022年，星期六 （1.428）式中HOT表示关于 和 的高阶项。用 表示 的增量，即 （1.429)第78页，共171页，编辑于2022年，星期六定义 的一次变分 为 的线性主部并令 ，得到 （1430）对上式积分号下第二项使用分部积分，整理后得 （1431）式中 任意，采用同第一节类似的分析方法，运用变分法的基本预备定理，可导出 （1.432）和 （1.433）第79页，共171页，编辑于2022年，星期六这两个方程正是变分问题的欧拉方程和横阶条件得出结论：条件：（1.434）同条件是等效的。再定义 的二次变分 为 的二次部分
41、，即 （1.435）将式(1435)同式(1216)相比较，可以看出它们之间的唯一差别仅在于式(1435)右边增加了一个比例系数12。由此可以得出结论：我们也可以利用 的符号来判断极值的性质。用变分方法来建立泛函取极值的必要条件和充分条件比微分方法要简明些。因此，在后面的讨论中我们都采用变分方法。第80页，共171页，编辑于2022年，星期六五、有约束条件的变分问题五、有约束条件的变分问题 前面几节讨论的变分问题都未涉及约束条件。但是，一切实际最优控制问题都存在着各式各样的约束条件，其中包括等式约束，也包括不等式约束。在这一节里，首先讨论有等式约束的变分问题；然后，简短地介绍有不等式约束的变分
42、问题。1有等式约束的变分问题 一切实际控制系统，其运动轨线都必须满足系统本身的运动方程式。这就是说，我们只能从满足系统运动方程式的函数的集合中来挑选使泛函 取极值的 。这就是一类有系统微分方程约束的条件极值问题。我们知道，在处理有等式约束的普通函数求极值问题时有一种比较方便的方法，叫做拉格郎乘子法。它通过拉格朗日乘子约束条件结合到原来求极值的函数，构成一个新函数。于是，在给定约束条件下求原来函致的条件极值问题就转化成求新函数的无条件极值问题。这个思想也可以用来处理有微分方程等式约束变分问题。首先，考虑只有两个自变量函数的泛函，然后再推广到有个自变量函数的情况。第81页，共171页，编辑于202
43、2年，星期六考虑性能泛函 （1.51）假设式中 和 固定，端点条件为 ，（1.52）函数 、满足下列微分方程：，（1.53）我们的任务是：在满足端点约束条件(152)和微分，方程(153)的函数的集合 、中，寻找出使性能泛函(1.51)取极值的函数。如前所述，泛函J取极值的必要条件它的一次变分等于取泛函(151)的一次变分，并令结果等于零，得到 （1.54）这是一个固定端点问题，式(154)中，第82页，共171页，编辑于2022年，星期六如果积分号内 和 彼此无关 且任意选取选，那么，我们就可以使用基本预备定理，得出积分号下两个方括弧内的函数别等于零。但是，这里 和 受方程(15，3)约束，
44、和 有关，只允许其中之一任意选取，因此不能使用基本预备本理。假设 和 是任意一组满足约束方程(153)的容许函数，取方程(1.53)的变分，得到 （1.55）将上式两边乘以待定函数因子(函数拉格郎乘子)，然后在到间积分，得（1.56）第83页，共171页，编辑于2022年，星期六将每个积分号下的第二项分别进行分部积分，并注意到 、都为零，得到 （1.57）将方程(157)的每一项加到方程(154)的相应项，得定义一个辅助泛函 （1.58）则上式可表示成 （1.59）由于 同 彼此相关，只有其中的一个（比如 ）可以任意选取，另一个(比如 )不 是 任 意 的，因 此 暂 还 不 能 对 式(15
45、9)直 接 使 用 基 本 预 备 定 理。那么，怎么办呢?我们知道，函数 中包含有待定函数 现在我们这样来选择 ，使它满足第84页，共171页，编辑于2022年，星期六 （1.510）经过这样选择的 ，就使方程(159)变成式中 是一任意函数，且 ，具备使用基本预备定理的条件，于是得到 （1.511）通过以上分析，可以得出如下结论：如果函数 、是在微分方程(153)约束下使泛函(151)取极值的极值函数，是有关的函数 拉格郎乘子，那么，、以及必须同时满足下列方程，即1）2）（1.512）3）第85页，共171页，编辑于2022年，星期六另一方面，如果利用辅助函数 构成一个新泛函，即那么，可以
46、看做3个自变量函数 、以及 的函数。利用前面几节的讨论结果，可得出使 取得极值的 、和 必须满足下列欧拉方程：1）2）（1.514）3）或第86页，共171页，编辑于2022年，星期六 容易看出，方程组(1512)与方程组(1514)完全相同。这个事实证明：如果通过函数拉格郎乘子 把约束方程(153)结合到原来求极值的泛函式(151)，构成一个如式(1513)所示的新泛函 ，那么，在微分方程（1.53）约束下求原泛函的条件取值问题转换成求新泛函的 无条件极值问题。上述结论虽然是针对只有两个自变量函数的泛函推导出来的但是进一步可以证明，这个结论也适用于具有多个自变量函数的泛函以及端点条件更为复杂
47、的情况。小小 结结总结以上讨论，得到如下结果：如果给出性能泛函 其中 是 维矢量，要求在矢量微分方程 (1516)的约束下求泛函J的极值，其中是m维矢量函数，且m n，那么，我们可以用一个m维矢量拉格郎乘子 将约束条件(1516)结合到泛函(1515)构成一个新泛函，即 (1517)第87页，共171页，编辑于2022年，星期六于是，在微分方程组(1516)约束下求泛函(15l 5)的条件极值问题就转化成求新泛函 (15l 7)的无条件极位问题。假设函数 ，使泛函取极值，那么，这n+m个函数必须满足下面n+m个方程：，（1.518）和 ，或写成矢量形式 和 在这里我们把辅助泛函 作为依赖于n+
48、m个自变量函数 ，的泛函来看待。方程(1520)同方程(1521)一起是它们的欧拉方程。为了说明有微分方程等式约束的变分问题的具体解法，下面分析几个例题。第88页，共171页，编辑于2022年，星期六例例151151 火箭在自由空间里的运动 用下列微分方程描述 图117是它的状态变量图。令 ，可以 建立状态方程 （1.522）其中 ，试求控制函数 ，把系统(1522)从初始状态 经过2秒钟转移到状态空间原点，即 且使性能指标 （1.523）取极小。第89页，共171页，编辑于2022年，星期六第90页，共171页，编辑于2022年，星期六解：解：选择矢量拉格郎乘子 ，把约束条件(1522)结合
49、到性能泛函(1。523)，构成辅助泛函 这里 式中控制函数 可视为泛函J的另一个自变量函数。假设 使取极值，那么，它们必须满足欧拉方程：联立求解上述欧拉方程，可得第91页，共171页，编辑于2022年，星期六利用边界条件可以算出4个积分常数于是得到最优控制和最优轨线，即 图118示出了最优控制函数 和最优轨线 。第92页，共171页，编辑于2022年，星期六例例1.52 1.52 给定线性系统状态方程初始条件为终端条件为 未规定试求控制 和 ，使性能泛函第93页，共171页，编辑于2022年，星期六解：利用拉格郎乘子 ，建立辅助泛函它的欧拉方程是 解欧拉方程，可得 ，第94页，共171页，编辑
50、于2022年，星期六利用边界条件可求出4个积分常数为于是得到最优控制 ，第95页，共171页，编辑于2022年，星期六例例1 15353 给定系统方程端点条件为试求控制 ，使性能泛函 ，未规定取极小。解：利用拉格郎乘子人，建立辅助泛函 它的欧拉方程是 由以上各式可得第96页，共171页，编辑于2022年，星期六由此可得由 ，可得由 ，可得 由 ，得第97页，共171页，编辑于2022年，星期六由 可得 联立求解上述方程，得 于是得到第98页，共171页，编辑于2022年，星期六2有不等式约束的变分问题在最优控制问题中，除存在微分方程等式约束外，还经常遇到不等式约束。比如，在姿态控制问题中，由于