4.1. Ejemplos elementales de sistemas dinámicos: números de Fibonacci, modelo de Malthus discreto
Los números de Fibonacci
Los números de Fibonacci son una secuencia de números enteros que aparecen en muchos fenómenos de la naturaleza. Esta secuencia se define de una manera muy simple: cada número, a partir del tercero, es la suma de los dos números anteriores. Es decir, la secuencia comienza así:
Formalmente, la sucesión de Fibonacci se expresa como:
Este patrón se repite para todos los números de la secuencia. Los números de Fibonacci aparecen, por ejemplo, en la distribución de las hojas de las plantas, las espirales de los caracoles, la proporción de las ramas de los árboles y muchos otros fenómenos naturales, por lo que no es gratuito cuestionarnos: ¿cómo se obtuvo dicha secuencia?
La sucesión de Fibonacci
En el año 1202, Leonardo de Pisa (1175–1250), mejor conocido como Fibonacci, publicó su libro Liber Abaci en el cual planteó el siguiente problema: Si se pone una pareja de conejos en un lugar rodeado por un muro, ¿cuántas parejas de conejos pueden salir de esa pareja en un año si se supone que cada mes cada pareja engendra una nueva pareja que, a partir del segundo mes, se vuelve fértil?
Para resolver el problema, Fibonacci supuso que cada pareja –a partir del segundo mes– daba a luz a una nueva pareja de conejos por mes, y que cada pareja está conformada por un macho y una hembra. Para modelar esto matemáticamente, se define como
donde
Además, el número de parejas de conejos de un mes de edad en n + 1 es igual al número de parejas de dos meses o más en el mes anterior, es decir
$y_{1,n+2}=y_{2,n+1}+y_{3,n+1}+…
Ahora, se supone que la pareja del inicio es adulta y no debe esperar dos meses para poder engendrar, por lo que las condiciones iniciales son
Por lo que para
entonces
Entonces podemos deducir que
luego y_1=1. Entonces las condiciones iniciales son y_0=y_1=1.
Para entender cómo crece esta población con el tiempo, tenemos que resolver la ecuación obtenida
Lo que sigue es encontrar la solución general de esta ecuación, y para eso utilizamos algo llamado la ecuación característica. La idea es proponer una solución de la forma:
donde C es una constante (que encontraremos más tarde), y
Ahora, sustituimos esta forma en la ecuación de diferencias para encontrar
Sustituyendo
Ahora, sustituyendo estas tres expresiones en
Como C es una constante no nula, podemos cancelarla en ambos lados de la ecuación, lo que nos deja con:
Ahora, podemos dividir toda la ecuación entre
Este es el resultado de la ecuación característica, y así es como queda simplificada la ecuación
Ahora, utilizando la fórmula general para resolver ecuaciones cuadráticas tenemos que:
En este caso, los coeficientes de la ecuación son:
- a = 1 (el coeficiente de
), - b = -1 (el coeficiente de
), - c = -1 (el término independiente).
Sustituyendo en la fórmula general, tenemos:
Simplificando:
Esto nos da dos soluciones para
Estas dos soluciones son los valores de
La solución negativa
Entonces, la solución general para la población de conejos en el mes n es:
donde
Usando las condiciones iniciales de la sucesión de Fibonacci para determinar las constantes
Sabemos que
Para n = 0:
Para n = 1:
Sustituyendo
Tenemos que las raíces de la ecuación característica para la sucesión de Fibonacci son:
Ahora, para calcular la diferencia
Para simplificar, primero agrupamos los términos de forma conveniente:
Luego, como
Por lo tanto,
Ahora que tenemos las constantes
Dado que
Ejemplo 1.
Calcula de los primeros 10 términos de la sucesión de Fibonacci:
Los primeros 10 términos son:
Ejemplo 2.
Calcula cuántas parejas de conejos habrá en 6 meses, usando la ecuación
El número de parejas de conejos después de seis meses es
Ahora, usando la fórmula
Ejercicio 1.
Calcula cuántas parejas de conejos habrá en 38 meses.
Respuesta modelo:
Usando la fórmula
En 38 meses habrá 39056979 conejos.
Ejercicio 2.
Calcula cuántas parejas de conejos habrá en 83 meses.
Respuesta modelo:
Usando la fórmula
En 83 meses habrá aproximadamente 9.9022 \times 10^{21} conejos.
Modelo de Malthus discreto
El modelo de Malthus describe cómo crece una población en condiciones ideales, es decir, cuando no hay limitaciones en los recursos disponibles (como alimentos, espacio o energía). Este modelo supone que la población crece de manera exponencial, lo que significa que el número de individuos en cada periodo de tiempo aumenta en función del tamaño de la población en el periodo anterior. El modelo se expresa como:
donde
En este modelo, no se tienen en cuenta las limitaciones de recursos, lo que significa que la población puede crecer indefinidamente sin restricciones. Este tipo de crecimiento es característico de poblaciones de microorganismos en cultivo, como las bacterias, cuando se encuentran en un ambiente con recursos abundantes y sin competencia.
Limitaciones del modelo
Aunque el modelo de Malthus proporciona una descripción útil del crecimiento rápido de poblaciones en condiciones ideales, no es realista para describir el comportamiento de poblaciones en ecosistemas naturales. En la realidad, los recursos son finitos. Esto significa que a medida que la población crece, los recursos disponibles (como comida y espacio) se vuelven limitados, lo que provoca que el crecimiento de la población disminuya.
En estos casos, el modelo de Malthus deja de ser aplicable, ya que no toma en cuenta los efectos de la competencia por recursos. Por esta razón, en la naturaleza, el crecimiento de las poblaciones se describe mejor mediante modelos más complejos, como el modelo logístico, que tiene en cuenta las restricciones de los recursos y permite predecir un crecimiento poblacional que eventualmente se estabiliza en un valor determinado.
Ejemplo 1.
Supón que una población inicial de bacterias es de 50 individuos, y la tasa de crecimiento es de 1.2 por mes. ¿Cuántos individuos habrá en la población al final de 6 meses?
Tenemos que el modelo original de Malthus es:
Aquí,
Entonces, en el mes 1 tenemos que
Luego, en el mes 2, sustituyendo
En el mes 3 la población será
Y así sucesivamente. De manera que para el mes n la población se calcula como
Por lo tanto, la población después de 6 meses será:
Ejercicio 1.
Si una población de 100 individuos tiene una tasa de crecimiento
Respuesta modelo:
Para r = 1.5:
• Después de 10 meses:
• Después de 20 meses:
• Después de 30 meses:
Para r = 1.2:
• Después de 10 meses:
• Después de 20 meses:
• Después de 30 meses:
Para r = 2.0:
• Después de 10 meses:
• Después de 20 meses:
• Después de 30 meses:
Comparación:
El modelo de Malthus muestra un crecimiento exponencial, entre más grande sea la tasa de crecimiento, mayor es también el aumento de la población en el tiempo. Se puede observar cómo las tasas r = 1.5 y r = 2.0 muestran un crecimiento mucho más rápido que r = 1.2.
Ejercicio 2.
Supón que una población de bacterias comienza con 500 individuos, pero su tasa de crecimiento cambia a lo largo del tiempo. Durante los primeros 3 meses, la tasa de crecimiento es 1.1, y durante los siguientes 3 meses, la tasa es 1.3. ¿Cuál será el tamaño de la población después de 6 meses?
Respuesta modelo:
Hemos de dividir el cálculo en dos partes:
- Primeros 3 meses con tasa de crecimiento 1.1:
- Próximos 3 meses con tasa de crecimiento 1.3:
Así, después de 6 meses, la población será aproximadamente 1462 individuos.
4.2. Órbitas: iteración y tipos de órbita
Definición de un sistema dinámico discreto
Cuando estudiamos sistemas dinámicos discretos, estamos interesados en cómo cambia el sistema con el paso del tiempo. Un sistema dinámico discreto es aquel cuya evolución se describe por medio de ecuaciones matemáticas que se actualizan en intervalos de tiempo específicos. Estos sistemas se «avanzan» en pasos discretos (por ejemplo, de un día a otro, de un año a otro, de una generación a otra).
Como menciona Torres Olin “Un sistema dinámico es un organismo cambiante, o propenso al cambio. Es un conjunto de elementos que describen la evolución de un sistema. […] Se define un sistema dinámico discreto cuando medimos el tiempo en forma tal que cualesquiera dos puntos en el tiempo
Introducción a las órbitas en sistemas dinámicos discretos
En el contexto de los sistemas dinámicos discretos, las órbitas se refieren a las trayectorias que siguen las soluciones de un sistema de ecuaciones a lo largo del tiempo, cuando el tiempo avanza en pasos discretos.
• Definición. Las órbitas son secuencias de puntos que representan el estado del sistema en diferentes momentos del tiempo. Cada uno de estos puntos se encuentra en el espacio de fases del sistema, que es el conjunto de todos los posibles estados del sistema en el que se puede encontrar.
Por ejemplo, si estamos estudiando la población de una especie en un ecosistema, cada punto de la órbita representará el número de individuos de esa población en un momento dado. A medida que el tiempo avanza, la población cambia, y por lo tanto, la órbita describe cómo esa población evoluciona con el tiempo.
• Relación con las ecuaciones recursivas. En un sistema dinámico discreto, la evolución del sistema se describe mediante ecuaciones recursivas. Esto significa que el valor del estado del sistema en el tiempo t+1 depende del valor en el tiempo t.
Un ejemplo sería un modelo de población descrito por la ecuación:
donde
• Importancia de las órbitas. Las órbitas permiten entender el comportamiento a largo plazo del sistema. A medida que el sistema evoluciona con el tiempo, las órbitas nos ayudan a identificar patrones. Según el tipo de órbita, el sistema puede estabilizarse en un valor constante, entrar en un ciclo periódico o volverse impredecible (lo que se conoce como caos).
¿Qué es la iteración en sistemas dinámicos discretos?
La iteración es el proceso mediante el cual repetimos una operación matemática para obtener el siguiente valor del sistema a partir del valor actual. En otras palabras, es aplicar repetidamente una función o ecuación para predecir el comportamiento futuro del sistema. Se trata de calcular el estado del sistema en el siguiente paso de tiempo, usando el valor actual como punto de partida. Cada «paso» o «iteración» produce un nuevo valor que se usa como base para la siguiente, es decir, cada vez que realizamos una iteración, tomamos el valor del sistema en el tiempo t y lo usamos para calcular el valor en el tiempo t + 1, y este valor calculado luego se convierte en el valor inicial para el siguiente paso
Ejemplo 1.
Consideremos un modelo de crecimiento poblacional donde la población en el tiempo t + 1 depende de la población en el tiempo t:
Supongamos una población inicial
• Para t = 0,
• Para t = 1,
• Para t = 2,
Este proceso continúa para muchos pasos, y la población va cambiando a medida que se realizan las iteraciones.
Tipos de órbitas
Las órbitas en los sistemas dinámicos pueden clasificarse en diferentes tipos, dependiendo de cómo se comporta el sistema con el paso del tiempo.
1. Órbitas periódicas: son aquellas en las que el sistema regresa al mismo estado después de un número finito de pasos, es decir, la secuencia de valores se repite cíclicamente.
Ejemplo: Si r = 2 y
Estas órbitas son estables y predecibles, ya que la población sigue un patrón constante que se repite con el tiempo.
2. Órbitas atractoras: son aquellas en las que el sistema tiende a estabilizarse en un valor constante o en un ciclo periódico a medida que pasa el tiempo. Estas órbitas son estables, y las pequeñas perturbaciones tienden a ser corregidas.
Ejemplo: En el modelo logístico de crecimiento poblacional, si los valores de r y K se eligen adecuadamente, la población puede estabilizarse en el valor P = K, que es la capacidad de carga del ambiente. Este valor se convierte en una órbita atractora, ya que la población tiende a estabilizarse alrededor de K con el paso del tiempo.
3. Órbitas caóticas: son aquellas en las que el comportamiento del sistema es extremadamente sensible a las condiciones iniciales. Pequeñas variaciones en el valor inicial del sistema pueden generar trayectorias completamente diferentes. Estos sistemas son impredecibles y no siguen un patrón repetitivo.
Ejemplo: Si r es demasiado grande (por ejemplo, r = 3.5), el sistema puede volverse caótico, donde las iteraciones no siguen un patrón predecible, sino que se distribuyen de forma irregular.