stringtranslate.com

Ecuaciones de Lotka-Volterra

Las ecuaciones de Lotka-Volterra , también conocidas como modelo depredador-presa de Lotka-Volterra , son un par de ecuaciones diferenciales no lineales de primer orden , frecuentemente utilizadas para describir la dinámica de sistemas biológicos en los que interactúan dos especies, una como depredadora y la otra como presa. Las poblaciones cambian a través del tiempo según el par de ecuaciones:

dónde

La solución de las ecuaciones diferenciales es determinista y continua . Esto, a su vez, implica que las generaciones tanto del depredador como de la presa se superponen continuamente. [1]

El sistema de ecuaciones de Lotka-Volterra es un ejemplo de un modelo de población de Kolmogorov (que no debe confundirse con las ecuaciones de Kolmogorov más conocidas ), [2] [3] [4] que es un marco más general que puede modelar la dinámica de los sistemas ecológicos con interacciones depredador-presa, competencia , enfermedades y mutualismo .

Interpretación biológica y supuestos del modelo

Se supone que las presas tienen un suministro ilimitado de alimentos y se reproducen exponencialmente, a menos que estén sujetas a depredación; este crecimiento exponencial está representado en la ecuación anterior por el término αx . Se supone que la tasa de depredación de la presa es proporcional a la tasa a la que los depredadores y la presa se encuentran; esto se representa arriba por βxy . Si x o y es cero, entonces no puede haber depredación. Con estos dos términos, la ecuación de la presa anterior puede interpretarse de la siguiente manera: la tasa de cambio de la población de la presa está dada por su propia tasa de crecimiento menos la tasa a la que es depredada.

El término δxy representa el crecimiento de la población de depredadores. (Obsérvese la similitud con la tasa de depredación; sin embargo, se utiliza una constante diferente, ya que la tasa a la que crece la población de depredadores no es necesariamente igual a la tasa a la que consume la presa). El término γy representa la tasa de pérdida de los depredadores debido a la muerte natural o la emigración; conduce a una disminución exponencial en ausencia de presas. Por lo tanto, la ecuación expresa que la tasa de cambio de la población del depredador depende de la tasa a la que consume presas, menos su tasa de muerte intrínseca.

El modelo depredador-presa de Lotka-Volterra hace una serie de suposiciones sobre el medio ambiente y la biología de las poblaciones de depredadores y presas: [5]

  1. La población de presas encuentra abundante alimento en todo momento.
  2. El suministro de alimentos de la población de depredadores depende enteramente del tamaño de la población de presas.
  3. La tasa de cambio de la población es proporcional a su tamaño.
  4. Durante el proceso, el entorno no cambia a favor de una especie y la adaptación genética es intrascendente.
  5. Los depredadores tienen un apetito ilimitado.
  6. Ambas poblaciones pueden describirse mediante una única variable, lo que equivale a suponer que las poblaciones no tienen una distribución espacial o por edad que contribuya a la dinámica.

Relevancia biológica del modelo

Cantidad de pieles de liebre de raquetas de nieve (amarillo, fondo) y de lince canadiense (línea negra, primer plano) vendidas a la Compañía de la Bahía de Hudson . Los linces canadienses se comen a las liebres de raquetas de nieve.

Es poco probable que ninguno de los supuestos anteriores se cumpla en el caso de las poblaciones naturales. No obstante, el modelo de Lotka-Volterra muestra dos propiedades importantes de las poblaciones de depredadores y presas, y estas propiedades suelen extenderse a las variantes del modelo en las que estos supuestos se relajan:

En primer lugar, la dinámica de las poblaciones de depredadores y presas tiende a oscilar. Se han observado fluctuaciones en el número de depredadores y presas en poblaciones naturales, como en los datos de la Compañía de la Bahía de Hudson sobre linces y liebres americanas [6] y en las poblaciones de alces y lobos del Parque Nacional Isle Royale [7] .

En segundo lugar, el equilibrio poblacional de este modelo tiene la propiedad de que la densidad de equilibrio de la presa (dada por ) depende de los parámetros del depredador, y la densidad de equilibrio del depredador (dada por ) de los parámetros de la presa. Esto tiene como consecuencia que un aumento en, por ejemplo, la tasa de crecimiento de la presa, , conduce a un aumento en la densidad de equilibrio del depredador, pero no en la densidad de equilibrio de la presa. Hacer que el entorno sea mejor para la presa beneficia al depredador, no a la presa (esto está relacionado con la paradoja de los pesticidas y con la paradoja del enriquecimiento ). Una demostración de este fenómeno la proporciona el aumento del porcentaje de peces depredadores capturados que había aumentado durante los años de la Primera Guerra Mundial (1914-18), cuando la tasa de crecimiento de la presa aumentó debido a un menor esfuerzo pesquero.

Otro ejemplo es la fertilización experimental del océano con hierro. En varios experimentos se disolvieron grandes cantidades de sales de hierro en el océano. La expectativa era que el hierro, que es un nutriente limitante para el fitoplancton, impulsaría el crecimiento del fitoplancton y que secuestraría dióxido de carbono de la atmósfera. La adición de hierro generalmente conduce a una floración corta en el fitoplancton, que es rápidamente consumido por otros organismos (como peces pequeños o zooplancton ) y limita el efecto del enriquecimiento principalmente al aumento de la densidad de depredadores, que a su vez limita el secuestro de carbono. Esto es como lo predicen las densidades de población de equilibrio del modelo depredador-presa de Lotka-Volterra, y es una característica que se traslada a modelos más elaborados en los que se relajan los supuestos restrictivos del modelo simple. [8]

Aplicaciones a la economía y al marketing

El modelo de Lotka-Volterra tiene aplicaciones adicionales en áreas como la economía [9] y el marketing. [10] [11] Puede utilizarse para describir la dinámica de un mercado con varios competidores, plataformas y productos complementarios, una economía colaborativa, etc. Existen situaciones en las que uno de los competidores expulsa a los demás competidores del mercado y otras situaciones en las que el mercado alcanza un equilibrio en el que cada empresa se estabiliza en su cuota de mercado. También es posible describir situaciones en las que hay cambios cíclicos en la industria o situaciones caóticas sin equilibrio y los cambios son frecuentes e impredecibles.

Historia

El modelo depredador-presa de Lotka-Volterra fue propuesto inicialmente por Alfred J. Lotka en la teoría de reacciones químicas autocatalíticas en 1910. [12] [13] Esta era efectivamente la ecuación logística , [14] originalmente derivada por Pierre François Verhulst . [15] En 1920 Lotka extendió el modelo, a través de Andrey Kolmogorov , a "sistemas orgánicos" usando una especie de planta y una especie de animal herbívoro como ejemplo [16] y en 1925 usó las ecuaciones para analizar las interacciones depredador-presa en su libro sobre biomatemáticas . [17] El mismo conjunto de ecuaciones fue publicado en 1926 por Vito Volterra , un matemático y físico, que se había interesado en la biología matemática . [13] [18] [19] La investigación de Volterra se inspiró en sus interacciones con el biólogo marino Umberto D'Ancona , que cortejaba a su hija en ese momento y que más tarde se convertiría en su yerno. D'Ancona estudió las capturas de peces en el mar Adriático y había notado que el porcentaje de peces depredadores capturados había aumentado durante los años de la Primera Guerra Mundial (1914-18). Esto lo desconcertó, ya que el esfuerzo pesquero se había reducido mucho durante los años de guerra y, como los peces presa son la captura preferida, uno esperaría intuitivamente que esto aumentara el porcentaje de peces presa. Volterra desarrolló su modelo para explicar la observación de D'Ancona y lo hizo independientemente de Alfred Lotka. Dio crédito al trabajo anterior de Lotka en su publicación, después de lo cual el modelo se conoció como el "modelo Lotka-Volterra". [20]

El modelo se amplió posteriormente para incluir el crecimiento de la presa dependiente de la densidad y una respuesta funcional de la forma desarrollada por CS Holling ; un modelo que se conoció como el modelo de Rosenzweig-MacArthur. [21] Tanto el modelo de Lotka-Volterra como el de Rosenzweig-MacArthur se han utilizado para explicar la dinámica de las poblaciones naturales de depredadores y presas.

A finales de la década de 1980, surgió una alternativa al modelo depredador-presa de Lotka-Volterra (y sus generalizaciones dependientes de la presa común): el modelo dependiente de la proporción o modelo de Arditi-Ginzburg . [22] La validez de los modelos dependientes de la presa o de la proporción ha sido muy debatida. [23]

Las ecuaciones de Lotka-Volterra tienen una larga historia de uso en la teoría económica ; su aplicación inicial se atribuye comúnmente a Richard Goodwin en 1965 [24] o 1967. [25] [26]

Soluciones a las ecuaciones

Las ecuaciones tienen soluciones periódicas . Estas soluciones no tienen una expresión sencilla en términos de las funciones trigonométricas habituales , aunque son bastante manejables. [27] [28] [29]

Si ninguno de los parámetros no negativos α , β , γ , δ se anula, tres de ellos pueden ser absorbidos en la normalización de variables para dejar sólo un parámetro: como la primera ecuación es homogénea en x , y la segunda en y , los parámetros β / α y δ / γ son absorbibles en las normalizaciones de y y x respectivamente, y γ en la normalización de t , de modo que sólo α / γ permanece arbitrario. Es el único parámetro que afecta a la naturaleza de las soluciones.

Una linealización de las ecuaciones produce una solución similar al movimiento armónico simple [30], con la población de depredadores 90° por detrás de la de presas en el ciclo.

Un ejemplo sencillo

Se menciona al margen el problema de la dinámica poblacional del conejo y el zorro.
Gráfico de espacio de fases para el problema depredador-presa para diversas condiciones iniciales de la población de depredadores.

Supongamos que hay dos especies de animales, un conejo (presa) y un zorro (depredador). Si las densidades iniciales son 10 conejos y 10 zorros por kilómetro cuadrado, se puede representar gráficamente la progresión de las dos especies a lo largo del tiempo; dados los parámetros de que las tasas de crecimiento y mortalidad de los conejos son 1,1 y 0,4, mientras que las de los zorros son 0,1 y 0,4, respectivamente. La elección del intervalo de tiempo es arbitraria.

También se pueden representar gráficamente soluciones paramétricamente como órbitas en el espacio de fases , sin representar el tiempo, pero con un eje representando el número de presas y el otro eje representando las densidades de depredadores para todos los tiempos.

Esto corresponde a eliminar el tiempo de las dos ecuaciones diferenciales anteriores para producir una única ecuación diferencial.

relacionando las variables x (depredador) e y (presa). Las soluciones de esta ecuación son curvas cerradas. Es susceptible de separación de variables : integración

produce la relación implícita

donde V es una cantidad constante que depende de las condiciones iniciales y se conserva en cada curva.

Un aparte: estos gráficos ilustran una limitación potencial grave en la aplicación como modelo biológico: para esta elección específica de parámetros, en cada ciclo, la población de conejos se reduce a números extremadamente bajos, pero se recupera (mientras que la población de zorros sigue siendo considerable en la densidad de conejos más baja). Sin embargo, en situaciones de la vida real, las fluctuaciones aleatorias de los números discretos de individuos podrían hacer que los conejos se extingan y, en consecuencia, también los zorros. Este problema de modelado se ha denominado el "problema del zorro atto", siendo un zorro atto un 10 −18 nocional . [31] [32] Una densidad de 10 −18 zorros por kilómetro cuadrado equivale a un promedio de aproximadamente 5×10 −10 zorros en la superficie de la Tierra, lo que en términos prácticos significa que los zorros están extintos.

Estructura hamiltoniana del sistema

Puesto que la cantidad se conserva a lo largo del tiempo, desempeña el papel de una función hamiltoniana del sistema. [33] Para ver esto, podemos definir el corchete de Poisson de la siguiente manera . Entonces , las ecuaciones de Hamilton se leen Las variables y no son canónicas, ya que . Sin embargo, utilizando transformaciones [34] y llegamos a una forma canónica de las ecuaciones de Hamilton que presenta el hamiltoniano : El corchete de Poisson para las variables canónicas ahora toma la forma estándar .

Diagrama de espacio de fases de otro ejemplo

Un ejemplo menos extremo cubre:

α = 2/3 , β = 4/3 , γ = 1 = δ . Suponga que x e y se cuantifican en miles cada uno. Los círculos representan las condiciones iniciales de la presa y el depredador desde x = y = 0,9 a 1,8, en pasos de 0,1. El punto fijo está en (1, 1/2).

Dinámica del sistema

En el sistema modelo, los depredadores prosperan cuando hay presas abundantes, pero, en última instancia, superan su oferta de alimentos y disminuyen. A medida que la población de depredadores es baja, la población de presas aumentará nuevamente. Esta dinámica continúa en un ciclo de población de crecimiento y disminución.

Equilibrio poblacional

El equilibrio poblacional ocurre en el modelo cuando ninguno de los niveles poblacionales cambia, es decir, cuando ambas derivadas son iguales a 0:

El sistema de ecuaciones anterior produce dos soluciones: y

Por lo tanto, hay dos equilibrios.

La primera solución representa efectivamente la extinción de ambas especies. Si ambas poblaciones están en 0, entonces seguirán así indefinidamente. La segunda solución representa un punto fijo en el que ambas poblaciones mantienen sus números actuales, distintos de cero, y, en el modelo simplificado, lo hacen indefinidamente. Los niveles de población en los que se alcanza este equilibrio dependen de los valores elegidos de los parámetros α , β , γ y δ .

Estabilidad de los puntos fijos

La estabilidad del punto fijo en el origen se puede determinar realizando una linealización utilizando derivadas parciales .

La matriz jacobiana del modelo depredador-presa es y se conoce como la matriz comunitaria .

Primer punto fijo (extinción)

Cuando se evalúa en el estado estable de (0, 0) , la matriz jacobiana J se convierte en

Los valores propios de esta matriz son

En el modelo, α y γ son siempre mayores que cero y, por lo tanto, el signo de los valores propios anteriores siempre será diferente. Por lo tanto, el punto fijo en el origen es un punto de silla .

La inestabilidad de este punto fijo es importante. Si fuera estable, poblaciones distintas de cero podrían ser atraídas hacia él y, como tal, la dinámica del sistema podría conducir a la extinción de ambas especies para muchos casos de niveles poblacionales iniciales. Sin embargo, como el punto fijo en el origen es un punto de silla y, por lo tanto, inestable, se deduce que la extinción de ambas especies es difícil en el modelo. (De hecho, esto solo podría ocurrir si las presas fueran erradicadas artificialmente por completo, lo que causaría que los depredadores murieran de hambre. Si los depredadores fueran erradicados, la población de presas crecería sin límite en este modelo simple). Las poblaciones de presas y depredadores pueden acercarse infinitesimalmente a cero y aún así recuperarse.

Segundo punto fijo (oscilaciones)

Evaluar J en el segundo punto fijo conduce a

Los valores propios de esta matriz son

Como los valores propios son puramente imaginarios y conjugados entre sí, este punto fijo debe ser un centro para órbitas cerradas en la vecindad local o una espiral atractiva o repulsiva. En sistemas conservativos, debe haber órbitas cerradas en la vecindad local de puntos fijos que existan en los mínimos y máximos de la cantidad conservada. La cantidad conservada se deduce anteriormente de que está en órbitas. Por lo tanto, las órbitas alrededor del punto fijo son cerradas y elípticas , por lo que las soluciones son periódicas, oscilando en una pequeña elipse alrededor del punto fijo, con una frecuencia y un período .

Como se ilustra en las oscilaciones circulantes en la figura anterior, las curvas de nivel son órbitas cerradas que rodean el punto fijo: los niveles de las poblaciones de depredadores y presas circulan y oscilan sin amortiguamiento alrededor del punto fijo con una frecuencia .

El valor de la constante de movimiento V , o, equivalentemente, K = exp(− V ) , , se puede encontrar para las órbitas cerradas cerca del punto fijo.

El aumento de K acerca la órbita cerrada al punto fijo. El valor máximo de la constante K se obtiene resolviendo el problema de optimización. El valor máximo de K se alcanza en el punto estacionario (fijo) y asciende a donde e es el número de Euler .

Véase también

Notas

  1. ^ Cooke, D.; Hiorns, RW; et al. (1981). La teoría matemática de la dinámica de las poblaciones biológicas . Vol. II. Academic Press.
  2. ^ Freedman, HI (1980). Modelos matemáticos deterministas en ecología de poblaciones . Marcel Dekker .
  3. ^ Brauer, F.; Castillo-Chavez, C. (2000). Modelos matemáticos en biología de poblaciones y epidemiología . Springer-Verlag .
  4. ^ Hoppensteadt, F. (2006). "Modelo depredador-presa". Scholarpedia . 1 (10): 1563. Bibcode :2006SchpJ...1.1563H. doi : 10.4249/scholarpedia.1563 .
  5. ^ "DINÁMICA DEPREDADOR-PRESA". www.tiem.utk.edu . Archivado desde el original el 15 de diciembre de 2012. Consultado el 9 de enero de 2018 .
  6. ^ Gilpin, ME (1973). "¿Las liebres comen linces?". American Naturalist . 107 (957): 727–730. doi :10.1086/282870. S2CID  84794121.
  7. ^ Jost, C.; Devulder, G.; Vucetich, JA; Peterson, R.; Arditi, R. (2005). "Los lobos de Isle Royale muestran una saciedad invariable en escala y una depredación dependiente de la densidad sobre los alces". J. Anim. Ecol . 74 (5): 809–816. doi :10.1111/j.1365-2656.2005.00977.x.
  8. ^ Pan, A.; Pourziaei, B.; Huang, H. (3 de junio de 2015). "Efecto de la fertilización con hierro del océano en la bomba biológica de carbono del fitoplancton". Avances en Matemáticas Aplicadas y Mecánica . 3 (1): 52–64. doi :10.4208/aamm.10-m1023. S2CID  124606355.
  9. ^ Prasolov, Alexander V. (2016). Algunos métodos y modelos cuantitativos en la teoría económica . Cuestiones económicas, problemas y perspectivas. Nueva York: Nova Publishers. ISBN 978-1-63484-937-1.
  10. ^ Hung, Hui-Chih; Chiu, Yu-Chih; Wu, Muh-Cherng (2017). "Un modelo Lotka-Volterra modificado para la difusión y sustitución de tecnologías de procesamiento DRAM multigeneracional". Problemas matemáticos en ingeniería . 2017 : 1–12. doi : 10.1155/2017/3038203 . ISSN  1024-123X.
  11. ^ Orbach, Yair (2022). Pronóstico de la dinámica del mercado y la tecnología . Israel: Ariel University Press. pp. 123–143. ISBN 978-965-7632-40-6.
  12. ^ Lotka, AJ (1910). "Contribución a la teoría de la reacción periódica". J. Phys. Chem. 14 (3): 271–274. doi :10.1021/j150111a004.
  13. ^ ab Goel, NS; et al. (1971). Sobre el modelo Volterra y otros modelos no lineales de poblaciones en interacción . Academic Press. ISBN 0-12-287450-1.
  14. ^ Berryman, AA (1992). "Los orígenes y la evolución de la teoría depredador-presa" (PDF) . Ecology . 73 (5): 1530–1535. doi :10.2307/1940005. JSTOR  1940005. Archivado desde el original (PDF) el 2010-05-31.
  15. ^ Verhulst, PH (1838). "Aviso sobre la ley que la población poursuit dans son accroissement". Correspondiente. Matemática y Física . 10 : 113–121.
  16. ^ Lotka, AJ (1920). "Nota analítica sobre ciertas relaciones rítmicas en sistemas orgánicos". Proc. Natl. Sci. USA 6 (7): 410–415. Bibcode :1920PNAS....6..410L. doi : 10.1073/pnas.6.7.410 . PMC 1084562 . PMID  16576509.  
  17. ^ Lotka, AJ (1925). Elementos de biología física . Williams y Wilkins .
  18. ^ Volterra, V. (1926). "Variazioni e fluttuazioni del numero d'individui in specie animali conviventi". Memoria. Acad. Lincei Roma . 2 : 31-113.
  19. ^ Volterra, V. (1931). "Variaciones y fluctuaciones del número de individuos en especies animales que viven juntas". En Chapman, RN (ed.). Ecología animal . McGraw–Hill .
  20. ^ Kingsland, S. (1995). Modelado de la naturaleza: episodios en la historia de la ecología de poblaciones . University of Chicago Press. ISBN 978-0-226-43728-6.
  21. ^ Rosenzweig, ML; MacArthur, RH (1963). "Representación gráfica y condiciones de estabilidad de las interacciones depredador-presa". American Naturalist . 97 (895): 209–223. doi :10.1086/282272. S2CID  84883526.
  22. ^ Arditi, R.; Ginzburg, LR (1989). "Coupling in predator-prey dynamics: ratio dependent" (PDF) . Journal of Theoretical Biology . 139 (3): 311–326. Bibcode :1989JThBi.139..311A. doi :10.1016/s0022-5193(89)80211-5. Archivado desde el original (PDF) el 2016-03-04 . Consultado el 2013-06-26 .
  23. ^ Abrams, PA; Ginzburg, LR (2000). "La naturaleza de la depredación: ¿dependiente de la presa, dependiente de la proporción o ninguna de las dos?". Tendencias en ecología y evolución . 15 (8): 337–341. doi :10.1016/s0169-5347(00)01908-x. PMID  10884706.
  24. ^ Gandolfo, G. (2008). "Giuseppe Palomba y las ecuaciones de Lotka-Volterra". Rendiconti Lincei . 19 (4): 347–357. doi :10.1007/s12210-008-0023-7. S2CID  140537163.
  25. ^ Goodwin, RM (1967). "Un ciclo de crecimiento" . En Feinstein, CH (ed.). Socialismo, capitalismo y crecimiento económico . Cambridge University Press .
  26. ^ Desai, M.; Ormerod, P. (1998). "Richard Goodwin: A Short Appreciation" (PDF) . The Economic Journal . 108 (450): 1431–1435. CiteSeerX 10.1.1.423.1705 . doi :10.1111/1468-0297.00350. Archivado desde el original (PDF) el 27 de septiembre de 2011 . Consultado el 22 de marzo de 2010 . 
  27. ^ Steiner, Antonio; Gander, Martín Jakob (1999). "Parametrische Lösungen der Räuber-Beute-Gleichungen im Vergleich". El Volterriano . 7 : 32–44.
  28. ^ Evans, CM; Findley, GL (1999). "Una nueva transformación para el problema de Lotka-Volterra". Journal of Mathematical Chemistry . 25 : 105–110. doi :10.1023/A:1019172114300. S2CID  36980176.
  29. ^ Leconte, M.; Masson, P.; Qi, L. (2022). "Oscilaciones de ciclo límite, tiempo de respuesta y la solución dependiente del tiempo para el modelo depredador-presa de Lotka-Volterra". Física de plasmas . 29 (2): 022302. arXiv : 2110.11557 . doi :10.1063/5.0076085. S2CID  239616189.
  30. ^ Tong, H. (1983). Modelos de umbral en análisis de series temporales no lineales . Springer–Verlag.
  31. ^ Lobry, Claude; Sari, Tewfik (2015). "Migraciones en el modelo Rosenzweig-MacArthur y el problema "atto-fox"" (PDF) . Arima . 20 : 95–125.
  32. ^ Mollison, D. (1991). "Dependencia de las velocidades de epidemia y población en parámetros básicos" (PDF) . Math. Biosci . 107 (2): 255–287. doi :10.1016/0025-5564(91)90009-8. PMID  1806118.
  33. ^ Nutku, I. (1990). "Estructura hamiltoniana de las ecuaciones de Lotka-Volterra". Physics Letters A . 145 (1): 27–28. Bibcode :1990PhLA..145...27N. doi :10.1016/0375-9601(90)90270-X. hdl : 11693/26204 . S2CID  121710034.
  34. ^ Baigent, Steve (2 de marzo de 2010). "Dinámica Lotka-Volterra: una introducción" (PDF) .

Lectura adicional

Enlaces externos