Mehrdimensionales Newton-Verf. < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) reagiert/warte auf Reaktion  | | Datum: | 12:46 Mo 14.06.2010 | | Autor: | leberkas2 |
| Aufgabe | Hallo,
hab folgendes Problem mit der Programmierung des Newton-Verfahrens in MATLAB. (nicht-lineare GLS)
In der Ausgabe sollen sämtliche Iterationsschritte mit Ergebnis angezeigt werden,
die man für's Ausrechnen der Nullstellen benötigt.
Bei mir wird aber nur das Endergibnis (x1=0,5; x2=0,5) angezeigt.
In meinem Beispiel werden genau 4 Schritte benötigt, um auf die Nullstellen zu kommen.
Vielleicht weiss jemand wie ich die Ausgabe aller Schritte in mein Verfahren implementiere...?
Hier seht ihr was ich bisher habe:
%%Nichtlineare Gleichungssysteme mit mehreren Variablen
%%Mehrdimensionales Newton-Verfahren
%%Für eine gegebene Funktion Funktion F(x,y) = [f1(x,y);f2(x,y)]
%%soll in Matlab das Newton-Verfahren implementiert werden.
%Gegeben sei:
%f1 = [mm] x^2+y^2+y-1=0
[/mm]
%f2 = [mm] x^2-y^2+x-y-2=0
[/mm]
%mit dem Startwert x0 = (0;0) |
HIER MEINE BISHERIGE LÖSUNG:
%Zur Vereinfachung werden die Variablen x,y in diesem Beispiel als x(1),x(2)
%angenommen. Aus der Ausgangsfunktion ergibt sich:
f1 = [mm] x(1)^2+x(2)^2+x(2)-1;
[/mm]
f2 = [mm] x(1)^2-x(2)^2+x(1)-x(2)-2;
[/mm]
N=20;
x=[0;0];
for i=1:N
[mm] F=[x(1)^2+x(2)^2+x(2)-1; x(1)^2-x(2)^2+x(1)-x(2)-2];
[/mm]
dF=[2*x(1)+2*x(2)+1; 2*x(1)-2*x(2)];
[mm] x=x-dF\F;
[/mm]
end
x
Vielen Dank schonmal falls Ihr mehr wisst;)
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
| |
|
| Status: |
(Antwort) fertig  | | Datum: | 13:04 Mo 14.06.2010 | | Autor: | leduart |
Hallo
ich seh nicht mal dass das das Newtonverfahren ist. Dein Ergebnis ist ne Nst von f1 nicht von f2,
wenn du keinen Ausgabebefehl gibst, seh ich nicht, wie du ne Ausgabe, auch des Endergebnisses kriegen sollst?
Gruss leduart
|
|
|
|
