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

最適化数学第 復習:制約なし最適化問題 4 回

N/A
N/A
Protected

Academic year: 2021

シェア "最適化数学第 復習:制約なし最適化問題 4 回"

Copied!
13
0
0

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

全文

(1)

最適化数学 第 4 回

[今回の項目]

1 復習:局所最適解の求め方

2 問題練習

(2)

復習:制約なし最適化問題

最小化 f(x, y) =x3−3xy+y3 制約 なし

x3−3xy+y3 のグラフ

(3)

復習: 1 次の最適性条件

[定理]

(

一次の最適性条件

)

f(x):多変数関数

が局所最適解ならば,

∇f(¯x) =0 (零ベクトル)

が成り立つ.

[定義]

pが,∇f(p) =0 を満たすとき,p f の停留点と呼ぶ.

(4)

復習:停留点

最小化 f(x, y) =x3−3xy+y3

x3−3xy+y3 のグラフ x3−3xy+y3 の勾配ベクトル

(5)

復習:ヘッセ行列を用いた判定法

停留点を局所最小解とすると,

¯

xで『グラフは局所的に下に窪んでいる』

←→ の近くで関数f は『局所的に凸関数である』

[定理]

( 2

次の最適性条件

)

1 (必要性)¯x が局所最小解

=⇒ ∇f(¯x) =0 かつ2f(¯x)が半正定値

2 (十分性)∇f(¯x) =0 かつ2f(¯x)が正定値

=⇒x¯ は局所最小解.

3 (否定)2f(¯x) が不定値のとき,¯x は局所最適解ではない.

局所最大解についても,それぞれ対応する箇所を半負定値,負定 値,極大値に置き換えたものが成り立つ.

(6)

復習: 2 次の最適性条件の幾何的イメージ

(十分性)∇f(¯x) =0 かつ2f(¯x)が正定値

=⇒x¯ は局所最小解.

(否定)2f(¯x) が不定値のとき,¯x は局所最適解ではない.

2f(x, y)

2f(x, y)

(7)

例題

f(x, y) =x3−3xy+y3 の局所最適値を求めよ.

[補足]

∇f(¯x) が正定値

⇐⇒ ∇f(¯x)の固有値がすべて正 f 2 変数の場合:

det

2f(¯x)

>0(行列式),かつ fxx(¯x)>0 =⇒ f は正定値

(8)

例題

f(x, y) =x3+y3−6x−6y の局所最適値を求めよ.

(解答例) ∇f(x, y) =

3x2−6 3y2−6

, ∇2f(x, y) =

6x 0 0 6y

より,

連立方程式

(3x2−6 = 0

3y2−6 = 0 を解くと,停留点 (x, y) = (±√ 2,±√

2)

(順不同)を得る.

(x, y) (√ 2,√

2) (−√ 2,−√

2) (±√ 2,∓√

2)

|∇2f(x, y)| + + −

fxx(x, y) + −

不定

f(x, y) −8√

2 8√

2 である.よって,(x, y) = (√

2,√

2)のとき,局所最小値 −8√ 2 (x, y) = (−√

2,−√

2) のとき,局所最大値 8√

2をとる.

(9)

練習問題

(1) f(x, y) =x3+ 8y3−12xy + 1 (2) f(x, y) =x3+ 3xy2−3x

(10)

練習問題

(3) f(x, y) =x4−4x2+ 8y2+ 4x2y

(4) f(x, y, z) =x3−y3+z2+ 3y2+ 3xy−3x−3y−3z+ 2

(11)

(1) f(x, y) =x3 + 8y3−12xy+ 1

-2 -1 0 1 2 3 4 5

-2 -1 0 1 2 3 4 -15 5 -10 -5 0 5 10 15 20 25 30

x**3 + 8*y**3 - 12*x*y + 1

-15 -10 -5 0 5 10 15 20 25 30

(12)

(2) f(x, y) =x3 + 3xy2−3x

-2 -1.5

-1 -0.5

0 0.5

1 1.5

2

-1.5 -2 -0.5 -1

0.5 0 1.5 1

-3 2 -2 -1 0 1 2 3

x**3 + 3*x*y**2 - 3*x

-3 -2 -1 0 1 2 3

(13)

(3) f(x, y) =x4 −4x2+ 8y2+ 4x2y

-3 -2

-1 0

1 2 3 -4

-3 -2

-1 0

1 2

-10 -5 0 5 10 15 20

x**4 - 4*x**2 + 8*y**2 + 4*x**2*y

x

y -10 -5 0 5 10 15 20

参照

関連したドキュメント

Optimal stochastic approximation algorithms for strongly convex stochastic composite optimization I: A generic algorithmic framework.. SIAM Journal on Optimization,

Dual averaging and proximal gradient descent for online alternating direction multiplier method. Stochastic dual coordinate ascent with alternating direction method

Hungarian Method Kuhn (1955) based on works of K ő nig and

of IEEE 51st Annual Symposium on Foundations of Computer Science (FOCS 2010), pp..

情報理工学研究科 情報・通信工学専攻. 2012/7/12

最大消滅部分空間問題 MVSP Maximum Vanishing Subspace Problem.. MVSP:

参考文献 Niv Buchbinder and Joseph (Seffi) Naor: The Design of Com- petitive Online Algorithms via a Primal-Dual Approach. Foundations and Trends® in Theoretical Computer

"A matroid generalization of the stable matching polytope." International Conference on Integer Programming and Combinatorial Optimization (IPCO 2001). "An extension of