最大值原理及其应用
内蒙古大学本科学年
第1页
最大值原理及其应用实例
[摘要]:最优控制中的最大值原理,是在目标泛函的最大化问题中得到最优控制的必要条件是使哈密顿函数达最大值而得名的。它被广泛应用于开放式捕鱼以及日常实际问题求最优策略的解决过程中,但是虽然它解决了古典变分法所遇到的困难,给出了最优控制问题解的必要条件,却绝非充分条件,所以在应用中具有一定局限性。
[关键词]:最大值原理;哈密顿函数;必要条件;最优控制;最优策略
正文:
1(最大值原理的内涵
设系统的状态方程为
xtx(),,xfxut,(,,)00
m控制u属于R中的某个有界闭集U,最优控制问题是求u?U,使得
tf
JuxtfLxutdt()((),)(,,),,,f,t0
最小。
假设f(x,u,t)的分量为f(x,u,t),并假设 i
,fxut(,,),Lxut(,,)i,fxut(,,),Lxut(,,)i,xfxut(,,),xLxut(,,)ji,t,ti,,,,,,
,,,,((),)((),)xttxttffff((),),,,(,1,2,,)xttijn,?,ff,,xtt()jff
***都是其自变量的连续函数。用u、x分别
示最优控制和最优轨迹,则u使J(u)取最小值的必要条件是:
** * 1.1存在协状态向量λ(t),它和x(t),u(t)一起满足正则方程
内蒙古大学本科学年论文 第2页
,xfxut,(,,),,,,H,,,,,,x,
* 1.2哈密顿函数作为u的函数在u=u(t)取最小值,即
,,,,,HxutHxut(,,,)min(,,,),,,uU,
1.3 正则方程的边界条件:
1.3.1 若x(t)=x是给定的,则边界条件为 ff
xtx(),xtx(),ff00
1.3.2 如果t给定,x(t)自由,那么边界条件为 ff
,,((),)xttff()t,,xtx(), f00,xf
1.3.3 如果t也是自由的,还要加一个条件: f
,,((),)xttffHxtuttt((),(),(),)0,,, ffff,tf以确定t。 f
1.3.4 如果要求x(t)落在m维流型S上, f
Nxtt((),)0,,1ff,Nxtt((),)0,,2ffS: ,?,
,Nxtt((),)0,rff,
那么边界条件为
xtx(),,00,T,,,((),)((),)xttNxttffff,()tv,,,f,,,xtxt()() ff,
,Nxtt((),)0,ff,T,vvv(,,)?1r,
如果t是自由的,再增加条件 f
T,,Nxttxtt((),)((),),ffffHxtutttv((),(),(),)0,,,, ffff,,ttff可以看出,上述的正则方程和边界条件与u无约束的情况用变分法导出的完全相同。
内蒙古大学本科学年论文 第3页
注:uU, 1 如果最优控制问题是求,使得目标函数J(u)最大,在最大值原理中,
*最优控制u应使哈密顿函数值最大,即
,,,,,HxutHxut(,,,)max(,,,),,, uU,
最大值原理就是在目标泛函的最大化问题中得到最优控制的必要条件是使哈密顿函数达最大值而得名的。
同样,由于目标函数的最小化问题中得到最优控制的必要条件是使哈密顿函数达最小值,也有人称这个原理为最小值原理。
注: 2哈密顿函数作为u的函数,在u=u*(t)处取最小值,因此式
,,,,,HxutHxut(,,,)max(,,,),,, uU,
*对区间[t,t]上的所有t成立,其含义是:对于由最优控制u=u*(t)引发的x(t)0 f
*和协状态λ(t),哈密顿函数作为u函数,在u=u*(t)处取最小值。
**在u为标量的情况,必要条件的含义是:在特定的时刻t和特定的x(t)、λ(t)曲线,哈密顿函数的右端可以看做u的函数,最优控制u*(t)使它取最小值
min()Hu uU,
当U为闭区间[a,b]时,允许的控制满足
aub,,
哈密顿函数仅作为u的函数,有如图所示的三种情况:
图1 关于哈密顿函数的示意图
第一种情况:哈密顿函数的最小值在区间[a,b]内的点达到。如果哈密顿函数关于u是可微的,则必要条件为
内蒙古大学本科学年论文 第4页
,H,0 ,u
其他两种情况是哈密顿函数的最小值在区间[a,b]的边界点达到,这时则必须用更
广泛的条件
min()Hu uU,
也就是用条件
,,,,,HxutHxut(,,,)min(,,,),,, uU,
来描述。
当x(t)是自由的时,应用最大值原理求最优策略的具体步骤如下: f
第1步:构造系统的哈密顿函数为
THxutLxutfxut(,,,)(,,)(,,),,,, 第2步:由
,,,,,HxutHxut(,,,)max(,,,)(,,,对于最大化)uU, ,,,,,HxutHxut(,,,)min(,,,)(,对于最小化),,uU,
导出决策变量u与状态变量x、协状态变量λ的关系,记为u=u(x,λ)。
第3步:写出以下正则(正规)方程组为
,xfxutxtx,,(,,)(),,00,,,((),)xtt,H ,ff,,,,,()t,,f,,,xxt()f,
**将u=u(x,λ)带入正则方程,解出x=x(t),λ=λ(t)。
**第4步:将x=x(t),λ=λ(t)代入u=u(x,λ)得到最优策略
**,,uuxttut,,((),())(),
如果决策问题还要求满足边界条件x(t)=x,则以这个边界条件取代正则方程组中的条ff
,,((),)xttff[1]()t,,件。 f,xf
2(最大值原理在开放式捕鱼中的应用
2.1.应用最大值原理建立模型
内蒙古大学本科学年论文 第5页
所谓开放式捕鱼是指任何渔船都可以任意捕捞。当然,要使捕渔业有长期利润,这就需要考虑怎样控制每年的捕鱼量才能有利于鱼的繁殖,使得在一定时间内鱼的产量最高,也就是使捕渔业赚得的经济收入总数为最大。为了讨论这个问题,首先建立如下数学模型:
设F(x)表示一个给定单种群鱼群的自然增长率,用h(t)表示收获率,x为鱼群的种群密度,t为时间,则鱼群的增长函数为:
dx,,Fxht()()t,0 (1) dt
赚得的经济收入总数为:
,,t,,JepCxFxxdt,,,[()][()] (2) ,0
其中,δ为经济效益指标(折扣率),p为所捕鱼的单位价格,c(x)为单位成本函数。
所以原问题是要求一个最优收获率h*(t)来使经济收总数J最大。文献用欧拉方程求解最优收获率h*(t)。
设
,,,,,,,epCxFxx[()][()]
那么由欧拉方程:
,,,,d, ,,,xdtx
则可得到方程:
'CxFx()()',,,Fx() (3) ,pCx()
若方程(3)有唯一解x=x*,则x*就为所求的最优种群密度(能使经济收入达到最大且
h。若用不会使鱼群灭绝的鱼群密度)表示最大可行收获率,则 max
*h,xx,max,***htFx()(), ,xx,
,*0xx,,
这种方法比较初等,应用范围很小,因为仅当控制变量h(t)无界,控制变量h(t)可以在m维空间中任取,才能在欧拉方程中取δh?0,而若控制变量h(t)有界时,欧拉方程则无法求解此最优控制问题,因而本文引入最大值原理来求解此问题。
内蒙古大学本科学年论文 第6页
2.2应用最大值原理来求解该问题
首先建立如下状态系统:
,xFxht,,()()状态方程
xx(0),初始条件 0
,,t,JepCxhtdt,,[()]()目标函数 ,0
由最大值原理,设λ(t)是伴随向量,则可确定哈密顿函数为:
,,tHxhtepCxhttFxht(,,,)[()]()()[()()],,,,,,
则协态方程为:
,H,,t'',,,()()()()teCxhtFx,,,, (4) ,x
耦合方程为:
,H,,t[()],,,,epCx ,h
由于控制变量h(t)在哈密顿函数H(x,h,λ,t)中线性的出现,说明该问题可以有一个奇异最优控制。而单凭最大值原理来求解奇异问题是不行的,下面引入广义的Legendre-Clebsch条件来对最大值原理进行补充。有了这个补充,奇异控制问题就迎刃而解了。
现在来看原问题:
,H,,t[()]0,,,,,epCx ,h
可得
,,t,,,epCx[()] (5)
H,0而 uu
dH,,,,,tt',,()[()](),,,,,,,epCxeCxx dth,
将(4)式代入上式可得
dH,,,,,,,ttt'''()[()]()[()()]()()(),,,,,,,,,epCxeCxFxhteCxhtFx dth,
,,,,tt'',,,,,,,,epCxeCxFxFx[()]()()()0
内蒙古大学本科学年论文 第7页
,,dH[()]0, ,,hdth
2dH,,,,,,,,,tttt2'',()[()]()()(),,,,,,,,epCxeCxxeCxFxe 2dth,
''''''',,t,,,,CxFxxeCxFxxFxFxx()()()()()(),,,,, (6)
将(4)式代入上式可得
2dH,,,,,,,ttt2''()[()]()[()()]()(),,,,,,,,epCxeCxFxhteCxFx 2dth,
,,,,,,ttt'''',,,,,,,,eCxFxFxhtepCxFxepCxFxFxht()()[()()][()]()[()]()[()()0,
(7)
2,,dH,,,,,,ttt''''(())()()()()[()](),,,,,eCxFxeCxFxepCxFx 2,,hdth
,,,,tt''',,,,eCxFxepCxFx()()(1)[()](),
只要δ<1,p>C(x),就有
2,,dH(())0, 2,,hdth
由(6)式可知:
'CxFx()()',,,Fx() (8) ,pCx()
由(7)式可知h(t)是一个关于C(x),F(x),F′(x),F″(x),δ的函数。只要由(8)式得出最优种群密度x*的表达式,就可以直接求出最优收获率h*的表达式。
2.3应用分析
下面以H.S.Mohring近年来提出的太平洋大比目鱼群的参数:
r=0.71,K=80.5×106kg
对Schaefer模型使用本文提出的最大值原理方法来研究一下关于捕鱼的最优收获问题。在最优控制问题中,只有某些简单的情况可以获得解析解,而绝大多数情况都只能获得数值解。在此,为了能得出x*和h*的解析式以便具体分析问题,不妨假设每单位重量所赚钱为常数α>0,即p-C(x)=α。以下为所给的状态系统:
内蒙古大学本科学年论文 第8页
xx,xrxhtxht,,,,,,(1)()0.71(1)()状态方程 6K80.510,
0.716,,,,xxht(80.510)() 680.510,
,,t,Jehtdt,,()目标函数 ,0
解:设λ(t)为伴随向量,则哈密顿函数为
0.71,,t6Hehttxxht,,,,,,,()()[(80.510)()] 6,80.510协态方程为:
,H0.710.716,,,,,,,,,,xxx(80.510)] 66,,,x80.51080.510
1.42,,,[0.71]x 6,80.510必要条件为:
,H,,,,tt,,,,0,,,,,ee (9) ,h
由于这是一个奇异控制问题,则要想使目标函数取最大值,还应满足
Legendre-Clebsch条件(此条件仅为必要条件)。
dH,,,t(),,,,,e dth,
1.42,,t,,,,,,,ex[0.71]0 (10) 6,80.510
,,dH[()]0, ,,hdth
2dH,1.421.42,,t2,,,,,,()[0.71],,,,exx 266dth,,,80.51080.510
1.420.711.42262,,,,,,ttt,,,,,,,,,,,,,eexxhtxe[(80.510)()][0.71] 666,,,80.51080.51080.510
,0
2,,dH1.42,,t,[()]0,,,e 26,,,hdth80.510由(9)式和(10)式可知:
1.42,,,,tt,,,exe,,,,[0.71]0 6,80.510
内蒙古大学本科学年论文 第9页
680.510(0.71),,,*,,x (12) 1.42
由(11)式可知:
220.71,,*6h,,,80.510 (13) 2.84
由(12)式可以看出,当δ=0时,产生一个最优种群密度x*=40.25×106kg,且产生了最大经济产量h*=14.29×106kg。随着δ的上升,最优种群密度x*下降,当δ>0.71,即折扣率δ大于内禀生长率r时,那么最优种群密度x*=0,在这种情况下,最优收获策略会使资源种群快速灭绝。所以,要使捕渔业获得长期利润,折扣率必须满足0?δ<0.71。
2.4
从过程上来看,最大值原理显得比较复杂,但它的应用范围很广,可以解决多种群鱼群,控制变量受限以及奇异最优控制问题。在开放式捕鱼中,人捕鱼会受到很多限制,
[2] 而欧拉方程又无法解决这些控制受限问题,因而最大值原理就比它适用得多。3(最大值原理应用实例举例
例3.1基金的最优管理问题
基金会达到一笔60万元的基金,现将这笔款存入银行,年利率为10%,该基金计划用80年,80年后要求只剩0.5万元用作处理该基金的结束事宜。根据基金会的需要,内年至多支取10万元作为某种奖金。我们的问题是制定该基金的最优管理策略,即每年支取多少元才能使基金会在80年中从银行取出的总金额最大。
解:用x(t)表示第t年存入银行的总钱数,u(t)第t年支取的钱数,则该问题的状态方程为
,xtrxtutr()()(),0.1,,,,
初值和终值分别为
xx(0)60,(80)0.5,,
控制约束条件为
内蒙古大学本科学年论文 第10页
5()10,,ut
目标泛函即性能指标为
80
Juutdt()(), ,0
从而,基金的最优管理问题就是求满足约束条件的u(t)使J(u)取最大值。用最
大值原理求解这个最优控制问题。
哈密顿函数为
Hxuurxurxxu(,,)()(1).,,,,,,,,,,
*根据最大值原理,u(t)应使哈密顿函数达到最大值,因此
5,10,,,,,,ut(), ,10,10,,,,,
状态方程和辅助和方程为
,xtrxuxx(),(0)60,(80)0.5,,,,,
,H,,,,,,,r ,h
于是
,rt,,()(0).te,
,rt,,(0)1,,1()1(0)0,,,,,,teut()如果,那么,由的表达式可治
,utt()10,080,,,,
,(0)1,,,(0),,()1,,这与实际不符,从而因此将由大于1单调下降到小于1,设则
最优管理策略为
5,0,,,t,,,ut(), ,10,80,,,t,,
于是由状态方程可得
10*0.1rttxtcecet,,,,,,,()100,80, r
x(80)0.5,,由边界条件得,
内蒙古大学本科学年论文 第11页
,,88cee,,,,(0.5100)99.5,
因此
10*0.180.1rttt,,xtcecet,,,,,,,,,,()10099.5100,80. r
x(0)60,,由状态方程和初始条件得
*0.1,txtet()1050,0,,,,,,
于是
0.1t,10100,e,0,,,t,,*xt(), ,,,80.1t,,t80.,,99.5100,,,
由连续性,有
0.180.1tt,,105099.5100,e,,,,
解得
50,,,10ln16.06. ,8,1099.5
故最优策略为
016.06,,,t5,,,ut(), ,16.0680,,,t10,,
即最优管理策略是:在16年以前每年支付5万元,16年以后每年支付10万元,共
支取720万元。
例3.2生产与消费问题
设x(t)表示总产值,u(t)表示用于扩大再生产的投资比例,0?u(t)?1。1-x(t)
表示用于消费的比例,假设产值的增长与投资额成正比,不妨设比例为1,则
,xtuxxx(),(0)0,,,, 0
问题是寻求在一个计划区间[0,T]内的控制u(t),使得消费总额最大。此最优控制
问题的状态方程为
,xtux(),,
xx(0)0,,, 0
xT()自由.
在u?[0,1]下,求使
内蒙古大学本科学年论文 第12页
T
Juxdt,,(1) ,0
,达到最大值的。 u
解: 哈密顿函数为
Huxuxxux,,,,,,(1)(1),,,
t*,0,,,xxxt()0,0(),,,xxtxe由于故,可见。根据最大值原理,u?[0,1]使H00
取最大值,因此,
,()10,t,,1,,,ut(), ,()10.t,,0,,,
,,,,()1t,()1t,t换接函数,换接点满足。
辅助方程为
H,,,,,,1(1)11,uuuuu,,,,,,,,,,,,, h,
,()0.T,
,,()t,()01T,,[,],Tut()0.,由的连续性和知。在末端的一个小区间内应取在上*,,()0T,,()tTt,,u=0,由,知,这个递减的线性函数只能一次达到值1,此时,,,1
*,,,()1,,,,T[1,]TT,,,,T1,从而上。 u,0
*T,1当,则u,0,即全部用于消费上,是一种短期行为。
**(,1),T,,()1t,T,1若,考虑换接的情形,在上u,0(考虑u,0,矛盾),由
Tt,,1,,(1)T,T,1,()te,,,,,,知,这个递减函数只在时取值为1,其余均大于1,故T,1,,0前不能再换接,从而,于是
[0,1],T,1,,,u, ,(1,],TT,0,,
[0,1]T,(1,]TT,即上全部用于投资,上全部用于消费。 4(最大值原理的局限性
内蒙古大学本科学年论文 第13页
最后需要指出的是,最大值原理虽然解决了古典变分法所遇到的困难,但是它也只给出了最优控制问题解的必要条件,而不是充分条件,所以由最大值原理所求的控制函数不一定是最优控制,因为有可能最优控制根本不存在。如果最优控制问题的解存在,但是从这方法得到的控制函数不止一个,就需要进行逐个检验,从中确定出最优解,如果该问题的实际物理背景有最优控制,而从最大值原理得到的解又只有一个,那么这个
[3] 解一定是最优控制。
[参考文献]:
[1] 王翼.经济系统的动态分析[J].北京:机械工业出版社, 2008:108-110.
吴斌.最大值原理在开放式捕鱼中的应用[J].西南工学院学报, 2001年01[2] 张红英, 荆海英,
期 :1-4.
[3] 吕显瑞,黄庆道编著.最优控制理论基础[M].科学出版社.