Ich versuche das folgende Gleichungssystem für die Variablen und zu lösen (alle anderen sind Konstanten):x 2
Ich kann sehen, dass ich dieses Gleichungssystem in eine einzige Gleichung einer einzigen Variablen umwandeln kann, indem ich die Gleichungen 1 und 2 für bzw. löse und sie in Gleichung 3 ersetze. Dabei kann ich Matlabs verwenden Befehl, um die Lösung zu finden. Unter Verwendung der Parameter , und fand ich die wahre Lösung als .fzero
Wenn ich jedoch Newtons Methode verwende, die auf das ursprüngliche 3-Variablen-3-Gleichungssystem angewendet wird, konvergieren die Iterationen nie zur Lösung, unabhängig davon, wie nahe ich an der wahren Lösung .
Zuerst vermutete ich einen Fehler bei der Implementierung der Newton-Methode. Nach mehrmaliger Überprüfung habe ich keinen Fehler gefunden. Dann habe ich versucht, eine anfängliche Vermutung , und siehe da: Der Jacobian ist singulär. Ich weiß, dass ein singulärer jacobianer die Konvergenzreihenfolge reduzieren kann, aber ich denke nicht, dass er notwendigerweise die Konvergenz zur wahren Lösung verhindert.
Meine Frage ist also, dass der Jacobi des Systems bei der wahren Lösung singulär ist:
Welche anderen Bedingungen sind erforderlich, um zu beweisen, dass die Newtonsche Methode nicht zur Wurzel konvergiert?
Würde eine Globalisierungsstrategie (z. B. Liniensuche) Konvergenz trotz des singulären Jacobi garantieren?