Lotka-Volterra predator-prey system - Ecuaciones del modelo Lotka-Volterra - Dinamica de sistemas

 

https://metasd.com/2011/06/lotka-volterra/

https://pybonacci.org/2015/01/05/ecuaciones-de-lotka-volterra-modelo-presa-depredador/

Lotka-Volterra predator-prey system

The Lotka-Volterra equations, which describe a predator-prey system, must be one of the more famous dynamic systems. There are many generalizations and applications outside of biology.

Wikipedia has a nice article, which I used as the basis for this simple model.

You can find an equilibrium at (Relative initial predators = 3, Relative initial prey = 0.5). An interesting experiment is to determine whether that point is stable.

Lotka-Volterra.vpm (runs in any Vensim version)

Lotka-Volterra-Alt.vpm (an alternative version that enlarges the diagram slightly and moves most of the graphs to a control panel on a second view)

---

Another version of the model, which requires an advanced version of Vensim or the Model Reader, uses a subscript (array) to randomize the initial conditions for the predator and prey levels. That allows you to see the effect of the parameters over a variety of random realizations in the predator-prey phase space.

Lotka-Volterra+Array.vpm (requires Vensim DSS or Pro, or the free Model Reader)

Update: I wrote a bit more on this topic here.

https://pybonacci.org/2015/01/05/ecuaciones-de-lotka-volterra-modelo-presa-depredador/

Saltar al contenido

• Inicio
• Foro
• Acerca de Pybonacci
• Contacto
• RunPy
• /

Ecuaciones de Lotka-Volterra: modelo presa depredador

• por Álex Sáez
• 2019-11-17
• Artículos
• 13 comentarios

Etiquetas:ecuaciones diferencialesEDOsmatplotlib.pyplotpython3

Introducción

Resulta intuitivo pensar que las poblaciones de un animal depredador y su presa están relacionadas de algún modo en el que si una aumenta, la otra lo hace también. Utilizaremos como ejemplo en este artículo un ecosistema aislado y formado por leones y cebras que viven en armonía, es decir, los leones comiéndose a las cebras. Imaginemos que por cualquier circunstancia, por ejemplo, por disponer de mayor cantidad de alimento, aumenta la población de cebras; los leones dispondrán de más alimento y su población aumentará, pero ¿qué ocurrirá a partir de este momento? Si la población de leones llega a ser demasiado grande para el número de cebras en nuestra sabana, podrían acabar con todas, provocando su propia extinción por inanición. Pero incluso si el festín no es tan grande como para comerse todas las cebras, pero sí para dejar una población muy mermada, probablemente los leones tendrán que pasar hambre una buena temporada y algunos de ellos morirán hasta que las cebras tengan tiempo suficiente para reproducirse y volver a ser pasto de los leones. ¿Cuántas cebras morirán en el atracón? ¿Cuánto tiempo pasarán los leones hambre? ¿Cuántos morirán?

Ecuaciones

Las ecuaciones de Lotka-Volterra son un modelo biomatemático que pretende responder a estas cuestiones prediciendo la dinámica de las poblaciones de presa y depredador bajo una serie de hipótesis:

• El ecosistema está aislado: no hay migración, no hay otras especies presentes, no hay plagas…
• La población de presas en ausencia de depredadores crece de manera exponencial: la velocidad de reproducción es proporcional al número de individuos. Las presas sólo mueren cuando son cazadas por el depredador.
• La población de depredadores en ausencia de presas decrece de manera exponencial.
• La población de depredadores afecta a la de presas haciéndola decrecer de forma proporcional al número de presas y depredadores (esto es como decir de forma proporcional al número de posibles encuentros entre presa y depredador).
• La población de presas afecta a la de depredadores también de manera proporcional al número de encuentros, pero con distinta constante de proporcionalidad (dependerá de cuanto sacien su hambre los depredadores al encontrar una presa).

Se trata de un sistema de dos ecuaciones diferenciales de primer orden, acopladas, autónomas y no lineales:
dxdt=αx–βxy
dydt=−γy+δyx
donde x es el número de presas (cebras en nuestro caso) e y es el número de depredadores (leones). Los parámetros son constantes positivas que representan:

• α: tasa de crecimiento de las presas.
• β: éxito en la caza del depredador.
• γ: tasa de decrecimiento de los depredadores.
• δ: éxito en la caza y cuánto alimenta cazar una presa al depredador.

Resolución

Resolveremos este sistema en Python usando la función odeint de scipy.integrate (puedes ver cómo funciona en el artículo El salto de Felix Baumgartner en Python) y representaremos los resultados con matplotlib. Para esto usaremos: Python 3.4, numpy 1.9.0, matplotlib 1.4.0, scipy 0.14.0.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

import numpy as np import matplotlib.pyplot as plt from scipy.integrate import odeint def df_dt(x, t, a, b, c, d):     """Función del sistema en forma canónica"""     dx = a * x[0] - b * x[0] * x[1]     dy = - c * x[1] + d * x[0] * x[1]     return np.array([dx, dy]) # Parámetros a = 0.1 b = 0.02 c = 0.3 d = 0.01 # Condiciones iniciales x0 = 40   # Presas y0 = 9    # Depredadores conds_iniciales = np.array([x0, y0]) # Condiciones para integración tf = 200 N = 800 t = np.linspace(0, tf, N) solucion = odeint(df_dt, conds_iniciales, t, args=(a, b, c, d))

Ahora que ya tenemos la solución, sólo tenemos que pintarla para ver cómo han ido las cosas en nuestra sabana virtual.

1 2	plt.plot(t, solucion[:, 0], label='presa') plt.plot(t, solucion[:, 1], label='depredador')

Vemos que se trata de una solución periódica en la que, como decíamos al principio, un aumento en la población de cebras va seguido del aumento del número de leones. Un gran número de depredadores merma la población de presas y a los pobres leones les toca pasar hambre una temporada. Otra forma interesante de visualizar estos datos es ver el número de presas en función del número de depredadores, en lugar de a lo largo del tiempo, es decir, podemos visualizar su mapa de fases. Podemos pintar también el campo de direcciones de nuestras ecuaciones usando la función quiver. El tamaño de las flechas se ha normalizado para que todas tengan la misma longitud y se ha usado un colormap para representar el módulo.

1 2 3 4 5 6 7 8 9 10 11

x_max = np.max(solucion[:,0]) * 1.05 y_max = np.max(solucion[:,1]) * 1.05 x = np.linspace(0, x_max, 25) y = np.linspace(0, y_max, 25) xx, yy = np.meshgrid(x, y) uu, vv = df_dt((xx, yy), 0, a, b, c, d) norm = np.sqrt(uu**2 + vv**2) uu = uu / norm vv = vv / norm plt.quiver(xx, yy, uu, vv, norm, cmap=plt.cm.gray) plt.plot(solucion[:, 0], solucion[:, 1])

Si nos fijamos en la línea azul, la coordenada $latex x$ en cada punto indica el número de presas y la coordenada $latex y$ el número de depredadores. La evolución a lo largo del tiempo que hemos representado antes, se obtiene al recorrer esta curva en sentido antihorario. Podemos ver también como el campo de direcciones nos señala la tendencia del sistema en cada situación. Por ejemplo, una flecha que apunta hacia arriba a la derecha indica que con ese número de cebras y leones en nuestra sabana, la tendencia será que aumenten ambos.

Llegados a este punto podemos preguntarnos qué habría ocurrido si el número inicial de cebras y leones hubiese sido otro. Como ya sabemos integrar ecuaciones diferenciales, bastaría con cambiar nuestra x0 e y0 y repetir el proceso (incluso podríamos hacer un widget interactivo). Sin embargo, se puede demostrar que a lo largo de las líneas del mapa de fases, como la que hemos pintado antes, se conserva la cantidad:

C=αlny–βy+γlnx−δx

Por tanto, pintando un contour de esta cantidad podemos obtener la solución para distintos valores iniciales del problema.

1 2 3 4 5 6 7 8 9

def C(x, y, a, b, c, d):     return a * np.log(y) - b * y + c * np.log(x) - d * x fig, ax = plt.subplots(1,2) ax[0].plot(solucion[:, 0], solucion[:, 1], lw=2, alpha=0.8) ax[0].scatter(c/d, a/b) levels = (0.5, 0.6, 0.7, 0.72, 0.73, 0.74, 0.75, 0.76, 0.77, 0.775, 0.78, 0.781) ax[0].contour(xx, yy, constant, levels, colors='blue', alpha=0.3) ax[1].plot(t, solucion[:, 0], label='presa') ax[1].plot(t, solucion[:, 1], label='depredador')

Vemos que estas curvas se van haciendo cada vez más y más pequeñas, hasta que, en nuestro caso, colapsarían en un punto en torno a (30,5). Se trata de un punto de equilibrio o punto crítico; si el sistema lo alcanzase, no evolucionaría y el número de cebras y leones sería constante en el tiempo. El otro punto crítico de nuestro sistema es el (0,0). Analizándolos matemáticamente se obtiene que:

• El punto crítico situado en (0,0) es un punto de silla. Al tratarse de un punto de equilibrio inestable la extinción de cualquiera de las dos especies en el modelo sólo puede conseguirse imponiendo la condición inicial nula.
• El punto crítico situado en (γ/δ,α/β) es un centro (en este caso los autovalores de la matriz del sistema linealizado son ambos imaginarios puros, por lo que a priori no se conoce su estabilidad).

Mejorando el modelo

Como se puede observar, este modelo tiene algunas deficiencias propias de su simplicidad y derivadas de las hipótesis bajo las que se ha formulado. Una modificación razonable es cambiar el modelo de crecimiento de las presas en ausencia de depredadores, suponiendo que en vez de aumentar de forma exponencial, lo hacen según una función logística. Esta curva crece de forma similar a una exponencial al principio, moderándose después y estabilizándose asintóticamente en un valor. Este modelo de crecimiento representa mejor las limitaciones en el número de presas debidas al medio (falta de alimento, territorio…). Llevando este modelo de crecimiento a las ecuaciones originales se tiene un nuevo sistema en el que interviene un parámetro más:

dxdt=(αx–rx2)–βxy dydt=−γy+δyx

Integrándolo como hemos hecho antes, el resultado obtenido es:

En este caso se puede observar como el comportamiento deja de ser periódico. El punto crítico que antes era un centro, se convierte en un atractor y la solución tiende a estabilizarse en un número fijo de presas y depredadores.

Referencias

Si tienes curiosidad sobre como seguir perfeccionando este modelo o cómo incluir otras especies, quizás quieras echar un vistazo a:

• Competitive Lotka–Volterra equations o The Predator-Prey Equations.

• Presentación ETSIINF-UPM.

• Apuntes Ecuaciones Diferenciales, Bartolomé Luque Serrano (ETSIA-UPM).
• Si te interesa ver cómo realizar la integración con diferentes métodos, puedes visitar Predator Prey Model – Bank Assignment of Numerical Mooc.

Actualización

En GitHub puedes encontrar el notebook con el que se ha hecho esta entrada. Si no tienes instalado el Notebook de IPython, puedes pinchar aquí para verlo con nbviewer (aunque no podrás jugar con los widgets interactivos).

ANTERIOR¿Quieres empezar con Python?

SIGUIENTEComprobando presencia de subconjuntos dentro de conjunto de forma eficiente