1 - Modelación Matemática
1 - Modelación Matemática
1 - Modelación Matemática
DEPARTAMENTO DE MATEMÁTICAS
MODELOS MATEMÁTICOS
DISCRETOS EN LA EMPRESA
GRADO EN ESTADÍSTICA Y EMPRESA
TEORÍA DE MODELOS DISCRETOS
2. Diagonalización de matrices 23
2.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.2. Valores y vectores propios . . . . . . . . . . . . . . . . . . . . . . . . 25
2.3. Diagonalización de matrices cuadradas. . . . . . . . . . . . . . . . . . 27
2.4. Cálculo de la potencia de una matriz diagonalizable . . . . . . . . . . 30
2.5. Métodos numéricos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.5.1. Método de las potencias . . . . . . . . . . . . . . . . . . . . . 32
2.5.2. El método de las potencias con cambio de escala . . . . . . . . 34
2.5.3. Deflación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3
ÍNDICE GENERAL
1
3
4 ÍNDICE GENERAL
INTRODUCCIÓN A LA
MODELIZACIÓN MATEMÁTICA
1.1. Introducción
Una de las herramientas más interesantes que actualmente disponemos para analizar
y predecir el comportamiento de un sistema económico es la construcción y posterior
simulación de un modelo matemático. Son muchas las razones que justifican la edad
de oro que hoy en dı́a vive la modelización matemática, pero debemos de destacar,
en primer lugar, el mejor conocimiento de los procesos económicos, y en segundo
lugar, el espectacular avance de los ordenadores y el software matemático.
Puesto que este trabajo es una introducción al estudio de los Modelos Matemáticos
en la Empresa, es conveniente comenzar esta primera sección precisando lo que en-
tendemos por un modelo matemático.
Con frecuencia la palabra modelo tiene distintas interpretaciones, nosotros la apli-
caremos en el sentido dado por el profesor Sixto Rı́os, [53]:
7
8 Tema 1 Introducción a la modelización matemática
Generalmente los métodos que se utilizan para estudiar un fenómeno económico son
la construcción de un modelo matemático o bien el uso del método cientı́fico, el cual
está basado en:
fenómeno económico.
Es evidente, que una de las ventajas del uso de los modelos matemáticos es su bajo
costo, si lo comparamos con los modelos fı́sicos. Por ejemplo, es mucho más barato
y rápido elaborar un modelo matemático que describa la evolución del número de
empresas que empezar con un determinado número de ellas y esperar cierto tiempo
para poder experimentar con ellas.
A continuación expondremos más detenidamente los pasos que debemos seguir para
construir un modelo matemático.
Tenemos que definir las variables y ver la manera en que están relacio-
nadas.
Debemos definir los parámetros del modelo, y asegurarnos de que cual-
quier otro parámetro es redundante.
ajustados por los valores, 100, 2 × 100, 22 × 100, 23 × 100, 24 × 100, · · · , que
corresponden a un crecimiento exponencial. Este último paso es el verdadera-
mente importante en el proceso de modelado.
Tabla 1.1
y(t) = N αt .
Predicción. Una vez que por la etapa anterior nos hemos asegurado de la
validez del modelo, pasamos a la etapa de predicción. Por ejemplo, en la si-
tuación que estamos analizando, si queremos obtener 3.200 células a partir de
400 células, necesitamos que pasen 3 perı́odos que equivalen a 60 minutos.
A pesar de la gran importancia que hoy en dı́a tienen los modelos matemáticos,
tenemos que tener en cuenta la siguiente observación. Por lo general, en los modelos
teóricos, se consideran sólo las relaciones cuantitativas entre las variables depen-
dientes e independientes del mismo, y entonces entran en juego las matemáticas.
Ahora bien, estos modelos describirán relaciones entre los organismos, pero nunca
pueden darnos el sentido económico del proceso. Por tanto, será imprescindible la
experimentación biológica.
Por último, un modelo matemático tiene que tener las siguientes cualidades:
Debe ser coherente, es decir, tiene que dar cuenta de todas las observaciones
anteriores y permitir prever el comportamiento futuro del fenómeno económico.
Y por último, debe ser flexible, en el sentido de que pueda ser cambiado y
adaptado a nuevas situaciones.
• Ecuaciones en diferencias.
• Teorı́a de bifurcaciones.
• Ecuaciones diferenciales (ordinarias y parciales).
• Análisis numérico.
• Procesos estocásticos.
Modelos mixtos:
16 Tema 1 Introducción a la modelización matemática
• Cadenas de Markov.
• Modelos de Leslie.
• Modelos de Lefkovitch.
Modelos en bioquı́mica.
Modelos en la genética.
Modelos en la epidemiologı́a.
“Los modelos son erróneos ... pero muchos de ellos son útiles ”
Entonces, ¿cómo pueden ser útiles si están equivocados? La respuesta a esta pre-
gunta puede ser que por la misma razón que en el pasado mapas erróneos, donde se
suponı́a que la tierra era plana y con distancias equivocadas, fueron muy útiles para
viajar.
A finales del siglo XIX Federico Engels se lamentaba de lo poco que estaban introdu-
cidas las Matemáticas en la vida real. Todo cambia a principios del siglo XX, cuando
Michaelis y Menten proponen un modelo bioquı́mico (que aún se utiliza hoy en dı́a),
para describir la catálisis enzimática. Dos años después, Lee presentó un modelo para
explicar los paradójicos efectos de las radiaciones sobre las células. Ahora, tenemos
que trasladarnos hasta mediados de siglo para encontrar otro ejemplo interesante.
Basándose en la propuesta de Galileo de establecer relaciones cuantitativas entre
magnitudes medibles, se intentaba encontrar un modelo que relacionase la intensi-
dad de un estı́mulo y la duración del mismo. El esfuerzo fue inútil, poniéndose de
manifiesto que para tener éxito en el modelado es importante atender no solamente
a la experimentación, sino también acertar en el tipo de relaciones cuantitativas a
estudiar. Además, a la hora de construir un modelo es fundamental saber separar la
información relevante que conocemos del problema de la que no lo es.
El contraste con esta última situación lo encontramos en el modelo de Hodking
y Huxley para la generación y transmisión del impulso nervioso. En este caso, se
proponı́an relaciones entre variables que fı́sicamente tenı́an sentido. Este modelo
construido en 1952 suele ponerse como ejemplo de modelo matemático aplicado a
la Biologı́a, de hecho, algunos autores piensan que juega un papel en la Neurologı́a
semejante a las ecuaciones de Maxwell en el estudio del Electromagnetismo, ya que
a través de él es posible explicar todas las propiedades experimentales conocidas
respecto a la generación y propagación del impulso nervioso. Al mismo tiempo, el
18 Tema 1 Introducción a la modelización matemática
modelo sugerı́a que la dinámica de muchos procesos económicos debı́a ser no lineal.
A partir de este momento, empieza la edad de oro para la construcción y posterior
interpretación de modelos matemáticos aplicados a la Biologı́a. En los años 60 se pu-
blicaron un gran número de trabajos, especialmente los relacionados con el sistema
nervioso, muchos de ellos con escaso interés práctico. El siguiente paso importan-
te se da en la década de los 70 cuando se descubre que las soluciones de sistemas
dinámicos presentaban un comportamiento caótico. Un ejemplo lo encontramos en
el modelo logı́stico de R. May, que supuso toda una revolución comparable al impac-
to causado por el modelo de Hodgkin y Huxley. La teorı́a del caos inmediatamente
entusiasmó a biólogos, fı́sicos y matemáticos, dedicados al estudio de los modelos
matemáticos.
De toda formas, muchas situaciones muy distintas, como pueden ser la actividad ce-
rebral, el electrocardiograma, la dinámica de poblaciones, el desarrollo embrionario,
la evolución de las enfermedades, son escenarios muy difı́ciles de modelar a través
de modelos elementales. Sin embargo, podemos realizar las simplificaciones conve-
nientes que expliquen parcialmente el comportamiento del sistema o bien aplicar
unas nuevas herramientas matemáticas, como es el uso de la geometrı́a fractal, para
explicar la variabilidad de la frecuencia del corazón.
La Econometrı́a, con una aportación importante de la Estadı́stica, es la ciencia que
se ocupa de la construcción de los modelos económicos; se inició con el profesor de
la Universidad de Oslo Ragnar Frisch. Con la ayuda de otros colaboradores crearon
en 1930 la Sociedad de Econometrı́a de los Estados Unidos, desde donde se edita
desde 1932 la revista Econometrı́a.
Un estado viene determinado por cada uno de los valores que puede tomar una
variable para cada uno de los valores del tiempo.
De una manera mas formal, un sistema dinámico está compuesto por la terna
(E, T, ϕ), siendo E el espacio de fase, T el conjunto del tiempo, y ϕ : T xE → E que
cumple:
(b) ∀x ∈ E; ϕ(0, x) = x
DIAGONALIZACIÓN DE MATRICES
2.1. Introducción
El álgebra matricial proporciona herramientas elementales para simplificar y resolver
problemas donde intervienen un número elevado de datos. El siguiente ejemplo es un
caso particular de la teorı́a de grafos, donde utilizamos la multiplicación de matrices
para modelar el proceso de dispersión de una enfermedad contagiosa.
EJEMPLO 2.1
Supongamos que tres personas han contraı́do una enfermedad contagiosa y que
este primer grupo tiene contacto con cinco personas de un segundo grupo. Estos
contactos, llamados contactos directos, pueden representarse mediante una matriz
A3×5 = (aij ) con aij = 1 si la i-ésima persona del primer grupo ha estado en contacto
con la j-ésima persona del segundo grupo, y aij = 0 si no han estado en contacto.
1 0 1 1 0
A= 0 1 1 0 1
1 0 0 1 0
Por ejemplo, el 1 que aparece en la fila 3 columna 4, significa que la tercera persona
del primer grupo (los infectados) ha estado en contacto con la cuarta persona del
segundo grupo.
Supongamos ahora que un tercer grupo de cuatro personas ha tenido diversos contac-
tos directos con individuos del segundo grupo. Esto también podemos representarlo
23
24 Tema 2 Diagonalización de matrices
Los contactos indirectos o de segundo orden entre las personas del primer y tercer
grupo se representan mediante una matriz C3×4 que es el producto AB
0 0 2 1
C = AB = 1 1 0 1
0 0 2 0
Se observa, que la tercera persona del último grupo tiene cuatro contactos indirectos,
mientras que la primera y la segunda sólo tienen uno.
que A⃗u = λ⃗u. Al vector ⃗u anterior se le conoce con el nombre de vector propio (o
autovector) de A asociado al valor propio λ.
EJEMPLO 2.2
Aplicaremos la definición anterior para encontrar los valores y vectores propios de la matriz
A asociada al endomorfismo (2.1).
3 2 −2 x x (3 − λ)x + 2y − 2z = 0
1 3 −1 y = λ y ⇒ x + (3 − λ)y − z = 0 (2.2)
1 2 0 z z x + 2y + (−λ)z = 0
Este sistema homogéneo debe tener solución distinta de la trivial (x, y, z) ̸= (0, 0, 0), por
tanto según el teorema de Rouché-Fröbenius el determinante de la matriz de los coeficientes
debe ser nulo
3−λ 2 −2
1 3 − λ −1 = 0 ⇒ λ3 − 6λ2 + 11λ − 6 = 0 .
1 2 −λ
Hacemos uso de la regla de Ruffini para resolver esta ecuación de tercer grado,
λ1 = 1 , λ2 = 2 , λ3 = 3 .
A continuación para cada uno de estos valores propios encontramos sus vectores propios
asociados. Si sustituimos λ1 = 1 en el sistema (2.2)
}
2x + 2y − 2z = 0
⇒ x = α, y = 0, z = α.
x + 2y − z = 0
S1 = {(α, 0, α) : α ̸= 0} ,
dicho subespacio tiene dimensión uno y está generado por el vector ⃗u1 = (1, 0, 1).
Del mismo modo, para λ2 = 2
}
x + 2y − 2z = 0
⇒ z =β, x = 0, y =β.
x+y−z = 0
S2 = {(0, β, β) : β ̸= 0} ,
Las raı́ces reales del polinomio caracterı́stico serán los autovalores o valores propios
de la matriz A. Llamaremos orden de multiplicidad de un autovalor λ a la multipli-
cidad de la raı́z λ del polinomio caracterı́stico.
Resumiendo, una vez obtenidos los valores propios λi de la matriz A resolviendo
la ecuación caracterı́stica (2.4), sustituimos dicho valores en el sistema homogéneo
(2.3) y calculamos los subespacios de vectores propios asociados. Es fácil probar que
las soluciones del sistema homogéneo constituyen un subespacio vectorial Si de IRn
cuya dimensión será siempre mayor o igual a uno.
Es conocido que el vector ⃗u1 respecto de la base B ′ tiene por coordenadas (1, 0, 0, · · · , 0).
Si calculamos la imagen de este vector,
d1 0 · · · 0 1 d1 1
0 d2 · · · 0 0 0 0
f (⃗u1 ) = D⃗u1 = .. .. .. .. = .. = d1 .. = d1⃗u1 .
. . . . . .
0 0 · · · dn 0 0 0
En general f (⃗ui ) = di⃗ui , i = 1, 2, · · · , n lo que nos indica que ⃗ui es un vector propio
asociado al valor propio di = λi , i = 1, 2, · · · , n.
Notemos que acabamos de probar que si el endomorfismo f es diagonalizable, enton-
ces es posible encontrar una base de IRn formada por los autovectores de f . Por otro
lado, si elegimos como base de IRn aquella que está formada por los autovectores de
f , entonces existirán n escalares λi , i = 1, 2, · · · , n tales que f (⃗ui ) = λi⃗ui . En este
caso, la matriz asociada al endomorfismo respecto de esta base es la matriz diagonal
λ1 · · · 0
D = ... . . . ...
0 ··· λn
LEMA 2.3.1 Autovectores asociados a autovalores distintos dos a dos son lineal-
mente independientes
Si estos dos vectores no son linealmente independientes, entonces ⃗vi = k⃗vj , lo que
implica que
f (⃗vi ) = f (k⃗vj ) ⇒ λi⃗vi = kλj ⃗vj = λj ⃗vi .
Pero al ser vectores no nulos, esta última igualdad implicarı́a que λi = λj , en contra
de lo supuesto.
Como consecuencia del lema, vectores no nulos pertenecientes a distintos subespacios
propios son linealmente independientes.
LEMA 2.3.2 La dimensión del subespacio propio asociado a un cierto valor propio
es como mucho igual al orden de multiplicidad del autovalor.
2.3 Diagonalización de matrices cuadradas. 29
EJEMPLO 2.3
Sea f : IR3 → IR3 un endomorfismo, cuya matriz respecto a una base B = {⃗v1 , ⃗v2 , ⃗v3 } es:
1 −1 −3
A= 2 4 3
0 0 3
cuyas soluciones son λ1 = 2 y λ2 = 3 como raı́z doble. Al no salir los tres autovalores
diferentes no podemos asegurar que la matriz A sea diagonalizable; tenemos que ver
si es posible encontrar una base del espacio vectorial R3 formada por autovectores
de f . Empezamos calculando el subespacio de autovectores asociado al autovalor
λ1 = 2,
−1 −1 −3 x 0 −x −y −3z = 0
2 2 3 y = 0 ⇒ 2x +2y +3z = 0
0 0 1 z 0 z = 0
El sistema tiene por soluciones z = 0, x = −y, dando lugar al subespacio
S2 = L(λ2 = 3) = {(α, −2α − 3β, β) : α, β ∈ IR} =< (1, −2, 0), (0, −3, 1) > .
∀n ∈ IN , An = CDn C −1 .
Al ser D diagonal
λn1 0 · · · 0
0 λn2 · · · 0
Dn = .. .. .. , ∀n ∈ IN,
. . .
0 0 ··· λnp
EJEMPLO 2.4
Dada la matriz
2 1 0
A = 0 1 −1 ,
0 2 4
para saber si es diagonalizable comenzamos resolviendo la ecuación caracterı́stica
EJEMPLO 2.5
Sea f : IR3 → IR3 un endomorfismo, cuya matriz respecto a una base B = {⃗v1 , ⃗v2 , ⃗v3 }
es:
1 −1 0
A= 2 4 0 .
0 0 3
Utilizamos el ordenador para encontrar los valores y vectores propios de f . Empe-
zamos introduciendo la matriz
A:={{1,-1,0},{2,4,0},{0,0,3}}
A continuación calculamos los valores propios:
Eigenvalues[A]
{2, 3, 3}
Como no existen tres valores propios distintos, de entrada no podemos afirmar que
la matriz A sea diagonalizable. Para ello es necesario conocer los vectores propios
de f
Eigenvectors[A]
32 Tema 2 Diagonalización de matrices
{{-1,1,0},{-1,2,0},{0,0,1}} .
Para ver si forman una base de IR3 calculamos su determinante
det[{{-1,1,0},{-1,2,0},{0,0,1}}]
-1
Como podemos ver los tres vectores son independientes y, por tanto, existe una base
de R3 formada por vectores propios de f . En consecuencia, la matriz A será diago-
nalizable.
entonces el valor propio dominante es el λ3 = −6. Observemos que si, por ejemplo,
λ1 = 6, entonces no existe un valor propio que sea dominante.
EJEMPLO 2.6
2.5 Métodos numéricos 33
Sea la matriz, ( )
−2 2
A= ,
−12 8
y tomemos como vector inicial uno cualquiera ⃗x(0) = (−1, 2)T . Ahora, calculemos
⃗x(1) = A⃗x(0) = γ1 (1, 4.6666)T
Como puede observarse, los diferentes valores de ⃗x(k) tiende a un múltiplo del vector
(1, 3)T cuando k es suficientemente grande. Éste será el vector propio dominante.
⟨⃗u1 , A⃗u1 ⟩ ⟨(1, 3)T , A(1, 3)T ⟩ ⟨(1, 3)T , (4, 12)T ⟩ 40
λ1 ≈ = = = = 4.
⟨⃗u1 , ⃗u1 ⟩ ⟨(1, 3)T , (1, 3)T ⟩ ⟨(1, 3)T , (1, 3)T ⟩ 10
Naturalmente, esta aproximación puede ser mejorada tomando en lugar del ⃗x(1) el
⃗x(10), pero recordemos que todo lo dicho es válido si λ1 es dominante.
La expresión ⃗x(k) = Ak ⃗x(0) ≈ λk1 c1⃗u1 , también nos permite encontrar un valor
aproximado de λ1 . En efecto,
⃗xj (k)
⃗xj (k) ≈ λ1⃗xj (k − 1), j = 1, 2, · · · , n ⇒ λ1 ≈ .
⃗xj (k − 1)
EJEMPLO 2.7
Aplicaremos el método de las potencias con cambio de escala, para encontrar el valor
propio y el vector propio dominante de la matriz
( )
−3 2
A= .
1 −1
2.5 Métodos numéricos 35
Empezamos tomando como vector inicial uno que nos permita simplificar las opera-
ciones, por ejemplo ⃗x(0) = (1, 1)T .
⃗x(1) = A⃗x(0) = (1, 0)T
Una posible estimación del vector propio dominante puede ser ⃗u1 ≈ (1, −0.71428)T .
Por otro lado, por definición de vector y valor propio
( )( ) ( )
−1 2 1 λ1
A⃗u1 = λ1 ⃗u1 ⇒ = ,
1 −1 −0.71428 −0.71428λ1
y podemos estimar el valor de λ1 ≈ −2.42856
Una de las preguntas que está presente en todo el estudio que estamos desarrollando
es la siguiente: ¿cuántos pasos debemos usar en el método de las potencias para
obtener una buena estimación?. No existe una respuesta rotunda a esta pregunta,
pero el análisis del error relativo puede aportar una solución parcial al problema.
DEFINICIÓN 2.5.2 Si ã es una estimación de la cantidad a, entonces el error
relativo se define como:
a − ã
a .
El error en porcentaje en la estimación se define como
a − ã
a × 100 % .
Supongamos que inicialmente el error relativo que estamos dispuestos a tolerar en el
valor propio dominante es ϵ. Si λ̃1 (i) es la estimación de λ1 en la i-ésima iteración,
entonces el método se terminarı́a cuando se cumpla:
λ − λ̃ (i)
1 1
< ϵ.
λ1
No obstante, aplicar la fórmula anterior tiene el gran inconveniente de que el valor
exacto λ1 no es conocido. En este caso, suele sustituirse la expresión anterior por
esta otra:
λ̃ (i) − λ̃ (i − 1)
1 1
< ϵ,
λ̃1 (i)
donde la cantidad que aparece a la izquierda se la conoce con el nombre de error
relativo estimado, y si se multiplica por 100 %, se llama error en porcentaje estimado.
36 Tema 2 Diagonalización de matrices
2.5.3. Deflación
A través del método de las potencias podemos estimar el valor propio dominante y
un vector propio asociado. Existen diferentes técnicas para conocer el resto de los
valores propios. De todas ellas nosotros utilizaremos la conocida con el nombre de
deflación. El método está basado en el siguiente teorema cuya demostración puede
verse en cualquier texto de métodos numéricos del Algebra Lineal.
TEOREMA 2.5.3 Sean λ1 , λ2 , · · · , λn , los valores propios de la matriz A. Supon-
gamos que λ1 es el valor propio dominante, ⃗u1 su vector propio asociado, y ⃗v un
vector tal que ⟨⃗u1 , ⃗v ⟩ = 1. Sea B la matriz definida como
B = A − λ1 ⃗u1 ⃗v T .
Entonces los valores propios de la matriz B son 0, λ2 , · · · , λn .
Se conoce con el nombre de deflación al proceso de encontrar el resto de los valores
propios a partir del valor propio dominante. Para ello se aplica el método de las
potencias con escala a la nueva matriz B.
EJEMPLO 2.8
Aplicaremos el método de las potencias con cambio de escala y la deflación para
encontrar todos los valores propios de la matriz
3 −0.5 1
A = −2 5 −2 .
−1 1.5 1
Si tomamos como vector inicial ⃗x(0) = (1, 1, 1)T , entonces ⃗x(1) = A⃗x(0) = (3.5, 1, 1.5)T
que dividiendo por 1.5 nos queda x⃗′ (1) = (2.3333, 0.6666, 1)T . De forma similar
⃗x(2) = Ax⃗′ (1) = (7.6666, −3.3336, −0.3334)T , o bien
x⃗′ (2) = (−0.2527, 9.998, 1)T .
Los valores siguientes se encuentran representados en la tabla siguiente. También
aparecen los cocientes entre la primera de las componentes de una iteración y la
anterior γk .
MODELOS BASADOS EN
ECUACIONES Y SISTEMAS EN
DIFERENCIAS
3.1. Introducción
En ocasiones, al construir un modelo matemático interesa elegir una variable que
tome valores discretos. Ası́ ocurre, por ejemplo, con el tiempo, ya que es común
realizar mediciones regulares a la hora de controlar un experimento. Estos datos
constituyen un conjunto finito, o infinito numerable, de valores de la variable in-
dependiente. Para este tipo de modelos determinı́sticos discretos, las herramientas
matemáticas más adecuadas para analizarlos son las ecuaciones en diferencias y los
sistemas en diferencias. El presente tema es una breve introducción a su estudio. Co-
menzaremos con los conceptos y definiciones básicas y nos centraremos en el estudio
de las ecuaciones en diferencias lineales de primer y segundo orden con coeficientes
constantes, ası́ como en los sistemas de ecuaciones diferenciales de primer orden con
coeficientes constantes.
A lo largo del capı́tulo llamaremos t a la variable independiente, y supondremos que
sólo toma los valores enteros t = 0, 1, 2, · · · . Generalmente, t representa el número
de generaciones (años, trimestres, meses, dı́as, · · · ) que han transcurrido desde un
momento inicial t = 0. Del mismo modo, {y0 , y1 , y2 , · · · } es una sucesión, donde yt
corresponde a un valor concreto de t.
39
40 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
Hoy en dı́a todavı́a se conservan copias de algunos de sus libros, como por ejemplo,
Liber abbaci (1202), Practica geometriae (1220), Flos (1225), y Liber quadratorum.
Un problema que se encuentra en la tercera sección de su libro Liber abbaci llevó a
la introducción de los números de Fibonacci y a la sucesión que lleva su nombre, y
es la razón por la que aún hoy en dı́a es tan recordado.
Al final del segundo mes, la hembra produce un nuevo par, ası́ que ahora
tenemos dos parejas de conejos en el campo.
Al final del tercer mes, la hembra inicial produce un segundo par, haciendo
que en el campo tengamos tres pares de conejos.
3.1 Introducción 41
Al final del cuarto mes, la hembra original ha producido otro nuevo par y la
hembra nacida dos meses antes produce su primer par, con lo que tendremos
5 pares.
Y asi sucesivamente......
una sucesión, que se inicia con 1 y 1 donde cada otro término es la suma de los dos
anteriores, y que recibe el nombre de la sucesión de Fibonacci.
1/1 = 1, 2/1 = 2, 3/2 = 1.5, 5/3 = 1.666, 8/5 = 1.6, 13/8 = 1.625, 21/13 = 1.61538,
42 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
Es fácil ver lo que sucede si dibujamos estos resultados en un gráfico (Figura 3.3),
el cociente tiende a un valor particular, el cual recibe el nombre de número áureo,
tiene un valor aproximado de 1.61804, y se representa por la letra griega ϕ (Phi).
Si la sucesión se inicia con otros dos números cualesquiera, {3, 7, 10, 17, 27, 71, · · · }
el cociente entre un término y el anterior en la sucesión tiende al número áureo ϕ.
También podemos hacer aparecer el número áureo a través de una cuestión estética,
que aparece en el libro La Divina proporción de Luca Pacioli. Si consideramos un
segmento y preguntamos cuál es la división más agradable en dos partes del mismo,
algunas personas pensarán que el punto medio es el más adecuado, otras en cambio
pensarán que la tercera o cuarta parte. La respuesta correcta no es ninguna de
ellas, ya que la división correcta es la conocida con el nombre de razón áurea o
divina proporción. Si el segmento es de longitud 1, entonces el segmento mayor tiene
longitud 0.618.... A un segmento dividido de esta forma decimos que está dividido
en la sección áurea. Pensemos que si u es la longitud del segmento, se tiene que
verificar:
u v
=
u+v u
si llamamos ϕ = u/v, entonces
u+v u 1
= ⇒ 1+ =ϕ ⇒ ϕ + 1 = ϕ2 ⇒ ϕ2 − ϕ − 1 = 0
u v ϕ
Simplificando, obtenemos la ecuación de segundo grado que tiene como raı́z el valor
√
1+ 5
ϕ= = 1.6180339887...
2
que como sabemos es el número áureo.
caléndulas granuladas 13; la achicoria 21; mientras que algunas margaritas pueden
encontrarse con 34, 55 y hasta 89 pétalos. Algunas especies son muy precisas sobre
el número de pétalos que tienen, (por ejemplo el botón de oro), pero otras tienen
pétalos que están muy cerca de ellos, pero que su media es justo un número de
Fibonacci.
Se puede comprobar que algunas semillas se disponen según espirales, una de ellas
a la izquierda y otra a la derecha. Si contamos la espiral de la derecha al final del
dibujo, nos encontramos con que hay 34. ¿Cuántas existen hacia el otro lado? Puede
comprobarse que estos dos números son vecinos en la sucesión de Fibonacci.
La pregunta inmediata que nos planteamos es ¿por qué sucede esto? La respuesta
se encuentra detrás de un problema de optimización. Si deseamos apilar objetos de
la “mejor manera posible”, la respuesta será que dependerá de la forma que posean.
Si estos objetos son cuadrados, la respuesta correcta será la figura 3.5 de la izquierda.
Sin embargo, si estos objetos son redondos la mejor disposición es la conocida con el
nombre de hexagonal (figura 3.5 derecha). Por otro lado, muchas plantas que poseen
un tallo alto tienen adheridas las hojas según un esquema bastante interesante, En
efecto, se cumple la llamada ley de la Filotaxis (ciencia que estudia el ordenamiento
de los elementos de una planta), “para cada especie de plantas el ángulo que forman
dos hojas consecutivas, llamado ángulo de divergencia, es constante”.
Pero, ¿por qué la naturaleza no utiliza alguna de éstas disposiciones? La mayorı́a
de las semillas son redondas, entonces ¿por qué no se disponen en forma hexagonal
las semillas del girasol? La razón es que aunque la simetrı́a hexagonal es la mejor
manera de empaquetar semillas circulares, esto no responde a la pregunta de por
qué las hojas se distribuyen alrededor del tallo o como se empaquetan las semillas
del girasol (las cuales son circulares porque es la forma que encierra máxima área
con una determinada longitud), cuando están creciendo en tamaño .
La naturaleza utiliza el mismo patrón para colocar las semillas del girasol, los pétalos
alrededor del borde de una flor y las hojas alrededor del tallo. Además, sigue siendo
eficiente cuando la planta continúa creciendo. Los botánicos han demostrado que
las plantas crecen desde un grupo diminuto de células dispuestas en el extremo de
cualquier planta en crecimiento, llamado el meristemo. Hay un meristemo separado
44 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
al final de cada rama o varita y allı́ es donde se forman las nuevas células, que crecen
en tamaño, pero las nuevas células solo se forman en esos puntos de crecimiento.
Pero debemos dejar claro que no siempre las sucesión de Fibonacci se encuentra
en el número de pétalos de cualquier flor, o en el número de semillas de plantas
del tipo del girasol, aunque si están cercanos a algunos términos de la sucesión de
Fibonacci. De todos modos, en un trabajo publicado por Jean en 1992, de 12700
observaciones correspondientes a 650 especies de plantas diferentes, la sucesión de
Fibonacci está presente en un 90 % de todas aquellas donde sus elementos están
dispuestos en formas de espirales.
Observemos que si notamos a yn al valor del término que ocupa el lugar n de la
sucesión de Fibonacci, entonces podemos modelizarla mediante la siguiente ecuación:
yt = y(t) = 2t + C,
EJEMPLO 3.1
y2 − 4y1 + 300 = 0 ,
lo cual nos indica que debemos saber otra medida, por ejemplo y1 , para poder
encontrar el resto de los valores. En las próximas secciones aprenderemos a resolver
este tipo de ecuaciones, y volveremos sobre (3.1).
EJEMPLO 3.2
Supongamos que una determinada población de insectos con 100 individuos, duplica
su número en cada generación, y que además, 10 nuevos individuos se incorporan
en cada generación procedente de otro lugar. Vamos a construir una ecuación en
diferencias que modele a esta situación y posteriormente la resolveremos.
Del enunciado se deduce,
yt = 110 × 2t − 10 .
De esta manera, yt queda definida por la ley de recurrencia anterior. Puede com-
probarse que yt es solución de la ecuación pedida y cumple las condiciones iniciales.
Además, es la única solución, ya que si wt es otra solución que cumple
TEOREMA 3.3.3 Si yt1 , yt2 son dos soluciones de (3.4), entonces k1 yt1 +k2 yt2 , con
k1 y k2 constantes, sigue siendo solución de (3.4).
A continuación veremos las condiciones bajo las cuales la combinación lineal de dos
soluciones particulares de la ecuación homogénea dan lugar a su solución general.
3.3 Ecuaciones lineales de segundo orden 49
y = k1 yt1 + k2 yt2 ,
cualesquiera que sean los valores de β1 y β2 , por hipótesis del teorema, el sistema es
compatible determinado. Pero por el Teorema 3.3.2 existe una única solución de la
ecuación homogénea que puede ser escrita como yt = k1 yt1 + k2 yt2 , pues basta tomar
β1 = y0 y β2 = y1 , y calcular α1 y α2 . Para finalizar asignamos los siguientes valores,
k1 = α1 y k2 = α2 .
A dos soluciones yt1 y yt2 cumpliendo las hipótesis del Teorema 3.3.2 le daremos
el nombre de sistema fundamental de soluciones. Siguiendo un razonamiento
similar al realizado en el Teorema 3.3.2, podemos demostrar el siguiente resultado.
TEOREMA 3.3.6 Si ytp es una solución particular de
a yt+2 + b yt+1 + c yt = 0 ,
es decir,
a λ2 + b λ + c = 0. (3.7)
50 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
yt1 = λt , yt2 = t λt ,
λ1 = α + i β , λ2 = α − i β ,
entonces
yt1 = λt1 , yt2 = λt2 ,
forman un sistema fundamental de soluciones. En este último caso, podemos
escribir la solución general de la ecuación homogénea de la siguiente manera,
yt = k1 λt1 + k2 λt2 ,
Al formar λt1 = ρt (cos tθ + i sen tθ) y λt2 = ρt (cos tθ − i sen tθ) un sistema
fundamental de soluciones, también lo será cualquier combinación lineal de
ellas, en particular
1
2
(λt1 + λt2 ) = ρt cos tθ
1
2i
(λt1 − λt2 ) = ρt sen tθ ,
la solución general será entonces
yt = k1 ρt cos tθ + k2 ρt sen tθ .
EJEMPLO 3.3
3.3 Ecuaciones lineales de segundo orden 51
yt = k1 (−1)t + k2 t (−1)t , k 1 , k2 ∈ R .
EJEMPLO 3.4
yt+2 = yt+1 + yt ; y0 = 1 ; y1 = 1
yt+2 − yt+1 − yt = 0 ; y0 = 1 ; y1 = 1
o bien (
√ )t+1 ( √ )t+1
1 1+ 5 1− 5
yt = √ −
5 2 2
simplificando
( √ )( √ )t+1 ( √ )( √ )t+1
1+ 5
2
1+ 5
2 − 1− 5
2
1− 5
2
ϕ = lı́m ( √ )t+1 ( √ )t+1
t→∞ 1+ 5
2 − 1− 5
2
√ √ ( b )t+1
Si hacemos a = 1+2 5 y b = 1−2 5 . Entonces, b < a y b
a < 1. En consecuencia, a
tiende a cero cuando t tiende a infinito.
El lı́mite anterior puede escribirse como,
( )t+1 √
a at+1 − b bt+1 a − b ab 1+ 5
ϕ = lı́m
t→∞ at+1 − bt+1
= lı́m ( )
t→∞ 1 − b t+1
=a=
2
a
AB AC 1 x
ϕ= = ⇒ = ⇒ x2 + x − 1 = 0
AC CB x 1−x
Resolviendo la ecuación de segundo grado,
√ √
−1 + 5 −1 − 5
x1 = ; x2 =
2 2
Es decir,
(√ ) √ √
1 1 2 2 5−1 2 5+2 5+1
ϕ= = √ =√ = (√ ) (√ )= =
x −1+ 5 5−1 5+1 5+1 5 − 1 2
2
O en este otro:
1 1
x=1+ 1 ⇒x=1+ ⇒ x2 = x + 1 ⇒ x2 − x − 1 = 0
1+ 1 x
1+ 1+..
al resolverlo obtenemos:
√ √
1+ 5 1− 5
x1 = =ϕ ; x2 =
2 2
yt = k1 yt1 + k2 yt2 , k1 , k2 ∈ R ,
(3.12)
+ [k2 (t + 1) − 2
k2 (t)] yt+2 .
o bien,
[ 1 1
] [ 2 ]
k1 (t) ayt+2 + byt+1 + cyt1 + k2 (t) ayt+2 2
+ byt+1 + cyt2
−q(t)
k1 (t + 1) − k1 (t) = 1
ayt+1 (λ2 − λ1 )
(3.14)
q(t)
k2 (t + 1) − k2 (t) = 2
ayt+1 (λ2 − λ1 )
EJEMPLO 3.5
con las condiciones iniciales y0 = 100 y y1 = 110. Para ello, empezamos encontrando
las raı́ces de la ecuación caracterı́stica
λ2 − 3λ + 2 = 0 ⇒ λ1 = 1, λ2 = 2 .
yth = k1 + k2 2t .
k1 (1) = k1 (0) − 10
k1 (2) = k1 (1) − 10 = k1 (0) − 2 × 10
.. ..
. .
k1 (t) = k1 (0) − 10t
que dan como solución k1 (0) = 90, k2 (0) = 10. La ecuación de los efectivos de la
población es:
yt = 80 − 10 t + 10 × 2t+1 , t = 0, 1, 2, · · ·
3.3 Ecuaciones lineales de segundo orden 57
Segundo método.
Para encontrar la solución general de una ecuación lineal completa de segundo orden
nos fijaremos en el término independiente q(t), y según sea, procederemos de una
manera u otra. Los casos más usuales que suelen presentarse son:
EJEMPLO 3.6
(a) 2yt+2 − yt = 2t
(b) yt+2 − 4yt+1 + 4yt = 3t2
(c) yt+2 − 4yt+1 + 4yt = 2t + 1.
1
2λ2 − 1 = 0 ⇒ λ = ±√ ,
2
(b) Para el segundo de los casos, es inmediato comprobar que λ = 2 es una raı́z
doble de la ecuación caracterı́stica. En consecuencia, la solución general de la
ecuación homogénea es
yth = k1 2t + k2 t 2t .
Al ser el término independiente un polinomio de segundo grado y el 1 no es
raı́z del polinomio caracterı́stico, probamos con una solución particular del
tipo
ytp = at2 + bt + c ,
llevando este valor en la ecuación en diferencias propuesta, e identificando
coeficientes, se obtiene a = 3, b = 12 y c = 24.
La solución buscada es
yt = k1 2t + k2 t 2t + 3t2 + 12t + 24 .
(c) Al coincidir la ecuación homogénea con la del caso anterior, lo único que tene-
mos que hacer es encontrar una solución particular. Para ello, buscamos una
solución particular de
yt+2 − 4yt+1 + 4yt = 2t ,
y otra de
yt+2 − 4yt+1 + 4yt = 1 .
Para la primera de ellas, al ser λ = 2 una raı́z doble, ensayamos con la solución
yt1 = kt2 2t .
yt = k1 2t + k2 t 2t + 81 t2 2t + 1 .
una breve introducción a las ecuaciones y a los sistemas en diferencias. Por este
motivo, solo abordaremos aquellos sistemas de ecuaciones en diferencias lineales y
de primer orden. Además, este tipo de sistemas son los que con más frecuencia se
presentan en las aplicaciones biológicas.
De entre este tipo de sistemas, el caso más elemental (aunque para casos más gene-
rales, el procedimiento a seguir es similar) consiste en dos ecuaciones y dos variables
{ 1
yt+1 = a11 yt1 + a12 yt2 + f1 (t)
2
yt+1 = a21 yt1 + a22 yt2 + f2 (t)
La clave para resolver este tipo de sistemas, es intentar expresarlo como una ecuación
en diferencias lineal de segundo orden con coeficientes constantes. En efecto, de la
primera de las ecuaciones
1 1 2
yt+2 = a11 yt+1 + a12 yt+1 + f1 (t + 1) , (3.15)
en la que sólo aparece un término a12 a22 yt2 , en el que no intervenga la función yt1 .
Despejando de la primera de las ecuaciones del sistema
1
a12 yt2 = yt+1 − a11 yt1 − f1 (t) ,
sustituyendo
1
yt+2 1
= a11 yt+1 1
+ a12 a21 yt1 + a22 (yt+1 − a11 yt1 − f1 (t)) + a12 f2 (t) + f1 (t + 1) ,
EJEMPLO 3.7
que es una ecuación en diferencias lineal de segundo orden con coeficientes cons-
tantes, que puede ser escrita
a − 5a + 5a = −4 ⇒ a = −4 ,
yt = −xt+1 + 3xt + 1
( √ )t+1 ( √ )t+1 ( √ )t
5+ 5 5− 5 5+ 5
= −k1 − k2 + 4 + 3k1
2 2 2
( √ )t
5− 5
+3k2 − 12 + 1 .
2
( √ ) ( √ )t ( √ ) ( √ )t
1− 5 5+ 5 1+ 5 5− 5
yt = k1 + k2 − 7.
2 2 2 2
( √ )t ( √ )t
√ 5+ 5 √ 5− 5
yt = (151 − 61 5) + (166 − 64 5) −7
2 2
62 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
EJERCICIOS PROPUESTOS
EJERCICIO 1
yt+2 − 2yt+1 + yt = 8
3.4 Sistemas de ecuaciones en diferencias 63
siendo x0 = 2 y x1 = 5.
3.5.1. Introducción
La teorı́a de sistemas dinámicos es una rama de las Matemáticas que se ocupa del
estudio del movimiento, y proporciona un lenguaje común para la Matemática, la
Biologı́a, la Economı́a, la Fı́sica, la Historia y la Literatura.
En la teorı́a de los sistemas dinámicos, un sistema se define como una colección
de elementos que continuamente interactuan para formar un conjunto unificado. A
las relaciones internas y las conexiones entre los componentes de un sistema se les
llama la estructura del sistema. Un ejemplo de un sistema es un ecosistema. La
estructura de un ecosistema está definida por las relaciones entre la población ani-
mal, nacimientos y muertes, cantidad de comida, y otras variables especı́ficas para
un ecosistema particular.
El término dinámico hace referencia al cambio a lo largo del tiempo. Si algo es
dinámico, es porque se está modificando constantemente. Un sistema dinámico es
aquel en el cual las variables se modifican para producir cambios a lo largo del tiem-
po. La manera por la cual los elementos o las variables de un sistema cambian con el
tiempo se denomina comportamiento del sistema. En el ejemplo del ecosistema,
el comportamiento está descrito por la dinámica que se produce como consecuencia
de los nacimientos y las muertes de la población. El comportamiento está expuesto a
la influencia de comida adicional, los depredadores, y al medio ambiente, los cuales
son todos elementos del sistema.
Los sistemas dinámicos también pueden usarse para analizar, como pequeños cam-
bios en una parte del sistema, pueden afectar al comportamiento del sistema com-
pleto. Si nos referimos al ejemplo del ecosistema, podemos analizar el impacto de la
sequı́a en el ecosistema o analizar el impacto de la eliminación de una determinada
especie animal en el comportamiento del ecosistema completo.
En relación con los sistemas dinámicos discretos, fue H. Poincaré en 1899 el primero
en utilizarlos al intentar simplificar un modelo continuo, pero ha sido en la década
de los cincuenta donde han sido estudiados y aplicados en problemas muy diversos.
En 1976 R. May, analizando el comportamiento de las ecuaciones en diferencias en
el modelo que lleva su nombre, observó que aún para el caso determinista, el modelo
podı́a presentar comportamientos “muy complicados”. En 1963 el meteorólogo E.
Lorentz descubre el caos matemático en un sistema dinámico continuo, presentando
a la comunidad cientı́fica el atractor que lleva su nombre. Poco después, en 1973,
M. Heron estudia el caso discreto, descubriendo un tipo de atractor muy parecido
al de Lorentz. Dos años después, Feigenbaum presentó por primera vez el diagra-
ma de bifurcación correspondiente al modelo logı́stico. Actualmente la teorı́a de las
bifurcaciones es un campo donde se investiga intensamente.
3.5 Sistemas dinámicos discretos 65
2.- La parábola logı́stica de May. En 1976 el biólogo Robert May formuló otra
ecuación para estudiar el crecimiento de una población de insectos en un eco-
sistema aislado, diferente de la de Malthus. May tuvo en cuenta los efectos
de saturación del ecosistema. Cuando la población se acerca al máximo po-
sible que el medio ambiente puede sustentar, entonces el parámetro α debe
disminuir, lo que equivale a considerar este parámetro función del número de
individuos. Con ello se llega a una ecuación de la forma
3.- Modelo matricial. Supongamos que una especie de aves que vive muchos
años, resulta capaz de reproducirse a partir del segundo año de vida y que,
por término medio, cada pareja de aves en edad reproductora crı́a anualmente
una nidada de la que sobreviven dos crı́as, una de cada sexo. Se supone que
66 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
a partir del segundo año todas las aves han emparejado. Se suelta una pareja
de aves en una región sin depredadores. ¿Cuál es la ley de evolución para la
población de aves¿.
Ante un problema de esta naturaleza, el primer paso consiste en seleccionar
las variables. Debido a las diferentes condiciones de reproducción, conviene
considerar dos segmentos en la población de aves: las de un año y las de dos o
más. Tomamos como variable x1 (k) el número de parejas adultas en el perı́odo
k. Debido a que existe siempre el mismo número de machos que de hembras,
la población de aves de un año puede también contarse por el número x2 (k)
de parejas que pueden formarse entre ellas. Se tienen entonces las siguientes
relaciones: {
x1 (k + 1) = x1 (k) + x2 (k)
x2 (k + 1) = x1 (k)
que pueden escribirse en forma matricial como
( ) ( )( )
x1 (k + 1) 1 1 x1 (k)
=
x2 (k + 1) 1 0 x2 (k)
xk+1 = f (xk ) , k = 0, 1, 2, · · ·
y en general
xk = f k (x0 ) .
La expresión f k (x), es la solución general o flujo de los sistemas dinámicos
discretos. Permite conocer el estado del sistema en cualquier instante a partir de
su posición inicial. El conjunto de valores
EJEMPLO 3.8
Puede comprobarse fácilmente que la ecuación de Malthus admite por solución gene-
ral la expresión
Φx (k) = αk x
El comportamiento de esta expresión es sencillo de comprender. Si x > 0 entonces
0 si 0<α<1
lı́m αk x = x si α=1
k→∞
+∞ si α>1
Menos sencillo resulta el problema de encontrar una fórmula explı́cita para el proble-
ma de las aves. Se obtiene
( ) ( )k ( )
x1 (k) 1 1 x1 (0)
=
x2 (k) 1 0 x2 (0)
Si calculamos las sucesivas potencias de la matriz, nos aparece un hecho curioso como
es la aparición en estas matrices de los célebres números de Fibonacci 1, 1, 2, 3, 5, 8,....
Supongamos una árida isla cerca de la costa de un rico continente. Estamos interesa-
dos en una especie particular de pájaros que anidan en estas islas. Desgraciadamente
el habitat de la isla es muy desfavorable ya que si los pájaros estuvieran aislados su
población disminuirı́a un 20 % cada año. Esta situación podemos reflejarla utilizando
el modelo de Malthus (exponencial)
xk+1 = 0.80xk , k = 0, 1, 2, · · · ,
donde xk es la población de pájaros en el tiempo k. Hay una colonia de pájaros en
el continente y cada año 1000 pájaros emigran a nuestra isla. Entonces, el cambio
de población en la isla puede ser descrito por el modelo
xk+1 = 0.80xk + 1000 , k = 0, 1, 2, · · · .
Observemos que el modelo es lineal en el sentido de que la función f (x) = 0.80x +
1000 representa a una lı́nea recta.
Ahora descubriremos y probaremos un teorema sobre sistemas dinámicos lineales
discretos, que corresponden al tipo
xk+1 = mxk + b , k = 0, 1, 2, · · · ,
donde m y b son constantes. Recordemos que los sistemas dinámicos discretos están
descritos por una ecuación de la forma
xk+1 = f (xk ) , k = 0, 1, 2, · · · .
En el caso particular de sistemas dinámicos discretos lineales, la función f es del tipo
f (x) = mx + b, y estamos interesados en puntos de equilibrio del sistema dinámico
discreto. Aquellos puntos x tales que f (x) = x. Estos puntos se llaman de equilibrio
porque si un término es uno de estos puntos, cada sucesión de términos siguientes
permanece en el mismo punto. De esta manera decimos que el sistema se encuentra
en equilibrio. Es inmediato comprobar el siguiente resultado.
RESULTADO 3.5.1 Si m no vale 1 entonces hay un único punto de equilibrio
b
x∗ =
1−m
Los modelos dinámicos discretos pueden comportarse de manera sorprendente. Al-
gunas veces una sucesión obtenida del sistema dinámico lineal discreto tiende direc-
tamente al punto de equilibrio. En otras ocasiones saltan alrededor de él, con saltos
cada vez más pequeños hasta tender al punto de equilibrio. O por el contrario los
saltos son cada vez más grandes y no tienden al punto de equilibrio.
Nuestro objetivo es formular y probar un teorema que nos determine cuando ocu-
rre cada una de estas clases de comportamiento. Comenzamos la construcción del
diagrama de Cobweb dibujando las gráficas
f (x) = mx + b , g(x) = x
3.5 Sistemas dinámicos discretos 71
1000
800
600
400
200
800
700
600
500
400
300
5 10 15 20
EJEMPLO 3.9
Vamos a calcular y dibujar x2 , x3 , · · · , para el modelo xk+1 = 0.80xk + 1000 y el
valor inicial x0 = 500.
El modelo anterior podemos escribirlo como
xk+1 = f (xk ) , k = 0, 1, 2, · · · ,
72 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
donde f (x) = 0.80x + 1000. Esta es una buena manera de representar a nuestro
modelo porque la función f (x) nos describe como la población, en cada año está de-
terminada por la población en el año anterior.
Los dos gráficos f (x) = 0.8x + 1000 y g(x) = x se cortan en el punto x∗ = 5000.
Este punto se llama punto de equilibrio ya que la población en los próximos años
será la misma que la población actual.
6000
5000
4000
3000
2000
1000
Hemos tenido ocasión de ver que en los sistemas dinámicos lineales discretos, el
punto de equilibrio
b
x∗ = ,
1−m
3.5 Sistemas dinámicos discretos 73
xk = f (xk−1 ) , f (x) = mx + b
lı́m |xk | = ∞ .
k→∞
x2 = mx1 + b
x3 = mx2 + b = m2 x1 + b(m + 1)
x4 = mx3 + b = m(m2 x1 + b(1 + m)) + b = m3 x1 + b(1 + m + m2 )
···
∑
k−2
xk = m x1 + b(1 + m + m + · · · + m ) = m x1 + b
k−1 2 k−2 k−1
mj .
j=0
lı́m mk−1 x1 = 0 .
k→∞
por ser la suma de los infinitos términos de una progresión geométrica de razón
|m| < 1. Por lo tanto,
b
lı́m xk = = x∗ .
k→∞ 1−m
∑
k−2
Por un razonamiento similar, si |m| > 1 se cumple que m k−1
x1 y mj no están
j=0
acotados cuando k → ∞.
74 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
lı́m xk = x∗ .
k→∞
EJEMPLO 3.10
f (x) = αx(1 − x) .
Por tanto, tendremos que resolver la ecuación f (x) = x, que tiene por solucio-
nes,
1
x∗ = 0, x∗ = 1 − .
α
3.5 Sistemas dinámicos discretos 75
Para clasificar estos puntos de equilibrio, tenemos que hacer uso del Teorema
3.5.3. La derivada de la función f (x) vale f ′ (x) = α − 2αx. Por tanto, el
primero de los puntos es asintóticamente estable si
( )
′ 1
f 1 − = |2 − α| < 1 ⇒ 1 < α < 3.
α
0.8
1
0.7
0.8
0.6
0.6
0.4
0.5
0.2
1 0.8
0.8
0.75
0.7
0.6
0.65
0.4
0.6
0.2 0.55
5 10 15 20
0.2 0.4 0.6 0.8 1
0.8
0.8
0.7
0.6
0.6
0.4
0.5
0.2
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
Tendremos ocasión de volver sobre este ejercicio cuando estudiemos los sistemas
caóticos.
mı́nimo entero k tal que f k (x∗ ) = x∗ , se llama orden del punto periódico.
En tal caso la órbita
3.5.7. Estabilidad
Un punto de equilibrio de un sistema dinámico representa un estado fijo del sistema.
Ahora bien, no todos los estados de equilibrio tienen la misma naturaleza.
EJEMPLO 3.11
0.4
1
0.2
0.5
5 10 15 20
-1 -0.5 0.5 1
-0.2
-0.5
-0.4
-1
Antes de la aparición de esta nueva teorı́a, sólo habı́a dos tipos de comportamientos
conocidos para un sistema dinámico: un estado fijo, donde los variables nunca cam-
bian, y el comportamiento periódico, donde el sistema está en un “circuito cerrado”
y se repite infinitamente. Las soluciones del sistema de Lorentz son definitivamente
ordenadas (siguen una espiral). Nunca se paran en un punto, ni se repiten, ni son
periódicas. A su representación gráfica se la conoce con el nombre Atractor de Lo-
rentz1 . Estas gráficas deben cumplir otra condición: no puede cortarse a sı́ misma
ya que, si ası́ fuese, habrı́a dos curvas diferentes a partir de ese punto de corte, lo
cual significarı́a dos realidades simultáneas y diferentes.
Una curva de estas caracterı́sticas no puede estar contenida en un plano, y por su-
puesto su dimensión es fraccionaria. Este tipo de atractores reciben el nombre de
atractores extraños, ya que su representación gráfica es un fractal. Queremos
insistir en la idea fundamental que encierra el concepto de atractor, como es la
siguiente: mientras es casi imposible predecir exactamente el estado futuro de un
sistema, es posible, y aún más, muchas veces es fácil modelar el comportamiento
general del sistema.
A continuación, resumimos algunos de los rasgos caracterı́sticos de los sistemas caóti-
cos.
Son muy sensitivos a las condiciones iniciales. Un cambio muy pequeño en los
datos iniciales dan lugar a resultados totalmente diferentes.
1
El atractor es la región del espacio hacia la cual convergen las trayectorias posibles dentro de
un sistema.
82 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
Parecen un desorden, o hechos al azar, pero no lo son, hay reglas que determi-
nan su comportamiento. Los sistemas hechos al azar no son caóticos.
Supongamos una órbita periódica del sistema dinámico discreto xt+1 = f (xt ). Si
tomamos como semilla, o condición inicial x0 , podemos elegir otro punto inicial
muy próximo al anterior y0 = x0 + δ0 , siendo |δ0 | la distancia (pequeña) entre estos
dos puntos de partida. A continuación iteramos la función n veces obteniéndose los
valores xn y yn . Podemos definir δn = yn − xn , esto es, la separación de los dos
puntos iniciales después de n itereaciones. Si ocurre que |δn | ≃ |δ0 |enλ , entonces este
valor de λ se le conoce con el nombre de exponente de Lyapunov, y despejandolo
de la expresión anterior nos permite dar su definición,
( )
1 δn
λ = ln
n δ0
Puede probarse que el exponente de Lyapunov de una órbita puede calcularse como,
1∑
n−1
λ = lı́m ln|f ′ (xk )|
n→∞ n
k=0
Tabla 3.1
84 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
La Tabla 3.1 muestra algunos de estos valores, y además los cocientes entre la
amplitud de un intervalo y el inmediatamente anterior, por ejemplo (3.5441 −
3.4495)/(3.4495 − 3). Lo llamativo de este hecho, es que los cocientes tienden al
número transcendente:
4.669201609110299067185532038204662016...
que se conoce con el nombre de constante de Feigembaun. Es un problema abierto
el estudiar la importancia que este número juega en la naturaleza. Se piensa que
puede tener un protagonismo similar al número e.
Las aplicaciones de la teorı́a del caos son múltiples y en campos muy diversos, en
Biologı́a, en Economı́a, en Medicina,... etc. Hasta ahora parecı́a que al estallar el caos
no serı́amos capaces de hacer nada, por ejemplo, si el avión empieza a moverse de
una manera extraña pensamos que la catástrofe es inevitable; o bien, si el corazón
empieza a latir rápidamente y sin ayuda inmediata puede ocurrir lo peor. En los
últimos años, en el campo de la Medicina, las investigaciones actuales, nos ofrecen
esperanzas de “domesticar” el caos. Edward Ott, Ceslo Grebogi (fı́sicos) y James A.
Yorke (matemático) han elaborado un algoritmo matemático que permite convertir
un determinado tipo de caos en un proceso periódico sencillo. La idea que encierra
el algoritmo, es que no es necesario comprender todo el movimiento caótico para
poderlo regular. Lo que tenemos que hacer es “mirar” continuamente a que dirección
tiende el proceso, y realizar perturbaciones pequeñas para volver a dirigirlo en el
camino deseado. Naturalmente aquı́ no se termina de vigilar el sistema, porque
después el caos aparecerá de nuevo. Por otro lado, el profesor A. Garfinkel de la
Universidad de California, ha conseguido transformar el movimiento caótico de un
corazón sacado a un conejo en un movimiento regular.
EJEMPLO 3.12
Supongamos
Cuando r > 1/4 las raı́ces son dos números complejos conjugados: λ1 = ρeiθ y
λ2 = ρe−iθ , donde
√ √
ρ = 14 + 4r−14 = r
√ √
4r−1/2
θ = arctan 1/2 = arctan 4r − 1
1.- Si 0 < r < 1/4, entonces λ1 y λ2 son dos números reales comprendidos es-
trictamente entre cero y uno. En consecuencia, si k tiende a infinito λk1 y λk2
tienden a cero.
Conclusión: El punto de equilibrio x∗ = 1 será estable. Además, después de
una pequeña perturbación la solución tiende de forma monótona al punto de
equilibrio. Puede verse gráficamente en la Figura 3.22 (izquierda).
86 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
2.- Si 1/4 < r, entonces λ1 y λ2 son dos números complejos conjugados, con
λ1 λ2 = |λ1 |2 = ρ2 = r. Además, si 1/4 < r < 1, entonces para que la solución
vk = Aλk1 + B λ̄k1 ,
sea un número real, debe suceder que B = Ā . Supongamos entonces que
A = α eiγ y B = α e−iγ , llevando estos valores en la solución
vk = Aλk1 + B λ̄k1 = α eiγ ρk eiθk + α e−iγ ρk e−iθk =
( )
αρk ei(γ+θk) + e−i(γ+θk) = 2αρk cos (γ + θk) ,
y por lo tanto,
√ cuando el parámetro r tiende a uno, entonces el ángulo θ tiende
hacia arctan 3 = π/3.
(b1) Cuando r sea mayor que uno, entonces |λ1 | > 1 y vk crece indefinidamente
al tender k hacia infinito. En consecuencia x∗ = 1 será un punto de equilibrio
inestable.
EJERCICIOS PROPUESTOS
EJERCICIO 2
xt+1 = (1 + r)xt − µ , t = 0, 1, 2, · · ·
xt+1 = 3xt (1 − xt ) , t = 0, 1, 2, 3, · · · .
λ 1−b
f (Nt ) = Nt+1 = N , α, b, λ > 0 , t = 0, 1, 2 · · · (3.25)
α t
10.- Encontrar y clasificar los puntos de equilibrio del siguiente modelo dis-
creto:
kxt
xt+1 = , k > a > 0,
a + xt
donde xt representa al número de individuos de una población en el tiem-
po t ¿‘Cuál será el comportamiento de la población a largo plazo, si
k = 30, a = 10 y x0 = 15 individuos?
1( )
Nt+1 = f (Nt ) = −Nt3 + Nt2 + 4Nt − 1
3
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.
Con estas rotundas palabras inicia el profesor Michael Barnsley su famoso libro
sobre fractales. En efecto, una vez que conocemos esta nueva geometrı́a es difı́cil
contemplar la naturaleza, y el mundo que nos rodea, con los mismos ojos. La geo-
metrı́a fractal, no sólo es interesante por la belleza de las imágenes que crea, sino que
en la actualidad es difı́cil encontrar una rama de conocimiento donde los fractales
no estén presentes.
El significado de la palabra Matemática es “lo que se aprende”, pero en el fondo de
su filosofı́a se encuentra la búsqueda de patrones. En un principio, nuestros antepasa-
dos observaban los fenómenos naturales, como los eclipses, con verdadero terror. No
podı́an entender cómo repentinamente el sol desaparecı́a del cielo y la más terrible
oscuridad se cernı́a sobre ellos. Tenı́an que buscar explicaciones y al no encontrarlas,
le atribuyeron a estos fenómenos naturales un carácter divino. Sin embargo, el po-
der del razonamiento de nuestra especie se impuso y algunos matemáticos griegos se
dieron cuenta que la naturaleza estaba regida por leyes que eran necesarias conocer.
La búsqueda de estos patrones ha sido lenta, larga y complicada. Se inició con Hi-
parco, continuó con el modelo heliocéntrico de Aristóteles y Ptolomeo, fue revisado
y mejorado por Copérnico y Galileo, posteriormente formalizado por Newton y mo-
dificado por Einstein. Lo importante es comprender que no es necesario recurrir a
los dioses para entender los fenómenos naturales. Nuestro universo, se rige por unas
leyes muy precisas y los cientı́ficos, los matemáticos entre ellos, se ocupan de descu-
brirlas.
En la larga búsqueda de estos patrones, los matemáticos griegos encontraron una
geometrı́a que representaba perfectamente al orden y lo regular . Apoyándose en
la herencia de las matemáticas egipcias y babilonias, los griegos construyeron los
cimientos de nuestras matemáticas actuales, con la invención de unas matemáticas
basadas en definiciones, axiomas, teoremas y demostraciones.
Desconocemos gran parte de la vida de Euclides, aunque la repercusión de su obra,
en especial Los Elementos, ha sido crucial para la humanidad. Con la invención de
la imprenta el libro de texto más reeditado ha sido la Biblia y a continuación Los
Elementos, con más de mil ediciones. Los Elementos son 13 libros donde se reco-
ge, estructura y organiza todo el saber conocido, especialmente en Geometrı́a y en
Teorı́a de Números.
3.6 Ejemplo de sistema dinámico: la Geometrı́a Fractal 91
Uno de sus rasgos más significativo de la geometrı́a euclı́dea es que a estos objetos se
le puede asignar un número, conocido como su dimensión topológica, que sólo puede
tomar valores enteros. De esta manera, un punto no tiene tamaño, su dimensión es
cero, una recta tiene dimensión uno, un plano dimensión dos, y ası́ sucesivamente.
Otra caracterı́stica de los objetos que representa la geometrı́a euclı́dea es que cuando
se amplı́an tienden a perder la forma que tienen inicialmente. Si se aumenta sucesi-
vamente una parte de una circunferencia al final adopta la forma de una recta y no
de una circunferencia.
Sin embargo, Galileo no estaba en lo cierto. Existen otros objetos y situaciones que
no son regulares o que son caóticos, ¿qué pasa con ellos¿, ¿cómo pueden ser repre-
sentados y manipulados¿Empecemos fijándonos en los sistemas caóticos. El caos es
una apasionante teorı́a, surgida en los últimos años, que ocupa a un elevado número
de cientı́ficos, con la esperanza de simplificar fenómenos complejos. Hasta finales del
siglo XIX se creı́a que todo nuestro entorno estaba regido por leyes determinista.
Los cientı́ficos estaban convencidos de que nuestro mundo era predecible, estable
92 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
“· · · no hemos de creer que las orillas del mar sean realmente defor-
mes por no tener la forma de un baluarte regular; que las montañas no
son exactamente como conos o pirámides, ni las estrellas están situadas
desmañadamente por no estar a una distancia uniforme · · · “
lápiz del papel (era continua) pero que estaba llena de irregularidades (no tiene
derivada en ninguno de sus puntos). Esto era algo que aborrecı́an el resto de sus
compañeros. Sin embargo, a pesar de lo extraño de esta función, lo sorprendente
es que su gráfica representa muy bien a fenómenos de la vida cotidiana, como por
ejemplo las cotizaciones en bolsa de las acciones de un banco español.
Estas nuevas extrañas estructuras fueron bautizadas por el gran matemático francés
Poincaré como “monstruos matemáticos”, curvas que estaban en contra de toda
intuición.
La comunidad matemática focalizó todo su esfuerzo en estudiar estas funciones llenas
de irregularidades, hasta tal punto que el famoso matemático Hermite dirigió una
carta a su buen amigo Stieljes diciendo, “ · · · abandono con espanto y horror esta
3.6 Ejemplo de sistema dinámico: la Geometrı́a Fractal 95
lamentable plaga de funciones sin derivada · · · ”. En este mismo sentido Poincaré co-
mentaba “ · · · cuando uno inventaba una nueva función, estaba en vista de alguna
meta conveniente; ahora se inventan intencionalmente para poner los defectos de los
razonamientos de nuestros padres”.
Hay que esperar hasta mediados del siglo XX con la aparición de Norbert Wiener
cuando, analizando el movimiento bowniano de las partı́culas, construye un modelo
matemático totalmente irregular. En este trabajo aparece por primera vez la palabra
caos. En la introducción del trabajo comenta:
“Las nubes no son esferas, las montañas no son conos, las costas no son
cı́rculos, como la corteza de un árbol no es plana ni un rayo viaja en
lı́nea recta, · · · La naturaleza no solamente exhibe un grado mayor sino
también un nivel diferente de complejidad ”
se unen los puntos medios de sus tres lados y se elimina el triángulo central. Tenemos
ahora tres triángulos y en cada uno de ellos se aplica el procedimiento descrito en el
triángulo inicial. Se reitera este proceso geométrico hasta el infinito, y el conjunto
teórico que se obtiene es el triángulo de Sierpinski.
Por supuesto que el mismo razonamiento puede aplicarse a un cuadrado, o a un
cubo, en lugar del triángulo. En este caso, el conjunto obtenido se llama la alfombra
de Sierpinski, y tiene la sorprendente particularidad de que el área (o el volumen en
el caso del cubo) de esta alfombra es cero. Es decir, una maravillosa alfombra que
nunca estarı́a sucia ya que no podrı́a acumular polvo, al carecer de área.
98 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
3.6.3. Autosemejanza
“ · · · un naturalista observa que una pulga tenı́a pequeñas pulgas en ella
y éstas a su vez otras más pequeñas que las mordı́an. Ası́ se procede hasta
el infinito.
De esta manera tan sencilla el escritor Jonathan Swift describe el concepto de au-
tosimilitud. La autosemejanza o autosimilitud es una de las caracterı́sticas de los
objetos fractales, por la cual una de estas imágenes puede ser descompuesta en tro-
zos más pequeños cada uno de los cuales es idéntico al objeto completo.
Existe una definición matemática más precisa que puede simplificarse diciendo que
la autosemejanza es una transformación que conserva la forma y las proporciones,
pero que puede variar en tamaño, posición y orientación. Lo importante es que ca-
da uno de los trozos se obtenga de la estructura completa a través de las mismas
transformaciones.
habı́a sido estudiada por muchos cientı́ficos en siglos anteriores, Leibniz entre ellos,
cuando propuso esta hermosa idea sobre la estructura del universo. Una gota de
agua contenı́a todo un universo, que a su vez contenı́a gotas de agua más pequeñas;
cada una de estas pequeñas gotas encerraba a su vez un universo que tenı́a en su
interior otras gotas de agua, todavı́a más pequeñas y ası́ sucesivamente.
En la naturaleza es fácil encontrarse con esta propiedad, por ejemplo la coliflor o el
brócoli la posee. Si tomamos un trozo podemos observar que ésta ramita es idéntica
a la coliflor completa. Esta misma experiencia podemos hacerla con la hoja de un
helecho, o en mucho más diferente como es el sistema vascular.
{z0 , z1 , z2 , z3 , · · · , zn , · · · , }
que recibe el nombre de órbita. Si todos los elementos de este conjunto de puntos
permanecen a menos de dos unidades del origen de coordenadas, entonces el punto
c, pertenecerá al conjunto de Mandelbrot, si la distancia es mayor de dos unidades c
no formará parte del conjunto. Es decir, el conjunto de Mandelbrot es el conjunto de
puntos cuyas órbitas generadas por la fórmula anterior nunca escapan de un cı́rculo
centrado en el origen de radio dos.
La genial idea de Mandelbrot, que pudo hacer porque la tecnologı́a estaba a punto,
fue la de pintar en color negro los puntos del conjunto de Mandelbrot y en blanco
los que estaban fuera. Después de un largo tiempo de espera el resultado obtenido
es el que puede apreciarse en la figura.
102 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
Por supuesto que una vez conocido el proceso, éste puede ser mejorado. Cuando los
números no pertenezcan al conjunto de Mandelbrot (el exterior del fractal) entonces
se le asignará un algoritmo de color, consistente en asociar colores diferentes a cada
punto, de acuerdo a ”la rapidez”del incremento de la órbita. Por ejemplo, si la
sucesión se incrementa lentamente, entonces la semilla se puede pintar de color azul,
por el contrario si crece más rápidamente, tendrá un color amarillo, naranja o verde
en función de la velocidad de este crecimiento. Los primeros que utilizaron esta
técnica fueron Peiten y Ritchter.
A mediados de los años 80, A. Douady y J.H. Hubbard probaron que este conjunto
y su complementario son conexos. Es decir, de un solo trozo, aunque al ampliar la
frontera del conjunto de la sensación de lo contrario. Debido a que siempre hemos
estado acostumbrado a lo regular a las estructuras representadas por la geometrı́a
2
hermosas imágenes. De entre ellos, uno de los más apropiados, por su sencillez de
manejo y potencia de resultados, es Ultra Fractal. El programa tiene la ventaja de
que no solo se puede construir bonitos fractales, sino que es posible realizar espec-
taculares animaciones de los mismos, ası́ como fondos fractales con movimiento.
Precisamente la belleza de las imágenes obtenidas junto con la sencillez de su crea-
ción ha acercado al mundo de las Matemáticas a un gran número de personas. En
boca de Mandelbrot “ayudan a tender puentes en el abismo que separa las cuestiones
matemáticas de la gente de la calle”.
Para poder medir el perı́metro de la costa se toma, en primer lugar, una unidad de
longitud. En este caso, el perı́metro se estima a partir del número de veces que es
necesario llevar el segmento para cubrir toda la costa. Es evidente, que al disminuir
el tamaño del segmento unidad, la longitud de la costa aumenta. Los objetos irre-
gulares carecen de medida exacta, puesto que la longitud dependerá de la unidad de
medida elegida.
Al no poderse caracterizar la frontera de un fractal por una longitud, es necesario
centrarse en un concepto diferente como es el de dimensión, procedente del latı́n di-
mensio, que significa medida. La idea es buscar un parámetro que, en cierto modo,
mida el grado de irregularidad del fractal, y eso será la dimensión fractal.
Desgraciadamente no hay una única definición de dimensión. Existe una definición
intuitiva de la dimensión euclı́dea que coincide con el número de valores reales que
son necesarios para situar a un punto en el espacio donde nos encontremos traba-
jando. Una lı́nea tiene dimensión uno, puesto que sólo se necesita un número real
para localizar cualquier punto de la lı́nea. Pensemos, por ejemplo, que cuando desea-
mos dar nuestra situación en una autopista aportamos sólo un número, su punto
kilométrico, aunque esta carretera no sea una lı́nea recta.
Otra definición intuitiva de dimensión euclı́dea hace referencia al número de direc-
ciones perpendiculares diferentes que se pueden tomar. En nuestro mundo cotidiano
contamos con tres direcciones: izquierda-derecha, adelante-atrás, y arriba-abajo; por
este motivo decimos que el espacio en que nos movemos es tridimensional.
Pero la aparición de los “monstruos matemáticos”, aquellas curvas que llenaban to-
do el espacio, desafiaba la idea intuitiva del concepto de dimensión. Este tipo de
objetos fractales fueron estudiados dentro de una nueva rama de las matemáticas
creada por Poincaré que recibe el nombre de Topologı́a. Una curiosa ciencia que
trata sobre la manera en que los objetos pueden ser deformados como si pareciesen
de goma. Un cı́rculo puede estirarse y convertirse en un triángulo, o en el copo de
nieve de Knoch. Desde el punto de vista topológico no es posible distinguir una lı́nea
recta de la curva de Hilbert, o una hoja de papel plana de aquella que se encuentra
completamente arrugada.
Definir el concepto de dimensión de manera rigurosa fue una tarea larga y compli-
cada que desafió a los mejores matemáticos de finales del siglo XIX y principios del
XX, entre los que se encontraban, Poincaré, Lebesgue, Brouwer, Cantor, Menger,
Peano, y Hilbert. Era inevitable que tan intenso trabajo diese como resultado un
número numeroso de definiciones, todas ellas relacionadas pero distintas: dimensión
topológica, dimensión por conteo por cajas, dimensión de autosimilitud, dimensión
euclı́dea,..., etc.
Por los comentarios anteriores queda claro que estudiar a fondo el concepto de
dimensión no es una tarea fácil, pero para el propósito que nos ocupa podemos
simplificarlo de la siguiente manera. Si tomamos un segmento de longitud uno y lo
dividimos en cinco partes iguales, es evidente que obtendremos 5 segmentos cada
3.6 Ejemplo de sistema dinámico: la Geometrı́a Fractal 105
La fórmula anterior se conoce con el nombre de ley de potencia, puesto que el número
de partes N (h) cambia como si fuese una potencia de 1/h. Con unos básicos cono-
cimientos matemáticos es posible despejar de esta expresión la dimensión, siendo su
valor
log(N (h))
D=−
h
. De esta manera, la dimensión fraccionaria o fractal del Conjunto de Cantor es
aproximadamente 0.63093, obtenido del cociente log(2)/log(3), y la dimensión de la
curva de Knoch log(2)/log(3) = 1.26186.
Tenemos que darnos cuenta de la importancia de estos resultados, ya que a través
de la dimensión fractal disponemos de un parámetro que caracteriza, de alguna ma-
nera, lo sinuoso e irregular de un objeto. Si el fractal es relativamente “suave” su
dimensión fractal estará próxima a su dimensión euclı́dea que es uno, por el con-
trario si “retuerce” hasta ocupar la forma de un cuadrado, su dimensión fractal
estará próxima a dos. La costa de Inglaterra tiene dimensión fractal 1.25, mientras
que la dimensión fractal de la costa de Australia es 1.13; en conclusión, la costa de
106 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
fórmula no lineal como el conjunto de Mandelbrot. Sin embargo, podemos hacer una
clasificación un poco más detallada de los mismos.
Fractales construidos a partir de un algoritmo de escape. Un ejemplo clásico
de ellos es el fractal de Mandelbrot obtenido, como sabemos, de la manera siguiente:
para cada punto del plano (la semilla), se calcula una serie de valores (su órbita),
mediante la repetición de una fórmula hasta que se cumpla una determinada con-
dición, y en ese momento al punto inicial se le asigna un color relacionado con el
número de repeticiones.
Por este trabajo recibió el gran premio de la Academia de las Ciencias. Durante
los años 20 disfrutó de gran fama, pero posteriormente su trabajo fue olvidado,
hasta que B. Mandelbrot lo puso de actualidad en 1970 con sus experimentos con
los ordenadores.
Otro tipo de fractales muy interesantes son los obtenidos por medio de un Sistema
de Funciones Iteradas (IFS). Es un método creado por M. Barnsley, basándose
en el principio de autosemejanza. En un fractal IFS siempre se puede encontrar una
parte de la figura que guarda una relación de semejanza con la figura completa.
Esa relación es a menudo muy difı́cil de apreciar, pero en el caso del helecho de la
figura siguiente es bastante clara: cualquier hoja es una réplica exacta de la figura
completa.
Existen fractales asociados con la teorı́a del caos que reciben el nombre de fracta-
les caóticos. Entre ellos, el más conocido es el atractor de Lorenz (izquierda) y el
diagrama de bifurcación (centro).
A continuación identificamos los lados 1 y 2 del dado con el punto 1, los lados 3
y 4 con el punto 2 y los lados 5 y 6 con el tercer punto. Ahora se elige un punto
cualquiera P de la hoja y lo desplazamos de la manera siguiente: lanzamos el dado
y movemos el punto P hasta la mitad de la distancia existente entre el punto P y
el vértice que indique el dado. Por ejemplo: si ha salido un 1 en el dado, entonces
situamos el punto P en la mitad de la distancia entre P y el vértice 1, si en el segundo
lanzamiento sale un 5, situamos el punto en la mitad de la última posición de P y
el vértice 3, y ası́ sucesivamente.
para incluir otros objetos de la naturaleza que adoptan esta forma, o el sistema
arterial y bronquial de un mamı́fero.
El proceso mediante el cual una rama principal se divide y subdivide sucesivamente
es el culpable de que el objeto tenga estructura fractal. La clave se encuentra en
la propiedad de autosemejanza por la cual se añade la misma estructura a escala
cada vez más pequeña. Pensemos, por ejemplo, en la red de afluentes de una cuenca
hidrográfica, donde el agua de lluvia tiene que ser desaguada, o en el sistema arte-
rial, donde la sangre tiene que llegar a todos los tejidos del cuerpo ocupando sólo
un 5 % del volumen total. Estos procesos pueden modelarse por un tipo de curva
(dimensión uno) que ocupa una superficie (dimensión dos), es decir, un fractal.
La regeneración de tejidos.
El gran problema para su tratamiento es que sólo puede ser diagnosticada con fia-
bilidad mediante una autopsia post-mortem o una biopsia, aunque existen criterios
no invasivos para diagnosticarla con aceptable certeza Sin embargo, su detección
precoz es muy importante puesto que se puede empezar el tratamiento más rápida-
mente, con el consiguiente aumento de la calidad de vida, y retraso en la aparición
de nuevos brotes. No hay ningún sı́ntoma tı́pico de la Esclerosis Múltiple que ayude
en el diagnóstico inicial. Es frecuente que el primer brote pase desapercibido sin que
la persona consulte con su médico.
Hasta ahora se habı́a descrito un único método capaz de detectar cambios en la
sustancia blanca aparentemente normal, la tasa de transferencia de magnetización,
el cual, además de ser costoso y difı́cil de realizar en la gran mayorı́a de los centros
hospitalarios, presenta resultados contradictorios en cuanto a su sensibilidad.
Se ha demostrado que, aunque están teniendo lugar procesos que contribuyen al
desarrollo de la enfermedad a nivel de la sustancia blanca cerebral, dicha sustancia
118 Tema 3 Modelos basados en ecuaciones y sistemas en diferencias
MODELOS DISCRETOS
MATRICIALES
4.1. Introducción
La dinámica de una población viene determinada por el número de nuevos nacimien-
tos y la probabilidad de morir que tienen los individuos que componen la población.
Por ello, es muy importante saber la estructura de edades de la población que es-
tamos estudiando. Es decir, cómo se reparten los individuos en las diferentes clases
de edad y lo que es más importante, conocer las probabilidades asociadas de su-
pervivencia, mortalidad y fecundidad. Generalmente esta información se refleja en
una tabla de vida, en la mayorı́a de los casos correspondiente a las hembras de la
población, ya que son las que contribuyen a la dinámica de la población en términos
de fecundidad.
El presente tema es una introducción al estudio de los modelos estructurados basa-
dos en el algebra matricial. Se inicia con un modelo probabilı́stico clásico como son
las cadenas de Markov y se centra fundamentalmente en el estudio de uno de los
modelos más conocidos como es el modelo de Leslie.
119
120 Tema 4 Modelos discretos matriciales
∑
m
Pij ≥ 0 , Pij = 1 , i, j ∈ {1, 2, · · · , m} . (4.1)
i=1
⃗
Con los números ptj formamos un vector columna X(t), y con los Pjk una matriz
cuadrada A. Ahora, las identidades (4.2) podemos expresarlas
⃗
X(t) ⃗ − 1) ,
= A X(t t = 1, 2, · · · n . (4.3)
La cadena de Markov homogénea está representada por el sistema de ecuaciones
lineales en diferencias anteriores. Observemos que por las condiciones (4.1) la matriz
A cumplirá:
La suma de los elementos de cada una de sus columnas vale la unidad. Sin
embargo, no ocurre lo mismo con la suma de los elementos de sus filas
Todos sus elementos Pjk son mayores o iguales que cero y menores o iguales
que uno.
A una matriz de transición de este tipo se la conoce con el nombre de matriz
estocástica. Aquellas matrices estocásticas donde la suma de los elementos de cada
fila sea la unidad reciben el nombre de doblemente estocásticas.
EJEMPLO 4.1
Para beber agua un animal puede ir a un lago o a un rı́o. Se sabe que no va al lago dos
dı́as seguidos y que si toma agua en el rı́o la probabilidad de que el dı́a siguiente beba
agua en cada uno de los sitios es la misma.
Estamos ante una cadena de Markov homogénea con dos estados, E1 que representa
al hecho de que el animal beba agua en el lago y E2 que beba agua en el rı́o.
La matriz de transición es ( )
0 1/2
A= .
1 1/2
EJEMPLO 4.2
Una tienda de animales que vende peces incluye una garantı́a por la que cualquier pez que
muera antes de cumplir tres meses se reemplaza de forma gratuita. Una vez que el pez se
ha reemplazado ya no queda cubierto por la garantı́a. Se sabe que:
2.- El 5 % de los peces que han cumplido un mes mueren durante el segundo.
3.- El 7 % de los peces que han cumplido dos meses mueren durante el tercer mes.
Podemos representar la situación anterior por medio de una cadena de Markov siendo
los estados, Ei := pez en el mes i = 1, 2, 3 de garantı́a, E4 := pez sin garantı́a por
haber sido repuesto, y E5 := pez sin garantı́a por tener más de 3 meses.
⃗
X(t) ⃗ − 1) ,
= A X(t t = 1, 2, 3, · · · ,
ya que uno de los objetivos que perseguimos al modelar una determinada situación
real es el de poder conocer su comportamiento a largo plazo. Puesto que la matriz
de transición A nos resuelve el problema de encontrar la ley de probabilidad, el
lı́mite de esta ley cuando el tiempo tiende a infinito, nos proporciona un método
para estudiar el comportamiento a largo plazo de la cadena de Markov.
∑
m
P [Ein /Ej0 ] = P [Ekn−1 /Ej0 ] P [Ein /Ekn−1 .Ej0 ] ,
k=1
Por la hipótesis (4.4) de inducción P [Ekn−1 /Ej0 ] = An−1 (k, j), y por la definición de
cadena de Markov P [Ein /Ekn−1 .Ej0 ] = P [Ein /Ekn−1 ] = Pik = A(i, k). Es decir,
∑
m
P [Ein /Ej0 ] = A(i, k) An−1 (k, j) ,
k=1
Las cadenas (a) y (b) de la Figura 4.3 no son regulares, ya que en el primer caso
sólo contiene un ciclo, y en el segundo el estado E4 no es accesible. Sin embargo,
la cadena (c) si es regular pues todos los estados son accesibles y además existen
los ciclos E2 E1 E3 E2 (impar) y E1 E2 E1 (par). En este último caso (c), la potencia
n-ésima de su matriz de transición es
0 1 0 1/3 1/3 1/3
A = 1/2 0 1/2 , At → 1/3 1/3 1/3
1/2 0 1/2 1/3 1/3 1/3
∑
4 ∑
4 ∑
4 ∑
4
aj1 = aj2 = aj3 = aj4 = 1 ,
j=1 j=1 j=1 j=1
y por tanto, todos los elementos de la primera fila del determinante anterior son
ceros, lo cual implica que este determinante es nulo, tal y como deseábamos probar.
TEOREMA 4.2.6 Si A es una matriz estocástica de orden n con todos sus ele-
mentos positivos (regular), entonces la sucesión de matrices An , n = 1, 2, · · · con-
verge hacia una matriz que tiene todas sus columnas iguales que coinciden con Π ⃗ tal
que:
∑
m
⃗ = 1.
Π(j)
j=1
La demostración de esta propiedad queda fuera del alcance del objetivo del curso
pero puede consultarse en la página 264 del libro “Modelos matemáticos y procesos
dinámicos ”de Santiago Pérez-Cacho y otros.
Como aplicación inmediata del Teorema 4.2.6 anterior, observemos que en el Ejemplo
4.1 la matriz de transición A de la cadena de Markov regular tiene como valor propio
λ1 = 1 , λ2 = −0.5 y autovectores U ⃗ 1 = (1, 2) , U
⃗ 2 = (−1, 1). En consecuencia,
⃗ = (1/3, 2/3) y por tanto
Π
( )
1/3 1/3
A →
t
,
2/3 2/3
cuando t → ∞.
126 Tema 4 Modelos discretos matriciales
EJEMPLO 4.3
Sea la matriz de transición correspondiente a seis estados
1 1/2 0 0 0 0
0 0 1/2 0 0 0
0 1/2 0 1/2 0 0
A= 0 0 1/2 0 1/2
.
0
0 0 0 1/2 0 0
0 0 0 0 1/2 1
Supongamos que en el momento inicial el sistema se encuentra en el estado E4 .
Figura 4.4.
Veamos como podemos pasar del estado inicial E4 al resto de los estados. Sabemos
que
⃗
X(0) = (0, 0, 0, 1, 0, 0)T .
Como puede apreciarse en la Figura 4.4, al cabo de un paso la probabilidad será,
⃗
X(1) = (0, 0, 1/2, 0, 1/2, 0)T ,
⃗
o bien X(1) ⃗
= A X(0). Del mismo gráfico deducimos que,
⃗
X(2) = (0, 1/4, 0, 1/2, 0, 1/4)T
⃗
X(3) = (1/8, 0, 3/8, 0, 1/4, 1/4)T
⃗
X(4) = (1/8, 3/16, 0, 5/16, 0, 3/8) .
O de forma matricial:
⃗
X(2) ⃗
= A X(1) , ⃗
X(3) ⃗
= A X(2) , ⃗
X(4) ⃗
= A X(3) .
4.2 Cadenas de Markov 127
A := {{1., 1/2, 0, 0, 0, 0}, {0, 0, 0.5, 0, 0, 0}, {0, 1/2, 0, 1/2, 0, 0}, {0, 0, 1/2, 0, 1/2, 0},
{0, 0, 0, 1/2, 0, 0}, {0, 0, 0, 0, 1/2, 1}}
MatrixPower[A, 200]
{{1., 0.8, 0.6, 0.4, 0.2, 0.}, {0., 1.07909*10−19 , 0., 1.746 *10−19 , 0., 0.}, {0., 0.,
2.82509*10−19 , 0., 1.746* 10−19 , 0.}, {0., 1.746*10−19 , 0., 2.82509*10−19 , 0., 0.}, {0.,
0., 1.746*10−19 , 0., 1.07909*10−19 , 0.}, {0., 0.2, 0.4, 0.6, 0.8, 1.}}
EJEMPLO 4.4
Figura 4.5.
Observemos que estamos ante una cadena de Markov cuyo diagrama de estados es
el siguiente:
128 Tema 4 Modelos discretos matriciales
es decir
⃗
X(1) ⃗
= AX(0) .
⃗
X(2) ⃗
= AX(1) ⃗
= A2 X(0) .
En general
⃗
X(t) ⃗
= At X(0) , t = 1, 2, · · · .
Ahora, estamos interesados en deducir este valor de una manera diferente. Obser-
vemos que si la matriz A fuese diagonal, entonces At serı́a muy fácil de encontrar,
bastarı́a elevar a t los elementos de la diagonal. Por esta razón, en primer lugar
procederemos a diagonalizar la matriz simétrica A.
Los valores propios de la matriz A son los siguientes:
−λ 2 1
2 3 2
|A − λ I| = 0 ⇒ 3 −λ 1
2
= 0,
1
3 −λ
1
3
desarrollando obtenemos la ecuación caracterı́stica
9 λ3 − 7 λ − 2 = 0 ,
cuyas soluciones son λ1 = 1, λ2 = −2/3, λ3 = −1/3. Por tanto, la matriz A es
diagonalizable siendo los subespacios propios asociados a estos valor propio
S1 =< (3, 3, 2) > , S2 =< (−1, 1, 0) > , S3 =< (−1, −1, 2) > .
En consecuencia, la matriz de paso C es,
3 −1 −1
C= 3 1 −1 .
2 0 2
Para encontrar At , actuamos de la manera siguiente
D = C −1 AC ⇒ A = CDC −1 ⇒ At = CDt C −1 ,
que en nuestro caso
−1
3 −1 −1 1 0 0 3 −1 −1
At = 3 1 −1 0 (−2/3)t 0 3 1 −1 .
2 0 2 0 0 (−1/3)t 2 0 2
Simplificando
( ( 2 )t ( 1 )t ) ( ( )t ( )t ) ( ( 1 )t)
1
3 + 4 −3 + −3 1
3 − 4 − 23 + − 31 8 (1 − − 3 )
3
8( ( )t ( )t )
8 ( ( )t ( )t ) ( 1 )t
At = 18 3 − 4 − 23 + − 13 1
3 + 4 − 23 + − 31 8( 1 − − 3
3
.
( ( 1 )t )
8 ( ( 1 )t ) ( 1 )t )
1
4 1 − −3 1
4 1 − −3 1
4 1 + 3 −3
Finalmente hacemos que t → ∞, entonces
3/8 3/8 3/8
At → 3/8 3/8 3/8 ,
1/4 1/4 1/4
y en consecuencia después de infinitas semanas la distribución de los ratones tiende
hacia
Primero = 83 X1 (0) + 38 X2 (0) + 38 X3 (0) = 38 (X1 (0) + X2 (0) + X3 (0)) = 38 Total
3 3 3 3 3
Segundo = 8 X1 (0) + 8 X2 (0) + 8 X3 (0) = 8 (X1 (0) + X2 (0) + X3 (0)) = 8 Total
1 1 1 1 1
Tercero = 4 X1 (0) + 4 X2 (0) + 4 X3 (0) = 4 (X1 (0) + X2 (0) + X3 (0)) = 4 Total
130 Tema 4 Modelos discretos matriciales
y en consecuencia si t → ∞,
3/8 3/8 3/8
At → 3/8 3/8 3/8
1/4 1/4 1/4
EJEMPLO 4.5
Supongamos que al realizar estudios climáticos en una determinada zona de nuestra pro-
vincia obtenemos los siguientes datos. Si un dı́a es caluroso, entonces la probabilidad de
que el dı́a siguiente sea también caluroso es 0.65, y 0.35 la probabilidad de que haga frı́o.
Por otro lado, si un dı́a es frı́o, entonces 0.7 es la probabilidad de que el dı́a siguiente siga
siendo frı́o y 0.3 de que sea un dı́a caluroso.
Si hoy hace frı́o, vamos a calcular la probabilidad de que pasado mañana haga frı́o.
Para encontrar la solución podemos utilizar el diagrama de la Figura 4.6. En él
observamos que la probabilidad pedida es:
Es decir, existe casi un 60 % de posibilidades de que si hoy hace frı́o pasado mañana
también lo siga haciendo. De forma similar, la probabilidad de que si hoy hace frı́o
pasado mañana sea un dı́a caluroso es
El ejemplo también puede resolverse utilizando las cadenas de Markov. Existen dos
estados E1 que representa al dı́a frı́o y E2 al dı́a caluroso, siendo la matriz estocástica
( )
0.70 0.35
A= .
0.30 0.65
⃗
Si hoy hace frı́o podemos representar esta situación por el vector X(0) = (1, 0)T .
⃗ ⃗ ⃗
El producto X(1) = AX(0) nos dará las probabilidades del próximo dı́a, y X(2) =
⃗
AX(1) 2 ⃗
= A X(0) las de pasado mañana
( )2 ( ) ( )
⃗ 0.70 0.35 1 0.7 ∗ 0.7 + 0.35 ∗ 0.3
X(2) = =
0.30 0.65 0 0.3 ∗ 0.7 + 0.65 + 0.3
Puede probarse fácilmente que La cadena de Markov es regular. Por tanto, para
realizar un estudio a largo plazo de los diferentes escenarios que pueden presentarse
debemos encontrar el valor de la matriz potencia An . Si utilizamos el programa
Mathematicar
Si hoy es un dı́a caluroso, ¿cuál es la probabilidad de que dentro de tres dı́as sea un
dı́a frı́o?.
132 Tema 4 Modelos discretos matriciales
EJEMPLO 4.6
Representemos por X1 (0) e X2 (0) a las poblaciones iniciales de conejos y zorros respecti-
vamente. Se sabe que el número de conejos en cualquier mes es la mitad de la población
de conejos del mes anterior y que el número de zorros en dicho mes es la suma de las
poblaciones de zorros mas la mitad de la de conejos en el mes anterior. Vamos a calcular
las poblaciones de zorros y conejos al cabo de “mucho”tiempo para estudiar la evolución
de las poblaciones a largo plazo.
Sean X1 (t) e X2 (t) las poblaciones de conejos y zorros al cabo de t meses. Del
enunciado del ejercicio se deduce
{
X1 (t + 1) = 0.5X1 (t)
, t = 0, 1, 2, · · ·
X2 (t + 1) = 0.5X1 (t) + X2 (t)
Si llamamos
( ) ( )
X1 (t) 0.5 0
⃗
X(t) = , A= , t = 0, 1, 2, · · · ,
X2 (t) 0.5 1
entonces
⃗
X(1) ⃗
= AX(0)
⃗ ⃗
X(2) = AX(1) ⃗
= A2 X(0)
.. ..
. .
⃗
X(t) ⃗ − 1) = At X(0)
= AX(t ⃗ .
Para completar el resto del ejercicio utilizamos el ordenador.
{0.5, 1}
P := Transpose[Eigenvectors[[A]]
Q := P.DiagonalMatrix[(0.5)k , 1].Inverse[P]
MatrixForm[Limit[Q, k → Infinity]]
( )
0 0
1 1
EJEMPLO 4.7
Supongamos que en una comunidad autónoma la población está dividida en cuatro clases,
E1 , E2 , E3 y E4 , ordenadas de mayor a menor de acuerdo con la riqueza. Una persona que
pertenece a una clase en un momento dado puede ascender, mantenerse o descender en el
siguiente con probabilidades dadas por la tabla:
E1 E2 E3 E4
E1 0.7 0.2 0.1 0
E2 0.2 0.4 0.1 0.3
E3 0.1 0.3 0.4 0.2
E4 0 0.1 0.4 0.5
⃗ k+1 = AX
X ⃗k
0.7 0.2 0.1 0 0.17 0.197
0.2 0.4 0.1 0.3 0.24 0.247
⃗ 2001 = AX
X ⃗ 2000 =
0.1
=
0.3 0.4 0.2 0.30 0.267
0 0.1 0.4 0.5 0.29 0.289
Del mismo modo,
0.7 0.2 0.1 0 0.197 0.214
0.2 0.4 0.1 0.3 0.247 0.251
⃗ 2002 = AX
X ⃗ 2001 =
0.1
=
0.3 0.4 0.2 0.267 0.258
0 0.1 0.4 0.5 0.289 0.276
⃗ 2000 = AX
X ⃗ 1999 ⇒ ⃗ 1999 = A−1 X
X ⃗ 2000
EJEMPLO 4.8
1.- Cada trabajador activo sólo tiene un hijo que sigue trabajando en el parque.
2.- La mitad de los hijos de los cientı́ficos lo son también, la cuarta parta pasa a ser
personal auxiliar especializado y el resto es personal colaborador no especializado.
3.- Los hijos del personal auxiliar se reparten entre las 3 categorı́as según los porcentajes
30 %, 40 %, 30 %
4.- Para los hijos de los colaboradores las proporciones de reparto entre las categorı́as
son 50 %, 25 % y 25 %.
Estamos ante una cadena de Markov donde el estado E1 representa a los cientı́ficos,
E2 al personal auxiliar, E3 al personal colaborador y la matriz de transición es A.
Es fácil ver que esta cadena es regular siendo el diagrama de estados el que aparece
dibujado en la Figura 4.7.
Para estudiar el comportamiento a largo plazo del modelo podemos calcular la matriz
potencia At cuando t → ∞. Un valor aproximado será
A:={{0.5,0.3,0.5},{0.25,0.4,0.25},{0.25,0.3,0.25}};
MatrixPower[A,500]
{{0.44176, 0.44176, 0.44176},{0.294118, 0.294118, 0.294118},
{0.264706, 0.264706 0.264706}}
Puesto que la cadena es regular, podemos utilizar el Teorema 4.2.6, para lo cual
necesitamos conocer los valores y vectores propios de la matriz A.
Eigenvalues[A]
{1., 0.15, −1.68812 10−17 }
Eigenvectors[A]
{{0.744438, 0.496292, 0.446663}, {0.784465, 0.496292, 0.446663},
{−0.707107, −3.18473 10−16 , 0.707107}}
La distribución estable vendrá dada por el vector propio asociado al valor propio 1
• el 44 % serán cientı́ficos,
• el 29 % serán personal auxiliar,
• el 27 % serán personal colaborador.
EJEMPLO 4.9
Supongamos que disponemos de una casa, un granero, un gato y un ratón. Los animales
pueden estar los dos en la casa, los dos en el granero o uno en el granero y otro en la casa.
Realizamos de forma sucesiva la siguiente experiencia:
Lanzamos dos monedas al aire, si salen dos caras cambiamos al ratón del lugar donde se
encuentre. Si salen una cara y una cruz, es el gato el que se cambia. Por último, si salen
dos cruces, entonces cambiamos al gato y al ratón del sitio donde se encuentran.
Si tenemos en cuenta las diferentes opciones para la casa, inmediatamente que-
dará también determinada las opciones para el granero. Los diferentes estados son:
Observemos que podemos modelizar la situación anterior por medio de una cadena
de Markov ya que la probabilidad Pij de pasar del estado Ej al Ei sólo depende
del i y del j. Por otro lado, como 1/4 es la probabilidad de sacar dos caras o dos
cruces y 1/2 la probabilidad de que salga una cara y una cruz, entonces la matriz
de transición para esta cadena es:
0 1/2 1/4 1/4
1/2 0 1/4 1/4
1/4 1/4 0 1/2
1/4 1/4 1/2 0
136 Tema 4 Modelos discretos matriciales
Para estudiar la evolución a largo plazo de esta cadena tenemos que ver en primer
lugar si es regular. Para ello al calcular
3/8 1/8 1/4 1/4
1/8 3/8 1/4 1/4
A2 =
1/4
1/4 3/8 1/8
1/4 1/4 1/8 3/8
(1, 1, 1, 1)T ,
Finalmente
0.25 0.25 0.25 0.25
0.25 0.25 0.25 0.25
At −→
0.25
cuando t → ∞.
0.25 0.25 0.25
0.25 0.25 0.25 0.25
⃗
X(0) = (1, 0, 0, 0)T ,
entonces
⃗
X(t) ⃗
= At X(0) = (0.25, .25, 0.25, 0.25)T ,
y es igual de probable que a largo plazo nos encontremos en cualquiera de los cuatro
estados posibles.
EJEMPLO 4.10
4.2 Cadenas de Markov 137
Figura 4.8.
En general
⃗x(t) = At ⃗x(0) , t = 1, 2, · · · .
Ahora, estamos interesados en deducir este valor de una manera diferente. Obser-
vemos que si la matriz A fuese diagonal, entonces At serı́a muy fácil de encontrar,
bastarı́a elevar a t los elementos de la diagonal. Por esta razón, en primer lugar
procederemos a diagonalizar la matriz simétrica A (recordemos que toda matriz
simétrica es diagonalizable) .
Los valores propios de la matriz A son los siguientes:
−λ 1 1
1 2 2
|A − λ I| = 0 ⇒ 2 −λ 1
2
= 0,
1
2 −λ
1
2
es decir,
1 Total
x1 (t) = 3 x1 (0) + 13 x2 (0) + 13 x3 (0) = 31 (x1 (0) + x2 (0) + x3 (0)) =
3
1 Total
x2 (t) = 3 x1 (0) + 13 x2 (0) + 13 x3 (0) = 31 (x1 (0) + x2 (0) + x3 (0)) =
3
1 1 Total
x3 (t) = 3 x1 (0) + 13 x2 (0) + 13 x3 (0) = (x1 (0) + x2 (0) + x3 (0)) =
3 3
EJEMPLO 4.11
Se dispone de dos urnas A y B y de dos bolas, una blanca y otra roja, ası́ como de
dos monedas. Las bolas pueden estar las dos en A, las dos en B, o una en A y otra
en B. Se realiza sucesivamente el siguiente experimento: se lanzan las dos monedas
al aire y si salen dos caras, se cambia la bola roja de urna; si se obtiene una cara y
una cruz, es la bola roja la que se cambia. Finalmente, si salen dos cruces se cambian
las dos bolas, blanca y roja, de la urna en la que se encuentre cada una a la otra
urna.¿Estamos ante una cadena de Markov?,¿cuál es la matriz de transición?
Evidentemente se trata de una cadena de Markov con los siguientes estados:
Estamos ante una cadena de Markov regular ya que la matriz A2 tiene a todos sus
elementos no nulos. En consecuencia, la distribución a largo plazo viene dada por el
vector propio (normalizado) asociado al valor propio estrictamente dominante λ = 1,
cuyo valor es (0.25, 0.25, 0.25, 0.25).
Conclusión: A largo plazo, existe la misma probabilidad (25 %) de que la urna A
se encuentre vacı́a, que contenga a la bola blanca, que contenga a la bola roja, o que
contenga a las dos bolas.
140 Tema 4 Modelos discretos matriciales
EJEMPLO 4.12
La tabla siguiente recoge la población humana entre los años 1800 y 1995 en miles de
millones de personas.
1800 1850 1870 1890 1910 1930 1950 1970 1975 1980 1985 1990 1995
0.91 1.13 1.30 1.49 1.70 2.02 2.51 3.62 3.97 4.41 4.84 5.29 5.75
Hay un amplio rango de modelos matemáticos todos ellos con un nivel más elevado de
dificultad. En este curso estudiaremos los más elementales, no obstante, a pesar de su
sencillez proporcionan un amplio número de aplicaciones.
Si representamos por yt al tamaño de una población en el tiempo t, existen cuatro procesos
que afectan al cambio de su tamaño, como son los nacimientos (N), las inmigraciones (I),
las muertes (M) y las emigraciones (E). Si suponemos el intervalo de tiempo [t, t + 1],
entonces el cambio de la población puede expresarse por medio de la siguiente ecuación
en diferencias:
yt+1 = yt + N + I − M − E
Los modelos que estudian el crecimiento de poblaciones independientemente de la densidad
de dichas poblaciones corresponden a los casos más simples, siendo las hipótesis más
generales que se suelen establecer las siguientes:
Todos los individuos son iguales.
EJEMPLO 4.13
Recuperemos el Ejemplo 4.12 relacionado con la población humana entre los años 1800 y
1995.
2 4 6 8 10 12
o bien,
y0 = e−0.0771347 = 0.925765 , r = e0.169038 = 1.18417 ,
y el modelo exponencial será ahora yt = 0.925765 ∗ 1.18417t
2 4 6 8 10 12
El modelo discreto exponencial es muy simplista, ya que estamos suponiendo que todos los
individuos de la población son iguales, pero es evidente que tanto la tasa de natalidad como
la de mortalidad dependen de la edad del individuo. Una manera de poder resolver este
inconveniente es dividir la población en clases de edades, y de esta manera asignar tasas
de fertilidad y natalidad a los individuos dependiendo de la clase a la que pertenezcan.
Observemos el parecido de las dos expresiones encontradas para los dos modelos estudiados
⃗
yt = rt y0 , en el modelo discreto exponencial, y X(t) ⃗
= At X(0), en el modelo discreto
matricial. No obstante, tenemos que hacer notar que ahora si podemos distinguir entre la
tasa de supervivencia de pájaros jóvenes y adultos.
EJEMPLO 4.14
Supongamos que en el modelo discreto matricial anterior, cada hembra adulta produce
por término medio cuatro hembras (lo cual indica al menos ocho huevos), la cuarta parte
de las hembras jóvenes sobreviven para llegar a adultas y las tres cuartas partes de las
hembras adultas sobreviven.
A la vista de estos datos, la matriz que representa al modelo es
( )
0 4
A=
0.25 0.75
Puede observarse que α, β ∈ [0, 1] y además α < β ya que no es tan probable que
sobrevivan más los pájaros jóvenes que los adultos.
144 Tema 4 Modelos discretos matriciales
En la Tabla 4.1 hemos escrito las razones X1 (t)/X2 (t) y Tt /Tt−1 del total de hembras
en los años sucesivos.
Notemos como la razón X1 (t)/X2 (t) se acerca a la constante 2.771 mientras que la
población total parece aumentar a una tasa constante del 44 % anual. Además, en
la Figura 4.11 se han representado las diferentes proporciones de hembras jóvenes
(en rojo) y de hembras adultas en función del tiempo, y se puede apreciar como a
partir de un determinado momento estas proporciones permanecen constantes.
0.8
0.7
0.6
0.5
0.4
0.3
5 10 15 20
¿Cuál es la razón teórica para que se cumplan las observaciones del ejemplo
anterior? Veamos que la respuesta a esta pregunta puede generalizarse y está muy rela-
cionada con el método de las potencias utilizado para estimar el valor propio dominante
de una matriz cuadrada.
Supongamos que U ⃗ 1, U
⃗ 2 sean los vectores propios asociados a los valores propios λ1 , λ2 de
la matriz A que representa al modelo discreto matricial. Entonces al formar estos vecto-
⃗
res una base del plano vectorial, podemos escribir X(0) = c1 U⃗ 1 + c2 U
⃗ 2 con unos valores
determinados de c1 y c2 . En consecuencia,
⃗
X(t) ⃗ − 1) = At X(0)
= A X(t ⃗ ⃗ 1 + c2 U
= At (c1 U ⃗ 2) ,
|A − λ I| = λ2 − βλ − γα = 0 ,
Cada grupo de edad cambiará por un factor λ1 cada año. Ası́, a la larga, la ecuación
⃗
X(t) ⃗
= At X(0) actúa igual que la ecuación yt = rt y0 . En un corto plazo (es decir, antes
de alcanzar la estabilidad) los números oscilan. La magnitud de esta oscilación depende
de la magnitud de λ2 /λ1 (que es negativa, con lo se explica la oscilación).
Para la matriz A con la que trabajamos en el ejemplo, obtenemos como valores propios:
λ1 ≈ 1.4436 ; λ2 ≈ −0.693 .
146 Tema 4 Modelos discretos matriciales
es decir si
1−β
γ> .
α
Como en nuestro caso β = 0.75, α = 0.25, resulta que γ > 0.25
0.25 = 1, con lo que la población
aumentará, cosa que ya sabı́amos (lo hace a un ritmo del 44 %).
Consideraciones finales:
1.- Para que podamos aplicar el modelo discreto matricial, es necesario que el medio am-
biente sea estable, es decir que los cambios ecológicos que se registren no modifiquen
las tasas de natalidad y mortalidad de los individuos de la población.
2.- En las hipótesis del modelo no hemos tenido en cuenta un factor muy importante
como es la densidad de la población. Es evidente que las las tasas de natalidad y
supervivencia varı́an con el tamaño de la población, ya que en caso contrario la
población crecerı́a de forma ilimitada y dominarı́a al resto de las especies.
Una clasificación muy general de los modelos matriciales aplicados al estudio de la evolu-
ción de las poblaciones en función de la matriz A de proyección, puede verse en la siguiente
figura:
CLASE EDAD
1 [0,E/4)
2 [E/4,E/2)
3 [E/2,3E/4)
4 [3E/4,E]
E 2E kE
t0 = 0; t1 = , t2 = ; · · · ; tk = ; ···
4 4 4
Los procesos de nacimiento y muerte entre dos tiempos consecutivos de observación se
pueden describir mediante los siguientes parámetros demográficos:
Al promedio del número de hijas que tiene una hembra durante el tiempo que
permanece en la clase de orden i, lo llamaremos ai con i = 1, 2, 3, 4
1.- ai ≥ 0, i = 1, 2, 3, 4 .
El caso bi = 0, no puede ocurrir ya que esto supondrı́a que ninguna hembra vivirı́a mas
allá de la clase i. También supondremos que hay al menos un ai > 0 lo que garantiza que
habrá nacimientos. A la clase donde ai > 0 la llamaremos clase fértil.
⃗
Si X(k) = (X1 (k), X2 (k), X3 (k), X4 (k))T representa al vector de distribución de las eda-
des en el tiempo tk , entonces el número de hembras de la primera clase en el tiempo tk
vendrá dado, únicamente por las nacidas entre los tiempos tk−1 y tk . Podemos escribir,
X1 (k) = a1 X1 (k − 1) + a2 X2 (k − 1) + a3 X3 (k − 1) + a4 X4 (k − 1) . (4.6)
4.3 Modelo de Leslie 149
EJEMPLO 4.15
Supongamos que las hembras de una población animal viven por término medio 25 años
y que esta población se divide en cinco clases de edades iguales con intervalos de 5 años.
Supongamos que la matriz de crecimiento de Leslie viene dada por
0 2 1 3 1
1 0 0 0 0
2 1
L= 0 3 02 0 0 .
0 0
3 0 0
3
0 0 0 4 0
Calculando
0 2 1 3 1 10 83
1 0 0 0 0 20 5
2
⃗
X(1) ⃗
= LX(0) =
0
1
3 0 0 0
5 =
20/3 .
0 0 2
0 0 10 10/3
3
3
0 0 0 4 0 8 15/2
simplificando
p(λ) = λ4 − a1 λ3 − a2 b1 λ2 − a3 b1 b2 λ − a4 b1 b2 b3 = 0 . (4.9)
A la vista de la expresión anterior, se justifica la introducción de una nueva función,
a1 a2 b1 a3 b1 b2 a4 b1 b2 b3
q(λ) = + 2 + + . (4.10)
λ λ λ3 λ4
Ahora, resolver la ecuación p(λ) = 0 es equivalente a resolver la equivalente q(λ) = 1. Un
rápido estudio de la función q(λ) nos permite deducir las siguientes propiedades:
4.3 Modelo de Leslie 151
Estas propiedades nos permiten esbozar una gráfica de la función q(λ), la cual puede verse
en la Figura 4.13.
Observemos que existe un único valor λ1 positivo, tal que q(λ1 ) = 1. Esto es, la matriz
de Leslie, L tiene un único valor propio λ1 positivo para el cual q(λ) = 1. Además, al ser
q ′ (λ1 ) ̸= 0 la raı́z será simple, o bien, su grado de multiplicidad es 1.
El paso siguiente será el de calcular el autovector propio asociado al valor propio λ1 . Por
definición, λ1 es aquel valor no nulo que cumple, LU⃗1 = λ1 U⃗1 , siendo U⃗1 el vector propio
asociado. Si X⃗ = (X1 , X2 , X3 , X4 )T , entonces
a1 − λ a 2 a3 a4 X1 0
b1 −λ 0 0 X2 0
(L − λ1 I) X
⃗ =0 ⇒ =
0 b2 −λ 0 X3 0
0 0 0 −λ X4 0
TEOREMA 4.3.1 Una matriz de Leslie L, tiene un único valor propio positivo λ1 . Este
⃗ 1 cuyas coordenadas son todas
valor propio es simple y tiene un vector propio asociado U
positivas
|λi | ≤ λ1
Para el estudio que estamos realizando se requiere que |λi | < λ1 para todos los valores
propios de L; en este caso, ya sabemos por el tema anterior, que λ1 será un valor propio
dominante de L.
EJEMPLO 4.16
Debemos observar que no todas las matrices de Leslie cumplen este requisito. En
1941 Harro Bernadelli publicó un trabajo en el Journal of the Burma Research Socie-
ty con el tı́tulo “Population Wawes”, donde observó un comportamiento periódico
en lugar de un comportamiento estable de la población. En concreto, propuso la
siguiente matriz:
0 0 6
L = 21 0 0 .
0 13 0
4.3 Modelo de Leslie 153
0 0 0 λk4
λi k
lı́m ( ) = 0, i = 2, 3, 4 .
k→∞ λ1
En consecuencia
1 0 0 0
1 ⃗ 0 0 0 0
lı́m X(k) = C
0
C −1 X(0)
⃗ .
k→∞ λk 1
0 0 0
0 0 0 0
Puede probarse, que el lado derecho de la igualdad, correspondiente al producto de las
matrices, coincide con
1 ⃗ ⃗1 ,
lı́m k X(k) = dU (4.11)
k→∞ λ1
⃗
X(k) ⃗ 1 ≃ λ1 X(k
⃗ 1 = λ1 dλk−1 U
= dλk1 U ⃗ − 1) .
1
EJEMPLO 4.17
3⃗
⃗
X(k) ≃ X(k − 1) .
2
Las hembras estarán distribuidas de acuerdo con la relación 1 : 1/3 : 1/18; lo cual
corresponde a una distribución del 72 % de hembras en la primera clase 24 % en la
segunda y 4 % en la tercera.
⃗
X(k) ⃗1
= dλk1 U
que da el vector de distribución de la población por edades, para valores grandes de tiempo.
Se presentan tres casos que dependen del valor propio positivo λ1 .
Este último caso es de especial interés ya que determina una población de crecimiento
cero.
Para cualquier distribución inicial de las edades, la población tiende a una distribución
⃗ 1 . Teniendo en cuenta que L U
en el limite que es algún múltiplo del vector propio U ⃗1 =
⃗1 = U
λ1 U ⃗ 1 , puede comprobarse que,
λ1 = 1 ⇐⇒ a1 + a2 b1 + a3 b1 b2 + · · · + an b1 b2 · · · bn−1 = 1
EJEMPLO 4.18
Supongamos que una población de animales hembras está dividida en dos clases de
edades. En cada perı́odo el 50 % de la primera pasa a la segunda. El número medio
de crı́as hembras de las de la primera clase es de 1 y de las de la segunda es 1.5.
Para construir el modelo de Leslie recordemos que conocemos por supervivencia el
porcentaje de hembras que sobreviven en un periodo de tiempo o más y por fertilidad
el número de hembras que por término medio tiene en un perı́odo de tiempo cada
una de las hembras de la población.
⃗
La expresión matricial del modelo de Leslie X(k) ⃗ − 1) es:
= LX(k
( ) ( )( )
X1 (k) 1 1.5 X1 (k − 1)
= .
X2 (k) 0.5 0 X2 (k − 1)
y ası́ sucesivamente.
Para estudiar su comportamiento en el lı́mite es necesario en primer lugar resolver
la ecuación caracterı́stica,
El vector propio correspondiente al valor propio positivo (que por el teorema estu-
diado será dominante) es,
( ) ( )
⃗1 = b1 T 1 T
U 1, = 1, .
λ1 3
Por tanto,
3⃗
⃗
X(k) ≃ X(k − 1) ,
2
lo cual indica, que para valores de k grandes, en cada perı́odo de tiempo aumentará el
número de hembras en un 50 % en cada una de las clases. Como además,
( )
⃗ = d( 3 )k 1 T
x(k) 1, .
2 3
EJEMPLO 4.19
Con ayuda del ordenador encontramos los valores propios de esta matriz,
Al ser λ1 = 1.05941 > 1 el valor propio estrictamente dominante, nos indica que la
población crece cada 15 años a un ritmo del 6 % (aproximadamente).
Si nos fijamos en el vector propio ⃗v1 asociado al valor propio λ1 ,
yt = y0 r t , t = 0, 1, 2, · · · ,
Definimos la fertilidad como el número medio de hembras que han nacido al finalizar
la primavera de una hembra con una edad x determinada, y la representaremos por
b(x).
4.4 Tablas de vida y modelo de Leslie 159
Por ejemplo b(5) = 3 significa que una hembra de 5 años tiene por término medio, al
finalizar la primavera, 3 hembras recién nacidas. La fertilidad será por tanto un número
positivo, que al expresar valores medios puede ser cero (el individuo de edad x no es fértil),
o bien un número decimal.
La Tabla 4.2 nos da una hipotética tabla de vida para un organismo que vive 4 años.
Tabla 4.2.
En la tabla anterior, comenzamos con 500 individuos, los cuales todos han fallecido al
iniciarse el quinto año.
S(x)
l(x) = .
S(0)
La representación gráfica de l(x) en función de x nos da una gráfica que se conoce con el
nombre de curva de supervivencia. Es bastante corriente utilizar la escala logarı́tmica, en
el eje de abscisas colocamos la edad del individuo x y en el de ordenadas ln(l(x)). Estas
curvas corresponden a algunos de los tipos que aparecen a la izquierda de la Figura 4.14.
Observemos en la Tabla 4.5 que l(1) = 0.8, es decir, el 80 % de la población inicial sobrevive
hasta llegar a la edad 1. Nuestra población corresponde al tipo I, como puede observarse
en la Figura 4.14 (derecha).
l(x + 1)
g(x) = .
l(x)
EJEMPLO 4.20
Supongamos que para x = 0 tenemos 100 peces en un acuario. Contamos la población una
vez al dı́a y obtenemos los siguientes datos:
x 0 1 2 3 4 5 6 7 8 9 10 11
Núm 100 85 72 61 52 44 37 31 26 22 19 16
x 12 13 14 15 16 17 18 19 20 21 22 23
Núm 14 12 10 8 7 6 5 4 3 3 2 2
De la tabla anterior obtenemos los diferentes valores de l(x) = S(x)/100, los cuales
se encuentran reflejados en la siguiente tabla:
x 0 1 2 3 4 5 6 7 8 9
l(x) 1 0.85 0.72 0.61 0.52 0.44 0.37 0.31 0.26 0.22
x 10 11 12 13 14 15 16 17 18 19
l(x) 0.19 0.16 0.14 0.12 0.10 0.08 0.07 0.06 0.05 0.04
x 20 21 22 23 - - - - - -
l(x) 0.03 0.03 0.02 0.02 - - - - - -
En nuestro ejemplo podemos ver en la Tabla 4.5 que la probabilidad de que un individuo
de edad 1 sobreviva y llegue a la edad 2 es de un 50 %.
Para poder estimar el valor r en (4.12) a partir de l(x) y b(x) es necesario encontrar, en
primer lugar, otros dos números como son la tasa neta de reproducción y el tiempo de
generación G.
Se define la tasa neta de reproducción R, como el número de individuos que por
término medio tiene una hembra durante toda su esperanza de vida. Es decir,
∑
k
R = l(0)b(0) + l(1)b(1) + · · · + l(k)b(k) = l(x)b(x) .
x=0
La matriz L sabemos que es la matriz de Leslie que tiene como primera fila los valores de
la natalidad y su subdiagonal principal son las probabilidades de supervivencia, el resto
de los elementos de la matriz son ceros. En la sección anterior hemos demostrado que para
una población con parámetros de nacimientos y muertes constantes, independientemente
de los valores iniciales, cuando ha transcurrido un “número adecuado”de generaciones el
porcentaje de individuos en cada una de las clases permanece constante, aunque el tamaño
total de la población crece exponencialmente.
EJEMPLO 4.21
Supongamos la siguiente tabla de vida para una población de caracoles:
Edad en años S(x) b(x)
0 500 0
1 400 2.5
2 40 3
3 0 0
∑
3
l(x)b(x)x
x=0 2.48
G= = = 1.107 años
∑3 2.24
l(x)b(x)
x=0
ln R
r= = 0.729 individuos/(individuos × año)
G
Para encontrar el valor exacto de la tasa de reproducción r utilizamos la ecuación
de Euler
∑3
1= e−rx l(x)b(x) ,
x=0
el valor r = 0.729 se encuentra por debajo del valor exacto. Probamos con diferentes
valores hasta llegar al r = 0.749.
164 Tema 4 Modelos discretos matriciales
Para construir el modelo de Leslie empezamos completando esta otra tabla, para
encontrar las tasas de natalidad ai , i = 1, 2, 3 y de supervivencia bi , i = 1, 2.
EJEMPLO 4.22
∑
4
l(x)b(x)x
x=0 4.3
G= = = 1.482 años
∑4 2.9
l(x)b(x)
x=0
ln R0
r= = 0.718 individuos/(individuos× año)
G
Supongamos que inicialmente la población de caracoles es de 200 en la primera clase,
0 en la segunda, 0 en la tercera, y 0 en la cuarta. Estamos interesados en construir
la matriz de Leslie para esta tabla de vida y proyectar la población “a largo plazo”.
Para ello elaboramos la tabla.
lo que nos permite proyectar la población para cualquier año. Por ejemplo, al cabo
de 5 años
⃗ (5) = LN
N ⃗ (4) = L5 N
⃗ (0) = (7613, 2804, 642, 75)T ,
Como antes, la diagonal principal representa la probabilidad de que una colonia perma-
nezca en la misma clase de tamaño. Los elementos de la subdiagonal principal representan
la probabilidad de que una colonia crezca y pase al tamaño siguiente. Sin embargo, aho-
ra existe la posibilidad de que parte de la colonia pueda fragmentarse (Pgm ) y pasar de
ser grande a ser mediana, o bien (Pmp ) pertenecer a las colonias pequeñas. Las colonias
pequeñas pueden agruparse y formar colonias medianas (Ppm ) o directamente colonias
grandes (Ppg ). Finalmente, observemos que la primera fila es la suma de dos términos, el
primero de ellos corresponde a la fecundidad, y el segundo a la transición de un estado
a otro. Puede probarse que para este tipo de modelos su comportamiento en el lı́mite es
exactamente igual al modelo de Leslie, es decir:
EJEMPLO 4.23
Es posible introducir esta nueva hipótesis en el estudio de los modelos matriciales. Por
ello, al construir un modelo debemos tener en cuenta las siguientes consideraciones:
¿La densidad depende sólo de la población total o por el contrario el efecto está dis-
tribuido sobre varias clases de edades?
p1 (N ) = p2 (N ) = 0.75
ps = 0.2
gs (N ) = 0.05
B1 (N ) = B3 (N ) = 100e−0.0001N
B2 (N ) = 200e−0.0001N
0.5 200
0.4
150
0.3
100
0.2
50
0.1
Como puede apreciarse, entre dos intervalos de tiempo, el 20 % de las semillas del semillero
sobrevivı́an pero seguı́an dormidas, mientras que 5 %de ellas germinaban. La proporción
de individuos que sobreviven p0 es inicialmente 0.75 cuando la densidad de la población
es muy baja, pero a medida que ésta aumenta, dicha proporción va disminuyendo.
Puede comprobarse que la población crece y además se mantienen constantes las propor-
ciones de cada una de las clases .
4.6 Modelos que dependen de la densidad 173
EJERCICIOS PROPUESTOS
EJERCICIO 3
1.- Supongamos un sistema con dos estados: un cazador dispara contra unos
animales. Existe el estado E1 de acierto, y el estado E2 de fallo. Se sabe
que si el cazador acierta, entonces en el segundo disparo tiene una proba-
bilidad de 3/4 de acertar. Si por el contrario en un disparo falla, entonces
la probabilidad de acertar en el siguiente es de 1/2.
2.- Se encierra a una rata en una caja dividida en compartimentos con puer-
tas que los comunican según se muestra en la figura. Cuando la rata sale
de un compartimento elige uno al azar.
3.- Las familias de un paı́s se clasifican según residan en áreas rurales, urba-
nas o suburbanas. Los estudios de movilidad demográfica estiman que, en
promedio, en el curso de un año, el 15 % de las familias urbanas cambian
174 Tema 4 Modelos discretos matriciales
8.- Una sala de cine decide programar las pelı́culas según el siguiente méto-
do: si una semana se proyectó una norteamericana, a la semana siguiente
se programará, dos de cada tres veces, una española, y una de cada tres
veces, una francesa. Si la pelı́cula programada fue francesa, dos de cada
tres veces será norteamericana y una de cada tres francesa. Finalmente,
si la pelı́cula programada fue española, la semana siguiente se progra-
mará española una de cada tres veces y norteamericana dos de cada tres
veces. Si inicialmente las cuotas de pantalla son el 50 % para el cine nor-
teamericano, el 35 % para el cine español, y el 15 % para el francés.
12.- Una población de ardillas está divida en tres clases de edades de dos años
de duración, a las que llamaremos jóvenes, medianas y adultas. La matriz
de Leslie viene definida de la siguiente manera: una hembra joven aporta
otra hembra y una mediana 24, además la cuarta parte de las jóvenes
sobreviven para llegar a medianas y el 50 % de las medianas se hacen
adultas.
15.- Sea una población de hembras dividida en tres clases de edades de 5 años
de duración. Su evolución está determinada por un modelo de Leslie
siendo su matriz,
1 a2 2
L = 1/2 0 0
0 2/3 0
¿Desaparecerá esta población a largo plazo?
Encontrar el valor de a2 para que cada 5 años la población aumente
en un 50 %
Para el valor de a2 anteriormente encontrado. Si a largo plazo el
número de hembras es de 800, ¿cuántas de ellas serán jóvenes?
178 Tema 4 Modelos discretos matriciales
BIBLIOGRAFÍA BÁSICA
179
180
[4] BURDEN, R.L.; FAIRES, J.D. Análisis Numérico, 2a ed. Grupo Editorial
Iberoamericano, S.A., 1996.
Es un libro interesante para introducirse en el Análisis Numérico, ya que es eminen-
temente práctico, prescindiendo en una gran parte de los casos de las demostraciones.
Cada tema se inicia con un ejemplo motivador del concepto a estudiar, posteriormen-
te se desarrolla el aspecto teórico completado con un número elevado de ejemplos y
finalmente, cuando es posible, se termina con los algoritmos adecuados para su im-
plantación en el ordenador.
Sus doce capı́tulos desarrollan gran parte del temario clásico de un curso de Análi-
sis Numérico. Es decir: Interpolación y Aproximación Polinomial. Diferenciación e
Integración Numérica. Problemas de Valor Inicial para Ecuaciones Diferenciales Or-
dinarias. Métodos Directos para la Solución de Sistemas Lineales. Técnicas Iterativas
en Álgebra de Matrices. Teorı́a de Aproximación. Aproximación de Valores Propios.
Soluciones Numéricas de Sistemas de Ecuaciones no Lineales. Problemas de Valores
en la Frontera para EDO. Soluciones Numéricas a Ecuaciones en Derivas Parciales.
Es el complemento adecuado para que los alumnos interesados puedan profundizar en
las cuestiones planteadas en clase, relacionadas con el cálculo aproximado de valores
y vectores propios, ası́ como con la resolución aproximada de ecuaciones y sistemas
de ecuaciones diferenciales ordinarias.
BIBLIOGRAFÍA COMPLEMENTARIA
MODELOS DISCRETOS
[6] ANTON, H. Álgebra Lineal, 3a ed. Limusa S.A., México, D.F., 1990.
[10] DEVANEY R.L. First Course in Chaotic Dynamical Systems: Theory and Experi-
ment. Addison - Wesley, 1992.
[11] DEVANEY R.L. Chaos, Fractals, and Dynamics: Computer Experiments in Mat-
hematics. Addison - Wesley, 1990.
[12] GROSSMAN, S.I. Álgebra lineal con aplicaciones, 4a ed. McGraw Hill Interameri-
cana de México S.A., México, D.F., 1991.
[13] HOLMGREN, R.A. A first course in discrete dynamical systems, 2a ed. Springer -
Verlag, New York, Inc., 1996.
ECUACIONES DIFERENCIALES
[17] BERMUDEZ, L., POCIELLO, E., RUÍZ, M.E.; VAREA, J. Ecuaciones dife-
renciales y en diferencias finitas. Ediciones Media, Sant Cugat del Vallés, 1995.
[26] QUESADA, J.M.; MOLINA, M.F.; SÁNCHEZ, F.T. Matemáticas II para In-
genierı́a Técnica Industrial. Los autores, Jaén, 2000.
CÁLCULO NUMÉRICO
[30] ISAACSON, E.; KELLER, H.B. Analysis of numerical methods. John Wiley and
Sons, 1966.
183
PAQUETES INFORMÁTICOS
[34] RAMÍREZ V., ET AL. Matemáticas con Mathematica, vol. I,II y III. Proyecto
Sur de Ediciones, Granada, 1995.
[36] BEGON, M.; MORTIMER, M.; THOMPSON, D.J. Population Ecology. A uni-
fied study of animals and plants, 3a . Blackwell Science, 2000.
[38] GIORDANO, F.R.; WEIR, M.D.; FOX, W.P. A first course in Mathematical
Modelling, 2a ed. Brooks/Cole, Pacific, California, 1997.
[39] GOTELLI, N.J. A primer of ecology. Sinauer Associates, Inc. Publishers, Sunder-
land, Massachusetts, 1995.
[40] DE GUZMAN, M.; ET al. Estructuras fractales y sus aplicaciones. Labor. Serie
Matemáticas. Barcelona, 1993.
[41] DOUCET; SLOEP Mathematical Modeling in the Life Sciencies. Ellis Horwood,
1992.
[44] HANNON, B.; RUTH, M. Modeling dynamic biological systems. Springer - Verlag,
New York, Inc., 1997.
[47] MAHAFFY, J.M. Modeling Mathematical. San Diego State University, USA, 2001.
184
[51] MURRAY, J.D. Mathematical Biology, 2a ed. Springer - Verlag, New York, 1993.
PAQUETES DE CÁLCULO
RECURSOS EN LA RED
http://smub.st-and.ac.uk/jason matthiopoulos/page5.html
http://www.epa.gov/epahome/models.htm
http://www.uni-klu.ac.at/ gossimit/links/bookmksd.htm
http://cisat.isciii.es/supercourse/
185
http://perso.wanadoo.fr/l.d.v.dujardin/ct/eng index.html
http://hypertextbook.com/chaos/links.shtml
http://www3.uakron.edu/biology/mitchell/ecology/ecolinks.html
http://rulbii.leidenuniv.nl/wwwkim/popdyn.html
http://www.gypsymoth.ento.vt.edu/ sharov/PopEcol/popecol.html
http://www-rohan.sdsu.edu/ jmahaffy/courses/s00/math121/
http://everest.ento.vt.edu/ sharov/3d/3dinsect.html