Cadenas de Markov
Cadenas de Markov
Cadenas de Markov
T I T U L O Cadenas de Marcov
Nombre de la Carrera
Ingeniería Industrial
noviembre de 2023
Av. Tecnológico #1, Col. Liberación C.P. 28976, Villa de Álvarez, Colima, Tel. 312
3129920 ext. 215 y 115 Correo electrónico: [email protected] tecnm.mx |
colima.tecnm.mx
Vv
1
ALEXIS GABRIEL RUVALCABA GALINDO
ÍNDICE
INTRODUCCIÓN…………………………………………………………………….3
CONCLUSIÓN………………………………………………………………………..26
BIBLIOGRAFIA……………………………………………………………………….26
2
ALEXIS GABRIEL RUVALCABA GALINDO
INTRODUCCIÓN
En las cadenas de Markov considera ciertos puntos discretos en el tiempo, para
toma valores de 1, 2, 3,…. Y la variable aleatoria que caracteriza al sistema. Y la
familia de las variables aleatorias forma un proceso llamado proceso estocástico.
Entonces las cadenas de Markov se usan para estudiar ciertos comportamientos de
largo y cortos plazos de sistemas estocásticos.
Las cadenas de Markov son modelos probabilísticos que se usan para predecir la
evolución y el comportamiento a corto y a largo plazo de determinados sistemas.
Ejemplos: reparto del mercado entre marcas; dinámica de las averías de máquinas
para decidir política de mantenimiento; evolución de una enfermedad.
Las cadenas de Markov se pueden aplicar en áreas como la educación,
comercialización, servicios
3
ALEXIS GABRIEL RUVALCABA GALINDO
4
Otro método para exhibir las probabilidades de transición es usar una matriz de
transición. . La matriz de transición para el ejemplo del diagrama de estados se muestra
en la tabla 4.1 .
5
Procesos estocásticos.
6
existencia de probabilidades de transición (de un paso) estacionarias también implica que,
para cada i, j y n (n = 0, 1, 2,...),
Ejemplo :
Una tienda de cámaras tiene en almacén un modelo especial de cámara que se puede
ordenar cada semana. Sean D1, D2, ...Dn las demandas de esta cámara durante la
primera semana, segunda semana, n-ésima semana, respectivamente. Se supone que las
Di son variables aleatorias independientes e idénticamente distribuidas que tienen una
distribución de probabilidad conocida. Sea X0 el número de cámaras que se tiene en el
momento de iniciar el proceso, X1 el número de cámaras que se tienen al final de la
semana uno, X2 el número de cámaras al final de la semana dos, etc. Suponga que X 0 = 3.
El sábado en la tarde la tienda hace un pedido que le entregan el lunes al momento
de abrir. La tienda usa la siguiente política ( s, S) para ordenar : si el número de cámaras
en inventario al final de la semana es menor que s = 1 (no hay cámaras en la tienda),
ordenar S = 3. De otra manera, no coloca la orden (si se cuenta con una o más cámaras
en el almacén, no se hace el pedido). Se supone que las ventas se pierden cuando la
demanda excede el inventario. Entonces, { Xt } para t = 0, 1, ..n es un proceso estocástico
de la forma que se acaba de describir. Los estados posibles del proceso son los enteros
0, 1, 2, 3, que representan el número posible de cámaras en inventario al final de la
semana.
7
Se observa que {Xi } ( en donde Xi es el número de cámaras en el almacén al final
de la semana i antes de recibir el pedido ), es una cadena de Markov. Se ve ahora cómo
obtener las probabilidades de transición (de un paso), es decir, los elementos de la matriz
de transición ( de un paso), Suponiendo que cada Dt tiene una distribución Poisson con
parámetro = 1.
P0 P P0 P0 p P1 P P1
P1
=P2 P P2 P2
P3 P P3 P3
0 1 2 3
0 1 – p(2) = . (1 – p(1)) – (1 – p(2)) = (1 – p(0)) – (1 – p(1)) = p(0) =
1 1 – p(0) = .080p(0) = 0.3680.184 0.3680 0.368 0
2 1 – p(1) = .632(1 – p(0)) – (1 – p(1)) = p(0) = 0.368 0
3 1 – p(2) = .264(1 – p(1)) – (1 – p(2)) 0.368= (1 – p(0)) – (1 – p(1)) = p(0) =
080 0.184 0.368 0.368
Probabilidad de transición estacionaria de 1 paso.
8
p
0 1 2 3
0 . 08 .18 . 36 .36
1 .63 .36 0 0
2 . 08 . 18 . 36 0
3 . . . .36
Las ecuaciones de Chapman-Kolmogorov proporcionan un método para calcular
estas probabilidades de transición de n pasos:
( ) ( ) ( )
∑
( )
∑ )( ( )y ( ) ( )
9
( )∑
Note que las ( ) son los elementos de la matriz P(2) ,pero también debe observarse
Ejemplo :
Retomando el ejemplo de las cámaras, para p(2) = p2, se tiene:
p * p = p2
0 1 2 3 0 1 2 3 0 1 2
3
0. . . .. . . .. .
. .
* =
108. 18. 360 36008. 18. 360 36024. 28. 300. 16. 263. 36. . 063. 36.
. 028. 25. 233. 23. 326. 36. 36. .26. 36. 36. .35. 31. 233. 09.
Así, dado que se tiene una cámara al final de una semana, la probabilidad de que
( ) no haya
cámaras en inventario dos semanas después es 0.283; es decir, .
De igual manera, dado que se tienen dos cámaras al final de una semana, la
probabilidad de que haya tres cámaras en el almacén dos semanas después es 0.097;
( )
esto es,
La matriz de transición de cuatro pasos también puede obtenerse como:
10
0 123012301230.... ........
11
De igual manera, dado que quedan dos cámaras en el almacén final de una
semana, se tiene una probabilidad de 0.171 de que haya tres cámaras en el
almacén 4
( )
semanas después; esto es,
Teorema
Ejemplo :
12
Compra 1 Refresco Refresco p = pAA pAB.90
=
.10
Refresco A 90A% 10B%pBA pBB.20
.80
Refresco B 20% 80
13
Entonces : .90 1+
.20
= .10 1+
1 1 2 .90 .
1 .80
*.20 1. =2
8=
De la condición 1 + 2= 1, 1= 1- 2. Al reemplazar esta expresión en la
Ejemplo :
Cierta compañía se especializa en el arrendamiento de camiones a personas que
desean realizar sus propias mudanzas. El gerente de distribución de la compañía está
considerando la posibilidad de aplicar un “cargo por traslado” para cubrir el costo de enviar
camiones desde áreas en las que hay sobrantes a otras en las que se necesitan. Antes de
decidir si se debe aplicar el cargo por traslado al costo de arrendamiento de los camiones
que se dirigen a áreas en las que hay sobrantes, desea determinar la proporción del
número total de camiones que, a largo plazo, acabarán en cada una de las áreas de renta.
Si las proporciones son aproximadamente las mismas, el cargo por traslado será
innecesario. Si no es así, el cargo dependerá de la proporción del total que termine en
cada región. El gerente de distribución ha dividido la parte del país que atiende la
compañía en tres regiones: norte, centro y sur. De registros previos se ha determinado la
siguiente información, referente a la proporción de camiones que se rentan en
determinada región y dónde se entregan:
14
Región Región donde se
Centr Sur
Nortrenta e20 % o30 % 50 %
e 30 % 40 % 30 %
Cent 20 % 40 % 40 %
donde se devuelveNort
M.C. José Alberto estrada Beltrán 8
En este momento 40% de los camiones están en el norte, 30% en el centro y 30%
en la región sur.
Dado el patrón de movimiento de los camiones, la compañía está interesada en
saber lo siguiente:
1. ¿ Qué proporción de camiones se encontrará en cada región después de un mes ?¿
Después de dos meses ?
2. Qué proporción de los camiones estará en cada región en un período “largo” ?
DIAGRAMA DE ÁRBOL
Probabilidad
Ubicación Ubicación Ubicación de cada
en el en el en el ubicación
mes 0 mes 1 mes 2 enel mes 2
0.2
0.3
Norte Centro 0.06
0.5
0.2
Norte 0.04
Sur 0.10
Norte 0.09
15
Norte0.3
0.3Centro
0.4
0.3
0.2
0.4
Sur
0.4
Centro 0.12
Sur 0.09
Centro 0.20
Sur 0.20
La probabilidad de estar en el norte en cada uno de los tres meses (mes 0, 1 y 2)
está dada por:
( 0.2 ) * ( 0.2 ) = 0.04
Para determinar la probabilidad de que un camión se encuentre en el norte después
de dos meses, se suman las tres probabilidades de encontrarlo en el norte:
p ( estar en el norte dado que se encontraba en el norte en el mes 0 ) = 0.04 + 0.09 +
0.10 = 0.23 p ( estar en el centro dado que se encontraba en el norte en el mes 0 ) =
0.06 + 0.12 +
0.20 = 0.38 p ( estar en el sur dado que se encontraba en el norte en el mes 0 ) = 0.10 +
0.09 + 0.20 = 0.39
9
Utilizando notación matricial, los cálculos para el mes 1 son:
N C S
N C S N 0 0. 0 N C S
1 0 0 * C 0 = 0. 0 0. 0. 0. S
0 0. 0
Vector de probabilidad Matriz de transición Vector de probabilidad
para comenzar en el de un mes después de un mes norte
Para el segundo mes, los cálculos son:
N C S
16
N C S N 0 0. 0. N C S
0. 0. 0. * C 0 =0. 0. 0.2 0.3 0.3 S
0 0. 0.
Vector de probabilidad Matriz de transición Vector de probabilidad
Después de un mes de un mes después de dos meses
N C S N C S N C S
N C S 0 0. 0 0 0 0 NN0.2 0.3
0.3 0. 0 0 0 0
0. 0. 0. * C 0 0. 0 S* 0 0 0 0 C =
S
17
ME 1 N C S MES 2
S
N C S N 0 0 0 N C S
0.2 0.3 0. C 0. 0. 0. 0.23 0.37 0.3
3 6 41 S 0. 0. 0. 6 7 87
A partir del mes 5 se observa estabilidad en los resultados, lo que significa que se
alcanzó el estado estable, es decir, ya se tienen definidas la tendencias a largo plazo,
pues ya no existe variabilidad en los resultados.
Con estos datos se puede dar respuesta a la segunda pregunta planteada. Esto
significa que a largo plazo, el 23.76% de todos los camiones se encontrará en el norte, el
37.62% estarán en la región central y el 38.61% en el sur.
Otro ejemplo:
18
Dado que los televidentes normales con bastante frecuencia tenderían a cambiar de
cadena al principio de la temporada con el objeto de ver programas nuevos o de volver a
ver programas antiguos, Goldman ha decidido esperar a que se estabilice la proporción de
televidentes que ven un programa determinado. Le ha pedido a Bill Washington, uno de
sus empleados, que estudie el periodo en el que cada cadena estará ofreciendo un
programa nuevo para que determine cuales serán las proporciones de televidentes finales.
Si Bill puede predecir estos valores, Mark Goldman estará en posibilidades de tomar una
decisión con respecto a un nuevo programa de la NBS, “Investigaciones en la
Universidad”, sin tener que esperar hasta que las preferencias de los televidentes se
vuelvan obvias a través de los datos de los “ratings” o recuentos de tasa de audiencia.
Bill supone que la selección de un televidente específico se ve influenciada más que
nada por el programa más reciente que ha observado en ese periodo y que las
proporciones finales en realidad son valores de estado estacionario. Para ello ha
elaborado la siguiente matriz de transición utilizando datos recopilados en años anteriores
y referentes a la forma en que los televidentes tienden a cambiar, de una cadena a otra,
semana a semana, para el tipo de programas que se considera. Dicha matriz es la
siguiente:
N C A
NB BS0. 0.B 0.B
CBS 0.2 0.4 0.4
ABC 0.3 0.3 0.4
S 2 2 6
En esta matriz, los valores son la fracción de televidentes que verán el programa de
cada cadena durante esta semana, dada la cadena que vieron la semana anterior.
Si la cadena NBS pretende proyectar el programa mencionado “Investigaciones en
la Universidad”, la cadena CBC el programa “Cuatro es una multitud” y la cadena ABS el
programa “Los demonios de Danny”, ¿ qué porcentaje de televidentes se calcula que
observará cada uno de ellos ?
Solución:
Sea S1 = % de televidentes que verán la cadena NBS en el momento del estado estable.
Sea S2 = % de televidentes que verán la cadena CBC en el momento del estado estable.
Sea S3 = % de televidentes que verán la cadena ABS en el momento del estado estable.
=
S1 S2 S1 S2 20. 0. 0. S2 = 0.4 S1 + 0.3
S2 + 0.2
SS30. 0. 0.S3 = 0.4 S1 + 0.4 S2 + 0.6
2 . .
Esto nos lleva a: 0.8 S1 + 0.3 S2 + 0.2 S3 = 0
19
0.4 S1 0.7 S2 + 0.2 S3 = 0
0.4 S1 + 0.4 S2 0.4 S3 = 0
Además: S1 + S2 + S3 = 1
Resolviendo este sistema de ecuaciones se tiene:
S1 = 1 S2 S3
0.8 (1 S2 S3) + 0.3 S2 + 0.2 S3 = 0
0.4 (1 S2 S3) 0.7 S2 + 0.2 S3 = 0
0.4 (1 S2 S3) + 0.4 S2 0.4 S3 = 0
Esto significa que cuando se alcance el estado estable, es decir, cuando los
televidentes tengan bien definidas sus preferencias por las cadenas de televisión, el
22.72% verán la cadena NBS, el 27.27% de televidentes verán la cadena CBC y el 50%
restante verán la cadena ABS.
Estos resultados tan bajos para las expectativas de la cadena NBS hacen que Mark
Goldman pueda tomar decisiones al respecto.
Uso de sofware
20
Para comenzar a armar un problema de este tipo es necesario ingresar
los campos:
• Título del problema (Problem Title)
• Número de estados (Number of States)
21
Para obtener las probabilidade estado estable, reemplace cualquiera de las ecuaciones en
(I) con 25(1)= 1 y resuelva las n nuevas ecuaciones sinultanease
La plantilla vacía representa una matriz con las relaciones entre los estados (State), sus
probabilidades iniciales (Initial Prob.) y el costo de cada uno de ellos (State Cost).
En el menú Resolver y analizar (Solve and Analyze) tenemos las opciones de Resolver los
estados completos (Solve Steady State) o mostrar el Proceso de Markov por pasos (Markov
Process Step)
22
Resolviendo el ejercicio paso a paso
Regresando a la matriz inicial y tomando la segunda opción del menú Resolver y analizar
(Solve and Analyze) tenemos una ventana que nos permite controlar las iteraciones del
proceso:
23
Podemos observar el Número de periodos procesados (The Number of Time Periods from
Initial).
Para el periodo dos (recuerde pulsar en NEXT PERIOD seguido del botón OK):
24
25
CONCLUSIÓN
REFERENCIAS
Taha, H. (2011). Investigación de operaciones. Editorial Pearson. § Hillier,
F. & Lieberman, G. (2015). Investigación de operaciones. McGraw-Hill
A.T. Bharucha-Reid (1960). “Elements Of The Theory of Markov Processes
Emilia-Saez-Cardenas. (2022). Cadenas de Markov Investigación de
Operaciones 2. dokumen.tips. https://dokumen.tips/documents/cadenas-de-
1
markov-investigacion-de-operaciones-2.html?page=1
And Their Applications”. McGraw Hill Series in Probability and Statistics.
Eugenia-Quijano. (2022). Investigación de Operaciones II UNIDAD 4
Cadenas de
Markov. dokumen.tips. https://dokumen.tips/documents/investigacion-de-
operaciones-ii-unidad-4-cadenas-de-markov.html?page=29