Teoria de Colas Tutorial
Teoria de Colas Tutorial
Teoria de Colas Tutorial
Teora de Colas.
Grupo ROGLE
Pgina 1 de 66
Teora de Colas
Contenido
1. Introduccin .............................................................................................................................. 5
Pgina 2 de 66
Teora de Colas
6. Problemas ................................................................................................................................ 44
Pgina 3 de 66
Teora de Colas
Pgina 4 de 66
Teora de Colas
1. Introduccin
Todos hemos experimentado en alguna ocasin la sensacin de estar perdiendo el tiempo al
esperar en una cola. El fenmeno de las colas nos parece natural: esperamos en el coche al estar en un
tapn, o un semforo mal regulado, o en un peaje; esperamos en el telfono a que nos atienda un
operador y en la cola de un supermercado para pagar....
Generalmente como clientes no queremos esperar, los gestores de los citados servicios no quieren
que esperemos.... Por qu hay que esperar?
La respuesta es casi siempre simple, en algn momento la capacidad de servicio ha sido (o es)
menor que la capacidad demandada. Esta limitacin se puede eliminar invirtiendo en elementos que
aumenten la capacidad. En estos casos la pregunta es: Compensa invertir?
La teora de colas intenta responder a estas preguntas utilizando mtodos matemticos analticos.
El trmino cliente se usa con un sentido general y no implica que sea un ser humano, puede
significar piezas esperando su turno para ser procesadas o una lista de trabajo esperando para imprimir
en una impresora en red.
clientes clientes
servicio
llegando servidos
clientes que
abandonan
Aunque la mayor parte de los sistemas se puedan representar como en la figura 1, debe quedar
claro que una representacin detallada exige definir un nmero elevado de parmetros y funciones.
Pgina 5 de 66
Teora de Colas
c) Disciplina de cola
Algunos autores incluyen una sptima caracterstica que es la poblacin de posibles clientes.
Tambin es posible que los clientes sean impacientes. Es decir, que lleguen a la cola y si es
demasiado larga se vayan, o que tras esperar mucho rato en la cola decidan abandonar.
Por ltimo es posible que el patrn de llegada vare con el tiempo. Si se mantiene constante le
llamamos estacionario, si por ejemplo vara con las horas del da es no-estacionario.
Los servidores pueden tener un tiempo de servicio variable, en cuyo caso hay que asociarle,
para definirlo, una funcin de probabilidad. Tambin pueden atender en lotes o de modo individual.
El tiempo de servicio tambin puede variar con el nmero de clientes en la cola, trabajando
ms rpido o ms lento, y en este caso se llama patrones de servicio dependientes. Al igual que el
Pgina 6 de 66
Teora de Colas
patrn de llegadas el patrn de servicio puede ser no-estacionario, variando con el tiempo
transcurrido.
La disciplina de cola es la manera en que los clientes se ordenan en el momento de ser servidos
de entre los de la cola. Cuando se piensa en colas se admite que la disciplina de cola normal es FIFO
(atender primero a quien lleg primero) Sin embargo en muchas colas es habitual el uso de la
disciplina LIFO (atender primero al ltimo). Tambin es posible encontrar reglas de secuencia con
prioridades, como por ejemplo secuenciar primero las tareas con menor duracin o segn tipos de
clientes.
En cualquier caso dos son las situaciones generales en las que trabajar. En la primera, llamada
en ingls preemptive, si un cliente llega a la cola con una orden de prioridad superior al cliente que
est siendo atendido, este se retira dando paso al ms importante. Dos nuevos subcasos aparecen: el
cliente retirado ha de volver a empezar, o el cliente retorna donde se haba quedado. La segunda
situacin es la denominada no-preemptive donde el cliente con mayor prioridad espera a que acabe
el que est siendo atendido.
En algunos sistemas existe una limitacin respecto al nmero de clientes que pueden esperar en
la cola. A estos casos se les denomina situaciones de cola finitas. Esta limitacin puede ser
considerada como una simplificacin en la modelizacin de la impaciencia de los clientes.
Es evidente que es preferible utilizar sistemas multiservidos con una nica lnea de espera para
todos que con una cola por servidor. Por tanto, cuando se habla de canales de servicio paralelos, se
habla generalmente de una cola que alimenta a varios servidores mientras que el caso de colas
independientes se asemeja a mltiples sistemas con slo un servidor.
Pgina 7 de 66
Teora de Colas
Se asume que en cualquiera de los dos casos, los mecanismos de servicio operan de manera
independiente.
Un sistema de colas puede ser unietapa o multietapa. En los sistemas multietapa el cliente
puede pasar por un nmero de etapas mayor que uno. Una peluquera es un sistema unietapa, salvo
que haya diferentes servicios (manicura, maquillaje) y cada uno de estos servicios sea desarrollado por
un servidor diferente.
En algunos sistemas multietapa se puede admitir la vuelta atrs o reciclado, esto es habitual
en sistemas productivos como controles de calidad y reprocesos.
2.1.7 Resumen
Las anteriores caractersticas bastan, de modo general, para describir cualquier proceso.
Evidentemente se puede encontrar una gran cantidad de problemas distintos y, por tanto, antes de
comenzar cualquier anlisis matemtico se debera describir adecuadamente el proceso atendiendo a
las anteriores caractersticas.
Una eleccin equivocada del modelo lleva a unos resultados errneos, y en muchos casos no
analizar adecuadamente nos puede llevar a pensar que el sistema no es posible de modelar.
2.2.1 Nomenclatura
Pgina 8 de 66
Teora de Colas
: Congestin de un sistema con parmetros: (,, c)
c
Con el paso del tiempo se ha implantado una notacin para representar los problemas de colas
que consta de 5 smbolos separados por barras.
A / B / X /Y / Z
Pgina 9 de 66
Teora de Colas
y significa que los clientes entran segn una distribucin exponencial, se sirven de manera
determinista con tres servidores sin limitacin de capacidad en el sistema y siguiendo una estrategia
FIFO de servicio.
La notacin anteriormente representada, por general, deja demasiados casos por resolver, pero
es suficiente para los casos ms importantes.
Disear eficientemente consiste, bsicamente, en definir un sistema cuyo coste (de diseo y de
operacin) se justifique por el servicio que da. Dicho servicio se puede evaluar mediante el coste de
no darlo. De este modo al disear se pretende minimizar unos supuestos costes totales.
A partir de los datos que nos suministra la teora de colas se puede obtener la informacin
necesaria para definir el nmero de asientos necesarios en una sala de espera, o la estructura de etapas
de un proceso de atencin al cliente.
Pgina 10 de 66
Teora de Colas
En cualquier caso, para poder tomar decisiones hacen falta datos que la teora de colas puede
dar en alguno de los siguientes tres aspectos:
Estos resultados son vlidos para cualquier problema de colas y por tanto sern utilizados en el
resto de desarrollo.
El nmero de clientes en el instante t, n(t), es el nmero de llegadas que han ocurrido hasta t
menos el nmero de servicios completados hasta t.
L En n p n
n 0
n c P
Lq E n q n
n c 1
Little, en su famosa frmula, establece una relacin entre la longitud de la cola y el tiempo de
espera:
L= W
Lq = Wq
El tiempo de estancia de un cliente en el sistema se relaciona con el tiempo de espera de un
cliente en la cola,
1
W Wq
Pgina 11 de 66
Teora de Colas
El nmero de clientes que por trmino medio se estn atendiendo en cualquier momento es:
r L Lq W Wq
En un sistema de un nico servidor:
L Lq n p n (n 1) p n p n 1 p 0
n 1 n 1 n 1
pb
c
El tiempo de estancia del cliente (i+1) en la cola es:
Wq(i ) S ( i ) T (i ) si Wq(i ) S ( i ) T (i ) 0
i 1
W q
0 si Wq(i ) S (i ) T (i ) 0
donde S(i) es el tiempo de servicio del cliente i, y T(i) es el tiempo que transcurre desde la
llegada del cliente y hasta la llegada del cliente (i+1)
Es mejor, sin embargo, utilizar una tcnica de recogida de informacin asociada a eventos.
Ejemplo
Pgina 12 de 66
Teora de Colas
12
clientes por unidad de tiempo
31
12
clientes por unidad de tiempo
30
40
El tiempo medio de estancia en la cola es de
12
70
El tiempo medio de estancia en el sistema es de
12
De aqu y a partir de la frmula de Little
Pgina 13 de 66
Teora de Colas
12 70 70
L W
31 12 31
12 40 40
Lq W q
31 12 31
En esta seccin se vern las caractersticas de una distribucin de Poisson y como se relacionan
con la distribucin exponencial. Posteriormente se analizan las ms importantes propiedades y algunas
generalizaciones al adoptar tal patrn de llegadas. Se cierra el apartado con argumentos que apoyan el
uso de la distribucin de Poisson. Adoptar la distribucin de Poisson implica que la probabilidad de
que lleguen n clientes en un intervalo de tiempo t es:
(t ) n t
p n (t ) e
n!
El tiempo entre llegadas se define, de este modo, como la probabilidad de que no llegue ningn
cliente:
p 0 ( t ) e t
siendo por tanto una distribucin exponencial.
El uso de este patrn de llegada (o de servicio) tiene, entre otras las siguientes propiedades:
*
Es habitual tambin admitir que el ritmo de atencin de cliente cuando el servidor est ocupado tiene
una distribucin de Poisson y la duracin de la atencin al cliente una distribucin exponencial.
Pgina 14 de 66
Teora de Colas
P2 La probabilidad de que una llegada ocurra entre el tiempo t y t+t es t+o(t), donde
o ( t )
es la tasa de llegada y o(t) cumple lim 0 . De hecho o(t) se podra entender
t o t
como la probabilidad de que llegue ms de uno.
P (t s )
(t s) (t s )
n
e t , s 0, t s
n
n!
P4 Si el nmero de llegadas sigue una distribucin de Poisson el tiempo entre llegadas sigue
una distribucin exponencial de media (1/) y al contrario
(t ) n t
Pn (t ) e Po (t ) e t
n!
P5 Si el proceso de llegada es Poisson, los tiempos de llegada son completamente aleatorios
con una funcin de probabilidad uniforme sobre el periodo analizado.
k!
f (t1 , t 2 ,..., t k / k llegadas en0, T )
Tk
P6 Para conocer los datos que definen un proceso de Poisson solo es necesario conocer el
nmero medio de llegadas
Pr T 1 / T t 0 Pr 0 T t1 t 0
a) Variabilidad de
t
(m(t )) n
Pn (t ) e m ( t ) , m(t ) (s )ds
n! o
b) Llegadas mltiples
Pgina 15 de 66
Teora de Colas
i 1
i
(t ) k ( k )
Pr N (t ) m e t
cm
k!
(k )
donde c m es la probabilidad de que k ocurrencias den un resultado total de m clientes.
Se dice que un proceso estocstico {X(t),t=0,1,...} es un proceso de Markov si, para cualquier
conjunto de n instantes t1<t2<...<tn, la distribucin de X(t) depende nicamente del valor de X(tn-1). Es
decir:
Una cola, con proceso de llegada Poisson-Exponencial de media , y con proceso de servicio
Poisson-Exponencial de media , se puede modelizar como una cadena de Markov continua, donde
en cada intervalo infinitesimal de tiempo puede ocurrir un nacimiento (llegada) o una muerte (salida)
Pr n n 1 en (t , t t ) n t o(t ) n0
Pr n n 1 en (t , t t ) n t o(t ) n 1
Pgina 16 de 66
Teora de Colas
n-1 n n+1
En situacin estacionaria, se puede decir que el balance de flujo alrededor del estado n debe
ser 0 (sino no sera estable). As las probabilidades de entrada en el estado n , deben ser iguales a la
probabilidad de las salidas:
n Pn n Pn n1 Pn 1 n 1 Pn1 n 0
En el origen
0 P0 1 P1
De las anteriores ecuaciones se puede extraer que:
n
i 1
Pn P0
i 1 i
y dado que
P 1
i 0
i
se puede calcular
1
P0 n
i 1
1
n 1 i 1 i
Aunque la resolucin de las anteriores ecuaciones parece complicada, no es estrictamente
necesario conocer cmo se puede resolver para poderlas aplicar. Slo en el caso de que nuestra
realidad no sea aplicable a un problema ya resuelto deberamos profundizar en los diferentes mtodos
que permiten resolver nuestro problema.
Pgina 17 de 66
Teora de Colas
a) Conocer todas las distribuciones para tener un conocimiento claro de las posibilidades de
caracterizacin una determinada distribucin.
c) Ser capaz de calcular la media y la desviacin tpica de un proceso en funcin de una muestra.
Si las posibles ocurrencias son un conjunto finito y uniforme de valores (e.g. el lanzamiento de un
dado perfecto) se conoce como variable Uniforme Discreta. Si la variable se da entre los valores a y
b enteros, la media de la distribucin es (a+b)/2 y la varianza es ((b-a+1)^2-1)/12
Tambin tiene un rango infinito la conocida como Poisson en la que se representan ocurrencias
para un conjunto grande e independiente de eventos distribuidos a lo largo del espacio o del tiempo La
distribucin tiene propiedades matemticas interesantes que la hacen muy utilizada. La media es y la
varianza es tambin .
Cuando las ocurrencias pueden tomar valores dentro de un rango continuo las distribuciones son
de tipo continuo. En teora de colas son especialmente adecuadas para representar intervalos de
tiempo.
Pgina 18 de 66
Teora de Colas
De hecho la distribucin Erlang es una parte de una clase ms amplia que son las distribuciones
gamma. Cada funcin gamma es definida por dos parmetros y . La media es y la varianza es
^2
La distribucin Normal
La distribucin logNormal
Continua-Uniforme
Exponencial
Erlang
Gamma
Weibull
Normal
LogNormal
Pgina 19 de 66
Teora de Colas
Pn 1 - n con
El nmero medio de clientes en la cola es:
L En n p n (1 ) n n (1 ) n n 1
n 0 n 0 n 0
Como
1
n ( )
( n ) 1 1
n 1
n
n0 (1 ) 2
De donde
L
1
L 1 2
W Lq Wq
( )
Pgina 20 de 66
Teora de Colas
N
L'q E q
N q 0
P (n X ) X
Se define r mientras que la tasa de ocupacin del sistema es
c
Cuando se consideran c servidores en paralelo, las tasas de llegada y de servicio pasan a ser:
a ( t ) e t
b ( t ) n e nt
donde
n 1 n c
n
c nc
La probabilidad de que haga n clientes en un sistema de este tipo es:
n
n P0 1 nc
Pn n!n
c n cc! n P0 nc
1
c 1 r n rc r
P0 1
n 0 n! c! (1 ) c
Pgina 21 de 66
Teora de Colas
rc
Lq P0
c!(1 ) 2
Lq
rc
Wq P0
c!(c )(1 ) 2
Y por tanto,
1 1 rc rc
W Wq P0 L r P0
c!(c )(1 ) 2 c!(1 )
2
Para facilitar el clculo de Lq se ha considerado interesante incluir el siguiente baco que relaciona
Pgina 22 de 66
Teora de Colas
n 1
n P0 1 n c c 1 r n r c 1 K c1
Pn n!n n! c! 1
1
P0 nc10 n c 1
cn K r r
c ncc! n P0 n! c! ( K c 1)
1
n0
P0 r c
Lq
c!(1 ) 2
1 K c 1 (1 )( K c 1) K c
L L 1
L Lq r (1 PK ) W Wq
(1 PK ) (1 PK )
Para facilitar el clculo de Lq se ha considerado interesante incluir los siguientes bacos que
relacionan el valor de con Lq para distintos valores de K-c.
Pgina 23 de 66
Teora de Colas
Pgina 24 de 66
Teora de Colas
/ n
Pn n!
c
( / ) i
i 0 i!
rc
Pc c! , r
c
ri
i 0 i!
Lo sorprendente de esta frmula es que es vlida, independientemente del tipo de distribucin del
servicio y por tanto es vlida para M/G/C/C
Pgina 25 de 66
Teora de Colas
L
L r (1 PC ) W
(1 PC )
r n e r
Pn n0 r
n!
1
L W
Y la probabilidad de que un cliente fuera del sistema entre en el periodo t y t+t es:
n t o(t )
n (0M n ) 0n M
n M
n cn 0 n c
nc
Pgina 26 de 66
Teora de Colas
r P
M
n
n
0 1 n c
Pn
c n! c! r
M
n nc
n
P0 cnM
c 1 M
n
P 1 1 P0 r n! r
M
n
n
n c
n
n 0 c c! n c
M
L nPn Lq L r ( M L )
n 1
L Lq
W Wq
( M L) ( M L)
1 n k
n 1
n k
Asumiendo que la llegada de clientes sigue una distribucin de Poisson de media , se puede
decir que:
n P 0n k
Pn 1k 10 n k 1 1
1 P0 n k 1
1 k k 1 1
1
1 1 1, 1
1 1 1
P0 1
k 1 1, 1
1
Pgina 27 de 66
Teora de Colas
L P0 1
1 (k 1) 1k k 1k 1 1k 1 k (k 1)
(1 ) 2
(1 ) 2
1
L Lq
Lq L (1 P0 ) W Wq
La literatura nicamente considera los dos primero tipos de clientes impacientes: los que no se
unen a la cola o las que la abandonan antes de tiempo
El cliente que no se une a la cola lo hace porque ya hay demasiados clientes antes. Si para
todos los clientes demasiados fuera la misma cantidad (k) el problema sera M/M/c/K. Sin embargo
lo normal es que k no sea constante para cada cliente. Por tanto la modelizacin es un poco diferente.
Se puede asumir para ello que el ratio de llegada , se ve afectado por una serie montona decreciente
tal que
n bn 0 bn 1 bn 1
En este caso
n
Pn P0 ( ) n bi 1
i 1
Se puede asumir que este tipo de clientes tiene una cierta probabilidad de irse si hay n clientes
en la cola r(n), donde
r(0)=r(1)=0
Por tanto
Pgina 28 de 66
Teora de Colas
n n
i 1 bi 1
Pn P0 P0 n
i 1 i i 1 r (i )
donde
n
bi 1
P0 1 n
n 1 i 1 r (i )
Pgina 29 de 66
Teora de Colas
3.9.1 M/G/1
Los clientes, en este modelo, siguen llegando con una distribucin de Poisson de media , pero
asumimos que son atendidos por un proceso ms general de duracin media 1/ y de desviacin tpica
.
En 1932 Pollaczek y Khintchine desarrollaron la frmula denominada P-K que permite evaluar
la longitud de la cola media.
2
2 2
L
2 1
De la anterior relacin se extrae directamente el tiempo de estancia en cola.
2
2 2
Wq
2 1
(Curry y Feldman, 2010) proponen una modificacin de la frmula que es bastante interesante
(adems de exacta) pues proporciona una relacin directa entre las colas M/M/1 y las colas M/G/1
permitiendo utilizar tablas ya conocidas.
1 2 z
Wq M G 1 Wq M M 1
2
3.9.2 G/G/1
Cuando la entrada tampoco sigue una distribucin exponencial se puede utilizar la aproximacin
de difusin Kingman para calcular el valor del tiempo de espera en cola en funcin de los coeficientes
Pgina 30 de 66
Teora de Colas
de variacin al cuadrado de la entrada y de la salida. De hecho esta es una aproximacin que adems
es siempre una cota superior.
2 e2 2 s2
Wq G G 1 Wq M M 1
2
2 e2 2 s2 1
Wq G G 1
2 1
Existen otras aproximaciones como la de Kraemer y Langenbach que mejoran la calidad del
resultado. En cualquier caso es interesante notar el efecto que tiene la variabilidad (ya sea a la entrada
o a la salida) en el tiempo de estancia en cola.
3.9.3 G/G/c
2 e2 2 s2 2c 2 1 1
Wq G G c
2 c 1
Las fuentes principales de variabilidad en el tiempo que un elemento ha de estar en el sistema son
las siguientes:
c) Disponibilidad de operarios.
Pgina 31 de 66
Teora de Colas
E[ Fi ]
a
E[ Fi ] E[ Ri ]
A partir de esta definicin Hopp and Spearman desarrollan la expresin para el tiempo de servicio
efectivo y su coeficiente de variacin asociado
E[Ta ]
E Te
a
C C Te C
2 2 2
1 C R a 1 a E[ R ]
2
i i
e s
E TS
A partir de estos valores es posible calcular los tiempos de espera utilizando la aproximacin al
problema M/G/1 ya expresada anteriormente con la siguiente formulacin.
Ca2 Ce2 e
Wq E Te
2 1 e
En esta frmula hay que destacar que la tasa de utilizacin efectiva viene afectada por la
disponibilidad del recurso tras eliminar el tiempo que est parado.
e
a
Hasta este momento siempre se ha considerado que slo un recurso (o conjunto de los mismos) era
limitador de la capacidad de la mquina. Podra ser que tuviramos que considerar que hay dos
recursos asociados a la utilizacin de la mquina: por ejemplo la mquina misma y un operario que le
haga la preparacin antes de empezar a ejecutar.
Desafortunadamente para este caso tan habitual no hay una solucin general y habra que recurrir a
diagramas de estado para calcular todos los parmetros bsicos.
Pgina 32 de 66
Teora de Colas
En ese caso la codificacin de cada uno de los estados es bsica para obtener e interpretar algn
resultado.
La propuesta de Curry y Feldman para un sistema con dos mquinas idnticas y un operador es
que cada estado lo representa una tripleta (n,i,j) donde n es el nmero de trabajos en el sistema, i y j
son el estado de cada una de las mquinas, pudiendo ser 0;s;p. 0 indica que la mquina est vaca y
parada, s indica que la mquina est sometida a un setup y p indica que la mquina est en
produccin.
4. Series y Redes
4.1 Introduccin
En este captulo se realiza una introduccin al tema de las redes de colas. Esta es un rea de
gran inters investigador y de aplicacin, con problemas muy complicados de plantear y de resolver.
Por este motivo se presenta nicamente una introduccin de los conceptos bsicos aunque su
aplicabilidad en el modelado de sistemas de fabricacin es ms que evidente.
Las Redes de Colas se pueden describir como un grupo de nodos (sean k), en el que cada nodo
representa una instalacin de servicio.
Dicha instalacin puede constar de ci servidores (i=,1...k) En el caso ms general los clientes
pueden entrar en cualquier nodo y, despus de moverse por la red, pueden salir en cualquier nodo.
Dentro de las Redes de Colas, se pondr especial inters en las denominadas Redes de
Jakcson. Estas tienen las siguientes caractersticas:
2. Los tiempos de servicio en cada nodo y son independientes y siguen una distribucin
negativa exponencial con parmetro i, que podra ser dependiente del estado
Si aadimos las caractersticas i y ri,0 para todo i estamos en el caso de las Redes de Jackson
cerradas Si no se da el caso anterior el problema se denomina de Redes de Jackson abiertas
Pgina 33 de 66
Teora de Colas
En cualquier caso las primeras son una variante reducida de las segundas.
Pgina 34 de 66
Teora de Colas
i 1
ri
0 i 1
1 j i 1 1 i k 1
rij 1 ik j0
0 en los dems casos
..
.
Figura 6: Sistema en serie.
Los clientes entran en el nodo 1 y salen en el nodo k, despus de pasar por cada uno de los
nodos.
Por tanto es necesario entender como salen los clientes de la primera etapa dada las
caractersticas de la entrada y de la etapa de servicio (i, i, ci)
Se puede demostrar que la salida de los clientes de un sistema M/M/c/ tienen una
distribucin idntica a la de la entrada, es decir Poisson con media . Por tanto una serie se compone
de k M/M/ci/ colas independientes, siempre que la entrada sea Poisson, el servicio sea exponencial
y no haya restricciones de capacidad.
Pn1,n2...nk=Pn1Pn2....Pnk
Pgina 35 de 66
Teora de Colas
c) De cada etapa i un cliente se mueve a otra etapa con probabilidad r ij, y al exterior con
probabilidad ri,0
i ri , 0
i
i
ri ,0
ri , j ri ,0
ri , 0
i i
ri , 0
El ratio de llegada i a cada etapa se obtiene mediante las denominadas ecuaciones de trfico
k
i i j r ji
j 1
de donde R
y por tanto ( I R ) 1
i
Definiendo i La probabilidad de que en el estado estacionario haya n i clientes en el
ci i
nodo 1, n2 en el nodo 2, etctera:
n
r i
Pn1n 2n 3... nk i Poi
a i ( ni )
donde
n
i ni ! ni ci ri i
ri
i
a (ni ) ni ci
ci ci ni ci
Po ,i / Poi a i ( ni )
1
concretamente si ci=1 i
Pgina 36 de 66
Teora de Colas
n n nk
Pn1n2 ...n k (1 1 ) 1 1 (1 2 ) 2 2 ...(1 k ) k
Cada etapa se comporta de modo independiente. Esto no significa que se comporte como
Poisson, aunque el comportamiento del valor medio permite considerar cada cola como una M/M/1
independiente. As pues:
i Li
Li Wi
1 i i
W
L i
Es una generalizacin bastante evidente que cada cliente tenga una matriz de ruta R(t), siendo
t=1...n el tipo de cliente.
Para abordar este problema en primer lugar hay que resolver las ecuaciones de trfico de
modo separado
( t ) (t ) ( I R (t ) ) 1
n
( t )
t 1
Todos los resultados anteriores son ahora aplicables. Adems la presencia media de un cliente
de tipo t se puede calcular como:
(t)
(t) i
Li (1) ( 2) ( n)
Li
i i ... i
Pgina 37 de 66
Teora de Colas
ri,j
k k
i i i j r ji j r ji j
j 1 j 1
Como estas ecuaciones son redundantes es posible asignar un valor cualquiera a cualquier i
(por ejemplo 1=1). Las probabilidades en el estado estacionario son:
1
Pn1n2 ...n k 1n1 2n 2 ... knk
G(N )
n1
donde G ( N ) 1 2n 2 ... knk
n1 n 2 ...nk N
1 k
in(i )
G( N )
Pn1n2 ... n k
i 1 a i ( n i )
Donde
ni ! ni ci
a i ( ni ) ni ci
ci ci ! ni ci
k
in (i )
G( N )
n1 n 2 ... nk N i 1 a i ( n i )
El mtodo anterior de analizar las redes de colas cerradas de Jackson tiene un coste
computacional elevado. Fundamentalmente debido a la coste de calcular G(N)
Pgina 38 de 66
Teora de Colas
El mtodo del Anlisis del Valor Medio (MVA), que se explica en este apartado, no requiere
calcular G(N) reduciendo de este modo el citado coste computacional. Este mtodo se basa en que la
frmula de Little es aplicable a travs de toda la red.
Sea
1 Li ( N 1)
Wi ( N )
i
donde
Wi(N)= Tiempo medio de espera en el nodo i para una red con N clientes
Li(N)=i(N)Wi(N)
Si pudiramos calcular i(N), podramos evaluar Wi y Li empezando desde N=0 hasta N=n de
modo recursivo.
j
vi v j r ji i
j 01
Dado que la anterior relacin es redundante porque la red es cerrada, establecemos vl=1 y
resolvemos para los otros nodos.
En general
i ( N )
Pi (n, N ) Pi (n 1, N 1)
i ( n) i
donde
ai ( n) j ci
i ( n)
a1 (n 1) j ci
Pgina 39 de 66
Teora de Colas
k
vi v j r ji asumiendo vl=1
j 1
c1 2
1
Paso3.1 Wi ( n) (1 Li (n 1) (ci 1 j ) pi ( j , n 1) i
ci i j 0
n
Paso 3.2 l ( n) k
con vl 1
v W ( n)
i 1
i i
Paso 3.3 i ( n) l ( n ) v i i, i l
i ( n)
Paso 3.5 Pi ( j, n) Pi ( j 1, n 1) i 1..k , j 1..n
i ( j ) i
Pgina 40 de 66
Teora de Colas
5. Simulacin
Algunos problemas no se pueden resolver mediante mtodos analticos. Entre otras razones
podran citarse la existencia de patrones no normalizados de entrada y de servicio, una gran
complejidad del sistema a modelar o la naturaleza de la disciplina de cola.
Adems , en ocasiones, los resultados analticos son para un estado estacionario que nunca se
alcanza, porque el sistema se interrumpe antes de abandonar el estado transitorio.
En estos casos el anlisis de las colas mediante simulacin puede ser una buena tcnica para
encontrar el resultado.
Hay que destacar, en cualquier caso que, si existen los modelos analticos, stos se deberan
utilizar. Aunque la simulacin permite resolver, o aproximar la resolucin, de muchos problemas
intratables no es la panacea dado que resolver mediante simulacin es equiparable a realizar una
experimentacin. Por tanto hay que utilizar todas las herramientas asociadas al diseo y anlisis de
experimentos: Recogida y Anlisis de Datos, realizacin de la experimentacin, anlisis y consistencia
de resultados, etc.
Otro de los defectos del uso de la simulacin frente a los mtodos analticos, se presenta cuando el
objetivo es el diseo de un sistema y no su evaluacin. En ese caso el anlisis por simulacin no
permite utilizar tcnicas de optimizacin convencionales, aunque hay que admitir que algunas
herramientas de simulacin incorporan tcnicas de optimizacin estocstica para resolver este tipo de
problemas
b) Simulacin.
Dado que estamos interesados en modelizar sistemas estocsticos, los datos de entrada deben
representar del modo ms fiable posible la realidad. En ocasiones se usan datos extrados de la misma
para reconocer la estructura de los datos de entrada. La ejecucin de la simulacin actualmente se
realiza mediante paquetes informticos avanzados (los denominados VIMS). El anlisis de resultados
tiene que ver con el clculo de la efectividad del sistema mediante las tcnicas estadsticas apropiadas.
Pgina 41 de 66
Teora de Colas
Adems la validacin del modelo es una exigencia que muchas veces se olvida al realizar modelos, q
consiste en comprobar que el sistema reacciona como lo hace la realidad.
Los dos mayores problemas en la modelizacin de los datos de entrada son la seleccin de la
familia de distribuciones estadsticas y una vez estimada la familia estimar los parmetros que definen
la funcin de los diferentes entradas.
En muchos casos los paquetes de simulacin suelen llevar una herramienta de ajuste estadstico.
Cuando esto no ocurra deberemos recurrir a las diferentes tcnicas estadsticas para definir tanto las
familias como los parmetros.
Por tanto para obtener conclusiones es necesario disear y ejecutar experimentos de una manera
lgica y comprehensiva.
Pgina 42 de 66
Teora de Colas
Los ms importantes paquetes de simulacin llevan incorporado herramientas que realizan estos
ejercicios mediante los cuales se pueden calcular los anteriormente citados intervalos de confianza.
El primer y fundamental paso en la validacin es verificar que el programa hace lo que est
previsto que haga. Otro paso es definir el grado de credibilidad, es decir hasta que punto los que van a
usar el modelo consideran que el mismo tiene una utilidad y representa la realidad en la media que
nos interesa. Para ello es necesario que los objetivos del estudio, las medidas de rendimiento y el nivel
de detalle debe pactarse y mantenerse en el nivel ms simple posible.
Cuando sea posible, los resultados de las simulaciones se deben comprobar con la realidad. Si esta
no estuviera disponible habra que intentar reproducir modelos tericos con soluciones conocidas
mediante mtodos analticos.
Pgina 43 de 66
Teora de Colas
6. Problemas
c) Cul es el tiempo medio que un cliente pasa en la biblioteca hasta que le han
atendido?
d) Cul es el tiempo medio que un cliente pasa en la cola esperando a que le atiendan?
Pgina 44 de 66
Teora de Colas
segn una distribucin de Poisson con una media de 60 por hora y el tiempo en que un cliente es servido
vara segn una distribucin exponencial de media 150 segundos. Asumiendo el estado permanente del
sistema, cul es el tiempo medio de estancia del cliente en el restaurante hasta que ha sido atendido? El
gerente de Panis ha credo ahora que es preferible una nica cola para distribuir finalmente a los tres
servidores y por tanto las barreras son eliminadas. cul es el tiempo de espera de este modo?
Pgina 45 de 66
Teora de Colas
a) cuntos camiones habr por trmino medio habr en cualquiera de los dos sistemas?
b) Cunto tiempo pasar cada camin en el taller en cualquiera de los dos sistemas?
c) El encargado estima que cada minuto que un camin pasa en el taller reduce los
beneficios en 2 euros. Se sabe que el sistema de dos camiones en paralelo tiene un
coste de un euro por minuto. Qu debera costar el segundo sistema para que no haya
diferencia econmica entre los dos?
Pgina 46 de 66
Teora de Colas
Pgina 47 de 66
Teora de Colas
Pgina 48 de 66
Teora de Colas
Pgina 49 de 66
Teora de Colas
b) Si se pone una caja ms que el mnimo requerido, cul es el tiempo medio de espera
en la cola? Cunta gente habr en cajas? cunta gente habr en todo el
supermercado?
Pgina 50 de 66
Teora de Colas
6.20 Mantemiento
Se desea que dos mquinas estn operativas el mximo tiempo posible. La mquina se rompe de
acuerdo con una exponencial de media . Una vez rota, una maquina tiene una probabilidad r12 de ser
reparada localmente por un responsable de mantenimiento que trabaja con una media de tiempo de 2.
Con probabilidad 1- r12 la mquina debe ser reparada por un especialista, que tambin trabaja segn
una exponencial de media 3. Despus de una reparacin local, existe una probabilidad r23 de que la
mquina requiera un servicio especial. Despus del servicio con el especialista la mquina siempre se
pone a trabajar. Se desea saber cmo se distribuyen los tiempos de estancia de las mquinas en
reparacin.
Pgina 51 de 66
Teora de Colas
abandona el restaurante. Se tarda 4 minutos en servir los rollitos y 5 minutos en servir el pollo frito, el
tiempo de servicio es exponencial. Cunta gente hay por trmino medio en el restaurante? Cul es el
tiempo medio de espera en cada ventanilla? Si alguien quiere los dos mens Cunto tiempo pasa en
el restaurante?
Pgina 52 de 66
Teora de Colas
g) Sobre el caso anterior Qu opinin le merece que vaya uno de los trabajadores de la
tercera seccin a ayudar al mecnico cuando haya dos o ms mquinas estropeadas?
Debe sustentar la opinin con datos, suponga para ello que el trabajador de la tercera
seccin se comporta como un mecnico ms, cuando trabaja como tal.
h) (1 punto) Cul sera en el caso anterior la probabilidad de que hubiera ms de una
mquina estropeada?
Resolucin
a) Cul es el nmero de puertas que habr en cada etapa, incluyendo las puertas en las
mquinas y las que estn siendo procesadas por los operarios?
En primer lugar el sistema dejara de ser una serie de colas convencional porque se
limita la capacidad de una de ellas. Esta sera una cola M/M2/7 Pero para saber cmo
afectara lo mejor es saber el porcentaje de veces que el almacn estara lleno
P(n=7)=4,1%. Por tanto durante un 3,8% de las ocasiones la primera etapa no podra
trabajar al estar bloqueado el sistema posterior, dado que la primera etapa trabaja al 95%,
la probabilidad de que el bloqueo del sistema afecte a la produccin total es alta.
c) Si en el sistema original el tiempo medio de entrega de una puerta es de 5 das.
Cuntas puertas hay?
70000 puertas al ao son 291,67 puertas al da. Por tanto habra que trabajar 2,6 horas
ms al da. El almacn anterior a la primera etapa aumentara en 71 puertas cada da, y se
reducira en la misma cantidad cada noche(o cuando se hagan las horas extra), pasando a
Pgina 53 de 66
Teora de Colas
Este es un sistema con fuente finita. El nmero de mquinas que estarn siendo
reparadas por trmino medio es de 0,67. Por tanto se tienen 2,33 mquinas trabajando, y
por tanto es ms que suficiente ya que 2,33*140>208.
0 1 2 3
f) (1 punto) Sobre el caso anterior Qu porcentaje de tiempo slo hay una mquina
trabajando? Y ninguna? Qu ocurre con los almacenes durante este tiempo que hay
menos de dos mquinas trabajando?
Slo hay una mquina trabajando cuando hay dos estropeadas P2=13,4%
A priori, si hay dos mquinas estropeadas, el ritmo al cual pasan puertas a la tercera
etapa es de 140 puertas/da. Con dos operarios es posible abastecer 192 puertas al da,
con lo que no pasara nada. Adems el porcentaje de veces que hay dos mquinas
estropeadas es muy bajo en el nuevo sistema.
h) (1 punto) Cul sera en el caso anterior la probabilidad de que hubiera ms de una
mquina estropeada?
Pgina 54 de 66
Teora de Colas
0 1 2 3
0,56536
P0 409
0,35335
P1 256
0,07361
P2 512
0,00766
P3 824
Pgina 55 de 66
Teora de Colas
La primera etapa se realiza simultneamente por dos equipos de trabajo, con un ritmo promedio
cada uno de ellos de una puerta cada 20 minutos. La segunda etapa la realiza un equipo de trabajo con
un tiempo de ciclo promedio de 11 minutos por puerta. La tercera etapa requiere del uso de otra
mquina con un tiempo de ciclo promedio de 10 minutos por puerta.
Por ltimo la cuarta etapa es de preparacin final. Como es un trabajo principalmente manual, que
realiza un nico operario, tiene un tiempo de ciclo de 18 minutos por unidad, y se dispone de tantos
trabajadores como se requieran, pues irn viniendo de otras secciones siempre que haya una puerta por
preparar.
a) Modele el problema segn teora de colas, estableciendo los parmetros bsicos, asumiendo
tiempos promedio exponenciales.
e) Cul es el tiempo promedio previsto de entrega de una puerta? Si le dicen que el tiempo de
entrega promedio es de 5 das. A qu puede ser debido?. Proponga un mecanismo de
correccin.
f) Cual ser el efecto sobre la cantidad de puertas en la primera etapa si en lugar de dos equipos de
trabajo con tiempos de ciclo como los citados se establece un nico equipo ms eficiente con un
tiempo de ciclo de 9 minutos por unidad.
Pgina 56 de 66
Teora de Colas
Solucin:
El problema propuesto es una serie de colas con una entrada =386 puertas/da.
a) La primera etapa es una cola M/M/2 con =386 puertas/da. y =22,5 puertas/da
La segunda etapa es una cola M/M/1 con =386 puertas/da. y =40,9 puertas/da
La tercera etapa es una cola M/M/1 con =386 puertas/da. y =45 puertas/da
La cuarta etapa es una cola M/M/ con =386 puertas/da. y =25 puertas/da
b) L1= 6,53 L2= 17 L3=6,07 L4= 1,54 LT=31,14 puertas
e) El tiempo que tardar ser el de salir la primera 0,806 das ms el que tardan en salir las 29
restantes. 0,806 + 29/ =0,806+0,751=1,557 das = 11,68 horas
f) Si en realidad tardan 5 das en salir, es porque en el sistema hay puertas de ms. La cantidad de
puertas que hay es L=(W-29/)=162 puertas de ms.
Pgina 57 de 66
Teora de Colas
6.25 Juguetes KP
Una empresa de transformados plsticos tiene 4 secciones (A,B,C,D). Los productos que fabrica se
pueden clasificar en 5 categoras, con demandas anuales diferentes.
Los productos de categora 1 tienen una demanda anual de alrededor de 500 unidades, y sus
especificaciones los harn circular por la seccin A, luego la seccin B y por ltimo la seccin C.
Los productos de categora 2 tienen una demanda anual de alrededor de 3000 unidades, y sus
especificaciones los harn circular por la seccin A, luego la seccin B y por ltimo la seccin D.
Los productos de categora 3 tienen una demanda anual de alrededor de 2000 unidades, y sus
especificaciones los harn circular por la seccin B y la seccin D.
Los productos de categora 4 tienen una demanda anual de alrededor de 2000 unidades, y sus
especificaciones los harn circular por la seccin A y la seccin C.
Los productos de categora 5 tienen una demanda anual de alrededor de 1000 unidades, y sus
especificaciones los harn circular por la seccin B y la seccin C.
Sabiendo que el ritmo de produccin por hora en una mquina de tipo A es de 2 unidades, el de B
es de 2 unidades, el de C es de 4 unidades y el de D es de 2 unidades por hora. Sabiendo que el ao
tiene 220 das laborables de 8 horas cada uno y asumiendo tiempos de servicio exponenciales.
Preguntas
e) Suponga que en la seccin B hacen falta dos mquinas. La experiencia con esas mquinas
indica que requieren un cierto mantenimiento con distribucin exponencial con media cada 5
das y que el tiempo que dura dicho mantenimiento se distribuye exponencialmente con
media de 1 da. Cuntas mquinas hacen falta para que haya al menos dos mquinas
funcionando el 90% del tiempo?.
Resolucin
Pgina 58 de 66
Teora de Colas
rAB=0,64
A B
rBC=0,23
rAC=0,36 rBD=0,77
C D
rC0=1 rD0=1
Nodo
A 5500 5500 3520
B 3000 6500 3520
C 0 3500 7040
D 0 5000 3520
El nmero de mquinas a ubicar en cada nodo son:
Nmero de
Nodo Mquinas
A 2
B 2
C 1
D 2
Para calcular el tiempo esperado medio en el sistema es necesario conocer la cola en cada nodo:
Nm
Nodo Mq L
A 5500 5500 3520 2 0.781 4.01
Pgina 59 de 66
Teora de Colas
L
i
i
4.01 12.51 0.99 2.86
0, 0017 aos 2,99horas
i
i 5500 6500 0 0
Parece evidente que se debe comenzar calculando con 3 mquinas, porque cada mquina est un
20% del tiempo en reparacin lo que hace imposible que M= 2 sea una respuesta vlida.
Con M=3 la probabilidad de que haya 0 1 mquina siendo reparadas, de acuerdo con las
frmulas, es
r P
M
n
n
0 1 n c
Pn
c n! c! r
M
n n c
n
P0 cnM
c 1 M
n
P 1 1 P0 r n! r
M
n
n
n c
n
n 0 c c!
n c
con M=3
p0 0.530
p1 0.318
p2 0.127
p3 0.025
Lo que indica que la probabilidad de que haya 0 1 mquina siendo reparada es 0.848
Con M=4
Pgina 60 de 66
Teora de Colas
P0 0.398
P1 0.319
P2 0.191
P3 0.076
P4 0.015
Lo que indica que slo en el 9.1% de los casos 3 4 mquinas estarn reparndose y por tanto se
cumplir la restriccin impuesta.
Pgina 61 de 66
Teora de Colas
Actualmente se reconocen 3 tipos de llamadas que se reciben en tres telfonos distintos. De tipo 1
se reciben 40 llamadas a la hora, de tipo 2 se reciben tambin 40 llamadas a la hora y de tipo 3 se
reciben en condiciones normales 20 llamadas a la hora.
El tiempo que se tarda en atender una llamada de tipo 1 es de 3 minutos igual que las llamadas de
tipo 2. Las llamadas de tipo 3 requieren una atencin en promedio de 5 minutos cada una.
Usted est diseando un nuevo sistema de atencin telefnica, que atendera a todos los clientes
con un nico nmero de telfono. Un sistema informtico discrimina el destino de la llamada
mediante una operacin que dura aproximadamente 20 segundos en promedio.
Una vez el sistema informtico decide el destino, tiene una probabilidad del 5% de equivocarse.
En ese caso el operador que recibe la llamada, enva sta al centro adecuado para que sea atendido.
Preguntas
a) Cuntos operadores pondra en la seccin que atiende llamadas de tipo 1?
b) Cuntos operadores pondra en la seccin que atiende llamadas de tipo 3?
c) Cuntas llamadas habra en promedio en la seccin que atiende llamadas de tipo 3 si pusieran
los operadores estrictamente necesarios?
d) Cuntas llamadas habra en promedio en la seccin que atiende llamadas de tipo 3 si pusieran
uno ms de los operadores estrictamente necesarios?
e) Cul es tiempo medio que un cliente estara en el sistema si pusiera los operadores
estrictamente necesarios?
f) Cul es tiempo medio que un cliente estara en el sistema si pusiera en cada seccin uno ms de
los operadores estrictamente necesarios?
g) Se le plantea una nueva alternativa. Consiste en hacer que todos los operadores atiendan todas
las llamadas, aunque en ese caso el tiempo de atencin de cada llamada es el doble del indicado
ms arriba para cada tipo. En este caso no hace falta discriminador. Cuntos operadores hacen
falta?
h) En el caso anterior, Cul es el tiempo medio de estancia en el sistema?
Si disponemos 12 operadores el tiempo medio de estancia en el sistema es de 0,25 horas.
Resolucin
a) Cuntos operadores pondra en la seccin que atiende llamadas de tipo 1?
El sistema es una red abierta de Jackson con las siguientes transiciones.
El Modelo representado es el siguiente.
Pgina 62 de 66
Teora de Colas
0 0,395 0,395 0, 21
0 0 0, 0333 0,0167
R
0 0, 0333 0,0167
0 0,025 0,025 1 (20)
rE,1=0,395
r1,2=0,0333
r2,1=0,0333
100
rE,2=0,395
E(180) 2 (20) r1,3=0,01666
r3,2=0,025
r1,2=0,0166 r3,2=0,025
rE,3=0,21
3 (12)
Pgina 63 de 66
Teora de Colas
No todo el mundo que entra en el Parque sigue el mismo camino, pero su inters radica en saber
cunto tardar usted en alcanzar el interior del Parque.
Entran aproximadamente 2400 coches por hora, para ser atendidos por 20 cajas en paralelo que
tardan en cobrar aproximadamente 29 segundos por cliente.
De todos los coches que entran un 18% van al aparcamiento VIP (al que usted no va ni ir y por
eso no sabe cmo funciona). El 82% restante va al aparcamiento convencional que, de un modo muy
eficiente es capaz de aparcar los coches, de uno en uno, con un tiempo de ciclo promedio de 3,5
segundos por coche.
Mediante otros medios (trenes y autobuses) se acercan junto con los clientes en coche particular
nuevos grupos de clientes.
Se calcula que en el parque que hemos elegido para depositar (y perder) nuestros ahorros, entran
aproximadamente 25.000 personas al da en 6000 grupos (los grupos son importantes porque compran
las entradas de modo conjunto). De los 6000 grupos slo 3500 grupos compran las entradas en taquilla
(los otros ya las compraron por agencia o llevan un pase de varios das comprado anteriormente).
Todos los clientes llegan aproximadamente en las 3 primeras horas de apertura del parque.
Pgina 64 de 66
Teora de Colas
Los que han de pagar tendrn que hacer 6 colas para pagar en 30 cajas (cada cola alimenta a 5
cajas). En cada caja tardarn en promedio 92 segundos en atenderles.
Tras pagar queda la ltima cola donde cada cliente pasa de modo individual, y pasan todos: los que
acaban de comprar y los que venan con ticket precomprado, por el detector de comidas y bebidas.
Estos son 12 carriles en paralelo, cada uno con su propia cola, que tardan 5 segundos en dejar pasar a
cada cliente.
a) Si cada coche mide 4 metros, cuantos metros de carretera hacen falta para que quepan en
promedio todos los coches que se pondrn en cola delante de las cajas de aparcamiento.
b) Cunto tiempo se tarda en hacer la cola para aparcar el coche, una vez haya pagado el
aparcamiento?
c) Cunto tiempo tardaremos en conseguir nuestra entrada desde que hemos aparcado el coche,
considerando que tarda 5 minutos desde que aparca hasta que llega a la cola?
d) Cunta gente habr como usted haciendo cola para pagar?
e) Cunto tiempo tardaremos en entrar en el parque desde que aparcamos?.
f) qu ocurrir en el sistema de cajas si el tiempo de atencin en la caja de compra de entradas el
tiempo medio de atencin es de 100 segundos?
g) qu repercusin tendra en la cola posterior dicha alteracin?
h) Con las condiciones de e, cul es la repercusin de aadir un carril adicional de venta de
entradas?
Pgina 65 de 66
Teora de Colas
El tiempo que en promedio tarda el transportador en cubrir un ciclo entero es de 25 segundos por
paleta transportada. Es decir en promedio se tarda 25 segundos en moverse a la mesa de rodillos en la
Se trata de definir la capacidad mnima que han de tener las 5 mesas de rodillos, para que el
sistema no se bloquee. Para ello
a) Definir el problema segn una red de colas de varios productos (suponer que las colas
inicialmente no tienen lmite en la capacidad). Definir la matriz de transicin para cada producto.
b) Calcular los valores de y para cada una de los servidores.
c) Cul ser la cola promedio en cada una de las mesas suponiendo que no tienen lmite de
capacidad?
d) Cmo afectara al sistema que el transportador eligiera para extraer, en cada ocasin, el
producto en cabecera de la mesa de rodillos con ms productos?
e) Cmo cree que afectar a la mquina que alimenta a la mesa 3 que sta tuviera una limitacin
de 10 unidades?. Y si fuera de 5 unidades?
Pgina 66 de 66