Matlab - metodo de euler

 
Vista:

metodo de euler

Publicado por nico (1 intervención) el 01/12/2006 19:16:48
Hola, vereis, el caso es que tengo que resolver unos problemas muy sencillitos y no tengo ninguna experiencia con MATLAB. No tengo tiempo y me estoy volviendo loco con tantos tutoriales (a los que diho sea de paso les podrían añadir algun que otro ejemplo más). Se trata de resolver unas ecuaciones diferenciales con el método de Euler. Hasta ahora he estado siguiendo el libro "Métodos numéricos con MATLAB" de John Mathews, y a pesar de que viene hasta el algoritmo, no hay manera de que a mi me salga. Me da un montón de errores que no entiendo, consulto la ayuda y entiendo todavía menos. La verdad es que estoy empezando a dudar si debería seguir pagando el alquiler de mi habitación o si sería más barato irme a vivir al zoo. Por favor ¿Alguien puede ayudarme?

Os describo los problemas:

-->El primero consiste en resolver un problema del tipo:

dx/dt=[(a*t+b)/(c-d*t)]*x donde a,b,c,d son reales

-->El segundo es un sistema de dos ecuaciones diferenciales acopladas:

dY/dt=F(t)-Y

dX/dt=Y+a*F(t)-b*X-d*X*F(t)

El problema fundamental es que no sé como definir las funciones, porque lo haga como lo haga nunca está bien. De todas formas si pudieseis decirme como se resuelve al menos una, o por lo menos darme un ejemplo lo agradecería.

Muchas gracias y un saludo

Nico
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
sin imagen de perfil

RE:metodo de euler

Publicado por Raúl Santiago (211 intervenciones) el 12/12/2006 10:52:51
Un ejemplo

alfa=1000;
lambda=5000;


t=0:0.0001:1;

tant=0;
sigma=0;
psi=1;
for i=1:length(t)

dsigma=-sqrt(sigma)+psi;
if dsigma>0
dpsi=-(1/lambda^2)*(alfa-C);
else
dpsi=-(1/lambda^2)*(alfa+C);
end
sigma=sigma+dsigma*(t(i)-tant);
psi=psi+dpsi*(t(i)-tant);
sigmag(i)=sigma;
psig(i)=psi;
tant=t(i);
end
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