Hintergrund
Ich löse eine Variante der Ornstein-Zernike- Gleichung aus der Flüssigkeitstheorie. Abstrakt kann das Problem als Lösung des Fixpunktproblems , wobei A ein integro-algebraischer Operator ist und c ( r ) die Lösungsfunktion ist (die OZ-Direktkorrelationsfunktion). Ich löse durch Picard-Iteration, wobei ich eine anfängliche Versuchslösung c 0 ( r ) bereitstelle und neue Versuchslösungen nach dem Schema c j + 1 = α ( wobei α einen einstellbaren Parameter steuertdass die Mischung aus c und A c in der nächsten Probelösung verwendet. Nehmen wir für diese Diskussion an, dass der Wert von α unwichtig ist. I wiederholenbis die Iteration konvergiertum innerhalb einer gewünschten Toleranz ε : Δ j + 1 ≡ ∫ d → r | c j + 1 ( r ) - c
Für einen weiten Bereich von Werten für konvergiert das obige Iterationsschema exponentiell schnell. Wenn ich jedoch λ verkleinere , erreiche ich schließlich ein Regime, in dem die Konvergenz nicht monoton ist (siehe Abbildung unten).
Schlüsselfrage
Hat die nicht-monotone Konvergenz bei iterativen Lösungen von Fixpunktproblemen eine besondere Bedeutung? Bedeutet das, dass mein iteratives Schema kurz vor der Instabilität steht? Sollte mich vor allem die nicht-monotone Konvergenz verdächtigen, dass die "konvergierte" Lösung keine gute Lösung für das Fixpunktproblem ist?
quelle