第1章最优控制问题

1第一章最优控制问题最优控制理论是现代控制理论中最早发展起来的分支之一.所谓控制就是人们用某种方法和手段去影响事件及其运动的进程和轨道，使之朝着有利于控制主体的方向发展.对于一个给定的受控系统，常常要求找到这样的控制函数，使得在它的作用下，系统从一个状态转移到为设计者希望的另一个状态，且使得系统的某种性能尽可能好.通常称这种控制问题为最优控制问题.最优控制理论主要讨论求解最优控制问题的方法和理论，包括最优控制的存在性和唯一性和最优控制应满足的必要条件等.最优控制理论始于20世纪50年代末，其主要标志是前苏联数学家庞特里亚金（L.C.Pontryagin）等人提出的“最大值原理”.最优控制理论在工矿企业、交通运输、电力工业、国防工业和国民经济管理等部门有着广泛的应用.本章我们通过几个具体实例介绍最优控制的基本问题和基本概念及其最优控制问题的数学描述.第一节最优控制实例下面列举几个简单但具有实际应用的例子，他们虽然来自完全不同的领域，但却反映了一个共同的问题-最优控制问题例1快速到达问题考虑一个机构（如车皮）W，其质量为m，沿着水平的轨道运动，不考虑空气的阻力和地面对车皮的摩擦力，把车皮看成一个沿着直线运动的质点，x(t)表示车皮在t时刻的位置，u(t)是施加在车皮上的外部控制力，假定车皮的初始位置和速度分别为00)0(,)0(yxxx,我们希望选择一个控制函数u(t)使车皮在最短时间内到达并静止在坐标原点，即到达坐标原点时速度为零.根据牛顿第二定律0),(ttuxm，(1.1)令211,xxxx，则（1.1）化为),(,221tuxmxx(1.2)其中)(),(21txtx分别表示车皮在t时刻的位置和速度，写成向量形式2,),()0(),(00TyxXtbuAXX(1.3)其中)()()(,/10,001021txtxtXmbA.由于技术上的原因，外部推力不可能要多大有多大，它在数量上是有界的，即Mtu)(，（1.4）其中M是正常数.101))((ttdttuJ（1.5）问题是寻找一个满足条件（1.4）的控制函数)(tu，把W由初态Tyx),(00转移到终态T)0,0(，且使性能指标（1.5）达到最小.任何能达到上述要求的控制函数，都称为最优控制.电梯的快速升降、轧钢机的快速控制和机械振动的快速消振问题都可以用上述问题阐述.例2国民收入的增长问题国民经济收入主要用于两个方面：扩大再生产的积累资金和满足人民生活需要的消费基金.我们的问题是如何安排积累和消费资金的比例使国民收入得到最快的增长.我们用x(t)表示t时刻的国民收入，y(t)表示用于积累基金的部分，)()()(txtytu称为积累率.我们的目的是寻求最优积累率u(t)，使国民收入x(t)增长最快.国民收入的增长率x取决于当时的收入总值x(t)和积累率u(t),即有),,()(uxtftx.（1.6）根据u(t)的实际含义，)(tu应满足1)(0tu.考虑一段时间T（5年或10年），使x(t)从初值0x达到尽可能大的x(T),即0)0(xx，（1.7）)(maxTx.（1.8）3问题归结为在（1.6）、（1.7）下求u(t)满足（1.8）。这等价于它的对偶问题：在固定端点条件1)(xTx下使T最小。记10)(,)0(xTxxx，（1.9）TdtuJ0)(.（1.10）国民收入的增长问题的一般提法是在条件（1.6）、（1.9）下求u(t)使性能指标（1.10）达到最小.例3登月艇的软着陆推力控制问题登月艇以最小能耗在月球表面进行软着陆的推力控制问题，可以用下列经过简化的问题表示.将登月艇视为一质点，用)(th代表艇距月球表面的距离，以远离月表面的方向为正方向，)(tv表示艇的速度，)(tu表示艇上火箭的推力，g表示月球上的重力加速度，k为一给定常数，)(tm为艇的质量，0h和0v分别表示艇在初始时刻的高度和速度，M表示艇在不装燃料时的质量，F代表燃料的初始质量.登月艇的运动方程为.,,kummugvvh（1.11）初始条件FMtmvtvhth)(,)(,)(00000.（1.12）终端条件0)(,0)(11tvth，（1.13）这是安全着陆所要求.作为控制函数的推力)(tu满足约束atu)(0，（1.14）其中a为艇的火箭所能达到的最大推力.要求燃料最少就是使10)()()(10ttdttmtmtm（1.15）4取最小，其中终止时刻1t是待定的.（15）等价于泛函10)()(ttdttuuJ.（1.16）登月艇的软着陆推力控制问题可以表述为：求满足（1.11）-（1.14）的)(tu使得性能指标泛函（1.16）取最小值.例4基金的最优管理问题某基金会得到一笔60万元的基金，现将这笔款存入银行，年利率为10%，该基金计划用80年，80年后要求只剩0.5万元用作处理该基金会的结束事宜.根据基金会的需要，每年至少支取5万元至多支取10万元作为某种奖金.我们的问题是制定该基金的最优管理策略，即每年支取多少元才能使基金会在80年中从银行取出的总金额最大。令x(t)表示第t年存入银行的总钱数，u(t)表示第t年支取的钱数.则,1.0),()()(rtutrxtx（1.17）x(0)=60，（1.18）x(80)=0.5，（1.19）根据基金会的需要，每年至少支取5万元至多支取10万元，因此,10)(5tu（1.20）基金会在80年中从银行取出的总金额为.)()(800dttuuJ（1.21）基金的最优管理问题是求满足（1.17）-（1.20）的u(t)使（1.21）中的J(u)取最大值.例5化学工程问题设)(,),(1txtxn表示反应器中n种物质在t时刻的数量（浓度），这n种物质在反应器中同时发生化学反应.设反应速率由微分方程组,,,1,)0()),(),(,,,(01nixxtptxxgdtdxiinii（1.22）来决定，其中)(),(tpt分别表示t时刻反应器内的温度和压力，我们能够在每个时刻控制温度和压力，而温度和压力满足约束条件,)(,)(2121ptppt（1.23）其中2121,,,pp都是常数，分别表示温度和压力可以达到的最小值和最大值.5假设反应已持续进行了一段时间T，T时刻的数量（浓度）为)(,),(1TxTxn.每种物质的价格为.,,1,niRci因此，最终产品的价值是niiiTxcpJ1).(),(（1.24）问题是寻求],0[T区间上的满足（1.22）和（23）的分段连续函数)(),(tpt，使得（1.24）中的),(pJ取最大值.例6生产计划问题设有若干台同样的机器，每台机器可以做两种工作，如果用于做第一种工作，每台每年可获利3万元，机器的损坏率为2/3；如果用于做第二种工作，每台每年可获利2.5万元，机器的损坏率为1/3。现考虑3年的生产周期，试确定如何安排生产计划可获得最大利润。设第k年可用机器的台数为)(kx，第k年分配做第一种工作的机器台数为)(ku，显然，第1k年可用机器的台数)1(kx满足状态方程)],()([32)(31)1(kukxkukx（1.25）)(ku满足约束条件),()(0kxku（1.26）第k年获得的利润为),(5.0)(5.2))()((5.2)(3)(kukxkukxkukR3年一共获得的利润为2020)](5.0)(5.2[)())((kkkukxkRkuJ，（1.27）生产计划问题是寻求满足状态方程（1.25）和约束条件（1.26）的)2,1,0)((*kku，使（1.27）中的目标))((kuJ达到最大.6第二节最优控制问题的数学描述前面所举的例子，虽然它们来自不同的领域，但是在数学上它们都有共同的表述。我们将这些最优控制问题作如下数学描述.2.1控制系统的状态方程控制系统的状态变量是指对事件及其运动起决定作用的量.控制系统的控制变量是指对事件及其运动起控制作用的量.控制系统的状态方程是指描述系统及其运动的方程，其中包含控制变量和状态变量.令nTnRxxx),,(1表示控制系统的状态变量，mTmRuuu),,(1表示控制系统的控制变量，It通常表示时间，Tnfff),,(1是IUX上有定义的n维向量函数.则控制系统的状态方程通常用一阶常微分方程组))(),(,(tutxtfx（2.1）来描述.当f不显含t时，称（1）为定常系统（或称为时不变系统）.当f关于x和u为线性时，称（1）为线性系统.这时方程可以写成utBxtAx)()(，（2.2）其中)(tA为n阶方阵，)(tB为n行m列矩阵。当A和B与时间t无关时，称（2.2）为线性定常系统或称为线性自治系统。2.2终止状态的目标集一般来说控制系统的初始时刻0t和初始状态)(0tx是给定的.但对控制系统的终止时刻1t和终端状态)(1tx来说，却因问题不同而有不同的要求，通常要求达到一个确定的目标集}0)),((,0)),((,)(:)({11211111ttxhttxhRtxtxAn.2.3容许控制函数集在实际问题中，控制变量通常是某种物理量，需要满足有界性等条件，满足这些条件的控制函数，称为容许控制函数，他们全体构成一个集合，称为容许控制函数集，记为}:)({并满足某些条件mRutu，书中，我们通常要求控制函数是分段连续的.72.4性能指标性能指标是指人们对某个控制过程及其结果作出评价的衡量尺度或标准.在数学上用泛函表示，主要有下面三种形式1）终端型性能指标也称麦耶（Mayer）型性能指标)),(()(11ttxuJ，2）积分型性能指标还称拉格郎日（Lagrange）型性能指标10))(),(,()(0ttdttutxtfuJ，3）混合型性能指标也叫包尔查（Bolza）型性能指标10))(),(,()),(()(011ttdttutxtfttxuJ.2.5最优控制问题的数学描述所谓最优控制问题就是寻求一容许控制)(tu，使系统的状态从给定的初值0x在终止时刻)(01tt转移到目标集A，并使性能指标)(uJ取最大值（或最小值）.若上述最优控制问题有解)(*tu，则)(*tu称为最优控制函数，相应的轨线)(*tx叫做最优轨线，而这时的性能指标叫做最优性能指标.2.6离散系统最优控制问题的数学表述在一些实际问题中，系统的状态变量和控制变量关于时间变量是离散的，这样的控制系统称为离散控制系统.令1001,,1,,))(,),(()(kkkkRkxkxkxnTn表示系统的状态变量，1,,1,,))(,),(()(1001kkkkRkukukumTm表示系统的控制变量.则离散控制系统的状态方程为差分方程,1,,1,)),(),(,()1(100kkkkkukxkkx其中，Tn),,(1.终止状态的目标集为},,2,1,0))(,(:)({111ljkxkNkxAj，容许控制集为})(:)({并满足某些条件mRkuku，性能指标为8)).(),(,()),(()(101110kukxkfkkxuJkkk离散系统的最优控制问题就是寻求一容许控制)(ku，使系统的状态从给定的初值00)(xkx在终止时刻)(01kk转移到目标集A，并使性能指标)(uJ取最大值（或最小值）.若上述最优控制问题有解)(*ku，则)(*ku称为最优控制函数，相应的轨线)(*kx叫做最优轨线，而这时的性能指标叫做最优性能指标.