Jacobimatrix bei Newtonverfahr < Matlab < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) beantwortet  | | Datum: | 03:32 Di 02.08.2011 | | Autor: | Hallo05 |
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
Hallo,
ich habe folgendes Problem. Ich habe eine Funtion $ f: [mm] \IR^{n} \rightarrow \IR^{n} [/mm] $, die als char übergeben wird (durch z.B. [mm] newton2('[x^2;y^3*x]',[1;1])). [/mm] Leider funktioniert dann der Befehl jacobian nicht mehr, da dort ein Vektor verlangt wird.
Nun hab ich versucht mittels eval(f) den char in einen Vektor umzuwandeln, aber dann ist die Variable 'x' unbekannt.
Mein Programmierversuch:
function [xa,u]=newton2(f,v);
g1=inline(f);
J=inline(jacobian(f)); oder J=inline(jacobian(eval(f)));
x0=v;
y=J(x0(1),x0(2))\ g1(x0(1),x0(2));
xa=x0-y;
r=norm(xa-x0);
while abs(r>0.000005)
x0=xa;
y=J(x0(1),x0(2))\ g1(x0(1),x0(2));
xa=x0-y;
r=norm(xa-x0);
end
u=g1(xa(1),xa(2));
Was könnte ich jetzt ändern, dass ich dort die Jacobimatrix herausbekomme?
Vielen Dank
|
|
| |
|
Hallo Hallo05,
![[willkommenmr]](/images/smileys/willkommenmr.png "[willkommenmr]")
> Ich habe diese Frage in keinem Forum auf anderen
> Internetseiten gestellt.
>
>
> Hallo,
> ich habe folgendes Problem. Ich habe eine Funtion [mm]f: \IR^{n} \rightarrow \IR^{n} [/mm],
> die als char übergeben wird (durch z.B.
> [mm]newton2('[x^2;y^3*x]',[1;1])).[/mm] Leider funktioniert dann der
> Befehl jacobian nicht mehr, da dort ein Vektor verlangt
> wird.
> Nun hab ich versucht mittels eval(f) den char in einen
> Vektor umzuwandeln, aber dann ist die Variable 'x'
> unbekannt.
>
>
> Mein Programmierversuch:
> function [xa,u]=newton2(f,v);
> g1=inline(f);
> J=inline(jacobian(f)); oder J=inline(jacobian(eval(f)));
> x0=v;
> y=J(x0(1),x0(2))\ g1(x0(1),x0(2));
> xa=x0-y;
> r=norm(xa-x0);
> while abs(r>0.000005)
> x0=xa;
> y=J(x0(1),x0(2))\ g1(x0(1),x0(2));
> xa=x0-y;
> r=norm(xa-x0);
> end
> u=g1(xa(1),xa(2));
>
> Was könnte ich jetzt ändern, dass ich dort die
> Jacobimatrix herausbekomme?
>
Die Syntax der Jacobimatrix ist doch: jacobian(f,v).
Siehe Dir dazu ein ![[]](/images/popup.gif) Beispiel an.
> Vielen Dank
Gruss
MathePower
|
|
|
|
