Utilizamos cookies propias y de terceros para mejorar la experiencia de navegación, y ofrecer contenidos y publicidad de interés.
Al continuar con la navegación entendemos que se acepta nuestra política de cookies.
Iniciar sesión Cerrar
Correo:
Contraseña:
Entrar
Recordar sesión en este navegador
Iniciar sesiónIniciar sesiónCrear cuentaCrear cuenta

Matlab - Mejor método de resolver sistema no lineal

  Imprimir  
Vista:

Mejor método de resolver sistema no lineal

Publicado por Lorenzo (1 intervención) el 29/05/2012 19:35:39
Hola, tengo un sistema no lineal de ecuaciones, para ello hago lo siguiente:

syms x y z t %Punto inicial para iterar X01=1 ; X02=1; X03=1; X04=1; X=[x y z t]; %Ahora tengo que poner el sistema f(1)=; f(2)=; f(3)=; f(4)=; b=[0 0 0 0]; TOLERANCIA=0.0000000000001; X0=[X01; X02; X03; X04]; % Es el punto inicial que se escoge para iterar; Jf=jacobian(f,X); t=1; j=1; while (t>=TOLERANCIA) && (j<100) A=double(subs(Jf,X,X0)); B=b'-subs(f,X,X0)'+A*X0; S=A\B; X0=S; t=norma(b-subs(f,X,X0)); j=j+1; end

En donde norma.m:

%X tiene que estar en columna function J=norma(X) J=sqrt(X*transpose(X));

Me gustaría saber, si existe un método mejor para resolver ecuaciones no lineales, por ejemplo, si hay alguno en el que no tenga que poner un punto inicial para iterar.

Un saludo y muchas gracias.
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder