21.02.2014 Aufrufe

Skript

Skript

Skript

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.

mit c < 1 2<br />

ausreichend ist, zumindest lokal um eine Lösung. Die Bedingung<br />

5.4 beschränkt die Art der nichtlinearität des Operators F und wird tangential<br />

cone condition genannt. Diese Bedingung ersetzt in gewisser Weise die stetige<br />

Invertierbarkeit des Operators F ′ (x † ), die bei der Konvergenzanalyse wohlgestellter<br />

Probleme eine zentraler Vorraussetzung ist. Conditions on the nonlinearity<br />

are not just artificial inventions, it can be shown that a slightly weaker<br />

condition than (5.4) is actually necessary for the convergence of the Landweber<br />

iteration (cf. [?]).<br />

5.2 Newton-artige Methoden<br />

Die grundlegende Ide von Newton-Methoden für nichtlineare Probleme wie (??)<br />

ist eine lokale Linearisierung. Daher beinhaltet jeder Schritt der Methode die<br />

Lösung des linearen Problems<br />

F ′ (x k )(x k+1 − x k ) = −(F (x k ) − y). (5.5)<br />

Da F ′ (x k ) jedoch im Fall schlecht-gestellter Probleme keine regulärer Operator<br />

ist, ist die Gleichung zur Berechnung von selbst ein schlecht-gestelltes, lineares<br />

Problem und x k+1 muss damit nicht wohl-definiert sein. Eine gebräuchliche<br />

Stategy to Konstruktion von Newton-Methoden für nichtlineare, schlechtgestellte<br />

Probleme ist daher Gleichung (5.5) mit einer Methode für lineare Problem<br />

zu regularisieren. Zum Beispiel ergibt Tikhonov-Regularisierung (wobei<br />

x k+1 − x k als die Unbekannte interpretiert wird)<br />

(F ′ (x k ) ∗ F ′ (x k ) + α k I)(x k+1 − x k ) = −F ′ (x k ) ∗ (F (x k ) − y), (5.6)<br />

die sogenannte Levenberg-Marquardt-Methode. Alternativ kann man die Tikhonov-<br />

Regularisierung auch mit anderen a-prioi Informationen versehen, die gebräuchlichste<br />

Form ist der Startwert x 0 in jedem Iterationsschritt:<br />

(F ′ (x k ) ∗ F ′ (x k )+α k I)(x k+1 −x k ) = −F ′ (x k ) ∗ (F (x k )−y)+α k (x 0 −x k ). (5.7)<br />

Diese Verfahren wird iterativ regularisierte Gauss-Newton-Methode genannt.<br />

Der zusätzliche Term erhöht in einigen Fällen die Stabilität der Methode. Wir<br />

bemerken, dass in beiden Fällen immernoch der Stop-Index der Regularierungsparameter<br />

ist, nicht α k ! Der Parameter α k kann so gewählt werden, dass er gegen<br />

Null konvergiert für k → ∞, um eine Überdämpfung der Newton-Methode<br />

zu vermeiden. Im fall der Gauss-Newton-Methode ist α k → 0 sogar eine notwendige<br />

Vorraussetzung für die Konvergenz der Methode, denn für α k → α ∞ ><br />

0 würde man die Lösung des Tikhonov-Regularisierten Problems approximieren<br />

und nicht die Lösung von (??).<br />

The convergence analysis of Newton-type method is rather involved and shall<br />

therefore be omitted here, we only mention that similar (sometimes weaker)<br />

nonlinearity conditions as for the nonlinear Landweber iteration are used.<br />

81

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

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!