27.01.2014 Aufrufe

Hesse-Matrix - imng

Hesse-Matrix - imng

Hesse-Matrix - imng

MEHR ANZEIGEN
WENIGER ANZEIGEN

Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.

YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.

<strong>Hesse</strong>-<strong>Matrix</strong><br />

Die quadratische Taylor-Approximation einer skalaren Funktion f in einem<br />

Punkt mit Koordinaten a = (a 1 , . . . , a n ) lässt sich in der Form<br />

f (x 1 , . . . , x n ) = f (a) + (grad f (a)) t (x − a)<br />

+ 1 2 (x − a)t H f (a)(x − a) + O(|(x − a)| 3 )<br />

schreiben, wobei die symmetrische <strong>Hesse</strong>-<strong>Matrix</strong><br />

⎛<br />

⎞<br />

∂ 1 ∂ 1 f (a) · · · ∂ 1 ∂ n f (a)<br />

⎜<br />

⎟<br />

H f (a) = ⎝ .<br />

. ⎠<br />

∂ n ∂ 1 f (a) · · · ∂ n ∂ n f (a)<br />

die zweiten Ableitungen enthält.<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 1-1


Beweis:<br />

Umschreiben der quadratischen Terme der Taylor-Entwicklung (n = 2) <br />

<strong>Hesse</strong>-<strong>Matrix</strong> 2-1


Beweis:<br />

Umschreiben der quadratischen Terme der Taylor-Entwicklung (n = 2) <br />

f (x, y) = f + f x (x − x 0 ) + f y (y − y 0 )<br />

+ 1 2! (f xx (x − x 0 ) 2 + 2f xy (x − x 0 ) (y − y 0 ) + f yy (y − y 0 ) 2 ) + R<br />

= f + f x (x − x 0 ) + f y (y − y 0 )<br />

+ 1 ( ) (<br />

2! ((x − x fxx f<br />

0) , (y − y 0 ))<br />

xy (x − x0 )<br />

f xy f yy (y − y 0 )<br />

wobei f und sämtliche partiellen Ableitungen im Punkt (x 0 , y 0 )<br />

ausgewertet werden<br />

)<br />

+ R ,<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 2-2


Restglied<br />

mit<br />

für ein θ ∈ [0, 1]<br />

R = ∑<br />

|α|=3<br />

1<br />

α! ∂α f (u, v) (x − x 0 ) α 1<br />

(y − y 0 ) α 2<br />

(u, v) = (1 − θ)(x 0 , y 0 ) + θ(x, y)<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 2-3


Restglied<br />

mit<br />

für ein θ ∈ [0, 1]<br />

Abschätzung<br />

R = ∑<br />

|α|=3<br />

1<br />

α! ∂α f (u, v) (x − x 0 ) α 1<br />

(y − y 0 ) α 2<br />

(u, v) = (1 − θ)(x 0 , y 0 ) + θ(x, y)<br />

R = O(|(x − x 0 , y − y 0 )| 3 )<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 2-4


Restglied<br />

mit<br />

für ein θ ∈ [0, 1]<br />

Abschätzung<br />

R = ∑<br />

|α|=3<br />

1<br />

α! ∂α f (u, v) (x − x 0 ) α 1<br />

(y − y 0 ) α 2<br />

(u, v) = (1 − θ)(x 0 , y 0 ) + θ(x, y)<br />

R = O(|(x − x 0 , y − y 0 )| 3 )<br />

Beweis im allgemeinen Fall (n ≥ 2) analog<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 2-5


Beispiel:<br />

f (x, y) = ln(x + 1/y)<br />

partielle Ableitungen<br />

1<br />

f x =<br />

x + 1/y , f y = − 1<br />

xy 2 + y ,<br />

1<br />

f xx = −<br />

(x + 1/y) 2 , f 1<br />

xy =<br />

(xy + 1) 2 , f yy = 2xy + 1<br />

(xy 2 + y) 2<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 3-1


Beispiel:<br />

f (x, y) = ln(x + 1/y)<br />

partielle Ableitungen<br />

1<br />

f x =<br />

x + 1/y , f y = − 1<br />

xy 2 + y ,<br />

1<br />

f xx = −<br />

(x + 1/y) 2 , f 1<br />

xy =<br />

(xy + 1) 2 , f yy = 2xy + 1<br />

(xy 2 + y) 2<br />

Gradient und <strong>Hesse</strong>-<strong>Matrix</strong> im Punkt (0, 1):<br />

( )<br />

( )<br />

1 −1 1<br />

grad f (0, 1) = , H f (0, 1) =<br />

−1<br />

1 1<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 3-2


die quadratische Taylor-Entwicklung<br />

( ) x<br />

f (x, y) = (1, −1)<br />

y − 1<br />

+ 1 ( ) ( −1 1 x<br />

2 (x, y − 1) 1 1 y − 1<br />

)<br />

+ O(|(x, y)| 3 )<br />

= x − y + 1 + 1 2<br />

(<br />

−x 2 + 2x(y − 1) + (y − 1) 2) + O(|(x, y − 1)| 3 )<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 3-3


Beispiel:<br />

quadratische Taylor-Entwicklung von<br />

f (x, y, z) = (xy) z<br />

im Punkt (1, 1, 1)<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 4-1


Beispiel:<br />

quadratische Taylor-Entwicklung von<br />

im Punkt (1, 1, 1)<br />

partielle Ableitungen<br />

f (x, y, z) = (xy) z<br />

f x = yz(xy) z−1 , f z = ln(xy)(xy) z ,<br />

f xx = y 2 z(z − 1)(xy) z−2 , f zz = (ln(xy)) 2 (xy) z ,<br />

f xy = z(xy) z−1 + xyz(z − 1)(xy) z−2 , f xz = y(xy) z−1 + yz ln(xy)(xy) z−1<br />

(f y , f yy und f yz durch Vertauschen der Variablen)<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 4-2


Gradient und <strong>Hesse</strong>-<strong>Matrix</strong>:<br />

⎛<br />

grad f (1, 1, 1) = ⎝<br />

1<br />

1<br />

0<br />

⎞<br />

quadratische Taylor-Entwicklung<br />

⎛<br />

⎠ , H f (1, 1, 1) = ⎝<br />

0 1 1<br />

1 0 1<br />

1 1 0<br />

⎞<br />

⎠<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 4-3


Gradient und <strong>Hesse</strong>-<strong>Matrix</strong>:<br />

⎛<br />

grad f (1, 1, 1) = ⎝<br />

1<br />

1<br />

0<br />

⎞<br />

⎛<br />

⎠ , H f (1, 1, 1) = ⎝<br />

quadratische Taylor-Entwicklung<br />

⎛ ⎞<br />

x − 1<br />

f (x, y, z) = 1 + (1, 1, 0) ⎝ y − 1 ⎠<br />

z − 1<br />

⎛<br />

+ 1 2 (x − 1, y − 1, z − 1) ⎝<br />

+ O(|(x, y, z)| 3 )<br />

0 1 1<br />

1 0 1<br />

1 1 0<br />

0 1 1<br />

1 0 1<br />

1 1 0<br />

⎞ ⎛<br />

⎠ ⎝<br />

= 1 + (x − 1) + (y − 1) + (x − 1)(y − 1)<br />

+(x − 1)(z − 1) + (y − 1)(z − 1)<br />

+ O(|(x − 1, y − 1, z − 1)| 3 )<br />

⎞<br />

⎠<br />

x − 1<br />

y − 1<br />

z − 1<br />

⎞<br />

⎠<br />

<strong>Hesse</strong>-<strong>Matrix</strong> 4-4

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!