九州大学学術情報リポジトリ
Kyushu University Institutional Repository
最適化問題と汎関数の極値問題
倉橋, 貴彦
長岡技術科学大学機械系機械情報・制御工学大講座
https://doi.org/10.15017/1470160
出版情報:九州大学情報基盤研究開発センター全国共同利用システム広報. 2 (1), pp.17-28, 2008-08.
九州大学情報統括本部広報委員会 バージョン:
権利関係:
∗
1
!"#$%&'()*&+, -./012345 64#, -./07!34589
:/2;<=4>?=@,9:ABC34DE!.FG()HAIJ342#KLMN%OP
.FQRS2TUVW4?2XNY4>Z[\%01G&'A+4.FQRSN,%O ],^_
`aAb[4?2Xc[@=%^_`aA.d234 eK!()fghi8XIJj=4>?
eK!. FQRS, 64HA. d234kl HRS!mnNY%9:A64ok2p$Jq
34?2NlHRSAkri!estu?2XNY4>vwN,%649:!es(Jj=D
okl HRSAtuDE!x y2j=4ti!]# $z{34>
Z[\% ok
f
X|kx
!esf (x) = x 2 − 2x + 5
2}[@=D~%okf
A.d234x
,%ok
f
x
!o34A)34?2Ndf dx (x) = 2x − 2
2#KX@=%?HX2G4 ?2Nok
f
X.d2G4x
X 4>?ZS~ %okf
A.d234|k,x = 1.0
2Gs%ok,
f (1.0) = 4.0
2G4>O% RS!"W4ok,%okA3|k/X
t
64#,A3 x
%y
%z
8|kok2G$#4?2X #>okA3|k/XG4|k ok2 G$#4/,ok2\=%oklHA E4RS!fp$,%okl HA E4RS2,G4tiX#@=4> oklHRSN,%okA 3|k|QaA. d234eK
G|kX E@=4
1
!fp%oklHRSN,%okA3ok|QaX.d2G4 e
KG okX E@=4>okA 3ok|Qa,|2\=%okl HRS,ok
|X2G4?2N%ok!lHA}[4okX E@=4>?eKGokl HRSA
tui,|i2\=4>vwN,%?oklHRS!A6$%|i!e 4o
klHRS!f 34ti!]#$z{A;K>
2
vN,oklHRSAtu|i!]# $z{34
[1]
>|i!e4N ,%%ok!lHA}[ 4okA]W4DE!% ¡¢JpDok
x(t)
2|Qaδx(t)
%£ L¤ Jk
ε
A#$okAJ (x(t) + εδx(t))
23 >¥¦1
§>?okJ (x(t) + εδx(t))
,¤|kε
¤k Hok2Gs
G(ε)
2¨W4 ¥©(1)
§>G(ε) = J (x(t) + εδx(t)) (1)
??N©
(1)
!"# $%£LJkε
ªA∆ε
2342%ª∆ε
!e4okª ,« ¬eK!¨W4>
G(ε + ∆ε) − G(ε) = J (x(t) + (ε + ∆ε)δx(t)) − J (x + εδx(t)) (2)
D%©
(2)
®,¯°±²A#42%©(3)
eK!j=4>G(ε + ∆ε) − G(ε) = ∆εδx(t) d
dx(t) J (x(t) + εδx(t)) + o( k ∆εδx(t) k )
∗
³´µ¶· ¸¹ ¸º
»¼ º
»½¾¿À Á
¸¹Ã Ä
E-mail: [email protected]
1
ÅÆÇÈÉÊËÌ ÍÎÏÐÑ
f(x)
ÇÒÑx
ÊÓÔÕÒÖ×df dx (x)
ØÙÚ(
Û)
ÜÔÕÝÞßÒÑx
àáâãäÕå1:
ok|Qa= dJ(x(t) + εδx(t))
dx(t) ∆εδx(t) + o( k ∆εδx(t) k ) (3)
?©
(3)
®Aª∆ε
N42©(4)
eK!Gs2
%ª∆ε
AsGu! W4?2Nl,©
(5)
eK!¨u?2XNY4>G(ε + ∆ε) − G(ε)
∆ε = J ′ (x(t) + εδx(t))δx(t) + 1
∆ε o( k ∆εδx(t) k ) (4)
∆ε→0 lim
G(ε + ∆ε) − G(ε)
∆ε = lim
∆ε→0 ( dJ(x(t) + εδx(t))
dx(t) δx(t) + 1
∆ε o( k ∆εδx(t) k ))
= dJ(x(t) + εδx(t))
dx(t) δx(t) (5)
??!%©
(5)
,JqN64DE%G(ε)
k2 p$ 3?2XNY% ©(6)
eK!¨u?2XNY4>
dG(ε)
dε = lim
∆ε → 0
G(ε + ∆ε) − G(ε)
∆ε = dJ(x(t) + εδx(t))
dx(t) δx(t) (6)
??N% ©
(6)
!" # $ε = 0
2342% ©(7)
A4?2XNY4>dG(ε) dε
ε=0 = dJ(x(t))
dx(t) δx(t) (7)
D%©
(1)
£LJkε
!f34A)p%ε = 0
2342©(8)
A4>dG(ε) dε
ε=0 = dJ(x(t) + εδx(t)) dε
ε=0 (8)
©
(7)
%(8)
es©(9)
oA4?2XNY4>dG(ε) dε
ε=0 = dJ(x(t) + εδx(t)) dε
ε=0 = dJ(x(t))
dx(t) δx(t) (9)
2
ÊÏ
J ′
ÎJ(δx(t))
Çδx(t)
ÊÜÌÞå©
(9)
,%64okA£LJkε
Npε = 0
2pD/,%dJ dx(t) (x(t)) δx(t)
2j=4?2ALMp$#4>?oA#$%lHRS!]# $z{34>
% ¡¢JpD
δx(t)
!fp$x(t) = x a (t)
2pD%©(9)
es%«¬eK!¨W 4>dG(ε)
dε = dJ(x a (t) + εδx(t))
dε = dJ(x a (t))
dx(t) δx(t) (10)
?~%okX
x(t) = x a (t)
Nl HA/ %]x a (t)
NN6=\%©(10)
,¤|kε
ok2p$ε = 0
Nl HA24>e$©(10)
es©(11)
eK!¨W 4>dG(0)
dε = dJ(x a (t) + εδx(t)) dε
ε=0 = dJ(x a (t))
dx(t) δx(t) = 0 (11)
©
(11)
!"#$δx(t)
,£LN64 ?2es%©(12)
XCs]>dJ (x a (t))
dx(t) = 0 (12)
«z{es%okX
x(t) = x a (t)
NlHA/ %]x a (t)
NN6=\%okXlHA24~,sc[XNY4>O ],ok
J
A£LJkε
!espε = 0
2pD!l HA]2#K?2N64>/KO]c[,%ok
J (x a (t))
kX2G4?2XlHA/] !G$#4?2N64>?=@?2!] #$% @tNY4 eK«¬z{
A;K>D%©
(11)
,okJ (x(t))
x a (t)
!"W4O|2\=%δJ
2 3?2!34>3
okl HRS!f 34ZS2p$%
P 0 (0, 0)
@P 1 (a, b)
N.A E4RS!]#$z{34 ¥¦
2
§[2]
>P 0 (0, 0)
@P 1 (a, b)
N.A E4RS,%.RS2/\=$#4>64A!"#$%£Ld
ds
,%x
%2:
.RSy
ddx
%dy
A#4 ?2NJes©(13)
eK!j=%ds dt
,©(14)
eK!¨u?2XNY4>ds = p dt 2 + dx 2 (13)
ds dt =
s
1 + dx 2 dt 2
= s
1 + dx(t) dt
2
(14)
??N%64AX.2G4RSAtuDE%©
(14)
!3dA64N pDHA1ok2p$Jq 34 ¥©
(15)
§>J (x(t)) = Z a
0
ds dt dt =
Z a 0
s
1 + dx(t) dt
2
dt (15)
??!JqpD1oklHRS!"#$%1ok
J
,|kx
okN6s%|kx
,|kt
ok2G$#4DE%oklHRS2G4> ??N%?oklHRSAtY%.
2G4ok
x(t)
A u>oklHRSAtuDE!%okx(t)
!fp$£LJkε
2|δx(t)
A[D okAc[4 ¥©(16)
§>J (x(t) + εδx(t)) = Z a
0
s
1 + dx(t)
dt + ε dδx(t) dt
2
dt
= Z a
0
s
1 + dx(t) dt
2
+2ε dx(t) dt
dδx(t) dt
+ε 2 dx(t) dt
2
dt (16)
%1ok
J
Xx(t)
NlHA/ %x(t)
NN6=\%©(16)
,¤|kε
ok2Gs%ε = 0
NlHA24>e$%@=D©(16)
!fp$%x(t)
!"W4O|δJ
)A342% ©
(17)
eK!G4>δJ = ∂J(x(t) + εδx(t))
∂ε
ε=0
= Z a
0
2 dx(t) dt dδx(t) dt +2ε dδx(t) dt 2 2
r
1 + dx(t) dt 2 +2ε dx(t) dt dx(t) dt +ε 2 dδx(t) dt 2 ε=0 dt
= Z a
0
dx(t)
dt
r
1 + dx(t) dt 2
dδx(t) dt
dt
= Z a
0
f (t) dδx(t) dt
dt (17)
??N%©
(17)
!fp$ A)p%|δx(t)
X ¡!"#$2G4?2Ac34 2% ©(18)
A4>Z a 0
f (t) dδx(t) dt
dt = Z a
0
d
dt (f (t)δx(t))dt − Z a
0
d
dt f (t) δx(t)dt
= f(a)δx(a) − f (0)δx(0) − Z a
0
d
dt f (t) δx(t)dt
= −
Z a 0
d
dt f (t) δx(t)dt = 0 (18)
D%£L|
δx
!fp$©(18)
X2G4DE!,%okf (t)
XJkA
2G4xyX64¥©
(19)
§>f (t) =
dx(t)
dt
r
1 + dx(t) dt 2
= A (19)
??N%©
(19)
!3 Atu?2Nt2G4okx(t)
A4?2XNY4>e$%©
(19)
es%. :!©(20)
A4>dx(t) dt = A
s
1 + dx(t) dt
2
dx(t) dt
2
= A 2 (1 + dx(t) dt
2
)
(1 − A 2 ) dx(t) dt
2
= A 2
dx(t) dt
2
= A 2 1 − A 2 dx(t)
dt =
s A 2
1 − A 2 = B (20)
D%©
(20)
es%dx(t) dt
,JkB
2G4xyX6s% ®A 34?2N©(21)
A4>dx(t) dt = B Z dx(t)
dt dt = Z
Bdt
x(t) = Bt + C (21)
??N% ¡
x
HesJkB
%C
AJ342©(22)
%©(23)
eK!G4>x(0) = 0 ...C = 0 (22)
x(a) = b ...B = b
a (23)
e$%©
(21)
!JkB
%C
A 34?2N.2G4okx(t)
,©(24)
eK!@=4>
x(t) = b
a t (24)
©
(24)
esP 0 (0, 0)
@P 1 (a, b)
NA.2G4okx(t)
,2
A N64 ?2X <4>
4
N,%oklHRS!fp$A;DX%vN,%okA 3okHX6
4 !es}[@=4RS!fp$%okA. 234okAIJ34RS!] #$z{3
4
[2]
>okA3ok!f34 , 3
2\=%?RS, VYoklHRS2G4> VYoklHRSZ2 p$%64
2
P 0
%P 1
N ¢JpDj
L
2P 0
%P 1
AN=D!esj=4A.234A E 4RSAc[ 4
4
¥¦
3
§>3:
8RS
P 0
%P 1
,¢Jj=$#4DE©(25)
! 3 /2%j"eXIE@=4>y(x 0 ) = y(x 1 ) = 0 (25)
??N%j
L
,©(26)
!es}[@=%S
HA1ok2Jq34 ¥©(27)
§>L = Z x 1
x 0
s
1 + dy(x) dx
2
dx (26)
S = Z x 1
x 0
y(x)dx (27)
??!%
L
"eP 0
%P 1
A!es =4S
A.!34oky(x)
A E4>A.234okA E4RS,%©
(25)
%(26)
/21ok ¥©(27)
§HA.234ok
y(x)
A E4RS2G4> AD3eK!okA.64#,. d!34RS!fp$,%O!±± Jki2#KiA#4>±±
JkiN,% !f34± ± Jk
5
A 34>1ok!fp
$±± Jk2 !esj=4AVpD34?2N% V
3
!"#Ü$%&äÕå
4
'Ç(ÉÎÊ)*(ÉÜ%&äÕå
5
+,-ÑÎÏ./ÇÀ0 !"#ÊÓÔÕ-Ñ
(
1
Ñ
)
Ü 2Í3ÌãäÏ+4 ×56Õ 'Ü7ã+,-Ñ 86ÕÌÎ+,1
Ñ9Ü%&äÍÌÕå
X2G48 Nj=$#4DE% !±± J kAW%1ok!
[4 ?2N@=41ok
6
lHRS,%.!JqpD1oklHRS281!G4>
1oklHRSAtu?2!es% A34okX 4>?RSN,%
1ok
(
©(27))
!f 34 ,©(26)
N64DE%±±JkAλ
2 p1okAC342©
(28)
eK!G4>?1oklHRSAtu?2N%A34ok
y(x)
A E4?2XNY 4>J ∗ (λ, y(x)) = Z x 1
x 0
y(x)dx + λ Z x 1
x 0
s
1 + dy(x) dx
2
dx − L (28)
??N%©
(28)
!fp$RS2%oky(x)
Aεδy(x)
!es342©(29)
A4>J ∗ (λ, y(x) + εδy(x))
= Z x 1
x 0
(y(x) + εδy(x))dx + λ Z x 1
x 0
s
1 + d(y(x) + εδy(x)) dx
2
dx − L
= Z x 1
x 0
(y(x) + εδy(x))dx + λ
Z x 1
x 0
s
1 + dy(x) dx
2
+2ε dy(x) dx
dδy(x) dx
+ε 2 dδy(x) dx
2
dx − L (29)
??N%©
(29)
es1okoky(x)
!"W 4O|δJ ∗
A)342©(30)
eK!G4>
δJ ∗ = ∂J ∗ (λ, y(x) + εδy(x))
∂ε
ε=0
= Z x 1
x 0
δy(x) + λ 1 2
2 dy(x) dx dδy(x) dx +2ε dδy(x) dx 2 r
1 + dy(x) dx 2 +2ε dy(x) dx dδy(x) dx +ε 2 dδy(x) dx 2
ε=0 dx
= Z x 1
x 0
δy(x) + λ dy(x)
dx
dδy(x)
dx
r
1 + dy(x) dx 2
dx = 0 (30)
D%©
(30)
!fp$A p%¡!"#$£Lokδy(x)
X2G4 Ac342%©
(31)
A4?2XNY4>Z x 1 x 0
dy(x)
dx
dδy(x)
dx
r
1 + dy(x) dx 2
dx = Z x 1
x 0
d dx
dy(x)
dx
r
1 + dy(x) dx 2
δy(x) dx − Z x 1
x 0
d dx
dy(x)
dx
r
1 + dy(x) dx 2
δy(x)dx
=
dy(x)
dx
r
1 + dy(x) dx 2 δy(x)
x=x 1 −
dy(x)
dx
r
1 + dy(x) dx 2 δy(x)
x=x 0
− Z x 1
x 0
d dx
dy(x)
dx
r
1 + dy(x) dx 2
δy(x)dx
= −
Z x 1
x 0
d dx
dy(x)
dx
r
1 + dy(x) dx 2
δy(x)dx (31)
6
ÐÑÜ$%&äÕå
©
(30)
!%©(31)
oA 342©(32)
eK!G4>δJ ∗ = ∂J ∗ (λ, y(x) + εδy(x))
∂ε
ε=0 =
Z x 1
x 0
1 − λ d dx
dy(x)
dx
r
1 + dy(x) dx 2
δy(x)dx = 0 (32)
??N%£L|
δy(x)
!fp$©(32)
X2G4DE! ,% ©(32)
¬ X2G4xyX64 ¥©
(33)
§>1 − λ d dx
dy(x)
dx
r
1 + dy(x) dx 2
= 0 (33)
D%©
(33)
AtuDE!dy(x) dx = tan θ
2342©(34)
eK!|34 ?2XNY4>d dx
dy(x)
dx
r
1 + dy dx 2 = 1
λ
d dx
tan θ
√ 1 + tan 2 θ = 1
λ
d dx
sin θ
cos θ 1 cos θ
= 1 λ d
dx
sin θ = 1
λ (34)
??N% ©
(34)
®!A 32©(35)
A4>Z d dx
sin θ dx = Z 1
λ dx sin θ = x
λ + A 1
λ sin θ = x + A 2 (35)
©
(35)
!"#$θ
!f34A)p342. :!©(36)
A4>d
dθ (λ sin θ) = d
dθ (x + A 2 )
λ cos θ = dx dθ
dx = λ cos θdθ (36)
??N%
dy(x)
dx = tan θ
2JpD,©(37)
eK!¨Yj=4>dy = tan θdx (37)
©
(37)
!fp$©(36)
Ap%342©(38)
eK!G4>dy = tan θ(λ cos θdθ)
dy = sin θ
cos θ λ cos θdθ dy = (sin θ)λdθ Z
dy = Z
(sin θ)λdθ
y = λ cos θ + B 1
λ cos θ = y + B 2 (38)
??N%©
(35)
%(38)
A= = p% Ap<42%L
A3ok,©(39)
eK!G4>
(x + A 2 ) 2 + (y + B 2 ) 2 = λ 2 × (sin 2 θ + cos 2 θ)
= λ 2 (39)
©
(39)
, Ap$"s%S
A.!34L
,Aλ
% HA(A 2 ,B 2 )
234 2G4>
5
v N,%oklHRS!"#$ X !es}[@=4~ !]#$z
{34
[2]
>ZS2 p$%t
!"W464Xx(t)
!es}[@=42 p%A
u(t)
!es34RSAc[ 47
¥¦4
§>|Qdx(t) dt
,%4:
)y7
ÇÎÏÊãÇ Ç!ÊË Ì Í$"ÇÜ 2Ï#$ÊÇ% & ÔÕ$ÇÜÔÕå
2!e$}[@=©
(40)
!esj=4234>dx(t)
dt = − ax(t) + u(t) (40)
©
(40)
!"#$a
,HA24kN 6s%dj#HA(J34~,¬X4 Xu%Y#HA(J34~,¬X4 X #?2ALM34>D%©
(40)
AtuDE!
t 0
2¡t f
!"W 4A©(41)
"e(42)
eK!(J34>x(t 0 ) = x 0 (41)
x(t f ) = 0 (42)
??!%©
(41)
"e(42)
t 0
2%¡t f
,t 0 = 0
%t f = ∞
2(J34>??N%A @¡N N./u !34?2Ac[4>
AAp$¬4DE%!o34 Ac[42% A4
Wb[$A¬4?2Xp#>e$%2!es1
okAJq 34>2HYjA134DE!%HA
pDHA)p7Np<D/A1ok234 ¥©
(43)
§8
>J (x(t), u(t)) = 1 2
Z t f
t 0
(x(t) 2 + u(t) 2 )dt (43)
??N% ¥©
(40)
%(41)
%(42)
§ADpGX@1ok ¥©(43)
§A.d234u(t)
A E4?2!G4>% VYoklHRS!fp$%±±JkiAF p±± JkA
λ(t)
2342%1ok,1ok"e es%©
(44)
eK!3?2XNY4>J ∗ (λ(t), x(t), u(t)) = 1 2
Z t f
t 0
(x(t) 2 + u(t) 2 )dt + Z t f
t 0
λ(t) dx(t)
dt + ax(t) − u(t) dt (44)
D%±± Jk
λ(t)
%x(t)
2u(t)
= =!fp$£LJk
ε 0
%ε 1
%ε 2
A p% ok!fp$ε 0 δλ(t)
%ε 1 δx(t)
%ε 2 δu(t)
A[42©(45)
A4
9
>
J ∗ (λ(t) + ε 0 δλ(t), x(t) + ε 1 δx(t), u(t) + ε 2 δu(t))
= 1
2 Z t f
t 0
((x(t) + ε 1 δx(t)) 2 + (u(t) + ε 2 δu(t)) 2 )dt +
Z t f
t 0
λ(t + ε 0 δλ(t)) d(x(t) + ε 1 δx(t))
dt + a(x(t) + ε 1 δx(t)) − (u(t) + ε 2 δu(t) dt (45)
??N%1ok
J ∗
Xokλ(t)
NlHA % ]λ(t)
N2342%©(45)
es1okok
λ(t)
!"W 4O |δJ 0 ∗
,©(46)
eK!G4>δJ 0 ∗ = ∂J ∗ (λ(t) + ε 0 δλ(t), x(t) + ε 1 δx(t), u(t) + ε 2 δu(t))
∂ε 0
ε
0 =0
= Z t f
t 0
δλ(t) dx(t)
dt + ax(t) − u(t) dt (46)
8
ÐÑØ ÔÕâÇ
ÇÜ
ÇÎ-5ãäÍËÏ5Ç# "!ÏÇ
Ø"ÞâÏÇ#$ÊÝ 1
Ñ %
2
&à'ÊÓ2Í(7Õ 'ÜÊßÕå'Ç'ÜØ),2Ï*Ê 1ÑÜ2Í %
1 2
&ØÊ(+,5ÐÑØ,-2Ï# ÊË ÌÍ 1
Ñà./0ä %
1
&ÜßÕÝÞÊ2ÍÌÕå9
53Ì Ç+,-ÑÎÐÑ5ß1ÒÑ56ÕâÏÒØ232 Ø"ÍÌßÌàÏ45
Î#$
t
ÊÐÔÕÐÑ56ÕâÏÒδλ(t)
Ø232 Ø"ÍÌÕå©
(46)
(40)
δλ(t)
>
!
(45)
"#$%ε 1
& '((47)
)*>∂J ∗ (λ(t) + ε 0 δλ(t), x(t) + ε 1 δx(t), u(t) + ε 2 δu(t))
∂ε 1
= Z t f
t 0
(x(t)δx(t) + ε 1 δx(t) 2 + λ(t) dδx(t)
dt + λ(t)aδx(t))dt (47)
+,-./0-."123 456 7%89#$:
δx(t 0 )
");
δx(t f )
>)<(47)
=>?@AB C ' (δx(t 0 )
");δx(t f )
DE(48)
)*>Z t f
t 0
λ(t) dδx(t) dt dt =
Z t f
t 0
d
dt (λ(t)δx(t))dt − Z t f
t 0
dλ(t)
dt δx(t) dt
= λ(t f )δx(t f ) − λ(t 0 )δx(t 0 ) − Z t f
t 0
dλ(t)
dt δx(t) dt
= −
Z t f
t 0
dλ(t)
dt δx(t) dt (48)
(48) (47)
FG(49)
)*>∂J ∗ (λ(t) + ε 0 δλ(t), x(t) + ε 1 δx(t), u(t) + ε 2 δu(t))
∂ε 1
= Z t f
t 0
(x(t)δx(t) + ε 1 δx(t) 2 − dλ(t)
dt δx(t) + λ(t)aδx(t)))dt (49)
H IJK:
J ∗
:x(t)
L MN Ox(t)
&PQ(49)
)RHIJK::
x(t)
"S?TδJ 1 ∗
(50)
)*>δJ 1 ∗ = ∂J ∗ (λ(t) + ε 0 δλ(t), x(t) + ε 1 δx(t), u(t) + ε 2 δu(t))
∂ε 1
ε
1 =0
= Z t f
t 0
(x(t)δx(t) − dλ(t)
dt δx(t) + λ(t)aδx(t)))dt
= Z t f
t 0
− dλ(t)
dt + aλ(t) + x(t) δx(t)dt = 0 (50)
(50)
" +,-. /0-. Uδx(t)
3V WR
(51)
XY− dλ(t)
dt + aλ(t) + x(t) = 0 (51)
(51)
)R123 45:Z[\ ]^_ `%a b(52)
)*Y
x(t) = dλ(t)
dt − aλ(t) (52)
cde
(45)
"#$%ε 2
& '((53)
)*Y∂J ∗ (λ(t) + ε 0 δλ(t), x(t) + ε 1 δx(t), u(t) + ε 2 δu(t))
∂ε 2
= Z t f
t 0
(δu(t)u(t) + ε 2 δu(t) 2 − λ(t)δu(t))dt (53)
H IJK:
J ∗
:u(t)
L MN Ou(t)
&PQ(53)
)RHIJK: :
u(t)
"S ?TδJ 2 ∗
(54)
)*YδJ 2 ∗ = ∂J ∗ (λ(t) + ε 0 δλ(t), x(t) + ε 1 δx(t), u(t) + ε 2 δu(t))
∂ε 2
ε
2 =0
= Z t f
t 0
(u(t)δu(t) − λ(t)δu(t))dt
= Z t f
t 0
(u(t) − λ(t))δu(t)dt = 0 (54)
(54)
"+,-. /0-. Uδu(t)
3
(55)
XYu(t) − λ(t) = 0 (55)
(55)
)R 245:Z[\ ]^_`%a b(56)
)*Y
u(t) = λ(t) (56)
(52)
(56)
"\ ]^_a)<89123 45 245JK:L X )*'(8912345-U
(40)
FG\]^_`%a[ X(57)
Yd
dt dλ(t)
dt − aλ(t) +a dλ(t)
dt − aλ(t) − λ(t) = 0 (57)
c
(57)
(58)
)*Yd 2 λ(t)
dt 2 − (a 2 + 1)λ(t) = 0 (58)
(58)
)<89\]^ _`%a :& [ T(59)
)*cYλ(t) = Ae √ 1+a 2 t + Be − √ 1+a 2 t (59)
\]^ _`%a:T
(58)
\ ]^_`%a123 45:(52)
FG(60)
XYx(t) = dλ(t)
dt − aλ(t)
= p 1 + a 2 (Ae √ 1+a 2 t − Be −
√ 1+a 2 t ) − a(Ae √ 1+a 2 t + Be −
√ 1+a 2 t ) (60)
(60)
%A
B
c9(41)
(42)
)<89+,-.");/0- ."S12345 %A
B
9Yc+ ,-.");/0-."S12456
(41)
(42)
)R(61)
(62)
XYx(t 0 ) = p 1 + a 2 (Ae √ 1+a 2 t 0 − Be −
√ 1+a 2 t 0 ) − a(Ae √ 1+a 2 t 0 + Be −
√ 1+a 2 t 0 ) = x 0 (61)
x(t f ) = p 1 + a 2 (Ae √ 1+a 2 t f − Be −
√ 1+a 2 t f ) − a(Ae √ 1+a 2 t f + Be −
√ 1+a 2 t f ) = 0 (62)
+ ,-.");/0-.
t 0 = 0
");t f = ∞
)R(61)
(62)
(63)
(64)
)*Y
x(0) = p 1 + a 2 (A − B) − a(A + B) = x 0 (63) x( ∞ ) = p 1 + a 2 (Ae ∞ ) − a(Ae ∞ )
= Ae ∞ ( p 1 + a 2 − a) = 0 (64)
c
(64)
"e ∞ 6 = 0
√
1 + a 2 − a 6 = 0
*)R%A
A = 0
Y(63)
A = 0
FG(65)
)*Yx 0 = − p 1 + a 2 B − aB
= − ( p 1 + a 2 + a)B (65)
(65)
)R%B
B = − √ 1+a x 0 2 +a
Y%A
B (60)
FG+,-.");/0-."S 45 DE12345
(66)
)*Yx(t) = p 1 + a 2 (Ae
√ 1+a 2 t
− Be −
√ 1+a 2 t ) − a(Ae
√ 1+a 2 t + Be −
√ 1+a 2 t )
= p 1 + a 2 x 0
√ 1 + a 2 + a e −
√ 1+a 2 t + a x 0
√ 1 + a 2 + a e −
√ 1+a 2 t
=
√ 1 + a 2 + a
√ 1 + a 2 + a x 0 e −
√ 1+a 2 t
= x 0 e −
√ 1+a 2 t (66)
)<
(66)
)Rc<1245x(t)
6(41)
FG JK: 4512345
x(t)
:[(67)
)*cYu(t) = dx(t)
dt + ax(t)
= − p 1 + a 2 x 0 e −
√ 1+a 2 t + ax 0 e −
√ 1+a 2 t
= ( − p 1 + a 2 + a)x 0 e − √ 1+a 2 t
= ( − p 1 + a 2 + a)x(t) (67)
5.1
!"#$%
(12)
)R&:L : %O '( )*Y !"#de12345-U
(68)
)R*+123"S+ ,-.
t 0
");/0-.t f
45(69)
(70)
)R*+Ydx(t)
dt = − ax(t) + u(t) (68)
x(t 0 ) = x 0 (69)
x(t f ) = 0 (70)
c,JK:-.de
(71)
)*%/\]^ _`%a%0 1 \ ]^_`%aλ(t)
)RJK: HI(72)
)*YJ (x(t), u(t)) = 1 2
Z t f t 0
(x(t) 2 + u(t) 2 )dt (71)
J ∗ (λ(t), x(t), u(t)) = 1 2
Z t f
t 0
(x(t) 2 + u(t) 2 )dt + Z t f
t 0
λ(t) dx(t)
dt + ax(t) − u(t) dt (72)
(11)
: bHIJK:?T(73)
)* YδJ ∗ =
Z t f t 0
∂J ∗
∂λ(t) δλ(t) + ∂J ∗
∂ dx(t) dt
δ dx(t) dt
+ ∂J ∗
∂x(t) δx(t) + ∂J ∗
∂u(t) δu(t) dt = 0 (73)
c,+,-.");/0 -."S45
x(t 0 )
");x(t f )
7%89δx(t 0 )
");
δx(t f )
Y(73)
? @ABC'( +,-.");/0-."S45 DE[
(74)
)*YZ t f
t 0
∂J ∗
∂ dx(t) dt
δ dx(t) dt
dt = Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt δx(t) dt − Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt
δx(t)dt
= ∂J ∗
∂ dx(t dt f ) δx(t f ) − ∂J ∗
∂ dx(t dt 0 ) δx(t 0 )
− Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt
δx(t)dt
= −
Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt
δx(t)dt (74)
)<
(74) (73)
FG(75)
)*YδJ ∗ =
Z t f
t 0
∂J ∗
∂λ(t) δλ(t) + − d dt
∂J ∗
∂ dx(t) dt
+ ∂J ∗
∂x(t)
δx(t) + ∂J ∗
∂u(t) δu(t) dt = 0 (75)
(75)
)RHIJK:?T,:&(76)
(77)
(78)
X Y∂J ∗
∂λ(t) = dx(t)
dt + ax(t) − u(t) = 0 (76)
− d dt
∂J ∗
∂ dx(t) dt
+ ∂J ∗
∂x(t) = − dλ(t)
dt + aλ(t) + x(t) = 0 (77)
∂J ∗
∂u(t) = u(t) − λ(t) = 0 (78)
(76)
12345 [ "R6(68)
)RYc,
(77)
\ ]^_`%a:&[ R\ ]^_`%a[ [9Yc,
(78)
HIJK: 245
u(t)
: Y -."S[ D+
(76)
6[ R(77)
(51)
"R(78)
(56)
Y
(76)
(77)
(78)
89,& [-.de % d X Y)<)*::C &:L#"
(11)
: b,[ + Y5.2
6[%,&[ T+ ,-."S6*+/0
-."S6 *+ Y.,"#"/0-."S456
(70)
*+9[ de Y"#" [B '( )*Y[
(74)
'(" /0-.
t f
" 457%89,δx(t f )
Y)<(74)
(79)
)*RBC)RX9,A /0-.AYZ t f
t 0
∂J ∗
∂ dx(t) dt
δ dx(t) dt
dt = Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt
δx(t) dt − Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt
δx(t)dt
= ∂J ∗
∂ dx(t dt f ) δx(t f ) − ∂J ∗
∂ dx(t dt 0 ) δx(t 0 )
− Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt
δx(t)dt
= ∂J ∗
∂ dx(t dt f ) δx(t f ) − Z t f
t 0
d dt
∂J ∗
∂ dx(t) dt
δx(t)dt (79)
(79)
)R"#"S(75)
(80)
)*YδJ ∗ =
Z t f
t 0
∂J ∗
∂λ(t) δλ(t) + − d dt
∂J ∗
∂ dx(t) dt
+ ∂J ∗
∂x(t)
δx(t) + ∂J ∗
∂u(t) δu(t) dt
+ ∂J ∗
∂ dx(t dt f )
δx(t f ) = 0 (80)
?T[ ,
(80)
AV W Rδx(t f )
(81)
9Y∂J ∗
∂ dx(t dt f )
= λ(t f ) = 0 (81)
(81)
\]^_`%a:/0 69 ['( 6Y6
!"#$%&'()*+,-./01. 234&5[6678,&:4L9#:O ;<'(=>Y"#8<
-.de#:";<2123?:)<<>45-U@6A
(82)
:)R*+2123:"S>+ ,-.
t 0
:"S>456A(83)
:)R*+ >Ydx(t)
dt = − Ax(t) + Bu(t) (82)
x(t 0 ) = x 0 (83)
:2)B
A
");B
123 :"S>:68,)B >Y c ,2JK:4A(84)
)*:%/82)BQ
2R
123:"S>JKCW56=C4:)<<89,DE)B6=
10
YF\F]^_`%a4 %:0 1 67 [:D=>F\ F]^_
`%a4
p(t)
6G82JK:46HI=>A(85)
)*:>YJ (x(t), u(t)) = 1
2 Z t f
t 0
(x(t) T Qx(t) + u(t) T Ru(t))dt (84)
10
GHIJKLMNO
1.0
PKLQRSTUNVWXYZ [\]^_U`a0.0
b1.0
U cdefghYiJ ∗ (p(t), x(t), u(t)) = 1 2
Z t f
t 0
(x(t) T Qx(t) + u(t) T Ru(t))dt +
Z t f t 0
p(t) T ( dx(t)
dt + Ax(t) − B u(t))dt (85)
2HIJK:4L96O67 8<2HIJK:4?T@5A
(86)
)*:>Y
δJ ∗ = Z t f
t 0
(δp(t)) T ∂J ∗
∂p(t) + ∂J ∗
∂ d x dt (t) T
δ dx(t) dt
+ ∂J ∗
∂x(t) T
δx(t)+ ∂J ∗
∂u(t) T
δu(t) dt = 0 (86)
c,2+ ,-.:"S>45
x(t 0 )
7%89<;> @5δx(t 0 )
>YA
(86)
?@A:D8<B5C 5'(6=>2+,-.:"S >45@5>6DE8[6 => 2A
(87)
)*:>YZ t f
t 0
∂J ∗
∂ d x dt (t) T
δ dx(t) dt
dt = Z t f
t 0
d dt
∂J ∗
∂ d x dt (t) T
δx(t) dt − Z t f
t 0
d dt
∂J ∗
∂ d x dt (t) T
δx(t)dt
= ∂J ∗
∂ d x dt (t f ) T
δx(t f ) − ∂J ∗
∂ d x dt (t 0 ) T
δx(t 0 )
− Z t f
t 0
d dt
∂J ∗
∂ d x dt (t) T
δx(t)dt
= ∂J ∗
∂ d x dt (t f ) T
δx(t f ) − Z t f
t 0
d dt
∂J ∗
∂ d x dt (t) T
δx(t)dt (87)
)<<2 A
(87)
6A(86)
:FG =>A(88)
)*:>YδJ ∗ =
Z t f
t 0
∂J ∗
∂p(t) δp(t) + − d dt
∂J ∗
∂ d x dt (t)
+ ∂J ∗
∂x(t)
δx(t) + ∂J ∗
∂u(t) δu(t) dt
+ ∂J ∗
∂ d x dt (t f ) δx(t f ) = 0 (88)
A
(88)
)RHIJK:4?T@5>,:2HIJK:4:4:D=>&5>=>A
(89)
2(90)
2(91)
2(92)
6X> >Y∂J ∗
∂p(t) = dx(t)
dt + Ax(t) − Bu(t) = 0 (89)
− d dt
∂J ∗
∂ d x dt (t)
+ ∂J ∗
∂x(t) = − dp(t)
dt + A T p(t) + Q T x(t) = 0 (90)
∂J ∗
∂u(t) = R T u(t) − B T p(t) = 0 (91)
∂J ∗
∂ d x dt (t f ) = p(t f ) = 0 (92)
2 2:"S>45
u(t)
612345x(t)
:4 =6D+>Y 245A
(91)
:)<<*+ 9 2F\ F]^_`%a4p(t)
:)R89<;> 2` :4) BP(t)
6G => 2:4x(t)
F\ F]^ _`%a4p(t)
:6=A(93)
6%=>Y
p(t) = P(t)x(t) (93)
c,2A
(93)
6A(90)
:FG =>A(94)
)*:>Y− d dt
∂J ∗
∂ d x dt (t)
+ ∂J ∗
∂x(t) = − dp(t)
dt + A T p(t) + Q T x(t)
= − d(P (t)x(t))
dt + A T P (t)x(t) + Q T x(t)
= − dP (t)
dt x(t) − P(t) dx(t)
dt + A T P (t)x(t) + Q T x(t) = 0 (94)
c,2A
(91)
)R 245A(95)
)*:>Yu(t) = (R T ) −1 B T p(t) = (R T ) −1 B T P (t)x(t) (95)
2 A
(94)
667 >A(89)
:FG=>2A(96)
)*:>Ydx(t)
dt + Ax(t) − Bu(t) = dx(t)
dt + Ax(t) − B(R T ) −1 B T P (t)x(t) = 0 (96)
A
(96)
6=>212345-U@A(97)
)*:*+ 9>Ydx(t)
dt = − (A − B(R T ) − 1 B T P (t))x(t) (97)
c,2A
(97)
6A(94)
:FG => ` :4 )BP (t)
:)>&5[2A(98)
6X>Y− dP (t)
dt x(t) − P (t) dx(t)
dt + A T P (t)x(t) + Q T x(t)
= − dP (t)
dt x(t) + P(t)(A − B(R T ) −1 B T P (t))x(t) + A T P(t)x(t) + Q T x(t)
= − dP (t)
dt + P (t)(A − B(R T ) −1 B T P(t)) + A T P (t) + Q T x(t) = 0 (98)
2
x(t) 6 = 0
)R2A(98)
3B5)B>Y)<<2JK:496=>)B
P (t)
::=>[A(99)
)*:>Y− dP (t)
dt + P (t)(A − B(R T ) −1 B T P (t)) + A T P(t) + Q T = 0 (99)
A
(99)
[9>Y 2F\F]^_` %a4:D=>/067p(t f ) = 0
*+ 9<"R2
x(t f )
#$9 > 2 A(93)
)RA(99)
:D =>678<2P(t f )
)B>67)R'(6)*YA(93)
:";<2`:4)BP(t)
:)>>Ac9<"R2T:[:D=>T6> 89<;>YA
(99)
6,:%Z6b82X 9,2 9,
P (t)
:)R 2:"S >45u(t)
A(100)
)*:'(89>Yu(t) = (R T ) −1 B T P (t)x(t) (100)
7
!"#$%&'()*+,-./01*+-."#)*:23 4&5[6678,&:4L9# 2[:" ;
<`:4) B`:4 )BC 89> A6&5[X 92T6
:> 89<;>Y 2&:4L9>67[6,=)*:2 '
(:)R2&:46=>:46> %T: 9>Y
=>
[3]
Y -.c [6R>2JK:4L9#:D8<2 67[F\F]^ _`a4C6JK:4:+,HIJK:4L 96>#dZ:R+,
267[>[,2J K:4:+<;>AJK:49:6*+
; JK:4L9#6HIJK:4L 9# 8<#;, >YD
+6 :2-."#:D8<2 :S >:4
u(t)
2l
l + 1
U >c )+,2
u(t) l+1 − u(t) l
::=>A6HIJK:4:+<L9#6; , 8<2JK:4L9#dZ:R >Y:2
l
@492 :;<TO!O"R8#$89>9 >Y 2HIJK:4:
u(t) l+1 − u(t) l
:%R89,A6&8,JK:4K
6A(101)
% :/ =>YK (p(t) l , x(t) l , u(t) l , u(t) l+1 ) = J ∗ (p(t) l , x(t) l , u(t) l )
+ 1
2 Z t f
t 0
(u(t) l+1 − u(t) l ) T W l (u(t) l+1 − u(t) l )dt (101)
:
W l
C 'F()*:%<<89,DE)B R2)B :; <6 +,8>,:-.8<*+ 9>'F()* >Yc,2/8,&:4K
2J ∗
6H IJK:4 :D8<./JK041>Y 2./2304
K
:D =>49#6,:2./23045 6@56 =>A
(102)
% :>7δK l =
Z t f t 0
(δp(t) l ) T ∂J ∗ l
∂p(t) l
+ ∂J ∗ l
∂ d x dt (t) l T
δ dx(t) dt
l
+ ∂J ∗l
∂x(t) l T
δx(t) l + ∂J ∗l
∂u(t) l T
δu(t) l dt = 0 (102)
-.de2A
(102)
:D8<289:;:<>=>x(t 0 )
?@1<;> @5δx(t 0 )
>7A A(102)
A5 BC:D8<D5E 5A6=>289:;:<>=>A@56FG8<H6IJ =>KL2 A
(103)
A% :>7Z t f
t 0
∂J ∗ l
∂ d x dt (t) l T
δ dx(t) l dt
dt = Z t f
t 0
d dt
∂J ∗ l
∂ d x dt (t) l T
δx(t) l dt − Z t f
t 0
d dt
∂J ∗ l
∂ d x dt (t) l T
δx(t) l dt
= ∂J ∗ l
∂ d x (t dt f ) l T
δx(t f ) l − ∂J ∗ l
∂ d x (t dt 0 ) l T
δx(t 0 ) l
− Z t f
t 0
d dt
∂J ∗ l
∂ d x dt (t) l T
δx(t) l dt
= ∂J ∗ l
∂ d x (t dt f ) l T
δx(t f ) l − Z t f
t 0
d dt
∂J ∗ l
∂ d x dt (t) l T
δx(t) l dt (103)
%M<2 A
(103)
6A(102)
:NO=>A(104)
A% :>7δJ ∗ =
Z t f t 0
∂J ∗ l
∂p(t) l δp(t) l + − d dt
∂J ∗ l
∂ d x dt (t) l
+ ∂J ∗ l
∂x(t) l
δx(t) l + ∂J ∗ l
∂u(t) l δu(t) l dt + ∂J ∗ l
∂ d x (t dt f ) l
δx(t f ) l = 0 (104)
A