动态最优化 徐高的笔记
T ∂F T dV (ε ) =∫ dt = ∫ Fy p (t ) + Fy′ p ′(t ) dt = 0 0 ∂ε 0 dε
[
]
(2.1.3)
又由分部积分法可得
∫
T
0
Fy′ p ′(t )dt = Fy′ p (t ) 0 − ∫ p (t )
T 0
[
]
T
T d d Fy′ dt = − ∫ p (t ) Fy′ dt 0 dt dt
F (t , y, y ′) = 0 ,则该积分收敛。
条件 II:
∫
∞
0
F (t , y, y ′)dt 中,如果被积函数具有形式 G (t , y, y ′)e − ρt ,其中 ρ > 0 且 G 有界,
则该积分收敛。 横截条件
[F − y ′F ]
y ′ t →∞
∆T + Fy′
[ ]
t →∞
∫ F dt = F
y
(2.1.7) (2.1.8)
Fy′y′ y ′′(t ) + Fyy′ y ′(t ) + Fty′ − Fy = 0
特殊情形下的欧拉方程
2
XG’s 动态最优化笔记
情形 I: F = F (t , y ′) ,解为
Fy′ = 常数
情形 II: F = F ( y , y ′) ,解为
s.t.
m
g (t , y1 ,L, y n ) ≤ c m
F = F + ∑ λi (t ) ci − g i
i =1 m
M
(
)
(2.5.6)
由欧拉方程
Fyj −
及互补松弛条件 可求得最优值。 等周问题
λi (t )(ci − g i ) = 0 , t ∈ [0, T ] , (i = 1,2,L, m)
Fyj −
不等式约束
d F y′ = 0 t ∈ [0, T ] dt j
T
( j = 1,2, L, n)
(2.5.4)
′ , L, y n ′ )dt V [ y1 ,L, y n ] = ∫ F (t , y1 ,L, y n , y1
0
(2.5.5)
g 1 (t , y1 ,L, y n ) ≤ c1
y ′ t =T
∆T + Fy′
[ ]
t =T
∆yT = 0
(2.2.10)
一般的横截条件(2.2.10)可以写成各种具体形式 特殊横截条件 情形 I:垂直终结线(固定时间水平问题) : 由 ∆T = 0 , F − y ′Fy′
[
]
t =T
∆T = 0 自动满足,可得(自然边界条件)
(2.2.11)
XG’s 动态最优化笔记
动态最优化1
徐 高2
1. 预备知识
1.1 目标泛函
标准问题:
V [ y ] = ∫ F [t , y (t ), y ′(t )]dt
T 0
(1.1.1)
迈耶问题(终结控制问题) :
V [ y ] = G[T , y (T )]
博尔扎问题:
(1.1.2)
V [ y ] = ∫ F [t , y (t ), y ′(t )]dt + G[T , F − λG , ( λ =常数)
得必要条件欧拉方程为
(2.5.10)
6
XG’s 动态最优化笔记
Fy −
d Fy′ = 0 dt
(2.5.11)
3. 最优控制理论
3.1 最大值原理
最优控制的最简单问题
max V = ∫ F (t , y, u )dt
0
T
(3.1.1)
s.t.
& = f (t , y, u ) y y (0) = A ,, y (T ) 自由(A,T 给定) u (t ) ∈ u ,对所有 t ∈ [0, T ]
∆yT = 0
(2.4.4)
两项中的任一项都必须单独趋于零。
5
XG’s 动态最优化笔记
2.5 约束问题
等式约束
′ , L, y n ′ )dt V [ y1 ,L, y n ] = ∫ F (t , y1 ,L, y n , y1
0
T
(2.5.1)
g 1 (t , y1 ,L , y n ) = c1
2. 变分法
2.1 变分法的基本问题
基本问题
max(min)V [ y ] = ∫ F [t , y (t ), y ′(t )]dt
T 0
(2.1.1)
s.t.
y (0) = A (A 给定) y (T ) = Z (T,Z 给定)
一阶条件(欧拉方程)的推导 推导的思路是将泛函极值的问题转化为一个一元函数极值的问题, 这样就可以用求驻点 的方法得到极值。 设 y (t ) 为最优函数,则满足约束的任意函数可以表示为 y (t ) + εp (t ) 。其中 p (t ) 是任
[欧拉—泊松方程]
(2.1.16)
2.2 可变端点的横截条件 max(min)V [ y ] = ∫ F [t , y (t ), y ′(t )]dt
T
0
(2.2.1)
s.t.
y (0) = A (A 给定)
y (T ) = yT (T, yT 给定)
一般横截条件的推导 推导的思想仍然是将泛函的极值问题转化为一个一元函数极值的问题。 不同的是需要将 终点也表示为 ε 的函数。
]
(2.2.13)
此式可以通过画一个图看出。详见蒋中一《动态最优化基础》76 页图 3.1 4
XG’s 动态最优化笔记
由于 ∆T 是任意的,可得横截条件为
[F + (φ ′ − y ′)F ]
y ′ t =T
=0
(2.2.14)
再加上 yT = φ (T ) 可确定曲线。 情形 IV:截断垂直(水平)终结线: 。做法是,先按照垂直终结线(水平终结线)方法 有终结约束 yT ≥ y min (或 T ≤ Tmax ) 求出最优曲线。检查是否符合约束,若是,则结束。否则按照固定终点问题 (T , y min ) (或
* * 假设 T 是已知的最优终结时间。则 T 邻近的任何 T 可以表示为
T = T * + ε∆T
其中, ε 为一个很小的数, ∆T 代表任意取定的 T 的小变动。则
(2.2.2)
3
XG’s 动态最优化笔记
dT dT (ε ) = = ∆T dε dε
式(2.2.1)的极值等价于以下函数的极值
(2.2.3)
0
T
(2.1.13) (2.1.14) (2.1.15)
Fyj −
T
d Fy′j = 0 ( j = 1,2, L, n) dt
[欧拉方程组]
V [ y ] = ∫ F (t , y, y ′, y ′′,L, y ( n ) )dt
0
Fy −
n d d2 n d Fy′ + 2 Fy′′ − L + (− 1) F (n) = 0 dt dt dt n y
d F y′ = 0 , t ∈ [0, T ] , ( j = 1,2, L, n) dt j
(2.5.7) (2.5.8)
′ , L, y n ′ )dt V [ y1 ,L, y n ] = ∫ F (t , y1 ,L, y n , y1
0
T
(2.5.9)
′ ,L, y n ′ )dt = k1 G 1 (t , y1 ,L, y n , y1 M s.t. T m ∫ G (t , y1 ,L, y n , y1′ ,L, y n′ )dt = k m
最大值原理 汉密尔顿函数4
H (t , y, u , λ ) ≡ F (t , y, u ) + λf (t , y, u )
(Tmax , yT ) )求解问题。
2.3 二阶条件
充分性定理:对于固定端点问题(2.1.1) ,如果 F (t , y, y ′) 关于 ( y , y ′) 是凹(凸)的, 那么欧拉方程对识别 V [ y ] 一个绝对最大值(最小值)是充分的。 对于可变端点问题, F (t , y, y ′) 的凹(凸)性,加上欧拉方程及横截条件,对识别 V [ y ] 的一个绝对最大值(最小值)是充分的
(2.2.5)
由
∫
T (ε )
0
T ∂F d dt = ∫ p (t ) Fy − Fy′ dt + [ Fy′ ]t =T p (T ) 0 ∂ε dt
(2.2.6)
得问题(2.2.1)的一阶条件为
∫
T
0
d p (t ) Fy − Fy′ dt + [ Fy′ ]t =T p (T ) + [ F ]t =T ∆T = 0 dt
(2.2.7)
由于 (2.2.7) 中的三项均含有任意元素, 故三项必须分别为零。 第一项可导出欧拉方程, 后二项导出横截条件。 又 则
∆yT = p(T ) + y ′(T )∆T 3 p(T ) = ∆yT − y ′(T )∆T
(2.2.8) (2.2.9)
故,一般的横截条件为
[F − y ′F ]
T 0
(1.1.3)
三种问题的转化:令 z (t ) ≡ G[t , y (t )] ,满足 z (0) = 0 ,则 可将迈耶问题与博尔扎问题均转化成为标准问题。
∫
T
0
z ′(t )dt = z (t ) 0 = G[T , y (T )] 。
T
1.2 求导一个定积分
I ( x) ≡ ∫ F (t , x)dt