• 検索結果がありません。

最適化問題と汎関数の極値問題

N/A
N/A
Protected

Academic year: 2022

シェア "最適化問題と汎関数の極値問題"

Copied!
23
0
0

読み込み中.... (全文を見る)

全文

(1)

九州大学学術情報リポジトリ

Kyushu University Institutional Repository

最適化問題と汎関数の極値問題

倉橋, 貴彦

長岡技術科学大学機械系機械情報・制御工学大講座

https://doi.org/10.15017/1470160

出版情報:九州大学情報基盤研究開発センター全国共同利用システム広報. 2 (1), pp.17-28, 2008-08.

九州大学情報統括本部広報委員会 バージョン:

権利関係:

(2)

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|k

x

!es

f (x) = x 2 − 2x + 5

2}[@=D~%ok

f

A.d234

x

,%

ok

f

x

!o34€A)34?2N

df dx (x) = 2x − 2

2#K‚Xƒ@=%?‚HX„2G

4 ?2Nok

f

X.d2G4

x

X…†4>?ZS~ %ok

f

A.d234|k,

x = 1.0

2

Gs%ok,

f (1.0) = 4.0

2G4>

O‡% RS!"W4ok,%okAˆ3|k‰/XŠ‹

t

64#,Œ‹Aˆ3 Ž

x

%

y

%

z

8|kok2G$#4?2X‘ #>okAˆ3|k‰/X’G4|k ok2 G$#4/,“ok2”\=%“oklHA…E4RS!fp$,%okl HA…E4R

S2,’G4tiX•#@=4> oklHRSN,%okAˆ 3|k|QaA. d234eK

G|kX…E@=4

1

!fp%“oklHRSN,%“okAˆ3ok|QaX.d2G4 e

KG okX…E@=4>“okAˆ 3ok|Qa,|2”\=%“okl HRS,“ok

|X„2G4?2N%“ok!lHA}[4okX…E@=4>?eKG“okl HRSA

tu‡i,|i2”\=4>vwN,%?“oklHRS!–—A6$%|i!e 4“o

klHRS!f 34ti!]#$z{A;K>

2

˜™š

v›N,“oklHRSAtu|i!]# $z{34

[1]

>|i!e4‚œN ,ž%“

ok!lHA}[ 4okAŸ]W4DE!% ¡¢JpDok

x(t)

2‰|Qa

δx(t)

%£ L¤ J

k

ε

A•#$“okA

J (x(t) + εδx(t))

2ˆ3 >¥¦

1

§>?“ok

J (x(t) + εδx(t))

,¤|k

ε

¤k Hok2Gs

G(ε)

2¨W4 ¥©

(1)

§>

G(ε) = J (x(t) + εδx(t)) (1)

??N©

(1)

!"# $%£LJk

ε

ªA

∆ε

2342%ª

∆ε

!e4“okª ,« ¬e

K!¨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

àáâãäÕå

(3)

1:

ok|Qa

= dJ(x(t) + εδx(t))

dx(t) ∆εδx(t) + o( k ∆εδx(t) k ) (3)

(3)

 ®Aª

∆ε

N42©

(4)

eK!Gs

2

∆ε

AsGu„! W4?2

Nl,©

(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)

Aƒ4?2XNY4>

dG(ε) dε

ε=0 = dJ(x(t))

dx(t) δx(t) (7)

†D%©

(1)

£LJk

ε

!f34€Ap%

ε = 0

2342©

(8)

Aƒ4>

dG(ε) dε

ε=0 = dJ(x(t) + εδx(t)) dε

ε=0 (8)

©

(7)

%

(8)

es©

(9)

oAƒ4?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)

ÊÜÌÞå

(4)

©

(9)

,%64“okA£LJk

ε

N€p

ε = 0

2pD/,%

dJ dx(t) (x(t)) δx(t)

2ˆj=4?2AL

Mp$#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)

N€N6=\%©

(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)

N€N6=\%“okXl

HA24~,sc[XNY4>O ],“ok

J

A£LJk

ε

!es€p

ε = 0

2p

!l HA]2#K?2N64>/KO]c[,%“ok

J (x a (t))

€kX„2G4?2Xl

HA/] !G$#4?2N64>?=@?2!] #$%  @tNY4 eK«¬z{

A;K>†D%©

(11)

,“ok

J (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…E

4RS!]#$z{34 ¥¦

2

§

[2]

P 0 (0, 0)

@—

P 1 (a, b)

†N.A…E4RS,%

.RS2/”\=$#4>64—‹A!"#$%£L€d‹

ds

,%

x

%

2:

.RS

(5)

y

‡ €d‹

dx

%

dy

A•#4 ?2N‡Jes©

(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%64—‹AX.2G4RSAtuDE%©

(14)

!3€d‹A64—

‹N 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

,|k

x

okN6s%|k

x

,|k

t

ok2G$#4DE%“oklHRS2G4> ??N%?“oklHRSAtY%.

2G4ok

x(t)

A u>†ž“oklHRSAtuDE!%ok

x(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

X

x(t)

NlHA/ %

x(t)

N€N6=\%©

(16)

,¤|k

ε

ok2Gs%

ε = 0

NlHA24>e$%ƒ@=D©

(16)

!fp$%

x(t)

!"W4O|

δJ

)A34

2% ©

(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$ Ap%|

δx(t)

X ¡—!"#$„2G4?2Ac34 2% ©

(18)

Aƒ4>

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)

(6)

†D%£L|

δx

!fp

(18)

2G4DE!,%ok

f (t)

XJk

A

2G4xyX64

¥©

(19)

§>

f (t) =

dx(t)

dt

r

1 + dx(t) dt 2

= A (19)

??N%©

(19)

!3€‡ ‚Atu?2Nt2G4ok

x(t)

Aƒ4?2XNY4>e$%

©

(19)

es%. :!©

(20)

Aƒ4>

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

,Jk

B

2G4xyX6s% ®A 34?2

(21)

Aƒ4>

dx(t) dt = B Z dx(t)

dt dt = Z

Bdt

x(t) = Bt + C (21)

??N%  ¡—

x

HesJk

B

%

C

AJ342©

(22)

(23)

eK!G4>

x(0) = 0 ...C = 0 (22)

x(a) = b ...B = b

a (23)

e$%©

(21)

!Jk

B

%

C

A 34?2N.2G4ok

x(t)

(24)

eK@

=4>

x(t) = b

a t (24)

©

(24)

es—

P 0 (0, 0)

@—

P 1 (a, b)

†NA.2G4ok

x(t)

,

2

—‹A N6

4 ?2X <4>

(7)

4

›N,%“oklHRS!fp$‚œA;DX%v›N,%“okAˆ 3okHX6

4 ‚!es}[@=4RS!fp$%“okA. 234okAIJ34RS!] #$z{3

4

[2]

>“okAˆ3ok!f34 ‚, ‚

3

2”\=%?RS, VY“ok

lHRS2G4> VY“oklHRSZ2 p$%64

2

—

P 0

%

P 1

N ¢JpD

j

L

2

P 0

%

P 1

AN†=D!esˆj=4A.234A…E 4RS

Ac[ 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

"e—

P 0

%

P 1

A!es †=4

S

A.!34ok

y(x)

A…E4>

A.234okA…E4RS,%©

(25)

%

(26)

‚/21ok ¥©

(27)

§

HA.234ok

y(x)

A…E4RS2G4> ‚AD3eK!okA.64#,. d

!34RS!fp$,%O!±± Jki2#K‡iA•#4>±±

JkiN,% ‚!f34± ± Jk

5

A 34>1ok!fp

$±± Jk2 ‚!esˆj=4AVpD34?2N% V

3

!"#Ü$%&äÕå

4

'Ç(ÉÎÊ)*(ÉÜ%&äÕå

5

+,-ÑÎÏ./ÇÀ0 !"#ÊÓÔÕ-Ñ

(

1

Ñ

)

Ü 2Í3ÌãäÏ+4 ×56Õ 'Ü7ã+,-Ñ 8ÌÎ+,

1

Ñ9Ü%&äÍÌÕå

(8)

X„2G48 ‚Nˆj=$#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 p

1okAC342©

(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%ok

y(x)

A

εδy(x)

!es342©

(29)

Aƒ4>

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)

es1okok

y(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)

2G4 Ac

342%©

(31)

Aƒ4?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

ÐÑÜ$%&äÕå

(9)

©

(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)

X„2G4DE! ,% ©

(32)

¬ X„2G4x

yX64 ¥©

(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)

Aƒ4>

Z d dx

sin θ dx = Z 1

λ dx sin θ = x

λ + A 1

λ sin θ = x + A 2 (35)

©

(35)

!"#$

θ

!f34€Ap342. :!©

(36)

Aƒ4>

d

dθ (λ sin θ) = d

dθ (x + A 2 )

λ cos θ = dx dθ

dx = λ cos θdθ (36)

??N%

dy(x)

dx = tan θ

2JpD‚,©

(37)

eK!¨Yˆj=4>

dy = tan θdx (37)

©

(37)

!fp

(36)

Ap%342©

(38)

eK!G4>

dy = tan θ(λ cos θdθ)

(10)

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

Aˆ3ok,©

(39)

eK!G4>

(x + A 2 ) 2 + (y + B 2 ) 2 = λ 2 × (sin 2 θ + cos 2 θ)

= λ 2 (39)

©

(39)

,‡ ‚Aˆp$"s%

S

A.!34

L

,A

λ

% ŽHA

(A 2 ,B 2 )

234‡ ‚2G4>

5

™

v ›N,%“oklHRS!"#$ ‚X€‡ ‚!es}[@=4~ !]#$z

{34

[2]

>ZS2 p$%Š‹

t

!"W464X

x(t)

!es}[@=42 p%

A

u(t)

!es34RSAc[ 4

7

¥¦

4

§>Š‹|Q

dx(t) dt

,%

4:

)y

7

ÇÎÏÊãÇ Ç!ÊË Ì Í$"ÇÜ 2Ï#$ÊÇ% & ÔÕ$ÇÜÔÕå

(11)

2!e$}[@=©

(40)

!esˆj=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>

AA•p$¬4DE%!o34 •Ac[42% A4

Wb[$A¬4?2X†p#>e$%2!es1

okAJq 34>2HYjA134DE!%HA

pDHA)p7Š‹Np<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.d234

u(t)

A…E4?2!G4% VY“oklHRS!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)

2

u(t)

‰= =!fp

LJk

ε 0

%

ε 1

%

ε 2

A p% ok!fp$

ε 0 δλ(t)

%

ε 1 δx(t)

%

ε 2 δu(t)

A[42©

(45)

A

ƒ4

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)

2342

(45)

es

1okok

λ(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 Ø"ÍÌÕå

(12)

©

(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 O

x(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 W

R

(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)

(13)

H IJK:

J

:

u(t)

L MN O

u(t)

&PQ

(53)

)RH

IJK: :

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)

XY

u(t) − λ(t) = 0 (55)

(55)

)R 245:Z[\ ]^_`%a b

(56)

)*

Y

u(t) = λ(t) (56)

(52)

(56)

"\ ]^_a)<89123 45 2

45JK:L X )*'(8912345-U

(40)

FG\]^_`%a[ X

(57)

Y

d

dt dλ(t)

dt − aλ(t) +a dλ(t)

dt − aλ(t) − λ(t) = 0 (57)

c

(57)

(58)

)*Y

d 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)

XY

x(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)

XY

x(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)

(14)

+ ,-.");/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)

)*Y

x 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)

)*Y

x(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<1245

x(t)

6

(41)

FG JK

: 4512345

x(t)

:[

(67)

)*cY

u(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*+Y

dx(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)

)*Y

J (x(t), u(t)) = 1 2

Z t f t 0

(x(t) 2 + u(t) 2 )dt (71)

(15)

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)

)*Y

Z 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)

)RY

c,

(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,[ + Y

(16)

5.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-.AY

Z 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 ['( 6Y

6

!"#$%&'()*+,-./01

. 234&5[6678,&:4L9#:O ;<'(=>Y"#8<

-.de#:";<2123?:)<<>45-U@6A

(82)

:)R*+2

123:"S>+ ,-.

t 0

:"S>456A

(83)

:)R*+ >Y

dx(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)B

Q

2

R

123:"S>JKCW56=C4:)<<

89,DE)B6=

10

YF\F]^_`%a4 %:0 1 67 [:D=>F\ F]^_

`%a4

p(t)

6G82JK:46HI=>A

(85)

)*:>Y

J (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`a

0.0

b

1.0

U cdefghYi

(17)

J (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 )

>Y

A

(86)

?@A:D8<B5C 5'(6=>2+,-.:"S >45@5>

6DE8[6 => 2A

(87)

)*:>Y

Z 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)

612345

x(t)

:4 =6D+>Y 2

45A

(91)

:)<<*+ 9 2F\ F]^_`%a4

p(t)

:)R89<;> 2` :4) B

P(t)

6G => 2:4

x(t)

F\ F]^ _`%a4

p(t)

:6=A

(93)

6

%=>Y

p(t) = P(t)x(t) (93)

(18)

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)

)*:>Y

u(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)

)*:>Y

dx(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>Y

dx(t)

dt = − (A − B(R T ) 1 B T P (t))x(t) (97)

c,2A

(97)

6A

(94)

:FG => ` :4 )B

P (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=>/067

p(t f ) = 0

*+ 9<"R2

x(t f )

#$9 > 2 A

(93)

)RA

(99)

:D =>678<2

P(t f )

)B>67)R'(6)*YA

(93)

:";<2`:4)B

P(t)

:)>

>Ac9<"R2T:[:D=>T6> 89<;>YA

(99)

6

,:%Z6b82X 9,2 9,

P (t)

:)R 2:"S >45

u(t)

A

(100)

)*:'(89>Y

u(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

(19)

=>

[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)

2

l

l + 1

U >c )+,2

u(t) l+1 − u(t) l

::=>A6HIJK:4:+<L9

#6; , 8<2JK:4L9#dZ:R >Y:2

l

@4

92 :;<TO!O"R8#$89>9 >Y 2HIJK:4:

u(t) l+1 − u(t) l

:%R89,A6&8,JK:4

K

6A

(101)

% :/ =>Y

K (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,&:4

K

2

J

6H IJ

K: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% :>7

Z 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

(104)

%PQR2304A56 @5>S:2QR2304AT04:D=>U5 >=>A

(105)

2

(106)

2

(107)

2

(108)

6V>KL >7

∂J l

∂p(t) l = dx(t) l

dt + Ax(t) l − Bu(t) l = 0 (105)

参照

関連したドキュメント

In this paper we study optimal control problems governed by fractional stochastic partial neutral functional integro-differential equations with infinite delay in Hilbert spaces..

initial functions are proved in the form of an integral maximum principle and conditions of transversality for nonlinear systems with a variable structure, delays and a

Optimal control problems for PDEs are most completely studied for the case in which the control functions occur either on the right-hand sides of the state equations, or the boundary

In the second part, using the Malliavin calculus approach, we deduce a general maximum principle for optimal control of general stochastic Volterra equations.. The result is applied

Their basic components are the representation of candidate solutions to the problem in a “genetic” form, the creation of an initial, usually random population of solutions,

Assume that Γ &gt; 3γ/2 and the control bound m is large enough such that the bang arc u m starting from the north pole intersects the singular arc z 0 γ/2δ, Then for the problem

ポートフォリオ最適化問題の改良代理制約法による対話型解法 仲川 勇二 関西大学 * 伊佐田 百合子 関西学院大学 井垣 伸子

This paper deals with the a design of an LPV controller with one scheduling parameter based on a simple nonlinear MR damper model, b design of a free-model controller based on