Matlab - lotka voterra

   
Vista:

lotka voterra

Publicado por fran (1 intervención) el 19/12/2008 20:46:29
Soy estudiante y necesito saber como resolver este sistema de ecuaciones de lotka-volterra, o algunas indicaciones para sacar la matriz: Y'(t)=AY(t)+f(t)


x'(t)=(a-bx(t)+cy(t))x(t) x(0)=x0
y'(t)=(d-ey(t)+fx(t))y(t) y(0)=y0

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

RE:lotka voterra matematicas

Publicado por angelica (1 intervención) el 07/01/2009 00:00:13
yo tambien soy estudiantey tengo eso mismo trabajo quiero saber si tienes idea como resolverlo.
gracias
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

RE:lotka voterra matematicas

Publicado por ramiro (54 intervenciones) el 08/01/2009 00:51:54
Angélica:
En los "demos" de matlab, a los que se puede acceder tecleando la palabra demos en la señal de inicio >> del matlab, siguiendo la secuencia:
matlab->numerics->numerical integration of differential equations
Aparece precisamente la ecuación de Lotka-Volterra como ejemplo, aunque hay que decir que ligeramente modificada, con respecto a la que ustedes plantean.
Les recomiendo seguir el ejemplo, haciendo las adecuaciones respectivas.
Por ejemplo, en el archivo:
function yp = lotka(t,y)
%LOTKA Lotka-Volterra predator-prey model.

% Copyright 1984-2002 The MathWorks, Inc.
% $Revision: 5.7 $ $Date: 2002/04/08 20:04:46 $

yp = diag([1 - .01*y(2), -1 + .02*y(1)])*y;
Modificar el nombre "lotka" a alguno similar, que podría ser "lotkag" o cualquier otro que sea más de su agrado.
Y la última linea cambiarla por:
yp=diag([a-b*y(1)+c*y(2),d-e*y(2)+f*y(1)]*y
En otras palabras, la matriz A que preguntan que cómo se definiría, quedaría así:
A=diag([a-b*y(1)+c*y(2),d-e*y(2)+f*y(1)])
Lo demás lo pueden copiar y pegar del demo.
Saludos y ojalá les sirva el tip.
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