Matlab - Sistema de ecuaciones diferenciales

   
Vista:

Sistema de ecuaciones diferenciales

Publicado por Cesar (2 intervenciones) el 06/01/2014 15:37:40
Hola, he tratado durante dias crear un programa con Rungge Kutta para solucioner un sistema de 3 ecuaciones y de segundo orden...Alguien tendra uno redactado parecido? Les coloco el que pude realizar que es solo de dos ecuaciones y de primer grado, espero poder partir de ahi.. Si tienen uno parecido bien sea de otro orden les agradeceria. Saludos

function [X1,Y1,T]=RK4Dosvariables
h=0.2;
n=2;
x=-1;
y=6;
%z=5
%xprima=1
%yprima=2
%zprima=3
Fx=inline('2*x+4*y','x','y');
Fy=inline('-x+6*y','x','y');
%Fz=inline(' ,x,y,z )

m1=Fx(x,y);
k1=Fy(x,y);;
m2=Fx(x+((m1*h)/2),y+((k1*h)/2))
k2=Fy(x+(m1*h)/2,y+(k1*h)/2)
m3=Fx(x+(m2*h)/2,y+(k2*h)/2)
k3=Fy(x+(m2*h)/2,y+(k2*h)/2)
m4=Fx(x+h*m3,y+h*k3)
k4=Fy(x+h*m3,y+h*k3)

%dx/dt
x1=x+(h/6)*(m1+2*m2+2*m3+m4)

%dy/dt
y1=y+(h/6)*(k1+2*k2+2*k3+k4)

x=x1;
y=y1;
T(n)=i+h;
n=n+1;
end
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
Imágen de perfil de JOSE JEREMIAS CABALLERO

Sistema de ecuaciones diferenciales

Publicado por JOSE JEREMIAS CABALLERO jjcc94@hotmail.com (3441 intervenciones) el 06/01/2014 16:25:26
una sugerencia.
Primero haz para una sola ecuación el método de Runge Kutta
luego para dos ecuaciones y observa que modificaciones haz hecho
luego para tres ecuaciones siguiendo el mismo análisis para dos ecuaciones.
...
luego generalizas para n ecuaciones.
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
0
Comentar