Fachbereich Mathematik - GSI
Fachbereich Mathematik - GSI
Fachbereich Mathematik - GSI
Sie wollen auch ein ePaper? Erhöhen Sie die Reichweite Ihrer Titel.
YUMPU macht aus Druck-PDFs automatisch weboptimierte ePaper, die Google liebt.
4.1 Numerische Minimierung mit Linesearch-Verfahren<br />
Üblicherweise wird im zweiten Punkt des allgemeinen Linesearch-Verfahrens als<br />
Abbruchkriterium die Stationarität an der Stelle Nk getestet, also:<br />
∇χ 2 ( Nk) = 0 . (4.5)<br />
Dieses Abbruchkriterium ist zum einem numerisch nicht realisierbar und zum anderen<br />
kann es bei restringierten Optimierungsproblemen nicht verwendet werden, da in<br />
diesem Fall Minima nicht notwendigerweise (4.5) erfüllen müssen. Dieser Fall kann<br />
zum Beispiel eintreten, wenn das Minimum auf dem Rand der zulässigen Menge liegt.<br />
Dies wurde bereits in Abschnitt 3.4 angesprochen. Um unnötig lange Rechenzeiten<br />
zu vermeiden, können folgende Abbruchkriterien verwendet werden:<br />
• Die relative Änderung der Zielfunktionswerte unterscheidet in aufeinanderfolgenden<br />
Iterationsschritten einen bestimmten Wert:<br />
|χ 2 ( Nk−1) − χ 2 ( Nk)|<br />
χ 2 ( Nk−1)<br />
< ɛ1 , χ 2 ( Nk−1) = 0 ∀ k . (4.6)<br />
• Wenn der Abstand zwischen zweier aufeinanderfolgender Schritte einen bestimmten<br />
Wert unterschreitet:<br />
|| Nk−1 − Nk|| < ɛ2 . (4.7)<br />
Als Norm wird irgendeine von den p-Normen für endlichdimensionale Vektorräume<br />
R n verwendet. In der Regel wird p = 2 (Euklidische Norm) oder p = ∞<br />
(Maximumsnorm) bevorzugt.<br />
• Ein Abbruchtest, ob mit einem Verfahren ein stationärer Punkt erreicht wurde,<br />
muss wie folgt implementiert werden:<br />
∇χ 2 ( <br />
< ɛ3 falls (−Nk)j < 0 ,<br />
Nk)j =<br />
(4.8)<br />
≥ 0 sonst .<br />
Hierbei handelt es sich um eine numerische Umsetzung der KKT-Bedingungen<br />
aus (3.36).<br />
• Wird keines der oberen Abbruchkriterien erfüllt, dann steigt das Verfahren<br />
nach einer festgelegten Anzahl an Iterationsschritten aus.<br />
Bemerkungen:<br />
• Für ɛ1, ɛ2 und ɛ3 werden üblicherweise Werte aus dem Intervall [10 −2 , 10 −8 ]<br />
gewählt.<br />
• Die ersten zwei Abbruchbedingungen können erst ab der zweiten Iteration des<br />
Linesearch-Verfahrens geprüft werden.<br />
51