Capı́tulo 6

APLICACIONES DE LOS SISTEMAS

DINÁMICOS DISCRETOS

6.1. Introducción
En este tema estudiaremos los casos más simples de crecimiento de poblaciones,
cuando la variable tiempo toma valores en un conjunto discreto, clasificados en
modelos independientes y dependientes de la densidad de la población.

DEFINICIÓN 6.1.1 Diremos que el crecimiento de una población es indepen-

diente de la densidad si las tasas de nacimiento y mortalidad no dependen del tamaño
de la población.

Recordemos que en el estudio de los modelos matriciales, ya hemos tenido ocasión de

analizar el comportamiento de ciertos modelos discretos y una breve introducción
a los modelos exponencial y logı́stico. Ahora, aplicaremos parte de los resultados
obtenidos en los temas anteriores y realizaremos un estudio más completo de algunos
de estos modelos.

6.2. Crecimiento independiente de la densidad de

la población
Comenzaremos analizando el modelo más simple de crecimiento de poblaciones de
una sola especie. Supondremos para empezar que:

La tasa de nacimientos es proporcional al número de individuos presentes.

79
80 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

La tasa de muertes es proporcional al número de individuos presentes.

Existen ciertos tipos de animales, como por ejemplo la mariposa Euphydrias editha,
que se reproduce una vez al año, poniendo sus huevos a primeros de Abril. Las
mariposas adultas vuelan durante un perı́odo corto de tiempo y entonces mueren.
Existen ratones que tienen crı́as solamente una vez al año en primavera, y que
viven alrededor de diez años. Para este tipo de especies, un modelo que suponga
que los nacimientos se dan continuamente y que las generaciones se superponen es
inapropiado.

Figura 6.1: Modelo discreto exponencial.

Mediremos el tiempo k en unidades de generación (un año, un mes, ...), y supondre-

mos que r es el número de individuos que nacen en la próxima generación a partir
de un individuo de la generación actual. Si xk simboliza al número de individuos de
la población en la generación k, entonces

xk+1 = r xk , k = 0, 1, 2, · · · .

Si x0 es el número inicial de individuos, de la expresión anterior se deduce

xk = x0 rk , k = 0, 1, 2, · · · , (6.1)

es decir, estamos ante un crecimiento exponencial o geométrico. El comportamiento

cualitativo de (6.1) está determinado por el valor de r y queda simbolizado en la
Figura 6.1.
Es evidente que este modelo representa a la población sólo en un intervalo corto de
tiempo, ya que el crecimiento es demasiado rápido. Además, este modelo basado en
la independencia de la densidad, no puede explicar la evolución de la mayorı́a de las
poblaciones que existen en la naturaleza.
Podemos preguntarnos por los valores reales, y no los teóricos, que se obtienen
6.2 Crecimiento independiente de la densidad de la población 81

del parámetro r en el laboratorio y en la naturaleza. En los experimentos en el

laboratorio puede encontrarse valores de r muy diferentes, dando lugar a crecimiento
muy rápido de poblaciones. Sin embargo, en la naturaleza este valor debe estar muy
cerca de uno, ya que en caso contrario la población desaparecerı́a o por el contrario
crecerı́a rápidamente.

Figura 6.2: Crecimiento de una población de pájaros.

La Figura 6.2 muestra la representación en escala logarı́tmica de una población de
pájaros de Gran Bretaña, desde el año 1955 al 1970. Observemos que al principio,
la población crece exponencialmente, pero después de algunos años, disminuye sus-
tancialmente. En la próxima sección trataremos de explicar este comportamiento.
La cuestión más importante de la dinámica de poblaciones es determinar las causas
y las consecuencias de la desviación del modelo exponencial.

EJEMPLO 6.1
El censo de los Estados Unidos se elabora cada diez años. En la Tabla 11.1. se
recogen los datos correspondientes al perı́odo 1790 - 2000.
La tasa de crecimiento en cada década se calcula dividiendo el censo correspondiente
al año superior entre el número de individuos en el año inferior. Por ejemplo, la tasa
de crecimiento en la década 1790 - 1800 es:
Población en 1800 5.308.483
= = 1.351 .
Población en 1790 3.929.214
El modelo matemático discreto más simple supone que la población en la próxima
década es igual a la población actual más la población actual por la tasa de creci-
miento medio, r, de la población. El modelo empieza con una población inicial, por
ejemplo, la correspondiente al año 1790. Para encontrar la población en la década
próxima, multiplicamos por (1+r). Con ello obtenemos una sucesión de poblaciones,
todas ellas encontradas a partir de la década anterior. Por ejemplo,
Población en 1800 = 1.349 × Población en 1790 = 5300510 ,
82 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

siendo 34.9 % la media de las tasas de crecimiento desde 1790 hasta 1860. Observemos
que existe una diferencia de aproximadamente 8000 individuos que equivale a un
error del 0.15 %. Podemos repetir el proceso anterior y encontrar las poblaciones
para las décadas 1810, 1820, ... , 1860, ya que en estos perı́odos la tasa de crecimiento
se mantiene razonablemente constante.

Tabla 6.1

La Tabla 6.2 muestra los datos obtenidos. En ella puede observarse que los errores
cometidos son pequeños hasta 1870, y además la población predicha por el modelo
es ligeramente superior a la población exacta, lo cual nos sugiere que durante el siglo
XIX bajó la tasa de nacimiento. Entre los años 1860 y 1870 tuvo lugar la guerra civil
americana, originando el brusco descenso en la tasa de crecimiento de la población
de Estados Unidos; además durante estos años aconteció la revolución industrial y
la sociedad pasó de ser mayoritariamente agrı́cola a una sociedad industrial con un
descenso significativo de los nacimientos.
Si continuamos usando el modelo anterior hasta 1920 o 1970 nos encontraremos con
una población predicha de 192365343 y 859382645 respectivamente, lo que supone
una estimación del 82 % y 323 % mayores que las reales. La conclusión que deducimos
es que el uso de este modelo de crecimiento está limitado a predecir la población
futura en años muy próximos, no se puede extrapolar a largo plazo.
Recordemos que el modelo matemático dado por

xk+1 = xk + rxk = (1 + r)xk , x0 = P (1790) = 3.929.214 , (6.2)

siendo r la tasa media de crecimiento, se conoce con el nombre de modelo de cre-

cimiento discreto exponencial o de Malthus. El modelo es un caso particular
de un sistema dinámico discreto o ecuación en diferencias. Las ecuaciones en dife-
rencias se usan con frecuencia en Ecologı́a, donde a menudo se puede determinar la
población de una especie o colección de especies, sabiendo la población en la gene-
ración anterior. El modelo de crecimiento malthusiano establece que la población
en la próxima generación es proporcional a la población de la generación actual. De
(6.2) se deduce inmediatamente

xk = (1 + r)k x0 , k = 1, 2, 3 · · · .
6.2 Crecimiento independiente de la densidad de la población 83

Tabla 6.2

A continuación modificaremos el modelo anterior para obligar a que la tasa de cre-

cimiento sea una función que dependa del tiempo. Hemos comprobado que la tasa
media de crecimiento que calculamos para las primeras décadas predice una pobla-
ción muy superior a la ofrecida por el censo. Para mejorar esta predicción, podemos
calcular para cada una de las décadas su tasa de crecimiento r y encontrar la recta
de regresión de todos estos datos.
Se pasa ası́ del modelo discreto autónomo xk+1 = f (xk ), al modelo discreto no
autónomo xk+1 = f (xk , tk ). La recta de regresión r(k) = 3.158 − 0.00155k ajusta a
la nube de puntos de las diferentes tasas de crecimiento. En este caso, la ecuación
en diferencia no autónoma será:

xk+1 = (1 + r(k))xk , (6.3)

siendo tk = 1790 + 10k, y k el número de décadas después de 1790.

Figura 6.3: Tasa de crecimiento para la población de EEUU.

La Figura 6.3 permite comparar los datos del censo con las diferentes proyecciones
que se obtienen al utilizar el modelo de crecimiento exponencial autónomo y no
autónomo (que no dependen/dependen del tiempo). Llamamos la atención sobre
84 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

el hecho de que si utilizamos (6.3) para encontrar la población en cada década, es

imprescindible conocer la población en la década anterior.

Figura 6.4: Modelos de crecimiento exponencial.

En la Tabla 6.3 se comparan numéricamente los datos reales con los obtenidos con
(6.3). El modelo (6.3) predice 278244477 individuos para el año 2000, cifra que se
encuentra ligeramente por debajo del valor real.

Tabla 6.3
6.3 Crecimiento dependiente de la densidad de población 85

6.2.1. Modelo discreto exponencial modificado

Hemos aplicado el modelo de crecimiento discreto exponencial para estudiar la evo-
lución de una población. Durante su aplicación, se ha considerado el sistema como
cerrado para poder trabajar con una tasa neta de crecimiento. Pero podemos modifi-
car dicho modelo para tener en cuenta el hecho de la inmigración y de la emigración.
Supongamos que una población xk crece de acuerdo al modelo discreto exponencial
y asumimos que el número de personas que entran y salen en cada intervalo de tiem-
po es constante (e − s = µ). Ahora, el crecimiento puede modelarse por la ecuación
en diferencias:
xk+1 = (1 + r)xk − µ , k = 0, 1, 2, · · · ,
donde r es la tasa de crecimiento. Conocidos estos datos y la población inicial x0
podemos encontrar una expresión general de xk . En efecto,

x1 = (1 + r)x0 − µ

x2 = (1 + r)x1 − µ = (1 + r) ((1 + r)x0 − µ) − µ =

(1 + r)2 x0 − ((1 + r) + 1) µ

x3 = (1 + r)3 x0 − ((1 + r)2 + (1 + r) + 1) µ

.. .. ..
. . .
( )
xk = (1 + r)k x0 − (1 + r)k−1 + (1 + r)k−2 + · · · + (1 + r) + 1 µ

Aplicando la fórmula que nos da la suma de un número finito de términos de una

progresión geométrica, se obtiene

(1 + r)k − 1
xk = (1 + r)k x0 − µ,
r

expresión más complicada que la correspondiente al modelo discreto exponencial

simple. Aunque en este caso concreto hemos podido encontrar una expresión para
xk en función de x0 , r y µ, tenemos que decir que en general este cálculo suele ser
complicado. Por esta razón, lo que se hace es estudiar el comportamiento cualitativo
del modelo, por ejemplo, a través de su diagrama de Cobweb.

6.3. Crecimiento dependiente de la densidad de

población
Ya hemos indicado que el análisis del modelo discreto exponencial y el sentido común,
nos dicen que este tipo de crecimiento no puede mantenerse durante mucho tiempo.
86 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

En todos los casos, llega un momento en que la población se regula. Se han pro-
puesto muchas hipótesis para explicar las causas que originan este autocontrol de la
población, entre otras:

Factores independientes de la densidad, como por ejemplo el clima.

La cantidad de comida disponible.

Problemas con su territorio o canibalismo.

Depredadores.

Parásitos o enfermedades.

De entre todos estos factores nosotros estudiaremos el segundo de ellos, es decir el

crecimiento dependerá de la densidad de la población, y por tanto, ésta se autoregula.
Un modelo clásico apropiado para describir poblaciones de animales (o plantas) que
viven un año, se reproducen y luego mueren, es de la forma:

xk+1 = f (xk ) , k = 0, 1, 2, · · · , (6.4)

donde f nos da el número de individuos para el próximo año en términos del número
de individuos actuales. Se han propuesto diferentes modelos, simplemente cambiando
la función f . Por ejemplo, en el estudio del caos se trabaja con el modelo de May
(1974) donde la función f es,

f (x) = cx (1 − x) .

6.3.1. El modelo de crecimiento discreto logı́stico

En 1913 T. Carlson estudió el crecimiento de un cultivo de levadura. La Tabla 6.4
muestra los datos recogidos en intervalos de una hora.

Tabla 6.4: Población de un cultivo de levadura

En ella se observa que la población no sigue un modelo de crecimiento discreto

exponencial, ya que a partir de cierto momento la población se estabiliza y no crece
exponencialmente. Es necesario que la función f (x), del sistema discreto dinámico
general xk+1 = f (xk ), ahora sea cuadrática en lugar de ser una ecuación lineal.
6.3 Crecimiento dependiente de la densidad de población 87

Este nuevo modelo se conoce con el nombre de modelo discreto logı́stico, y viene
expresado por ( xk )
xk+1 = xk + rxk 1 − , k = 0, 1, 2, · · · . (6.5)
M
Observemos que para valores pequeños de la población 1 − xMk ≈ 1 y el modelo
coincide con el exponencial. Sin embargo, para valores de la población xk ≈ M
entonces xk+1 ≈ xk . El parámetro M recibe el nombre de capacidad de carga de
la población.

Figura 6.5: Modelo para un cultivo de levadura.

El comportamiento de (6.5) es bastante más complicado que (6.2). No existe una

solución exacta de este sistema dinámico discreto. El ecologı́sta Robert May (1974)
estudió dicha ecuación para diferentes poblaciones y descubrió que podı́a presentar
dinámicas muy diferentes. Este hecho lo pusimos de manifiesto al analizar el caos
matemático, ya que (6.5) puede ser escrita como xk+1 = µxk (1 − xk ).
A continuación aplicaremos este modelo para estudiar la evolución del cultivo de
levadura.

EJEMPLO 6.2

En la Figura 6.5 hemos dibujado xk+1 como función de xk . Por ejemplo, los dos
primeros puntos son (9.6 , 18.3) y (18.3 , 29). Posteriormente utilizando el programa
Mathematicar se ha encontrado la parábola que pasa por el origen y = ax − bx2
que mejor ajusta a estos datos, obteniéndose

xk+1 = 1.5612xk − 0.000861x2k .

Podemos utilizar un programa de simulación, como por ejemplo POPULUSr , y ob-

tendrı́amos la Figura 6.5 De forma cualitativa podemos ver que inicialmente se
produce un crecimiento exponencial y que posteriormente la población se estabiliza
alrededor de 650 que es la capacidad de carga del modelo.
88 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

Figura 6.6: Simulación del modelo.

Observemos también que el punto de inflexión está situado en la mitad de la capaci-

dad de carga, que corresponde a un tiempo entre las 9 y 10 horas. En este momento
se produce el máximo crecimiento de la población.

6.3.2. Generalización del modelo discreto logı́stico

La mayorı́a de otros modelos comparten los rasgos cualitativos observados en el
modelo de May. Si representamos en el eje de abscisas la población en el tiempo k,
y en el eje de ordenadas la población en el perı́odo siguiente xk+1 , en gran parte de
ellos se obtiene una curva del tipo representado en la Figura 6.7.

Figura 6.7: Representación de los puntos (xk , xk+1 )

Observemos que esta curva tiene un único máximo. Cuando el nivel de la población es
pequeño, entonces aumenta en función de la población actual, pero cuando el número
de individuos es elevado, los mecanismos propios relacionados con la densidad de la
población (competición, por ejemplo) reducen su nivel en los próximos años.
6.3 Crecimiento dependiente de la densidad de población 89

De entre los modelos más citados en el estudio de dinámica de poblaciones, se

encuentran: ( x )
f (x) = x 1 + x(1 − ) ,
k
x
f (x) = x er(1− k ) ,

λx
f (x) =
(1 + αx)β
En una de las prácticas del Laboratorio Matemático, realizamos un estudio intensivo
del segundo de los modelos, conocido con el nombre de modelo de Ricker (1954).
Para los otros dos casos, se puede hacer un tratamiento similar.

EJEMPLO 6.3

Un modelo matemático dependiente de la densidad de la población y alternativo al

modelo logı́stico de May, ha sido propuesto por Gilpin y Ayala (1973), y se expresa
como: ( ( )α )
xk
xk+1 = f (xk ) = r xk 1 − , k = 0 , 1, 2, · · · (6.6)
β
donde α es un parámetro positivo que depende del organismo en cuestión.
El punto de equilibrio no nulo de este modelo se obtiene resolviendo la ecuación
( ( )α )
x
f (x) = x ⇒ r x 1 − =x
β

cuyo valor es
( )1
∗ r−1 α
x =β .
r

Para estudiar la estabilidad del modelo primero debemos derivar la función f (x).
Una vez simplificada se obtiene
( ( )α ( )α )
′ x x
f (x) = r 1 − −α .
β β

Luego
( ( )1)
r − 1 α
f ′ (x∗ ) = f ′ β = 1 − αr + α .
r

Este punto de equilibrio será estable cuando |f ′ (x∗ )| < 1, lo cual ocurre cuando
1 < r < 1 + α2 .
90 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

En ciertas ocasiones, como por ejemplo en el modelo logı́stico de May

f (x) = c x (1 − x/M ),

si el nivel de la población es demasiado bajo, entonces el número de individuos tiende

a largo plazo al punto de equilibrio x∗ = 0 y la población desaparece. Este fenómeno
es conocido en ecologı́a con el nombre de Efecto Allen. Muchas poblaciones bio-
lógicas que presentan este efecto, decrecen en su tamaño si el número de individuos
se encuentran por debajo de cierto nivel crı́tico xc . La región donde xk < xc es
conocida con el nombre de zona de depredación.
Podemos modificar el modelo anterior, para tener en cuenta este hecho, de la manera
siguiente: ( x)
f (x) = c x 1 − (x − a) , a > 0 .
M

6.4. Ejemplo de modelo discreto para la pesca

En los últimos años los modelos discretos han sido muy utilizados en el diseño
de estrategias para la pesca. Se ha demostrado que son muy útiles para evaluar
diversas tácticas de capturas de peces con un doble objetivo, en primer lugar para
maximizar los beneficios y en segundo lugar para realizar una explotación de recursos
mantenidos en el tiempo. El modelo que vamos a estudiar también puede ser aplicado
a cualquier otro tipo de recurso renovable.
Supongamos que la densidad de la población en ausencia de capturas viene dada
por
xk+1 = f (xk ) , k = 0, 1, 2, · · · .
Si suponemos que ϵ(k) es la captura realizada en la población en el tiempo k, la cual
es la que genera la población en el tiempo k + 1, entonces el modelo que estudia la
dinámica de la población viene dado por:

xk+1 = f (xk ) − ϵ(k) , k = 0, 1, 2, · · · . (6.7)

Las dos preguntas que debemos contestar son:

¿Cuál es el máximo rendimiento biológico sostenible YM ?
¿Cuál es el máximo rendimiento económico EM ?
Si encontramos los puntos de equilibrio de (6.7), deducimos que

x∗ = f (x∗ ) − ϵ∗ ⇒ ϵ∗ = f (x∗ ) − x∗ .

Si el máximo rendimiento sostenible del punto de equilibrio YM se alcanza cuando

x∗ toma el valor xM , entonces su valor podemos encontrarlo haciendo
∂ϵ∗
=0 ⇒ f ′ (x∗ ) = 1 .
∂x∗
6.4 Ejemplo de modelo discreto para la pesca 91

El valor de YM será
YM = f (xM ) − xM (6.8)
y esta situación sólo es interesante cuando YM ≥ 0.
Una estrategia podrı́a ser mantener la población de peces en estos niveles con el
objetivo de hacer máxima la captura YM . Pero como es difı́cil tener un conocimiento
exacto de la población actual de peces, entonces este método puede ser difı́cil llevarlo
a la práctica. Por esta razón, es más interesante formular el problema de optimización
en términos de capturas y esfuerzos.
Supongamos que el esfuerzo para capturar un pez, de una población x, es ax, donde
a es el parámetro de captura (que es independiente de la densidad x). Entonces el
esfuerzo para reducir x en 1 unidad es 1/(ax) y f (x) en 1 unidad es 1/(af (x)). De
esta manera, el esfuerzo EM para obtener la captura YM = f (xM ) − xM es

∑
f (xM )
EM = (axi )−1 .
xi =xM

Frecuentemente los valores de este sumatorio son de tal manera que se pueden apro-
ximar por la siguiente integral
∫ ( )
1 f (xM ) 1 1 f (xM )
EM ≈ dx = ln . (6.9)
a xM x a xM

Las ecuaciones (6.8) y (6.9) nos dan la relación de YM , EM en función de x.

EJEMPLO 6.4

Para terminar, aplicamos estos resultados a un modelo concreto, conocido como

disco de Holling, que viene definido por:

βxk
xk+1 = , 0<α<β.
α + xk

En primer lugar encontramos el valor de xM resolviendo 1 = f ′ (xM ). Es decir,

( )′
βxM αβ √ (√ √ )
1= = ⇒ x M = α β − α .
α + xM (α + xM )2
Si sustituimos en las ecuaciones (6.8) y (6.9), nos da
βxM
YM = − xM
α + xM
( )
1 β
EM = ln .
a α + xM
92 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

En este ejemplo, podemos eliminar entre las dos expresiones xM y obtener una
relación explı́cita entre YM y EM ,
( )( )
YM = βe−cEM − α ecEM − 1 .

6.5. Ejemplo de modelo discreto para la economı́a.

Modelo de la telaraña.
Es un modelo elemental que simula el comportamiento de un bien en el mercado,
sujeto a las variaciones de la oferta y de la demanda.
La empresa que ofrece el bien, cambiará su oferta , Ot , durante el perı́odo t propor-
cionalmente a la variación del precio del bien en el perı́odo anterior t − 1. Pensemos
que si el precio del bien en el perı́odo anterior ha aumentado, entonces incremen-
tará la oferta en el perı́odo siguiente t. Por el contrario, si el precio disminuye en el
perı́odo t − 1, entonces la empresa ofertará menos cantidad en el perı́odo siguiente t,
intentando contrarestar la tendencia a la baja y evitando la disminución de ingresos.
Por lo tanto, si el precio del bien sube (baja) en el perı́odo t − 1, entonces la oferta
del mismo bien en el perı́odo siguiente t sube (baja) proporcionalmente,
Ot = aPt−1 ; a > 0; t ∈ IN
Además, debemos añadir a la ecuación anterior un sumando constante b que se
interpreta como la acción de las fuerzas independientes de la variación del precio del
bien, que también tiene influencia en la variación de la oferta del bien en el mercado.
Ot = aPt−1 + b ; a > 0; t ∈ IN (6.10)
Vamos a suponer que el incremento de la demanda en el perı́odo t de dicho bien,
Dt , varı́a proporcionalmente al aumento de su precio en dicho perı́odo. Es decir, el
consumo del bien crecerá si su precio disminuye y recı́procamente. En consecuencia,
la constante de proporcionalidad c debe ser negativa. Además, y al igual que con la
oferta, incluiremos un término constante d.
Dt = cPt + d ; c < 0; t ∈ IN (6.11)
Por último, supondremos que la dinámica del mercado hace que la oferta y la de-
manda del bien tiendan a coincidir en cada perı́odo t,
O t = Dt ; ∀t ∈ IN ⇒ aPt−1 + b = cPt + d ⇒ cPt − aPt−1 + (d − b) = 0
que puede expresarse como el sistema dinámico discreto:
a b−d
Pt = f (Pt−1 ) = Pt−1 + ; t ∈ IN (6.12)
c c
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 93

EJEMPLO 6.5

Dadas las siguientes funciones Dt = 100 − 2Pt y Ot = −20 + 3Pt−1 , hallar el valor
de equilibrio del precio y comprobar si es estable o inestable. Suponer que el valor
inicial es P0 = 25, y calcular los valores numéricos de Pt hasta t = 4.

Incluir solución

EJERCICIO 30 Sea yt el número de individuos de una determinada

especie de animales en el tiempo t. Sabiendo que su evolución sigue una
relación de la forma:
3 1
yt+2 = yt+1 − yt , t = 0, 1, 2, · · · ,
2 2
probar que la población se estabiliza a largo plazo.

La ecuación en diferencias anterior es homogénea ya que puede ser escrita como,

2yt+2 − 3yt+1 + yt = 0 ,

y tiene como ecuación caracterı́stica

1
2λ2 − 3λ + 1 = 0 ⇒ λ1 = 1, λ2 = ,
2
luego la solución general

( 1 )t k2
y t = k1 + k2 2 = k1 + , k1 , k2 ∈ IR
2t

Si tomamos lı́mites cuando t tiende a infinito se obtiene de manera inmediata que

yt → k1 .

EJERCICIO 31 Resolver la ecuación en diferencias de orden dos

yt+2 + yt = 1 + t .

La solución general se construye a partir de una solución particular de la ecuación

completa y la solución general de la ecuación homogénea asociada. Empezamos, por
tanto, encontrando las raı́ces del polinomio caracterı́stico

λ2 + 1 = 0 ⇒ λ = ±i .
94 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

Es decir, dos números complejos conjugados de módulo 1 y argumento π/2. La

solución general de la ecuación homogénea es:
(π ) (π )
yth = k1 cos 2 t + k2 sen 2 t , k1 , k2 ∈ IR

Para buscar una solución particular de la ecuación completa, observamos que el

término independiente 1 + t, es un polinomio de primer grado. Ensayamos con la
solución ytp = a + bt. Al imponer que sea solución de la ecuación en diferencias, se
obtiene,

a + b(t + 2) + a + bt = 1 + t ⇒ 2a + 2b = 1, 2b = 1 ⇒ a = 0, b = 1/2 ,

luego, la solución particular buscada es ytp = 1/2 t. La solución general de la ecuación

completa será
(π ) (π ) 1
yt = k1 cos 2 t + k2 sen 2 t + 2 t , k1 , k2 ∈ R

EJERCICIO 32 Resolver la siguiente ecuación en diferencias lineal de

coeficientes constantes

yt+2 − 6yt+1 + 5yt = 3t .

Para encontrar la solución yt de la ecuación completa empezamos buscando yth , que

es la solución general de la homogénea

yt+2 − 6yt+1 + 5yt = 0.

Al ser las raı́ces del polinomio caracterı́stico λ1 = 1 y λ2 = 5,

yth = k1 + k2 5t .

Para conocer una solución particular de la ecuación completa nos fijamos en el

término independiente 3t y ensayamos la solución ytp = c3t . Sustituyendo en la
ecuación inicial y simplificando

1
c3t+2 − 6c3t+1 + 5c3t = 3t ⇒ c=− .
4

La solución general vendrá dada por

yt = k1 + k2 5t − 14 3t , k1 , k2 ∈ IR
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 95

EJERCICIO 33 [Modelo de la telaraña.] Para el ajuste dinámico de

un bien en el mercado, se suele utilizar un modelo discreto, que se
fundamenta en las siguientes hipótesis:

La oferta del bien depende del precio del perı́odo anterior (la pro-
ducción del producto se decide teniendo en cuenta el precio en
ese momento, pero tarda en realizarse un perı́odo de tiempo, por
ejemplo en los productos agrı́colas),

St = −c + dPt−1 , c, d > 0

La demanda en cada perı́odo depende del precio del bien en el

mismo perı́odo de tiempo

Dt = a − bPt , a, b > 0

La condición de equilibrio será Dt = St , siendo el valor del precio

inicial P0 .

Analizar el comportamiento del modelo.

De la tercera de las hipótesis que nos da la condición de equilibrio, obtenemos la

siguiente ecuación en diferencias

d a+c d a+c
Pt = − Pt−1 + ⇒ Pt+1 = − Pt + ,
b b b b

que para resolverla, damos al tiempo los valores t = 0, 1, 2, · · ·

d a+c
P1 = − P0 +
b b
( )
d a+c d d a+c a+c
P2 = − P1 + =− − P0 + +
b b b b b b
( ) ( )
d 2 d a+c a+c
= − P0 − +
b b b b
(( )2 ( ) )
d a+c d d d a+c a+c a+c
P3 = − P2 + =− − P0 − + +
b b b b b b b b

( ) ( ) ( ) ( )( )
d 3 d 2 a+c d a+c a+c
= − P0 + − + − +
b b b b b b
.. ..
. .
96 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

.. ..
. .
( )t ( )t−1 ( ) ( )( )
d d a+c d a+c a+c
Pt = − P0 + − + ··· + − +
b b b b b b
( )t ( ) [( ) ( ) ( ) ]
d a+c d t−1 d t−2 d
= − P0 + − + − + ··· + − +1
b b b b b

Los sumandos que se encuentran dentro del corchete son la suma1 de t términos de
una progresión geométrica de razón − db , cuyo valor es

1 − (−d/b)t
,
1 − (−d/b)

si sustituimos este valor en Pt y simplificamos convenientemente,

( ) [ ( )] ( ) [ ]
d t a+c d t d t a+c a+c
Pt = − P0 + 1− − = − P0 − + .
b b+d b b b+d b+d

Llamando

a+c
Pe := ,
b+d

que se conoce con el nombre de precio teórico de equilibrio para las funciones de
oferta y demanda dada. Sustituyendo

( )
d t
Pt = (P0 − Pe ) − + Pe .
b

Si suponemos que d < b, entonces d/b < 1 y Pt tiende a Pe cuando t tiende a infinito.
Es decir, las fuerzas del mercado harán que el precio del producto tienda al precio
de equilibrio teórico. Esta situación queda reflejada en la figura siguiente:

1
La suma de n términos de una progresión geométrica de razón r vale

1 − rn
Sn = a 1
1−r
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 97

Figura 6.8: Modelo de telaraña.

Para un precio inicial P0 , los productores ofrecen en el perı́odo la oferta S1 , pero la

demanda cubrirá dicha oferta a un precio diferente P1 . A este último precio la oferta
del perı́odo siguiente será S2 , que será cubierta por la demanda a un precio P2 , y
ası́ sucesivamente.

EJERCICIO 34 Encontrar la raı́z de la ecuación x3 − x − 1 = 0 utilizando

el método del punto fijo.

Recordemos que para aplicar el método debemos escribir la ecuación en la forma

x = g(x) y obtener una sucesión xk+1 = g(xk ) partiendo de un valor x0 ∈ (a, b).
Observemos que si aplicamos el Teorema de Bolzano a la función φ(x) = x3 − x − 1
en el intervalo [1, 2] nos aseguramos que la ecuación x3 − x − 1 = 0 tiene una raı́z
en el intervalo (1, 2),
Por otro lado, si escribimos la ecuación como x3 − 1 = x y consideramos la función
g(x) = x3 − x, podemos tomar como valor inicial o semilla un número entre 1 y 2,
por ejemplo x0 = 1.5. Al ser g ′ (x) = 3x2 , en cualquier entorno de 1.5 se cumple
|g ′ (x)| > 1. En consecuencia, el método no es convergente.
1
También es posible escribir la ecuación de esta otra manera x = (x + 1) 3 , y ahora
1
considerar otra función h(x) = (x + 1) 3 con derivada h′ (x) ≈ 0.165 en un entorno
de x0 = 1.85. Es decir, el método del punto fijo es convergente.
Para encontrar la raı́z utilizamos el software Mathemathicar

h[x ] := (x + 1)1/3
FixedPointList[h, 1.85, 14]
98 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

{1.85, 1.417799, 1.342167, 1.328024, 1.325345, 1.3248371, 1.324740,

1.324722, 1.324718, 1.324718, 1.324717, 1.324717, 1.324717, 1.324717}

Figura 6.9: Diagrama de Cobweb.

EJERCICIO 35 El modelo formal a tiempo discreto, que describe la

convivencia de dos especies con funciones de efectivos xt e yt , con medi-
das mensuales, es el siguiente:
{
xt+1 = xt + yt − 12
2t
yt+1 = yt + 28t .

Al principio x0 = 40, y0 = 31. ¿En qué situación está el ecosistema al cabo

de 4 meses?

Al no depender la segunda de las ecuaciones de xt , empezamos resolviéndola. Su

ecuación homogénea asociada es λ − 1 = 0, que tiene por raı́z λ = 1, dando lugar a
la siguiente solución general de la ecuación homogénea yt = k1 .
Para encontrar la solución general de la ecuación completa, ensayamos la solución
particular yt = a/2t . Sustituyendo
a a 8
= + t ⇒ a = −16 ,
2t+1 2 t 2
y la solución general de la ecuación completa es:

16
y t = k1 −
2t

Si sustituimos este valor en la primera de las ecuaciones del sistema

16 12 28
xt+1 = xt + k1 − − t ⇒ xt+1 − xt = − + k1 .
2t 2 2t
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 99

La ecuación homogénea tiene a λ = 1 como raı́z de su ecuación caracterı́stica asocia-

da, por tanto, y = k2 será su solución general. Para encontrar una solución particular
de la solución completa nos fijamos en el término independiente −28/2t + k1 , que
como podemos ver está formado por dos términos −28/2t y la constante k1 (un
polinomio de grado cero). Al ser λ = 1 raı́z de la ecuación caracterı́stica, debemos
tener en cuenta la observación realizada en la teorı́a, y tenemos que ensayar con un
polinomio de un orden mayor. En resumen, debemos probar con xt = a/2t + bt + c.
Sustituyendo en la ecuación, se obtiene
a (a ) 28
+ b(t + 1) + c − + bt + c = − t + k1 .
2t+1 2t 2
Simplificando
( )
a 1 28
− + b = − t + k1 ⇒ a = 56, b = k1 .
2t 2 2
En consecuencia, la solución general de la ecuación completa es

56
xt = k2 + k1 t + , k1 , k2 ∈ IR
2t

Si estamos interesados en encontrar la solución particular para los valores x0 = 40 e

y0 = 31, debemos sustituir en las soluciones generales encontradas
40 = k2 + 56 ⇒ k2 = −16
31 = k1 − 16 ⇒ k1 = 47
La solución particular que cumple las condiciones iniciales es:

7
xt = 47t − 16 + t−3
2
1
yt = 47 − t−4
2

Un método alternativo para resolver el ejercicio es el siguiente.

Dando los valores t = 1, 2, 3, · · · , se obtiene
y1 = 31 + 8 210
( ) ( )
y2 = 31 + 8 210 + 8 211 = 31 + 8 210 + 1
21
(1 1
) (1 )
y3 = 31 + 8 20
+ 21
+ 8 212 = 31 + 8 20
+ 1
21
+ 1
22

..
. ( )
yt = 31 + 8 210 + 1
21
+ 1
22
+ ··· + 1
2t−1
.
Ahora, utilizando la fórmula que nos da la suma de t términos de una progresión
geométrica de razón 1/2,
( )t
1 − 12
yt = 31 + 8 ,
1 − 12
100 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

y simplificando
( )t
1
yt = 47 − 16
2

Para encontrar xt , sustituimos el valor de yt en la primera de las ecuaciones

( )t ( )t
1 1
xt+1 = xt + 47 − 16 = xt + 47 − 28 ,
2 2
y resulta ser del mismo tipo a la anterior,
( )0
x1 = 40 + 47 − 28 12
( ( )0 ) ( )1 [( ) ( )1 ]
0
x2 = 40 + 47 − 28 12 + 47 − 28 12 = 40 + 2 × 47 − 28 21 + 12
[( ) ( 1 )1 ] ( 1 )2
1 0
x3 = 40 + 2 × 47 − 28 2 + 2 + 47 − 28 2

[( ) ( 1 )1 ( 1 )2 ]
1 0
= 40 + 3 × 47 − 28 2 + 2 + 2
..
. [( )
0 ( )1 ( )t−1 ]
xt = 40 + t × 47 − 28 21 + 12 + · · · + 12 .

Sumando los t términos de esta progresión geométrica

( )t
1
1−
2
xt = 40 + 47t − 28 ,
1
1−
2
que una vez simplificada
( )t
1
xt = −16 + 47t + 56
2

Para saber la situación de las poblaciones al cabo de los 4 años, sustituimos en las
soluciones encontradas t = 4,
x4 ≈ 176, y4 = 46 .

EJERCICIO 36 Dos especies que conviven en un mismo territorio si-

guen un crecimiento descrito por el sistema de ecuaciones en diferencias
siguiente: {
xt+1 = 5xt − 2yt + t
yt+1 = 4xt − yt + 3
donde el tiempo t está medido en años. Si inicialmente el número de
individuos de cada especie es x0 = 130 e y0 = 250, resolver el sistema y
analizar el comportamiento a la larga de las dos especies.
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 101

Comenzamos el ejercicio convirtiendo el sistema en una ecuación en diferencias li-

neal de segundo orden con coeficientes constantes. Para ello, de la segunda de las
ecuaciones deducimos
yt+2 = 4xt+1 − yt+1 + 3 .
Ahora, sustituimos xt+1 de la primera de las ecuaciones en la expresión anterior

yt+2 = 4 (5xt − 2yt + t) − yt+1 + 3 = 20xt − 8yt + 4t − yt+1 + 3 .

Por último, sustituimos el valor xt de la segunda de las ecuaciones del sistema, y

simplificamos
yt+2 − 4yt+1 + 3yt = 4t − 12 .

Para resolverla, empezamos encontrando la solución general de su ecuación ho-

mogénea asociada

λ2 − 4λ + 3 = 0 ⇒ λ1 = 1, λ2 = 3 ;

la solución buscada es:

yth = k1 + k2 3t .

A continuación necesitamos una solución particular de la ecuación completa. Al ser

el término independiente un polinomio de primer grado y λ = 1 raı́z del polinomio
caracterı́stico, ensayamos la solución ytp = at2 + bt + c. Si sustituimos en la ecuación
y simplificamos

(−4a)t − 2b = 4t − 12 ⇒ a = −1, b=6 ⇒ ytp = −t2 + 6t .

La solución general de la ecuación completa es:

yt = k1 + k2 3t − t2 + 6t , k1 , k2 ∈ IR

Para encontrar xt , despejamos de la segunda de las ecuaciones y sustituimos el valor

de yt
xt = 1
4 (yt+1 + yt − 3)
( )
= 1
4 k1 + k2 3t+1 − (t + 1)2 + 6(t + 1) + k1 + k2 3t − t2 + 6t − 3
( )
= 1
4 2k1 + 4k2 3t − t2 − 10t + 2 .
Es decir,
1 1 5 1
xt = k1 + k2 3t − t2 − t + , k1 , k2 ∈ IR
2 4 2 2

Finalizamos el ejercicio encontrando la solución particular del sistema correspon-

diente a las condiciones iniciales, x0 = 130 e y0 = 250. Sustituyendo en las expresio-
nes de xt e yt obtenemos el sistema

 130 = 12 k1 + k2 + 12

250 = k1 + k2 ;
102 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

cuya solución nos proporciona los valores k1 = 241 y k2 = 9. En consecuencia,

241 1 5
xt = + 9 × 3 t − t2 − t
2 4 2

yt = 241 + 9 × 3t − t2 + 6t

Como podemos apreciar, si en las expresiones anteriores tomamos lı́mites cuando t

tiende a infinito, nos encontramos con xt → ∞ e yt → ∞.

EJERCICIO 37 Dos especies admiten el siguiente modelo de coexisten-

cia: {
xt+1 = 4xt + 6yt − 3t
yt+1 = −2xt − 4yt + 3t .
Obtener las expresiones de las funciones de efectivos de las dos especies
que satisfacen las condiciones iniciales: x0 = 21, y0 = 150, expresando por
separado el caso t par del caso impar.

Utilizando el método de reducción en el sistema anterior, obtenemos.

2xt+1 + 3yt+1 = 2xt + 3t . (6.13)

Ahora, aumentamos un paso en la primera de las ecuaciones y despejamos yt+1 ,

1( )
yt+1 = xt+2 − 4xt+1 + 3 3t .
6
Sustituyendo en (6.13) y simplificando se obtiene,

xt+2 − 4xt = −3t . (6.14)

Es fácil comprobar que λ1 = 2 y λ2 = −2 son las raı́ces de la ecuación caracterı́stica.

Por lo tanto, la solución general de la homogénea es

xht = k1 2t + k2 (−2)t .

El término independiente sugiere una solución particular del tipo xpt = c 3t . Susti-
tuimos en (6.14)

1 1
c 3t+2 − 4c 3t = −3t ⇒ c=− ⇒ xpt = − 3t .
5 5
La solución general de la ecuación completa es:

1 t
xt = k1 2t + k2 (−2)t − 3 , k1 , k2 ∈ IR
5
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 103

Para encontrar el valor correspondiente de yt despejamos su valor en la primera de

las ecuaciones
1( )
yt = xt+1 − 4xt + 3t .
6
Al conocer el valor de xt podemos sustituir
( )
1 1 t+1 1 t
yt = (k1 2 t+1
+ k2 (−2) t+1
− 3 ) − 4(k1 2 + k2 (−2) − 3 ) + 3 .
t t t
6 5 5
Finalmente, simplificando se llega a

1 1
yt = − k1 2t − k2 (−2)t + 3t , k1 , k2 ∈ IR
3 5

Para determinar la solución particular, es necesario tener en cuenta las condiciones

iniciales x0 = 21, y0 = 150. De aquı́ determinarı́amos las constantes k1 = 256.5 y
k2 = −235.3. Para finalizar, observemos que si t es par (−2)t = 2t , y en el caso
impar (−2)t = −2t .

EJERCICIO 38 Resolver el sistema en diferencias

{
xt+1 = −3xt + 6yt + e−t
yt+1 = −xt + 2yt − e−t ,

donde xt e yt representan los efectivos de dos especies animales y el tiem-

po t, considerado como variable discreta, se mide en años. Comprobar
que en este método, los efectivos iniciales de una de las especies han de
necesariamente condicionar los de la otra. Tomar, por ejemplo, x0 = 35
y calcular y0 . Analizar las posibilidades de extinción.

El método que utilizamos está basado en convertir el sistema anterior en una única
ecuación en diferencias que dependa de una sola variable. Para ello, sumamos a la
primera ecuación del sistema la segunda multiplicada por (-3)

xt+1 − 3yt+1 = 4e−t . (6.15)

Necesitamos hacer desaparecer yt+1 de está ecuación, y esto lo conseguimos sustitu-

yendo t por t + 1 en la primera de las ecuaciones del sistema
1( )
xt+2 = −3xt+1 + 6yt+1 + e−t−1 ⇒ yt+1 = xt+2 + 3xt+1 − e−1 e−t ,
6
sustituyendo en (6.15)
xt+2 + xt+1 = e−t (e−1 − 8) .

Esta ecuación tiene como solución general de la ecuación homogénea asociada

xht = k1 (−1)t ,
104 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

y como consecuencia de la forma del término independiente, probamos con la solución

particular xpt = Ae−t . Es fácil obtener el valor A = 16.72. Por tanto,

xt = k1 (−1)t + 16.72e−t , k1 ∈ IR

En la primera de las ecuaciones del sistema, despejamos yt y sustituimos el valor

encontrado de xt

1( ) 1
yt = xt+1 + 3xt − e−t ⇒ yt = k1 (−1)t + 9.22e−t
6 3

Observemos que en las soluciones del sistema sólo aparece una constante (k1 ), esto
obliga a que los valores iniciales de cada una de las especies tengan que estar relacio-
nados. Como sabemos que x0 = 37, sustituimos y determinamos el valor k1 = 20.28,
lo que nos permite saber el valor inicial de la segunda de las especies

1 20.28
y0 = 20.28(−1)0 + 9.22e0 = + 9.22 = 16 .
3 3

La solución particular pedida es:

xt = 20.28(−1)t + 16.72e−t

yt = 6.76(−1)t + 9.22e−t

Analicemos el problema de la extinción. Para la primera de la especie xt ,

observamos que si t es par al ser e−t y (−1)t positivos, es imposible que xt se anule.
Si consideramos un año impar

16.72
e−t = < 1,
28.28

pero si despejamos t tenemos que tomar logaritmos y aparecerá una cantidad ne-
gativa, la cual no tiene sentido biológico. En consecuencia, la primera de la especie
nunca desaparecerá.
Si repetimos el mismo razonamiento para la segunda de las especies, deducimos que
si el tiempo es positivo yt ̸= 0 . Supongamos que t es impar

9.22
e−t = >1 ⇒ t ≈ 0.3 .
6.76

Si somos estrictos, esta solución al no ser entera no deberı́amos considerarla, pero

podemos interpretarla diciendo que la extinción de la segunda especie se producirá a
los 0.3 años, o bien a los 109 dı́as.
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 105

EJERCICIO 39 Dos especies, una depredadora x y otra presa y, se re-

producen de manera que, en solitario, sus poblaciones se duplicarı́an
y cuadruplicarı́an cada año, respectivamente. La presencia de depreda-
dores produce el efecto de disminuir cada año la población de presas
en cuatro veces el efectivo de los depredadores existentes al comien-
zo del año, y la de presas hace aumentar la especie depredadora en k
veces el efectivo de las presas existentes al comienzo del año. Movimien-
tos migratorios suman cada año 20 individuos a la especie depredadora
procedente de otra región del ecosistema.

1.- Describir la evolución cuantitativa de estas especies mediante un

sistema de ecuaciones en diferencias de primer orden.

2.- Determinar el valor de k, sabiendo que la ecuación caracterı́stica

de la ecuación en diferencias de segundo orden que satisface la
población presa tiene una raı́z doble igual a 3

3.- Si inicialmente las poblaciones depredadora y presa constan de 7 y

80 individuos, respectivamente, se desea saber el número de indivi-
duos que componen cada una de ellas al cabo de 6 años y si alguna
se extingue a tiempo finito.

El sistema de ecuaciones en diferencias cuando las poblaciones están en solitario es:

{
xt+1 = 2xt
yt+1 = 4yt ,

donde el tiempo se encuentra expresado en años.

Al poner en contacto ambas especies, el sistema anterior se transforma en
{
xt+1 = 2xt + kyt + 20
yt+1 = 4yt − 4xt .

Para el segundo de los apartados, resolveremos el sistema anterior. Comenzamos

aumentando un paso en la segunda de las ecuaciones, y sustituyendo el valor de
xt+1 dado en la primera,

yt+2 = 4yt+1 − 4xt+1 = 4yt+1 − 4(2xt + kyt + 20) = 4yt+1 − 8xt − 4kyt − 80 .

Despejamos en la segunda ecuación del sistema xt y sustituimos

8
yt+2 = 4yt+1 − (4yt − yt+1 ) − 4kyt − 80 ,
4
que da lugar a la siguiente ecuación en diferencias

yt+2 − 6yt+1 + 4(2 + k)yt = −80 . (6.16)

106 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

Su ecuación caracterı́stica λ2 − 6λ + 4(2 + k) = 0, tiene por raı́ces

√
3 ± 9 − 4(2 + k) ,

y al ser el 3 una raı́z doble, entonces k = 1/4. El sistema nos quedará


 xt+1 = 2xt + 14 yt + 20

yt+1 = 4yt − 4xt ,

El tercer apartado consiste en resolver el sistema con las condiciones iniciales x0 = 7,

y0 = 80. Sabemos que la solución general de la ecuación homogénea vale

yth = (k1 + k2 t)3t .

Para encontrar una solución particular de la solución completa, probamos con yt =

A. Sustituimos en (6.16)

A − 6A + 9A = −80 ⇒ A = −20 .

La solución general buscada es

yt = (k1 + tk2 )3t − 20 , k1 , k2 ∈ IR

que nos permite, sustituyendo en

1
xt = − yt+1 + yt
4
escribir
xt = 1
4 (k1 + (t − 3)k2 ) 3t − 15 , k1 , k2 ∈ IR

Si hacemos que t = 0 
 7= 1
4 (k1 − 3k2 ) − 15

80 = k1 − 20
cuya solución es: k1 = 100, k2 = 4. Es decir,

xt = (22 + t)3t − 15

yt = (100 + 4t)3t − 20

Ahora podemos encontrar el número de presas y depredadores al cabo de 6 años

x(6) = 20397 , y(6) = 90376 .

Como se puede observar,

xt → +∞ , yt → +∞ ,
cuando t → ∞, por lo que ninguna de las dos especies desaparecerá.
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 107

EJERCICIO 40 Supongamos que la función oferta y la función deman-

da, de un animal exótico, vienen dadas por:

S(p) = 1000p − 400 D(p) = 5000 − 500p

donde p denota el precio del animal. Supongamos que el cambio del

precio viene descrito por

pt+1 = pt + α(D(pt ) − S(pt )) , α ∈ IR+ , t = 0, 1, 2, · · · .

Demostrar que este modelo es lineal y encontrar el punto de equilibrio.

Basta sustituir los valores de la oferta y de la demanda en la ecuación en diferencias

pt+1 = pt + α(D(pt ) − S(pt )) = pt + α(5000 − 500pt − 1000pt + 400)

= (1 − 1500α)pt + 5400α .

Estamos ante un modelo discreto lineal, siendo f (x) = (1 − 1500α)x + 5400α.

El punto de equilibrio se encuentra resolviendo la ecuación f (x∗ ) = x∗ , cuyo valor
es x∗ = 18/5.
Para saber si es un punto de equilibrio estable o inestable, nos fijamos en la pendiente
de la recta,

|1 − 1500α| < 1 ⇒ −1 < 1 − 1500α < 1 ,

en consecuencia, para que el punto de equilibrio sea estable tiene que ocurrir

1
0<α< .
750

Observemos que en el momento que el precio del animal corresponde al punto de

equilibrio,
18
S(18/5) = 1000 × − 400 = 3200
5
18
D(18/5) = 5000 − 500 × = 3200 ,
5

la oferta y la demanda coinciden.

108 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

EJERCICIO 41 Para dos poblaciones de un mismo tipo de bacterias,

que crecen independientemente una de la otra, obtenemos los siguientes
datos:

t 0 1 2 3 4 5 6 7 8 9 10
N1 (t) 10 11 10 20 25 60 110 140 165 175 185
N2 (t) 30 50 75 110 145 170 180 185 180 180 180

Dibujar N1 (t) y N2 (t) en función del tiempo t. ¿Cuál de estas dos pobla-
ciones se parece más a la ecuación logı́stica?.

Si utilizamos el software Mathematicar obtenemos la siguiente representación gráfi-

ca

175

150

125

100

75

50

25

2 4 6 8 10

La curva en color verde tiende a 180 cuando aumentamos el valor del tiempo, tiene
forma en S, y además su punto de inflexión se encuentra hacia la mitad de la
capacidad de carga 180/2 = 90. En consecuencia N2 (t) es la más parecida a la
ecuación logı́stica.

EJERCICIO 42 Sea r = 0.69 y K = 100. Dibujar en el plano los puntos

(N (t), N (t + 1)) correspondientes a las siguientes ecuaciones:

N (t + 1) = N (t)er (6.17)

N (t)er K
N (t + 1) = (6.18)
N (t)(er − 1) + K

Empezamos la resolución del ejercicio con la ecuación (6.18) que podemos reescribirla
200N (t)
N (t + 1) = ,
N (t) + 100
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 109

El resto lo resolveremos con el software Mathematicar . Lo iniciamos escribiendo las

funciones
f[x ] := 200 ∗ x/(x + 100)
g[x ] := x
continuamos encontrando 20 términos de la órbita correspondiente al valor x0 = 5,
iters = NestList[f, 5., 20]
cuyos valores son:
{ 5., 9.5238, 17.3913, 29.6296, 45.7142, 62.7450, 77.1084, 87.07482, 93.0909, 96.4218,
98.1783, 99.0807, 99.5382, 99.7686, 99.8841, 99.9420, 99.9710, 99.9855, 99., 99.9963,
99.9981 }
Se observa que los valores de esta población tienden al punto de equilibrio 100.

En ciertas ocasiones, es frecuente enfrentar las representaciones gráficas de N (t + 1)

en función de N (t), con la de N (t).
fg = Plot[{x, f[x]}, {x, 0, 150}, PlotStyle− > {RGBColor[1, 0, 0],
RGBColor[0, 0, 1]}, DisplayFunction− > Identity]

gi = ListPlot[Partition[Flatten[Transpose[{iters, iters}]], 2, 1],

PlotJoined− > True, DisplayFunction− > Identity]
Show[fg, gi, AspectRatio− > 1, DisplayFunction− > $DisplayFunction]
ListPlot[iters, PlotJoined− > True]

Figura 6.10:

Como puede apreciarse, la población sigue un modelo logı́stico discreto.

La ecuación (6.17) observamos que es lineal en el plano (N (t), N (t + 1)); pasa por
el origen de coordenadas y tiene de pendiente er . Para su estudio, realizamos un
análisis similar al caso anterior. En este caso los primeros términos de la órbita son:
{ 5. 9.9685, 19.8745, 39.6241, 78.9992, 157.5019, 314.0141, 626.0548, 1248.1751,
2488.5062, 4961.3735, 9891.5675, 19720.9719, 39318.0080 }
que experimentan un crecimiento exponencial, como se pone de manifiesto en las
siguientes gráficas
110 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

Figura 6.11.

EJERCICIO 43 Dada la ecuación logı́stica del crecimiento

K
N (t) = ( )
K−N (0)
1+ N (0)
e−rt

Expresar N (t + 1) en función de N (t)

Empezamos resolviendo el ejercicio calculando N (t + 1) de la expresión N (t),

k k
N (t + 1) = = ( )
k−N (0) −rt −r k−N (0) −rt
1+ N (0) e e e−r er + N (0) e

k k
= ( )= ( k−N (0) −rt
)
k−N (0) −rt
e−r e −1+1+
r
N (0) e ke−r er −1
+
1+ N (0)
e
k k

1 er
= ( er −1 ) = er −1 ,
e−r k + e−r N (t)−1 k + 1
N (t)

que una vez simplificada, nos da como solución

er kN (t)
N (t + 1) =
k + N (t)(er − 1)
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 111

EJERCICIO 44 [Modelo de Varley, Gradwell y Hassell]. Muchas po-

blaciones de insectos se rigen por el siguiente modelo
λ 1−b
f (Nt ) = Nt+1 = N , α, b, λ > 0 , t = 0, 1, 2 · · · (6.19)
α t
donde λ representa a la tasa reproductiva (λ > 1) y Nt−b /α es la fra-
cción de la población que sobreviven desde la infancia a la edad adulta
reproductiva. Podemos expresar (6.19) como
( )
1 −b
Nt+1 = N (λNt ) ,
α t

es decir, Nt+1 será igual a la fracción de insectos que sobreviven en la

generación t + 1 por el número de insectos que nacen de la generación t.
Estudiar los puntos de equilibrio del modelo.

Tenemos que estudiar el sistema dinámico discreto Nt+1 = f (Nt ), siendo

λ 1−b
f (x) = x .
α
Sus puntos de equilibrio se obtienen resolviendo la ecuación
λ 1−b λx
f (x) = x ⇒ x =x ⇒ = x,
α αxb

es decir, x∗1 = 0, x∗2 = (λ/α)1/b . Para poderlos clasificar es necesario conocer el valor
de la derivada de f (x) en cada uno de estos puntos. Al ser f ′ (x) = αλ (1 − b)x−b
tenemos
f ′ (x∗2 ) = 1 − b ⇒ |1 − b| < 1 ⇒ 0 < b < 2 .
El punto de equilibrio (λ/α)1/b es estable siempre que 0 < b < 2. En caso contrario
b > 2, el modelo tiene en (λ/α)1/b un punto de equilibrio inestable.
En el primero de los puntos no existe f ′ (x∗1 ), y por lo tanto no podemos seguir el
procedimiento anterior. No obstante su análisis a nivel biológico no es interesante
pues indicarı́a que inicialmente no existen individuos en la población.

EJERCICIO 45 Un modelo discreto frecuentemente utilizado en

dinámica de poblaciones, consiste en la ecuación
Nt
Nt+1 = f (Nt ) = Nt er(1− k ) r , k ∈ IR+ , t = 0, 1, 2, · · · .

Encontrar y analizar los puntos de equilibrio.

112 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

Procedemos de forma idéntica al ejercicio anterior, pero ahora con la función f (x) =
x
xer(1− k ) , obteniéndose de forma inmediata los puntos de equilibrio x∗1 = 0 y x∗2 = k.
La derivada de la función f (x) es
( xr )
f ′ (x) = er(1− k ) 1 −
x
,
k
que nos permite clasificar el punto de equilibrio que es más interesante. Puesto que
f ′ (k) = 1 − r, entonces el x∗2 = k será estable si 0 < r < 2.

EJERCICIO 46 Analizar los puntos de equilibrio del modelo discreto

no lineal siguiente:
kNt
Nt+1 = f (Nt ) = b, k > 0 , t = 0, 1, 2, · · · ,
b + Nt

Es inmediato comprobar que este modelo presenta en x∗1 = k − b un punto de

equilibrio estable, si k > b.

EJERCICIO 47 Calcular la posición de los puntos fijos y de los puntos

2-periódicos en el modelo logı́stico Nt+1 = 3.3Nt (1 − Nt ), t = 0, 1, 2, · · · .

Para resolver gráficamente el ejercicio, necesitaremos representar las funciones g(x) =

x, f (x) = 3.3x(1 − x), f 2 (x) = f (f (x)) y encontrar los puntos de corte. En nuestro
caso
f 2 (x) = 3.32 x(1 − x)(1 − 3.3x + 3.3x2 ) .

Figura 6.12: Puntos de equilibrio y 2-periódicos de f (x) = 3.3x(1 − x).

6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 113

Los puntos de equilibrio son xA = 0, xB = 0.6969.. y los puntos 2-periódicos xc =

0.47492701.., xD = 0.8236032.

EJERCICIO 48 Supongamos que la tasa de crecimiento de una pobla-

ción P satisface
g(P ) = 0.03P (1 − P/600) .
El modelo discreto de crecimiento logı́stico viene dado por

Pt+1 = Pt + g(Pt ) , t = 0, 1, 2, · · · .

1.- Encontrar la población cuando g(P ) es cero y cuando tiene un máxi-

mo (vértice).

2.- Calcular P1 , P2 , P3 para un valor inicial P0 = 100. Encontrar los

puntos de equilibrio.

1.- Es evidente que g(P ) se anula para los valores P = 0 y P = 600. Al ser una parábola
que corta al eje de abscisas en los puntos 0 y 600, su vértice estará situado en P = 300
siendo g(300) = 4.5. Es decir el valor máximo de la parábola será el punto (300, 4.5).

2.- Para encontrar las poblaciones pedidas solamente tendremos que sustituir los valores
adecuados en el modelo presentado. De esta manera,

P1 = P0 + 0.03P0 (1 − P0 /600) = 102.5 .

De forma similar P2 = 105.05, P3 = 107.65.

EJERCICIO 49 La población de China en 1980 era de 985 millones,

y el censo de 1990 mostró que la población habı́a crecido hasta 1.137
millones. Suponiendo que la población crece según la siguiente ley de
crecimiento discreto exponencial

Pt+1 = (1 + r)Pt , t = 0, 1, 2, · · · , (6.20)

donde t es el número de décadas después de 1980 y Pt la población t

décadas después de 1980.

1.- Encontrar la constante r de crecimiento y predecir la población

para el año 2000.

2.- Encontrar el tiempo necesario para que se duplique la población.

114 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

1.- Llevando los datos en (6.20),

P (1990) = P1 = (1 + r)P (1980) = (1 + r)P0 ⇒ 1.137 = (1 + r)985 ,

y el valor buscado es r = 0.1543.

2.- La segunda parte del ejercicio se deduce de la ecuación

ln 2
Pt = 1.1543t P0 = 2P0 ⇒ t= = 4.83 décadas ,
ln 1.1543
se necesitan 48.3 años para que la población se duplique.

EJERCICIO 50 Un invertebrado vive en un lago que está afectado por

el efecto de la contaminación que penetra lentamente en el ecosistema.
La dinámica poblacional para este invertebrado viene dada por el si-
guiente modelo de crecimiento exponencial no autónomo

Pn+1 = (1 + k(tn ))Pn , P0 = 40.000 , (6.21)

donde tn = n es el número de dı́as desde la medida inicial de la población

y k(t) = 0.08 − 0.01t es la tasa de crecimiento, que claramente decrece con
el tiempo.

1.- Encontrar la población para este organismo en los próximos 5 dı́as.

2.- Cuando la tasa de crecimiento cae a cero, la población alcanza su

máximo. Encontrar cuando ocurre y el tamaño de la población en
este momento.

3.- Encontrar cuando se da el nivel máximo de polución, lo que obliga

a la extinción de la especie.

1.- La población para el primer dı́a es P1 = 1.08 × 40000 = 43200, ya que k(t0 ) = 1.08.
Para el resto de los dı́as la solución es P2 = 46224, P3 = 48997, P4 = 51447, P5 =
53505.

2.- La tasa de crecimiento es cero al cabo de los t = 8 dı́as. La población será de

P8 = 1.01P7 = 1.01 × 1.02 × 1.03 × 53505 = 56775.

3.- Para encontrar la respuesta del tercer apartado, tenemos que determinar cuando el
factor 1 + k(tn ) se anula. En nuestro caso,

1 + 0.08 − 0.01t = 0 ⇒ t = 108 .

Esto ocurre al cabo de los 108 dı́as que es la cota superior teórica para la extinción
de la especie.
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 115

EJERCICIO 51 Para una determinada especie animal se considera que

P es la proporción de individuos que como máximo pueden pertenecer
a un hábitat concreto (esto es, P = 0 indica la ausencia y P = 1 indica
que no puede haber más individuos). Además, se ha comprobado que
las tasas de fertilidad y mortalidad son, respectivamente, las siguientes:
3a a a
f (P ) = (1 − P ) ; m(P ) = 1 − + P,
8 8 8
siendo a un número real comprendido entre 2 y 6.

1.- Comprobar que la ecuación en diferencias yt+1 = a2 yt (1−yt ), modeliza

a la dinámica de la población.

2.- Encontrar los puntos de equilibrio del modelo y estudia su estabi-

lidad.

1.- Para un año cualquiera t, tenemos que

yt+1 = yt + (fertilidad - mortalidad)yt

es decir
(( ) ( )
3a a a ) ayt
yt+1 = yt + (1 − yt ) − 1 − + yt yt = (1 − yt ) .
8 8 8 2

2.- La función que define al modelo es f (x) = a2 x(1 − x). Para encontrar los puntos fijos
o de equilibrio resolvemos la ecuación no lineal f (x) = x, cuyas soluciones son:
a−2
x∗1 = 0 ; x∗2 = .
a

Para analizar la estabilidad de los puntos fijos, en primer lugar, calculamos la deri-
vada,
a
f ′ (x) = − a x ,
2
y sustituimos cada uno de los puntos,
a
f ′ (0) =
2
Para que este punto sea estable tiene que ocurrir que |f ′ (0)| < 1, lo que obliga a
que el valor del parámetro a sea menor que dos, y esto no puede ocurrir ya que el
enunciado indica que 2 < a < 6.
Para el segundo punto,
a−2 4−a 4−a
|f ′ (x∗2 )| = |f ′ ( )| = | |<1 ⇒ −1 < <1 ⇒ 2 < a < 6.
a a a
En consecuencia, el punto de equilibrio x∗1 es inestable y el x∗2 es estable.
116 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

EJERCICIO 52 Dado el sistema dinámico discreto xt+1 = f (xt ) = x3t − x2t + 1.

Analizar su comportamiento a largo plazo.

En primer lugar realizaremos el estudio cuantitavivo encontrando los puntos de

equilibrio del modelo y su clasificación.

F (x) = x ⇒ x3 − x2 + 1 = x ⇒ x∗1 = 1 (doble) , x∗2 = −1

si calculamos la derivada f ′ (x) = 3x2 − 2x y sustituimos en los puntos de equilibrio,

|f ′ (1)| = 1

estamos ante un punto de equilibrio x∗1 = 1 dudoso. Por otro lado,

|f ′ (−1)| = 5 > 1

en este caso el punto de equilibrio x∗2 = −1 es inestable.

Para poder clasificar el punto de equilibrio dudoso, encontraremos 25 términos de

la órbita de f correspondiente a la semilla x0 = 1.2
NestList[p,1.2,25]

{1.2, 1.288, 1.47778, 2.04338, 5.35651, 125.998, 1.98439 × 106 , 7.81414 × 1018 , 4.77137 × 1056 ,
1.08625 × 10170 , 1.281696405368169 × 10510 , 2.105501227265631 × 101530 ,
9.33397206326970 × 104590 , 8.1320396988399 × 1013772 , 5.3777235179070 × 1041318 ,
1.5552328140679 × 10123956 , 3.761717978576 × 10371868 , 5.32302735792 × 101115605 ,
1.50825958784 × 103346817 , 3.43105977494 × 1010040451 , 4.0391022997 × 1030121354 ,
6.589531805 × 1090364063 , 2.861301849 × 10271092191 , 2.342561635 × 10813276574 ,
1.285502955 × 102439829723 , 2.12431657 × 107319489169 }

que podemos observar tiende a infinito. Por otro lado, la órbita de f correspondiente
a la semilla x0 = 0.2 es:
NestList[p,0.2,25]

{0.2, 0.968, 0.970015, 0.971786, 0.973356, 0.974757, 0.976015, 0.977152, 0.978184,

0.979126, 0.979988, 0.980781, 0.981513, 0.98219, 0.982819, 0.983404, 0.98395, 0.984461,
0.98494, 0.985391, 0.985814, 0.986214, 0.986591, 0.986949, 0.987287, 0.987608}

que se estabiliza en el punto de equilibrio x∗1 = 1. En consecuencia, este punto de

equilibrio es un nodo.
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 117

0.99

0.98

0.97

0.96

5 10 15 20 25 30

Figura 6.13: Órbita para la semilla x0 = 0.2

Como ejercicio idéntico al anterior, se propone el estudio del siguiente sistema

dinámico: √
xt+1 = f (xt ) = 4xt − 3
que presenta un punto de equilibrio estable en 3, y un punto de equilibrio inestable
en 1.
118 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

EJERCICIOS PROPUESTOS

EJERCICIO 53

1.- Calcular k para que la ecuación en diferencias de segundo orden:

xt+2 − 2kxt+1 + (k + 1)xt = t2 + 3 + 3t

tenga a 1 y a 3 como raı́ces de la ecuación caracterı́stica.

Resolver, a continuación, la ecuación completa con las condiciones inicia-
les: x0 = 10 , x1 = 18.

2.- El crecimiento de una especie viene descrito por la siguiente ecuación en

diferencias:
xt+2 − 4xt = −3t+1 , t = 0, 1, 2, 3, · · ·
donde xt representa a la cantidad de animales en el año t. Determinar
el número de animales al finalizar un año cualquiera ”t”, sabiendo que
inicialmente hay 10 y que transcurrido un año su número es de 20

3.- El incremento de la población de una determinada especie animal en un

año es la mitad del incremento del año anterior, si no intervienen factores
externos. La población inicial es de 950 individuos y de 975 al finalizar el
primer año.

Escribir la ecuación en diferencias que modeliza a la situación plan-

teada.
Determinar la cantidad de individuos de dicha especie al finalizar
un año cualquiera ”t”.
Estudiar el comportamiento a “largo plazo ”de la población.

4.- Sea xt el número de individuos de una determinada especie de animales

en el tiempo t. Se sabe que año tras año sobreviven la tercera parte de
los animales y además se incorporan 200 a la población.

Construir un modelo discreto lineal para la situación planteada.

Calcular los seis primeros términos de las órbitas correspondientes
a las semillas:
x0 = 90 , x0 = 600 .

Construir los diagramas de Cobweb del apartado anterior, e inter-

pretar biológicamente los resultados obtenidos.
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 119

5.- Se sabe que la evolución de una población de una determinada especie

de peces viene dada por el sistema dinámico discreto lineal

yt+1 = f (yt )

Encuentra una función f de tal manera que la población a largo

plazo se estabilice en 10 individuos independientemente del número
inicial de peces.
Comprueba el resultado anterior por medio del diagrama de Cobweb
con los valores iniciales: y0 = 1 peces y y0 = 20 peces.

6.- La siguiente ecuación en diferencias describe la evolución de una pobla-

ción en años sucesivos,
10er xt
xt+1 = , t = 0, 1, 2, 3, ··· .
10 + (er − 1)xt
Demuestra que si el parámetro r es positivo, entonces la población
se estabiliza en 10 individuos, mientras que si r es negativo, la po-
blación desaparecerá.

7.- La siguiente ecuación en diferencias describe la población de insectos en

un manglar en años sucesivos,

xt+1 = αxt e−xt , t = 0, 1, 2, · · ·

siendo α un parámetro positivo y xt el número de insectos en el año t

¿Cuál debe ser el valor de α para que el punto de equilibrio no trivial sea
estable?

8.- Si sobre una población no influyen factores que modifiquen el crecimiento,

se observa que,
( )t
1
5xt+2 − 6xt+1 + xt = , t = 0, 1, 2, · · ·
5
siendo xt el número de individuos en el año t. Encontrar el número de
individuos en el quinto año, sabiendo que inicialmente eran 10 y al año
siguiente 20.

9.- La siguiente ecuación en diferencias representa la dinámica de una po-

blación
αyt
yt+1 = f (yt ) =
10 + yt
donde yt representa al número de individuos en el tiempo t.

9.a.- Encuentra el valor del parámetro α para que uno de sus puntos de
equilibrio sea 40.
9.b.- Con el valor calculado para α, encontrar y clasificar los puntos de
equilibrio del modelo
120 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

9.c.- ¿Qué puedes decir de la evolución “a largo plazo”de esta población?

10.- Contestar de forma razonada a las siguientes cuestiones:

10.a.- Sea la ecuación en diferencias yt+2 − 2yt+1 + yt = 0, donde yt repre-

senta a la cantidad de individuos en el año t. Si el número inicial de
individuos es 2 y al cabo de un año es 5, ¿cuál será el valor de la
población al cabo de 10 años?
10.b.- Encontrar la solución general de la ecuación en diferencias yt+2 −
2yt+1 + yt = 8.

11.- La dinámica de una determinada especie responde a la siguiente ecuación

en diferencias:
xt+1 = 2xt e1−xt ; t ∈ IN
Responder, de forma razonada, a las siguientes cuestiones:

¿Cuáles son los puntos fijos (puntos de equilibrio)?

¿Cómo son los puntos fijos con respecto a la estabilidad?

12.- Considera la ecuación en diferencias:

xt−1 = xt ext −2

Encuentra los puntos de equilibrio de la ecuación y estudia su estabilidad.

13.- Supongamos los siguientes modelos matemáticos utilizados en dinámica

de poblaciones de aves:

Modelo 1: Pt+1 = λPt ; con λ > 0

Modelo 2: Pt+1 = (λPt + a)Pt ; con λ < 0 ; a > 0

La información de la que se dispone es la siguiente:

El primer año de observación se contaron 500 aves y el segundo

1000.
La población no se espera que se extinga ni que crezca ilimitada-
mente.
Si llegara a observarse 4000 aves es de esperar que ese número no
varı́e ya en el transcurso del tiempo (Es decir, P=4000 se mantendrı́a
constante).

Justifica de manera razonada, si los dos modelos planteados cumplen los

datos biológicos o no. Además, investigando sobre esta población se ha
descubierto que la población constante de 4000 aves es estable. Justifica
si los modelos matemáticos propuestos cumplen esta nueva condición.
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 121

14.- Dada la ecuación en diferencias,

1( 3 )
xt+1 = −xt + 4x2t − 3xt + 2
2
Estudiar la estabilidad de sus puntos de equilibrio.
Suponiendo que xt represente al tamaño de una población, discutir
su evolución a largo plazo, dependiendo de los valores iniciales x0

15.- Sea yt el número de individuos de una población en el año t. La evolución

de dicha población viene dada por la siguiente ecuación en diferencias:

yt+2 + byt+1 + cyt = 2t

Si la solución general de la ecuación homogénea asociada es yth =

c1 2t + c2 3t , encontrar el valor de los parámetros a y b.
Si inicialmente el número de individuos era de y0 = 4 y al año si-
guiente era de 9, encontrar el número de individuos para un año
cualquiera t.

16.- Para una especie se ha comprobado que las tasas de fertilidad y mortali-
dad vienen dadas por,

1 1
f (P ) = , m(P ) = 1 − ,
2 1+P

respectivamente.

Determinar la ecuación en diferencias que rige la dinámica de dicha

población.
Calcular y clasificar los puntos de equilibrio del modelo.
Realizar una interpretación del comportamiento de la especie a largo
plazo, a partir de los resultados obtenidos en los apartados anterio-
res.

17.- Se considera la serie de pesos siguiente:

p37 = 1.35 Kg. ; p38 = 2.90 Kg. ; p39 = 1.74 Kg. ; p40 = 3.06 Kg. ; p41 = 1.35 Kg.

que se ajusta a una ecuación en diferencias del tipo xt+1 = xt (a−xt ), siendo
a un parámetro a determinar.

Estimar el valor de a.
Determinar los puntos de equilibrio del modelo y estudiar su esta-
bilidad.
122 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

18.- Considera la siguiente versión del modelo logı́stico discreto de crecimien-

to:
xt+1 = f (xt ) = 3.2xt (1 − 0.25xt )
Encuentra, mediante una simulación con ordenador, las primeras 100 ite-
raciones de la órbita de x0 y dibújala, suponiendo que x0 = 2.75, que
x0 = 2.5 y finalmente con x0 = 2. ¿Qué tipo de equilibrio tiene el mode-
los?¿Existen 2-ciclos? Justifica las respuestas.

19.- Una población de palomas parte de 1000 ejemplares. Se reproduce de tal

manera que la población en cada año es el doble que la del año anterior
más cuarenta y cinco cuartos de la de hace dos años. Además cada año se
extraen 20 individuos para su estudio. ¿Cuál es la población de la colonia
en un año cualquiera t? ¿cuántos individuos hay después de cinco años?

20.- La siguiente ecuación en diferencias describe la población de ardillas en

años sucesivos,

xt+1 = x3t − 3x2t − 3xt + a , t = 0, 1, 2, 3, ··· .

siendo a un parámetros positivo y xt el número de ardillas en el año t .

Encuentra el valor del parámetro a sabiendo que existe un punto de

equilibrio en x∗ = 2
Clasificar los puntos de equilibrios que tienen sentido biológico para
conocer el comportamiento a largo plazo de la población.

21.- Sea el modelo discreto logı́stico

Nt+1 = 3.3Nt (1 − Nt ) ,

donde Nt representa al número de individuos de la población en el perı́odo

t. Clasifica el punto de equilibrio no trivial, y comprueba el resultado
haciendo uso del diagrama de Cobweb.

22.- Calcular y clasificar los puntos de equilibrio del siguiente modelo discreto,
con r > 0 y k > 0,
( ( ))
Nt
f (Nt ) = Nt+1 = Nt 1 + r 1 − , t = 0, 1, 2, 3, .
k

23.- La siguiente ecuación en diferencias:

α xt
xt+1 = , α,β > 0, xt ≥ 0 ,
1 + β xt
fue propuesta por Kaplan & Glais en 1995 y juega un papel muy importante
en análisis de modelos no lineales genéticos y en redes neuronales.

23.a.- Encontrar y analizar los puntos de equilibrio

6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 123

23.b.- Sea α = β = 1. Dibujar de forma aproximada el diagrama en telaraña

(cobweb) tomando como semilla x0 = 4.

24.- Si sobre una población no influyen factores que modifiquen el crecimiento,

se observa que,
( )t
1 1
(yt+2 − yt+1 ) − (yt+1 − yt ) = ,
3 3

siendo yt el número de individuos en el tiempo t.

Explicar el significado “biológico”de la ecuación anterior

¿Crecerá la población a largo plazo?.

25.- Indicar en cada una de las siguientes ecuaciones si es lineal o no lineal.

Si es lineal determinar la solución; si es no lineal encontrar y analizar el
tipo de puntos de equilibrio.

xt = (1 − α)xt−1 + βxt α , β ∈ IR+

xt
xt+1 =
1 + xt
xt+1 = xt e−axt a ∈ IR+
(xt+1 − α)2 = α2 (x2t − 2xt + 1) α ∈ IR+
k
xt+1 = k1 , k2 , k ∈ IR+
k1 + k2 /xt

26.- Se sabe que la evolución de una población de una determinada especie

de peces viene dada por el sistema dinámico discreto lineal

yt+1 = f (yt )

26.a.- Encuentra una función f de tal manera que la población a largo

plazo se estabilice en 10 individuos independientemente del número
inicial de peces.
26.b.- Comprueba el resultado anterior por medio del diagrama de Cobweb
con los valores iniciales: y0 = 1 peces y y0 = 20 peces.

27.- Sea yt el número de individuos de una población en el año t. La evolución

de dicha población viene dada por la siguiente ecuación en diferencias:
yt+2 + ayt+1 + byt = 3

Si la solución general de la ecuación homogénea asociada es yth =

C1 1t + C2 2t , encontrar el valor de los parámetros a y b.
Si inicialmente el número de individuos era de y0 = 2 y al año si-
guiente su número era y1 = 4, encontrar la población después de 4
año.
124 Capı́tulo 6 Aplicaciones de los sistemas dinámicos discretos

28.- Si sobre una población no influyen factores que modifiquen el crecimiento,

se observa que,

yt+2 − 4yt+1 + 3yt = 3t , t = 0, 1, 2, 3 · · · ,

siendo yt el número de individuos en el año t. Encontrar el número de

individuos en el tercer año, sabiendo que inicialmente eran 2 y al año
siguiente 8 individuos.

29.- La evolución de una población xt viene determinada por el siguiente mo-

delo discreto exponencial con inmigración y emigración,

xt+1 = (1 + r)xt − µ , t = 0, 1, 2, · · ·

siendo el parámetro positivo µ la diferencia entre el número de personas

que entran y las que salen, el parámetro r la tasa de crecimiento de la
población, y x0 el número inicial de individuos.

Estudiar el comportamiento a largo plazo del modelo según los di-

ferentes valores del parámetro r.
Comprueba el resultado anterior por medio del diagrama de Cob-
web, para r = 0.2 y µ = 10.

30.- Dos especies conviven de acuerdo con el siguiente modelo discreto:

{
xt+1 = −xt + yt + 3t
yt+1 = 4xt + 2yt + 3t ,

donde el tiempo t se encuentra expresado en años. Hallar sus funciones

de efectivos.

31.- Dos especies que conviven en un mismo territorio, evolucionan del modo
descrito por el sistema de ecuaciones en diferencias siguiente:
{
xt+1 = 7xt − 2yt − t + 2
yt+1 = 6xt − yt + 5t ,

donde el tiempo t se encuentra expresado en años.

Si, inicialmente, el número de individuos de cada especie es x0 =

70, y0 = 251, resolver el sistema para obtener los efectivos en función
de t y analizar el comportamiento a la larga de las dos especies.
Comprobar que al cabo de un año se ha extinguido la primera es-
pecie y analizar el comportamiento a posteriori de la segunda.

32.- El crecimiento de dos especies que coexisten viene descrito por el sistema
de ecuaciones en diferencias siguiente:
{
xt+1 = (9 − a)xt + (4 − b)yt + 5t ; x0 = 10
yt+1 = (b − 2)xt + 5yt + 5t ; y0 = 45
6.5 Ejemplo de modelo discreto para la economı́a. Modelo de la telaraña. 125

donde el tiempo t se encuentra expresado en años. Si la ecuación ho-

mogénea, asociada a la ecuación en diferencias de segundo orden, que
resulta de eliminar yt en el sistema es:

xt+2 − 8xt+1 + 15xt = 0.

Resolver todos los sistemas que cumplan las condiciones anteriores y

analizar el comportamiento a la larga de las dos especies.