Teoria Colas
Teoria Colas
Teoria Colas
TEORÍA DE COLAS
PROCESOS DE LLEGADA POISSONIANOS.
PROCESOS DE NACIMIENTO Y MUERTE
MODELO M/M/1
MODELOS M/M/s M/M/s// M/M/s//k
MODELOS M/M M/M///k
MODELOS M/M/1/c M/M/s/c M/M/c/c
MODELOS NO EXPONENCIALES
MODELOS DETERMINISTAS
SIMULACIÓN DE SISTEMAS
TEORÍA DE COLAS
Elementos básicos de los sistemas de colas.
Patrones de servicio y llegadas.
Análisis de sistemas de colas simples.
Procesos de llegada poissonianos.
Procesos de nacimiento y muerte
Portal Estadística Aplicada: TEORÍA DE COLAS ‐ PROCESOS DE NACIMIENTO Y MUERTE 2
INTRODUCCIÓN
TEORÍA DE COLAS
Una Cola se presenta con frecuencia cuando se solicita un servicio por parte de una
serie de clientes y tanto el servicio como los clientes son de tipo probabilístico.
La primera aplicación de teoría de Colas se debe al matemático danés Erlang sobre
conversaciones telefónicas en 1909 para el cálculo de tamaño de las centralitas.
La Teoría de Colas es una disciplina de Investigación Operativa que se encarga de
poponer modelos para el manejo eficiente de Líneas de Espera.
Una Línea de Espera es una hilera formada por uno o varios clientes que aguardan para
recibir un servicio. Los clientes pueden ser personas, objetos, máquinas que requieren
un mantenimiento, contenedores de mercancías para ser embarcados, elementos de
inventario para ser utilizados, etc.
Una Línea de Espera se forma por un desequilibrio temporal entre la demanda de un
servicio y la capacidad del sistema para gestionarlo.
Los Modelos de Líneas de Espera son muy útiles para determinar cómo operar un
sistema de colas de la manera más eficaz, permiten encontrar un balance adecuado
entre el costo de servicio y la cantidad de espera: Proporcionar demasiada capacidad
de servicio para operar el sistema implica costos excesivos. De otra parte, si no se
cuenta con suficiente capacidad de servicio surgen esperas excesivas con
desafortunadas consecuencias.
El término 'cliente' se utiliza en un sentido general, pudiendo ser una persona, piezas
esperando su turno para ser procesadas, una lista de trabajo esperando para ser
impresas en una impresora de red, etc.
Básicamente la mayoría de los modelos de colas consiste: Los clientes que requieren
un servicio se generan en el tiempo en una fuente de llegada, después entran al
sistema y se unen a una cola.
(a) Cuando hay varios canales en paralelo es conveniente mantener una cola única.
(b) Cuando hay tiempos de servicios muy dispares para los diferentes clientes que
forman la cola conviene establecer canales separados.
(c) Si la cola aumenta hasta cierto limite conviene aumentar la capacidad de los
canales.
Disciplina de cola: Es la regla en el orden que se van a seleccionar los clientes que se
encuentran a la espera de ser atendidos en la cola, existen varias reglas, entre las más
comunes se pueden encontrar:
♠ FIFO (first in first out): Se atiende al cliente en el orden que llegan a la cola, el
primero en llegar será el primero en ser atendido. En los modelos básicos de colas se
supone como normal la disciplina de primero en entrar, primero en salir, a menos que
se establezca de otra manera.
♠ LIFO (last in first out): Consiste en atender primero al que ha llegado de último,
también se le conoce como 'pila'.
♠ RSS (random selection of service): Se selecciona a los clientes de una cola de forma
aleatoria, con algún procedimiento de prioridad o algún otra preclasificación.
♠ Processor Sharing: Todos los clientes experimentan con eficacia el mismo retraso,
ya que comparten entre todos los clientes de la cola la capacidad del sistema
atendiendo a todos por igual.
Capacidad del sistema: Es el número máximo de clientes que pueden estar dentro del
sistema haciendo cola antes de ser atendidos para recibir el servicio, al igual que la
fuente de llegada este número puede ser finito o infinito.
COLA
COSTE DE ESPERA
⇒ SERVICIO
COSTE DE SERVICIO
COSTE DE ESPERA: Coste del tiempo que se produce cuando los clientes tienen que
esperar, como un activo que no está siendo utilizado repercute negativamente,
afectando a la imagen de la empresa o sistema, por lo que se considera de forma
económica.
EQUILIBRIO DE COSTES
♦ Los datos que definen un proceso de Poisson vienen dados por el número medio de
llegadas.
( m(t) )
n
− m(t)
♦ λ puede variar con el tiempo, con lo que pn (t) = e donde
n!
t
m(t) = ∫ λ(s) ds
0
⎢⎣ ⎝ β ⎠ ⎝β ⎠ ⎥⎦
PROCESOS DE POISSON
Si los tiempos entre llegadas/servicios de clientes se distribuyen según una exponencial
Exp(λ ) , el número de llegadas/servicios de clientes hasta un cierto tiempo es un
proceso de Poisson.
⎧1 − e − λ x x≥0
∫
x
La función de distribución: F(x) = f(t) dt = ⎨
−∞ ⎩ 0 x<0
1 1
X ∼ Exp(λ ) → E(X) = V(X) =
λ λ2
La distribución de probabilidad del tiempo que falta para que ocurra el evento es
siempre la misma independientemente del tiempo que haya pasado.
En la mayor parte de los modelos elementales de colas es común que las entradas y
salidas del sistema ocurran según un proceso de nacimiento y muerte. El proceso
explica cómo varía el estado del sistema N(t) al aumentar t.
En este contexto,
N(t) estado del sistema en tiempo t ≡ Número de cliente en el sistema
Nacimiento ≡ Llegada de clientes al sistema
Muerte ≡ Salida de clientes una vez servidos
• La distribución del tiempo que falta para la llegada es exponencial Exp(λ n ) con
n = 0, 1, 2, ... , siendo λ n la tasa de llegada de clientes al sistema.
• La distribución del tiempo que falta para la salida es exponencial Exp(μn ) con
n = 0, 1, 2, ... , siendo μ n la tasa de salida de clientes del sistema.
• Hay independencia entre el tiempo hasta próxima llegada y tiempo hasta próxima
salida.
⎧n → n + 1 un nacimiento
La transición del estado será: ⎨
⎩ n → n − 1 una muerte
Tomando estos supuestos, el proceso es un tipo especial de cadena de Markov de
tiempo continuo.
Los parámetros λ n y μ n son tasas medias en la distribución exponencial, en ocasiones
estos valores son constantes ∀μ
Por ser el sistema estacionario, la tasa media de llegada es igual a la tasa media de
salida para cualquier estado n , es decir, λ n − 1 pn − 1 + μ n + 1 pn + 1 = λ n pn + μ n pn
λ0 λ1
n = 1 : λ 0 p 0 + μ 2 p 2 = λ 1 p1 + μ 1 p1 → p 2 = p0
μ1 μ 2
λ 0 λ1 λ2
n = 2 : λ 1 p1 + μ 3 p 3 = λ 2 p 2 + μ 2 p 2 → p 3 = p0
μ1 μ2 μ3
λ 0 λ 1 λ 2 ... λ n − 1 ∞
por recurrencia, pn =
μ1 μ 2 μ 3 ... μn
p0 ∑p
n=0
n =1
∞ λn −1 λ2 λ1 λ 0 1 1
∑ μn μ 3 μ 2 μ1
p0 = 1 → p0 = ∞ λn −1 λ1 λ 0
= ∞ λn−1 λ1 λ 0
n=0
∑ μn μ 2 μ1
1+ ∑ μn μ 2 μ1
n=0 n=1
El procedimiento para resolver estas ecuaciones no es otro que despejar todas las
variables en términos de una de ellas, siendo la más conveniente p 0 .
λn − 1 λn − 2 λ1 λ 0
Cn =
μn μn − 1 μ2 μ1
Para n = 0 se define Cn = 1
Estos resultados son de estado estable pues al aplicar el limite cuando t tiende a
En un puesto del mercado se compran y venden mentiras, los clientes entran según
un proceso de Poisson de tasa 4 clientes/minuto y la permanencia en el puesto es un
tiempo exponencial con un tiempo medio de 10 minutos. Suponiendo que la capacidad
del puesto es infinita para recibir a clientes, se solicita:
a) Probabilidades de equilibrio
b) Número medio de clientes en el puesto del mercado
Solución:
λ0 4
n = 0 : μ 1 p1 = λ 0 p 0 → p1 = p0 = p 0 = 40 p 0
μ1 1 / 10
λ0 λ1 4x 4 1 1
p2 = p0 = p 0 = 402 p 0 = 402 p 0
μ1 μ 2 1 2
x
2 2!
10 10
λ0 λ1 λ2 4 x 4x 4 1 1
p3 = p0 = p 0 = 403 p 0 = 403 p 0
μ1 μ 2 μ 3 1 2 3
x x
6 3!
10 10 10
1 n
por recurrencia, pn = 40 p 0 ∀ n ≥ 1
n!
∞ ∞
1 1 1
∑ pn = ∑ n! 40n p0 = 1 → p0 = ∞
1
= ∞
1
= e− 40
n=0 n=0
∑ n! 40n 1 + ∑ n! 40n
n=0 n=1
x x2 xn
Se considera el desarrollo de Taylor e = 1 + + +
x
+
1! 2! n!
1 n 40n − 40
Finalmente, pn = 40 p 0 = e n ≥ 0 → distribución de Poisson P(λ = 40)
n! n!
b) El número medio de clientes en el mercado:
∞ ∞
⎛ 40n −40 ⎞ −40 ∞ 40n ∞
40n − 1
L = ∑ n pn = ∑ n ⎜ e ⎟=e ∑ n = 40 e−40 ∑ =
n=0 n=1 ⎝ n! ⎠ n=1 n (n − 1)! n=1 (n − 1)!
∞
40n − 1 ∞
40k
= 40 e −40
∑ (n − 1)!
= 40 e ∑
−40
= 40 e−40 e40 = 40
n=1 k = 0 k!
En una oficina bancaria rural un empleado atiende a los clientes que llegan según
porceso de Poisson de parámetro λ . Los clientes son reacios a esperar demasiado
tiempo, de modo que si hay k personas esperando a ser atendidas con probabilidad
q(k) se quedan en otro caso se van. Una vez que entran en la cola esperan
independientemente según su paciencia o abandonan un tiempo exponencial de tasa
φ . Por otra parte, el tiempo de servicio del empleado es exponencial de tasa μ .
Se quiere modelar el número de personas en la cola como un proceso de nacimiento y
muerte, calulando tasas y probabilidades de la variable estado.
Solución:
X(t) ≡ Número de personas en la cola de espera
Los clientes salen atendidos (muertes) con distribución de Poisson de tasa μ o por
impaciencia con distribución de Poisson de tasa φ
Suponiendo que la variable X(t) se comporta como una cadena de Markov de tiempo
continuo, las probabilidades de transición se rigen por la ecuación:
Se supone que existe una distribución estacionaria, basta con imponer que las tasas de
muerte sean mayores que las taas de nacimiento.
Como los cambios son a estados contiguos, se puede modelar X(t) como un proceso
de Poisson de nacimiento y muerte.
⎧μ + φ k k > 0
Tasa de muerte: μ (k) = ⎨
⎩ 0 k=0
∏
q(k)
por recurrencia, pn = λ n
p0 n≥1
k=0
μ+kθ
1
donde p 0 =
∞⎛ n n−1
q(k) ⎞
1+ ∑ ∏⎜λ
⎜
n=1 ⎝ k=0
⎟
μ + k θ ⎟⎠
Solución:
υ j p j = ∑ pk υk pkj
k≠j
Como los cambios de estado corresponden solo a los artículos (j) y (j + 1) las variables
X 1 (t) y X 2 (t) corresponden a procesos de nacimiento y muerte.
Variable X 1 (t) :
⎧λ 1 X<N ⎧0 j=0
Tasa de nacimiento λ j = ⎨ Tasa de muerte μ j = ⎨
⎩0 X≥N ⎩μ1 j≥1
λ 0 1 ⎛ λ1 ⎞ 1
N = 0 : μ1 p11 = λ 0 p10 → p11 = p0 = ⎜ ⎟ p0
μ1 ⎝ μ1 ⎠
2
λ λ ⎛λ ⎞
N = 1 : λ 0 p + μ2 p = λ 1 p + μ1 p
1 1 1 1
→ p = 0 1 p10 = ⎜ 1 ⎟ p10
1
μ1 μ 2 ⎝ μ1 ⎠
0 2 1 1 2
3
λ λ λ ⎛λ ⎞
N = 2 : λ 1 p + μ3 p = λ 2 p + μ2 p
1 1 1 1
→ p = 0 1 2 p10 = ⎜ 1 ⎟ p10
1
μ1 μ 2 μ 3 ⎝ μ1 ⎠
1 3 2 2 3
N
⎛λ ⎞ 1
por recurrencia, p = ⎜ 1 ⎟ p10 → p10 =
1
⎝ μ1 ⎠
N k
⎛λ ⎞∞
1+ ∑⎜ 1 ⎟
k = 1 ⎝ μ1 ⎠
Variable X 2 (t) :
⎧λ 2 X<N ⎧0 j=0
Tasa de nacimiento λ j = ⎨ Tasa de muerte μ j = ⎨
⎩0 X≥N ⎩μ 2 j≥1
2
λ λ ⎛λ ⎞
N = 1 : λ 0 p20 + μ 2 p22 = λ 1 p12 + μ 1 p12 → p22 = 0 1 p20 = ⎜ 2 ⎟ p20
μ1 μ 2 ⎝ μ2 ⎠
3
λ λ λ ⎛λ ⎞
N = 2 : λ 1 p12 + μ 3 p23 = λ 2 p22 + μ 2 p22 → p23 = 0 1 2 p20 = ⎜ 2 ⎟ p20
μ1 μ 2 μ 3 ⎝ μ2 ⎠
N
⎛λ ⎞ 1
por recurrencia, pN2 = ⎜ 2 ⎟ p20 → p20 =
⎝ μ2 ⎠
k
⎛λ ⎞
∞
1+ ∑⎜ 2 ⎟
k = 1 ⎝ μ2 ⎠
k N−k
N
⎛λ ⎞ ⎛λ ⎞
N
En consecuencia, p = 1 − ∑ p p 1 2
N−k = 1 − p p ∑⎜ 1 ⎟ ⎜ 2 ⎟
1 2
k =0 ⎝ μ1 ⎠ ⎝ μ 2 ⎠
k 0 0
k =0
Solución:
Es una cadena finita por lo que existen probabilidades estacionarias, dado que es
irreductible.
⎧ i −1
⎪ ∏
⎪ k=0
(β + k λ )
p0 0<i≤N
⎪⎪ i! μ i
p i = ⎨ N− 1
∏
⎪
⎪
(β + k λ ) λ i − N
⎪ k=0 p0 N < i ≤ 3N
⎪⎩ i (N − 1)! μ
i
N ∏ (β + k λ ) 3N ∏ (β + k λ ) λ i − N
∑ ∑
k=0 k=0
+
i=1 i! μ i
i=N+1 i (N − 1)! μ i
i
⎛1⎞
Las ecuciones del estado estacionario son: pi = ⎜ ⎟ p 0
⎝2⎠
1 1
k 1−
⎛1⎞
9 9
1 1
∑ pi = ∑ ⎜ ⎟ p 0 = 1 → p 0 = 9
k= 0 ⎝ 2 ⎠ ⎛1⎞
k = 9
⎛1⎞ ⎛1⎞
=
⎛1⎞
2 =
10
2
⎛1⎞
10
k=0
∑
k=0
⎜ ⎟
⎝ 2 ⎠
1−⎜ ⎟ ⎜ ⎟ 1−⎜ ⎟
⎝ 2 ⎠ ⎝ 2 ⎠ ⎝ 2 ⎠
1−⎜ ⎟
⎝2⎠
1
1−
2
El término 'cliente' se utiliza en un sentido general, pudiendo ser una persona, piezas
esperando su turno para ser procesadas, una lista de trabajo esperando para ser
impresas en una impresora de red, etc.
Básicamente la mayoría de los modelos de colas consiste: Los clientes que requieren
un servicio se generan en el tiempo en una fuente de llegada, después entran al
sistema y se unen a una cola.
En determinado momento se selecciona a un cliente de la cola para proporcionarle el
servicio mediante alguna regla conocida como disciplina de la cola. Se lleva a cabo el
servicio que el cliente requiere mediante un Mecanismo de servicio, y después el
cliente sale del sistema de colas.
• En el sistema se puede actuar en las siguientes características:
(a) Ley que rige las llegadas.
(b) Disciplina de la cola.
(c) Ley que determina el servicio (elección entre tipo y número de canales).
(a) Cuando hay varios canales en paralelo es conveniente mantener una cola única.
(b) Cuando hay tiempos de servicios muy dispares para los diferentes clientes que
forman la cola conviene establecer canales separados.
(c) Si la cola aumenta hasta cierto limite conviene aumentar la capacidad de los
canales.
(λ t)n − λt
La distribución de llegadas: pn (t) = e , n = 0, 1, 2, clientes
n!
pn (t) ≡ Probabilidad de que n clientes estén en el sistema en el tiempo t
PROCESOS DE POISSON
Si los tiempos entre llegadas/servicios de clientes se distribuyen según una exponencial
Exp(λ ) , el número de llegadas/servicios de clientes hasta un cierto tiempo es un
proceso de Poisson.
⎧1 − e − λ x x≥0
∫
x
La función de distribución: F(x) = f(t) dt = ⎨
−∞ ⎩ 0 x<0
1 1
X ∼ Exp(λ ) → E(X) = V(X) =
λ λ2
Disciplina de cola: Es la regla en el orden que se van a seleccionar los clientes que se
encuentran a la espera de ser atendidos en la cola, existen varias reglas, entre las más
comunes se pueden encontrar:
♠ FIFO (first in first out): Se atiende al cliente en el orden que llegan a la cola, el
primero en llegar será el primero en ser atendido. En los modelos básicos de colas se
supone como normal la disciplina de primero en entrar, primero en salir, a menos que
se establezca de otra manera.
♠ LIFO (last in first out): Consiste en atender primero al que ha llegado de último,
también se le conoce como 'pila'.
♠ RSS (random selection of service): Se selecciona a los clientes de una cola de forma
aleatoria, con algún procedimiento de prioridad o algún otra preclasificación.
Capacidad del sistema: Es el número máximo de clientes que pueden estar dentro del
sistema haciendo cola antes de ser atendidos para recibir el servicio, al igual que la
fuente de llegada este número puede ser finito o infinito.
MECANISMO DE SERVICIO
λ
us = ≡ Utilización promedio del sistema.
μ
λ
ρ= (factor de utilización) ≡ Congestión de un sistema.
s. μ
El factor de utilización ρ da una idea de la capacidad del sistema que es utilizada por
los clientes entrantes.
ρ < 1 → Tasa de servicio > Tasa de llegada de clientes
ρ > 1 → Tasa de llegada de clientes > Tasa de servicio → La cola crece con el tiempo
estado: L S = ∑n . p
n=0
n
David Kendall introdujo en 1953 una notación que permite describir las colas y mostrar
sus características pudiendo clasificar los diferentes tipos de colas por medio de
iniciales. De este modo, se tiene:
A / S / c/K /N/D
K Capacidad del sistema, es decir, el número máximo de clientes que puede haber en el
sistema. Cuando se trata de una cola infinita el parámetro se puede omitir.
N Cualquier tipo de disciplina de la cola (FIFO, LIFO, RSS, etc), se puede omitir el parámetro
en caso de ser FIFO.
Un único servidor s = 1
MEDIDAS DE RENDIMIENTO: FÓRMULA DE LITTLE
λ
El Factor de utilización en el caso de un servidor ρ = coincide con la probabilidad de que
μ
λ
un cliente nuevo tenga que esperar para ser servido ρ = p =
μ
Sí ρ < 1 → El sistema es estable. En otro caso, es inestable.
λ
Probabilidad de que ningún cliente se encuentre en el sistema de colas: p 0 = 1 −
μ
λ
Probabilidad de que n clientes se encuentren en el sistema de colas: pn = 1 −
μ
(tiempo promedio en sistema = tiempo promedio en cola + tiempo de servicio)
1
Tiempo de servicio:
μ
1
Tiempo promedio de estancia en el sistema: Ws =
μ −λ
λ ρ
Número promedio de clientes en el sistema: L s = λ .Ws = =
μ−λ 1−ρ
RELACIONES SISTEMA: Wq = ρ . Ws Lq = ρ . Ls
1 λ 1 λ+ μ− λ 1 ⎛ 1 ⎞
Wq + = + = = = Ws ⎜ Ws = Wq + μ ⎟
μ μ . (μ − λ ) μ μ . (μ − λ ) μ − λ ⎝ ⎠
λ λ λ2
Número promedio de clientes en cola: L q = ρ . L s = x = (L q = λ . Wq )
μ μ − λ μ . (μ − λ )
Lq
Número promedio de clientes en la cola para un sistema ocupado: L b =
p
Wq
Tiempo promedio de clientes en la cola para un sistema ocupado: Wb =
p
Número de clientes servidos: L = L s + L q
∞
Número promedio de clientes en ela cola: L q = ∑
n=2
(n − 1) . p n
Customer population (Tamaño de la población de clientes): Aparece por defecto M, indicando que
es infinita. En caso de fuente limitada se pone el tamaño de la población.
Busy server cost per hour: Coste del servidor ocupado ≡ c q + c s
Idle server sost per hour: Coste del servidor desocupado ≡ c s .s
Customer waiting cost per hour: Coste de espera de los clientes ≡ c q
Customer being served cost per hour: Coste de los clientes siendo servidos
Cost of customer being balked: Coste por la pérdida de clientes, en el caso que la cola sea finita
Unit queue capacity cost: Coste unitario de capacidad de cada unidad de cola
Solución:
λ2 0,752
Lq = = = 2,25 clientes/minuto
μ .(μ − λ ) 1. (1 − 0,75)
λ 0,75
Ls = = = 3 clientes/minuto
μ − λ 1 − 0,75
Solución:
⎛ λ ⎞⎛ λ ⎞
Cuando n = 1 → p1 (t) = P(L s = 1) = ⎜ 1 − ⎟ ⎜ ⎟ = (1 − 0,75) x 0,75 = 0,1875
⎝ μ⎠⎝ μ⎠
n
⎛ λ ⎞⎛ λ ⎞
→ p 4 (t) = P(L s = 4) = ( 1 − 0,75 ) ( 0,75 ) = 0,079
4
d) pn (t) = P(L s = n) = ⎜ 1 − ⎟ ⎜ ⎟
⎝ μ⎠⎝ μ⎠
Solución:
Existe un 51,2% de posibilidad de encontrar una cola con más de dos cleintes.
P(Ws > t) = e−μ (1 − ρ ) t → P(Ws > 25) = e−0 ,25 (1 − 0,8) 25 = e− 1,25 = 0,286
P(Wq > t) = ρ e−μ (1 − ρ ) t → P(Wq > 25) = 0,8 e−0,25 (1 − 0 ,8) 25 = 0,8 e− 1,25 = 0,23
Solución:
Lq λ
El tiempo promedio de espera en cola: Wq = =
λ μ . (μ − λ )
λ 0,125
Wq = = = 8,33 minutos
μ . (μ − λ) 0,2 ( 0,2 − 0,125 )
Solución:
En los dos casos se trata de colas M/M/1, puesto que tanto los tiempos entre llegadas
como los tiempos de servicio son variables aleatorias con distribución exponencial.
• En la situación actual, la tasa de llegadas es λ = 1 / 1 = 1 aviones al día, y la tasa de
servicio es μ = 1 / (1 / 2) = 2 aviones al día.
λ 1
Factor de utilización: ρ = = = 0,5
μ 2
λ 1
Número promedio de aviones en el sistema: L s = = = 1 avión
μ−λ 2−1
λ2 1 1
Número promedio de aviones en cola: L q = = = avión
μ . (μ − λ ) 2 . (2 − 1) 2
1 1
Tiempo promedio de estancia en el sistema: Ws = = = 1 día
μ−λ 2−1
λ 1 1
Tiempo promedio de espera en cola: Wq = = = día
μ . (μ − λ ) 2 . (2 − 1) 2
1 1
Tiempo promedio de estancia en el sistema: Ws = = = 4 días
μ − λ 0,5 − 0,25
λ 0,25
Tiempo promedio de espera en cola: Wq = = = 2 días
μ . (μ − λ ) 0,5 . (0,5 − 0,25)
CONFIGURACIONES DE MANTENIMIENTO:
Actual: L s = 1 avión L q = 1 / 2 avión Ws = 1 día Wq = 1 / 2 día
Propuesta: L s = 1 avión L q = 1 / 2 avión Ws = 4 días Wq = 2 días
Con la propuesta de cambio se observa que cada vez que un avión vaya a ser revisado
pasará en el sistema el cuádruple del tiempo que pasaba con el sistema anterior, pero
como cada avión va a ir con frecuencia cuatro veces menor, el tiempo perdido en el
taller a largo plazo va a ser igual.
En este caso, la decisión entre las dos configuraciones se toma en función de los costes
de operación.
Solución:
λ 2
Factor de utilización: ρ = = = 0,6667
μ 3
λ 2
Número promedio de clientes en el sistema: L s = = = 2 clientes
μ−λ 3−2
λ2 4
Número promedio de clientes en cola: L q = = = 1,333 clientes
μ . (μ − λ ) 3 . (3 − 2)
El modelo supone que los tiempos entre llegadas y los tiempos de servicio son
variables aleatorias distribuidas exponencialmente, la disciplina es FIFO y la población
es infinita.
Se diferencia respecto al modelo M/M/1 en que el número de servidores s puede ser
cualquier número natural tal que s ≥ 1 . Cuando el número de servidores es mayor que
1, las expresiones anteriores no son tan sencillas.
En esta línea, en la tasa de servicio μ n hay que distinguir dos casos:
⎧μn = n μ cuando n ≤ s
μn = ⎨ = min(n μ , s μ)
⎩ μ n = s μ cuando n > s
μ ≡ tasa media de servicio de todos los servidores en conjunto
s μ ≡ tasa máxima de servicio para s servidores
El siguiente diagrama de tasas (cadena de Markov del modelo M/M/s) representa los
posibles estados del sistema y las transiciones entre ellos.
⎧ 1 ⎛ λ ⎞n
⎪ ⎜ ⎟ n = 1, 2, ... , s − 1
λ 0 λ 1 ... λ n − 1 ⎪ n! ⎝ μ ⎠
cn = =⎨ n−s
μ1 μ2 ... μn s
⎪1 ⎛λ⎞ ⎛ λ ⎞
⎪ s! ⎜ μ ⎟ ⎜ s μ ⎟ n = s, s + 1, ...
⎩ ⎝ ⎠ ⎝ ⎠
⎧ 1 ⎛ λ ⎞n
⎪ ⎜ ⎟ p0 n = 1, 2, ... , s − 1
⎪ n! ⎝ μ ⎠
pn = c n p 0 = ⎨ s n−s
⎪1 ⎛λ⎞ ⎛ λ ⎞
⎪ s! ⎜ μ ⎟ ⎜ s μ ⎟ p 0 n = s, s + 1, ...
⎩ ⎝ ⎠ ⎝ ⎠
∑p = ∑c p
1 1
n n 0 = 1 → p0 = = s−1 s−1 n−s
=
n ∞
(λ / μ) ⎛ λ ⎞
∑c ∑ ( λ / μ) ( λ / μ)
n n s
n=0 n=0
n=0
n
n=0
n!
+∑n=s
n!
+ ∑
n=0
⎜ ⎟
s! ⎝ s . μ ⎠
1 1
= s−1 s
= s−1
(λ / μ)n 1 ⎛ λ ⎞ 1 (λ / μ)n (λ / μ)s ⎛ s . μ ⎞
∑
n=0
n!
+
s! ⎜⎝ μ ⎟⎠ 1 − ρ ∑
n=0
n!
+
s!
⎜ ⎟
⎝ s .μ − λ ⎠
MEDIDAS DE RENDIMIENTO:
1
Tiempo en el servicio:
μ
λ
Utilización promedio del sistema: us =
μ
λ
Factor de utilización (factor de carga/ intensidad tráfico): ρ =
s. μ
∞
Número promedio de clientes en cola: L q = ∑
n=s+1
(n − s) . pn
s
1 ⎛λ⎞ ρ ( λ / μ) s λ μ
Número promedio de clientes en cola: L q = ⎜ ⎟ p0 = p0
s! ⎝ μ ⎠ (1 − ρ)2 (s − 1)! (s μ − λ)2
L
Tiempo promedio de espera en cola: Wq = q
λ
L q = λ Wq ( )
λ
Número promedio de clientes en sistema: L s = L q + ( L s = λ Ws )
μ
1 ⎛ Ls ⎞
Tiempo promedio de estancia en el sistema: Ws = Wq + ⎜ Ws = ⎟
μ ⎝ λ⎠
L
Número promedio de clientes en la cola para un sistema ocupado: L b = q p w =
pw n≥s
pn ∑
Wq
Tiempo promedio de clientes en la cola para un sistema ocupado: Wb =
pw
pw = ∑p
n≥s
n
A medida que se añaden servidores al sistema las fórmulas van siendo más complicadas,
en especial para el cálculo de probabilidades.
Se asume que la probablidad de la función de tiempos de servicio es una exponencial
negativa de parámetro μ n.
⎛ ⎛λ⎞
s ⎞
⎜ ⎜ μ ⎟ ⎛ 1 − e − μ t (s − 1 − λ / μ ) ⎞ ⎟
−μt ⎜ ⎟
P(W > t) = e ⎜ 1 + ⎝ ⎠ ⎜ ⎟⎟
s! (1 − ρ) ⎝ s − 1 − λ / μ ⎠
⎜ ⎟
⎜ ⎟
⎝ ⎠
− μ t (s − 1 − λ / μ )
1−e
Cuando s − 1 − λ / μ = 0 se utiliza = μ .t
s−1−λ/μ
Customer population (Tamaño de la población de clientes): Aparece por defecto M, indicando que
es infinita. En caso de fuente limitada se pone el tamaño de la población.
Busy server cost per hour: Coste del servidor ocupado ≡ c q + c s
Idle server sost per hour: Coste del servidor desocupado ≡ c s .s
Customer waiting cost per hour: Coste de espera de los clientes ≡ c q
Customer being served cost per hour: Coste de los clientes siendo servidos
Cost of customer being balked: Coste por la pérdida de clientes, en el caso que la cola sea finita
Unit queue capacity cost: Coste unitario de capacidad de cada unidad de cola
Solución:
1
(80 / 50)n (80 / 50)2 ⎛ 100 ⎞
∑
n=0
n!
+
2!
x ⎜ ⎟=
⎝ 100 − 80 ⎠
⎡ (80 / 50)0 (80 / 50) ⎤ (80 / 50)2
=⎢ + ⎥+ x 5 = ( 1 + 1,6 ) + 6,4 = 9
⎣ 0! 1! ⎦ 2!
o bien
1 1
p0 = s−1 s
= 1 n 2
= 0,1111
( λ / μ)
n
1 ⎛λ⎞ 1 (80 / 50) 1 ⎛ 80 ⎞ ⎛ 1 ⎞
∑ n!
+
s! ⎜⎝ μ ⎟⎠ 1 − ρ
∑
n=0
n!
+ ⎜ ⎟ ⎜
2 ⎝ 50 ⎠ ⎝ 1 − 0,8 ⎟⎠
n=0
( λ / μ) s λ μ (80 / 50)2 x 80 x 50
Lq = p0 = 0,111 = 2,8444 clientes
(s − 1)! (s μ − λ )2 1! (2 x 50 − 80)2
L s 4,44
o bien, Ws = = = 0,0556 horas = 3,33 minutos
λ 80
(λ / μ)n
c) El porcentaje de tiempo que determinado operario esté libre: pn = p0
n!
Lq Wq
∑p
2,8444 0,0356
pw = n = 1 − 0,29 = 0,71 ⎯⎯
→ Lb = = =4 Wb = = = 0,0500
n≥2
pw 0,71 pw 0,71
Solución:
El servicio del ambulatorio está utilizado un 87%, esto es, pasa ocioso el 13% del
tiempo, sistema estable al ser ρ = 0,87 < 1
( λ / μ) s λ μ
b) Número promedio de pacientes en la sala de espera: L q = p0
(s − 1)! (s μ − λ )2
1 1 1
p0 = s−1 s
= 2
= = 0,033
(λ / μ)n 1 ⎛ λ ⎞ 1 2,6n 1 ⎛ 1 ⎞ 29,51
∑ n!
+
s! ⎜⎝ μ ⎟⎠ 1 − ρ
∑
n=0
n!
+ 2,63 ⎜
3! ⎝ 0,13 ⎟
⎠
n=0
2 2
2,6n 1 3 ⎛ 1 ⎞ 2,6n
∑
n=0
n!
+ 2,6 ⎜
3! ⎝ 0,13 ⎟=
⎠ n=0
n! ∑
+ 22,53 = 1 + 2,6 + 3,38 + 22,53 = 29,51
( λ / μ) s λ μ 2,63 x 12 x 4,615
con lo que, L q = p0 = x 0,033 = 4,72 pacientes
(s − 1)! (s μ − λ )2 2! x 1,8452
Lq 4,72
Tiempo promedio de espera en cola: Wq = = = 0,39 horas
λ 12
1 1
Tiempo promedio estancia en el sistema: Ws = Wq + = 0,39 + = 0,60 horas
μ 4,615
Solución:
0,5
Tasa de servicio por empleado ≡ μ = 4 x = 1 camiones/hora
2
Utilización promedio del sistema: us = λ / μ = 3 / 1 = 3
λ 3
Utilización promedio de las cuatro plataformas: ρ = = = 0,75
s. μ 4 x 1
3 3
3n 1 4 ⎛ 1 ⎞ 3n
∑
n=0
n!
+ 3
4! ⎜ 0,25 ⎟ =
⎝
∑
⎠ n = 0 n!
+13,5 = 1 + 3 + 4,5 + 4,5 + 13,5 = 26,5
( λ / μ) s λ μ 34 x 3 x 1
Lq = p0 = x 0,03774 = 1,5284
(s − 1)! (s μ − λ )2 3! (4 − 3)2
L 1,5284
Tiempo promedio de espera de camiones cola: Wq = q = = 0,5094 hora
λ 3
Tiempo promedio de estancia de camiones en el sistema:
1 1
Ws = Wq + = 0,5094 + = 1,5094 horas
μ 1
Costes: c q = 60 , c s = 40 , c q + c s = 100
La sala de urgencias de un hospital tiene una tasa media de 3 pacientes a la hora, siguiendo una
distribución de Poisson. La sala cuenta con dos enfermeras que invierten un promedio de 15 minutos
por paciente, según una distribución exponencial. Para evitar la cola de espera surgen dos opiniones:
El jefe de sala solicita una enferma más, la dirección del hospital plantea que en ocasiones las dos
enfermeras están ociosas y considera que se debe reducir a una la cantidad de enfermeras.
Se sabe que una enfermera cobre 10 euros/hora y se ha valorado que por cada hora que un paciente
permanece en la sala el coste es de 5 euros/hora.
¿Cuántas enfermeras debe tener la sala para minimizar el coste total del sistema?
s = 2 servidores
λ = 3 pacientes/hora
μ = 4 . 15 = 60 pacientes/hora
Costes: c q = 5 , c s = 10
9. Número medio de clientes en la cola cuando el sistema esté lleno: L b = 0,6 clientes
12. Tiempo medio de estancia de clientes en la cola cuando el sistema está lleno: Wb = 0,2000 hora
13. Probabilidad de que no haya clientes en el sistema o probabilidad de que todos los servidores
estén ociosos: p0 = 45,4545%
14. Probabilidad de que un cliente llegue al sistema y tenga que esperar, equivalente a la
probabilidad de que esté ocupado el sistema: P(T > 0) = 20,4545%
15. Número medio de clientes que abandonan la cola por hora (para el caso de cola finita), en este
caso como la cola es infinita es cero.
18. Costo total de la espera de los clientes por hora: CQ = c q . L q = 5 . 0,1227 = 0,6136 euros/hora
22. Costo total del sistema por hora: CTS = c q . L s + c s . s = 5 . 0,8727 + 10 . 2 = 24,3635 euros/hora
60
Tasa de servicio por operario: μ = = 40 clientes/hora
1,5
λ 30
Utilización promedio del sistema: us = = = 0,75
μ 40
λ 30
Factor de utilización o congestión del sistema: ρ = = = 0,375
s. μ 2 x 40
( λ / μ) s λ μ 0,752 x 30 x 40
Lq = p 0 = 0,454545 = 0,1227 clientes
(s − 1)! (s. μ − λ)2 (2 x 40 − 30)2
λ
Promedio de clientes en el sistema: L s = L q + = 0,1227 + 0,75 = 0,8727 clientes
μ
Lq 0,1227
b) Tiempo medio de espera en cola: Wq = = = 0,0041 horas
λ 30
1 1
Ws = Wq + = 0,0041 + = 0,0291 horas
μ 40
(30 / 40)
n = 1 : p1 = p 0 = 0,75 x 0,4545 = 0,3409
1
Lq
∑
0,1227
pw = pn = 1 − 0,7955 = 0,2045 ⎯⎯
→ Lb = = = 0,6
n≥2
pw 0,2045
Wq 0,0041
Tiempo promedio de clientes en la cola para un sistema ocupado: Wb = = = 0,0200
pw 0,2045
Solución:
6,51 ⎛ 1⎞
Tiempo medio de respuesta: Ws = = 0,1085 horas ⎜ Ws = Wq + μ ⎟
60 ⎝ ⎠
Número medio de clientes en la cola: L q = 0,17 clientes
Luego,
Lq 0,17
Tiempo medio de espera en cola: Wq = = = 0,0085 horas
λ 20
1 1
b) Ws = Wq + → = Ws − Wq = 0,1085 − 0,0085 = 0,1 → μ = 10 clientes/hora
μ μ
λ 20
Factor de utilización del sistema: ρ = = = 0,5
s μ 4 x 10
λ 20
Con dos comerciales menos ρ = = = 1 el sistema se vuelve inestable.
s μ 2 x 10
λ 20
c) Número medio de comerciales ocupados ≡ = =2
μ 10
Una llamada queda en espera cuando todos los comerciales están ocupados.
P(N ≥ 4) = 1 − P(N < 4) = 1 − p 0 − p1 − p2 − p3
Por ser el sistema estacionario, la tasa media de llegada es igual a la tasa media de
salida para cualquier estado n , es decir, λ n − 1 pn − 1 + μn + 1 pn + 1
2
p0 =
19
20 2 4
10p1 = 20p 0 → p1 = p0 = 2 x =
10 19 19
10p = 20 p 4
20p2 + 20p 0 = 10p1 + 20p1 ⎯⎯⎯⎯⎯
1 0
→ 20p2 = 20p1 → p2 = p1 =
19
20 p = 20 p 20 8
30p3 + 20p1 = 20p2 + 20p2 ⎯⎯⎯⎯⎯
1 2
→ 30p3 = 20p2 → p3 = p2 =
30 57
2 4 4 8 1
P(N ≥ 4) = 1 − P(N < 4) = 1 − p 0 − p1 − p2 − p3 = 1 − − − − =
19 19 19 57 3
1 1 2 4 1 2 2 4 1 3 2 8
p1 = (2) = p2 = (2) = p3 = (2) =
1! 19 19 2! 19 19 3! 19 57
Fuente finita de tamaño k. Los clientes una vez servidos vuelven a la fuente.
En este caso (sistema cerrado) λ es la tasa de retorno de un cliente, no la tasa de
llegadas de los clientes al sistema.
La tasa de retorno es entonces el número de servicios solicitados por unidad de tiempo
y por un cliente.
n
1 k! ⎛ λ ⎞
p0 = pn = p0 si 1 ≤ n ≤ k
k
k! ⎛ λ ⎞
n
(k − n)! ⎜⎝ μ ⎟⎠
∑
n=0
(k − n)! ⎜⎝ μ ⎟⎠
k k n k n
k! ⎛ λ ⎞ k! ⎛ λ ⎞
∑ ∑ ∑
1
pn = ⎜ ⎟ po = po ⎜ μ ⎟ = 1 → p0 =
(k − n)! ⎝ μ ⎠ (k − n)! ⎝ ⎠
k
k! ⎛ λ ⎞
n
n=0 n=0 n=0
∑
n=0
(k − n)! ⎜⎝ μ ⎟⎠
k
Número promedio de clientes en el sistema: L s = ∑n . p
n=1
n
k
Número promedio de clientes en en la cola: L q = ∑ (n − 1) . p
n=2
n
MEDIDAS DE RENDIMIENTO:
Solución:
1 1 1
p0 = k n
= 10 n
= = 0,53796
k! ⎛ λ ⎞ ⎛ 0,005 ⎞ 1,85886
∑ ∑ (10 − n)! ⎜⎝
10!
n=0
(k − n)! ⎜⎝ μ ⎟⎠ n=0
0,1 ⎟⎠
10 n 10
10! ⎛ 0,005 ⎞
∑ ∑
10!
⎜ ⎟ = 0,05n = 1 + 10 x 0,05 + 90 x 0,052 + 720 x 0,053 +
n=0
(10 − n)! ⎝ 0,1 ⎠ n = 0 (10 − n)!
+ 5040 x 0,054 + 30240 x 0, 055 + 151200 x 0, 056 + 604800 x 0,057 + 1814400 x 0,058 +
Tiempo promedio de espera de los robots en la cola para ser atendidos por el
encargado de mantenimiento:
Lq 0,2972
Wq = = = 6,4330 horas
λ (k − L s ) 0,005(10 − 0,76)
1 1
o bien, Ws = Wq + = 6,4330 + = 16,4330 horas
μ 0,1
Número promedio de clientes en la cola Lq
0,2972
para un sistema ocupado: Lb = = = 0,6433
p 0,4620
Tiempo promedio de clientes en la cola W 6,4330
para un sistema ocupado: Wb = q = = 13,9230
p 0,4620
Costes: c q = 10 , c s = 30 , c q + c s = 40
Solución:
1 1 1
p0 = k n
= 6 n
= = 0,4845
k! ⎛ λ ⎞ 6! ⎛ 0,05 ⎞ 2,06392
∑ (k − n)! ⎜⎝ μ ⎟⎠ ∑
n=0
(6 − n)! ⎜⎝ 0,5 ⎟⎠
n=0
6 n 6
6! ⎛ 0,05 ⎞
∑ ∑
6!
⎜ ⎟ = 0,1n = 1 + 6 x 0,1 + 30 x 0,12 + 120 x 0,13 +
n=0
(6 − n)! ⎝ 0,5 ⎠ n = 0 (6 − n)!
p1 p2 p3 p4 p5 p6
0,290708942 0,145354471 0,058141788 0,017442537 0,003488507 0,000348851
1 1
o bien, Ws = Wq + = 1,278 + = 3,279 horas
μ 0,5
En el modelo con población finita los clientes alternan entre estar dentro y fuera del
sistema. Por analogía con el modelo M/M/s se supone que el tiempo que pasa cada
cliente fuera del sistema es una variable aleatoria exponencial Exp(λ ) .
Cuando n clientes están dentro, (k − n) clientes están fuera, y por tanto la distribución
de probabilidad del tiempo que falta para la próxima llegada al sistema es el
mínimo de (k − n) variables exponenciales independientes de parámetro
λ (k − n) . De este modo,
⎧ λ (k − n) 0 ≤ n ≤ N ⎧n μ 1 ≤ n ≤ s − 1
λn = ⎨ μn = ⎨
⎩ 0 n>N ⎩s μ n≥s
Es un sistema cerrado de cola, con fuente finita de tamaño (k). Los clientes una vez
servidos vuelven a la fuente. La tasa de retorno de un cliente es λ
⎧ k! ⎛λ⎞
n
⎪ ⎜ ⎟ p0 n<s
⎪ (k − n)! n! ⎝ μ ⎠
pn = ⎨ n
⎪ k! ⎛λ⎞
⎪ (k − n)! sn − s s! ⎜ μ ⎟ p 0 n≥s
⎩ ⎝ ⎠
k s−1 n k n
⎛λ⎞ ⎛λ⎞
∑ ∑ ∑
k! k!
pn = ⎜ μ ⎟ p0 + n−s ⎜ μ ⎟ p0 = 1
n=0 n=0
(k − n)! n! ⎝ ⎠ n=s
(k − n)! s s! ⎝ ⎠
k ⎡s − 1 ⎛λ⎞
n k
⎛λ⎞ ⎤
n
∑ ∑ ∑
k! k!
pn = p 0 ⎢ ⎜μ⎟ + n−s ⎜μ⎟ ⎥ = 1
⎢ (k − n)! n! ⎝ ⎠ (k − n)! s s! ⎝ ⎠ ⎥
n=0 ⎣n = 0 n=s ⎦
1
p0 = s−1 n k n
⎛λ⎞ ⎛λ⎞
∑ ∑
k! k!
⎜μ⎟ + ⎜ ⎟
n=0
(k − n)! n! ⎝ ⎠ n=s
(k − n)! sn − s s! ⎝ μ ⎠
λ
Ls = Lq + → λ = μ (L s − L q )
μ
k −1
Tasa efectiva de llegada al sistema: λ = λ ∑ (n − s) . p
n=0
n
MEDIDAS DE RENDIMIENTO:
λ
Utilización efectiva del sistema: ρ =
s.μ
Ls 1
Tiempo promedio de clientes en el sistema: Ws = Ws = Wq +
λ μ
Lq
Tiempo promedio de clientes en la cola: Wq =
λ
Lq
Número promedio de clientes en la cola para un sistema ocupado: L b =
pw
pw = ∑p
n≥s
n
Wq
Tiempo promedio de clientes en la cola para un sistema ocupado: Wb =
pw
pw = ∑p
n≥s
n
Solución:
Es un modelo de cola M / M / 2 / /10 de un sistema cerrado, las k = 10 aeronaves constituyen la
población de clientes, verificándose las demás condiciones.
1
Tasa de llegadas: λ = = 0,005 averías/hora
200
1
Tasa de servicio: μ = = 0,1 mantenimiento/hora
10
Probabilidad de que ninguna aeronave se encuentre en el sistema de colas:
1 1
p0 = s−1
= = 0,6099
⎛λ⎞
n k
⎛λ⎞
n
1,5 + 0,1396
∑ (k − n)! n! ⎝⎜ μ ⎠⎟ + ∑ (k − n)! s
k! k!
n−s ⎜ ⎟
n=0 n=s
s! ⎝ μ ⎠
s −1 n 1 n
⎛λ⎞ ⎛ 0,005 ⎞
∑ ∑
k! 10!
⎜ ⎟ = ⎜ 0,1 ⎟ = 1 + 0,5 = 1,5
n= 0
(k − n)! n! ⎝μ⎠ n= 0
(10 − n)! n! ⎝ ⎠
k n 10 n
⎛λ⎞ ⎛ 0,005 ⎞
∑ ∑
k! 10!
⎜ ⎟ = ⎜ 0,1 ⎟ = 45 x 0,05 + 180 x 0,05 + 630 x 0,05
2 3 4
n−s n−2
n=s (k − n)! s s! ⎝ μ ⎠ n=2 (10 − n)! 2 2! ⎝ ⎠
+ 1890 x 0,055 + 4725 x 0,056 + 9450 x 0,057 + 14175 x 0,058 + 14175 x 0,059 + 70875 x 0,0510 =
= 0,1396
k
Número promedio de clientes en la sistema: L s = ∑n . p
n=1
n
k
Ls = ∑n . p
n=1
n = 1 x 0,3050 + 2 x 0,0686 + 3 x 0,0137 + 4 x 0,024 + 5 x 0,0004 = 0,4951
10
Lq = ∑ (n − 2) . p
n=3
n = 1 x 0,0137 + 2 x 0,0024 + 3 x 0,0004 = 0,0198
λ 0,047525
O bien, L q = L s − = 0,4951 − = 0,0198
μ 0,1
k 10
pw = ∑ p = ∑p
n≥ s
n
n≥
n = 1 − 0,9149 = 0,0851
Ls 0,4951
Tiempo promedio de clientes en el sistema: Ws = = = 10,4169
λ 0,0475
Lq 0,0198
Tiempo promedio de clientes en la cola: Wq = = = 0,4169
λ 0,0475
Lq 0,0198
Número promedio de clientes en la cola para un sistema ocupado: L b = = = 0,2327
pw 0,0851
Wq 0,4169
Tiempo promedio de clientes en la cola para un sistema ocupado: Wb = = = 4,8957
pw 0,0851
Coste total en espera: CQ = c q . L q = 10 x 0,0198 = 0,198
Modelo M / M ‐ Modelo M / M / / / k
Modelo M / M / 1 / c ‐ Modelo M / M / s / c
Modelo M / M / c / c
Formalización de sistemas de espera (Colas)
A / B / C / D / E / F
A: Ley de llegadas al sistema. Las llegadas corresponden a un proceso markoviano. Si la tasa de llegadas
es constante a lo largo del tiempo, los tiempos entre llegadas siguen una ley exponencial y el número
de llegadas por unidad de tiempo una ley de Poisson. De ahí que indistintamente a las llegadas
constantes se pueden referir con los términos exponenciales y poissonianos. Cuando las llegadas son
infinitas se denotan por M. El parámetro más utilizado para la tasa de llegadas es λ
B: El parámetro más utilizado para la tasa de servicio es μ (número medio de unidades por unidad de
tiempo que es capaz de entrar en un canal o servidor). Cuando las llegadas al sistema son infinitas se
denotan por M
F: Las llegadas del centro emisor pueden ser infinitas o finitas (k, un número natural). Que el centro
emisor sea infinito significa que no se modifica por el hecho de que algunas unidades se encuentren
en el sistema, por lo que la ley de llegadas es independiente del número de llegadas que contenga el
sistema. Por el contrario, si el centro emisor es finito, con k unidades inicialmente, sus características
dependen del número de elementos en el sistema.
La formalización general de un sistema de espera (Cola) con disciplina FIFO es: M / M / s / c / k
Dejar en blanco un símbolo de la formalización se interpreta que es infinito (M)
En este sentido,
M / M → M / M / ∞ / ∞ / ∞
M / M / / / k → M/M/∞/∞/k
M / M / s / c → M / M / s / c / ∞ (k = ∞)
M / M / s / s Erlang B → M / M / s / c / ∞ (c = s , k = ∞)
1
Tiempo promedio de clientes en el sistema: Ws =
μ
Tiempo promedio de clientes en la cola: Wq = 0
a) Cada vez que un empleado deja la empresa, su puesto entra en la cola de puestos
vacantes. Suponiendo que se inicia de inmediato la búsqueda de un empleado, el modelo
de cola es de una cantidad infinita de servidores. Se trata, pues, de una cola M / M o una
cola (M / M / ∞ / ∞ / ∞)
1 1
Tiempo promedio de empleados en el sistema: Ws = = = 4 meses
μ 0,25
Tiempo promedio de empleados en la cola: Wq = 0
− 13,6
b) Probabilidad de ningún empleado en el sistema: p 0 = e = 0,0000012
n λ
1 ⎛λ⎞ −μ 1
Probabilidad de n empleados en el sistema: pn = ⎜ ⎟ e = 13,6n e − 13,6
n! ⎝ μ ⎠ n!
1 1
p1 = 13,6 e − 13,6 = 0,0000169 p2 = 13,62 e − 13,6 = 0,0001147
1! 2!
c) Para que la empresa tenga siempre 100 empleados trabajando en cualquier momento,
tendrá que tener:
100 + ρ = 100 + 13,6 114 empleados
Tasa de llegadas: λ
Tasa de servicio: μ
λ
Intensidad de tráfico: ρ =
μ
n
k! ⎛λ⎞
(k − n)! n! ⎜⎝ μ ⎟⎠
Probabilidad de n clientes en sistema: pn = k
n = 0 , 1, 2 , ... , k
⎛ λ⎞
⎜1 + μ ⎟
⎝ ⎠
k −1
Tasa efectiva de servicio: λ = λ ∑np
n=1
n
k
⎛ λ ⎞
Número promedio de clientes en el sistema: L s = ∑
n=1
n pn Ls = k ⎜ ⎟
⎝ λ+μ⎠
Número promedio de clientes en la cola: L q = 0
Ls
Tiempo promedio de clientes en el sistema: Ws =
λ
Tiempo promedio de clientes en la cola: Wq = 0
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 101
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 102
Sistemas de colas:
M/M/1/ c / ≡ M/M/1/ c / ∞
M/M/ s / c ≡ M/M/ s / c / ∞
M/M/ c / c ≡ M/M/ c / c / ∞ (s = c , Erlang B)
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 103
MODELO DE COLAS M / M / 1 / c ≡ M / M / 1 / c / ∞
Este tipo de sistemas de colas se caracteriza por tener una cola finita, como indica la cuarta
inicial de la notación de Kendall.
El modelo M / M / 1 / c es aquel en el que un servidor atiende todas las peticiones.
El sistema es aquel en la que los clientes que llegan dejan la cola a partir de una
determinada longitud ya que no están dispuestos a soportar una larga espera.
El número máximo de clientes en el sistema en estos modelos se encuentran limitado a c ,
que coincide con la suma del número de servidores y el tamaño de la cola, por lo que la
capacidad de la cola es (c − 1)
λ
El factor de saturación (intensidad del tráfico) : ρ = determina como varían las
μ
probabilidades pn de que haya n clientes en el sistema.
Sí ρ < 1 → La oferta de demanda supera a la oferta de servicio.
Sí ρ = 1 → Todos los estados son equiprobables.
Sí ρ > 1 → La oferta de servicio supera a la oferta de demanda.
(1 − ρ) . ρn 1−ρ
Probabilidades del estado: pn = n = 0 , 1, 2, ... , c → p =
1 − ρc + 1 1 − ρc + 1
0
pn = ρn . po n = 1, 2, ... , c
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 104
c
Número promedio de clientes en el sistema: L s = ∑n . p
n=1
n
c
Número medio de clientes en cola: L q = ∑ (n − 1) . pn ⎡⎣L s = L q + (1 − p 0 ) ⎤⎦
n=2
c −1
Tasa efectiva de llegada: λ = λ . ∑p
n=0
n
Ls ⎛ 1⎞
Tiempo promedio de estancia en el sistema: Ws = W
⎜ s = W +
μ ⎟⎠
q
λ ⎝
Lq
Tiempo promedio de espera en la cola: Wq =
λ
Lq
Número promedio de clientes en la cola para un sistema ocupado: L b =
1 − po
Wq
Tiempo promedio de clientes en la cola para un sistema ocupado: Wb =
1 − po
MEDIDAS DE RENDIMIENTO:
⎡ (1 − ρ) . ρc ⎤
Tasa efectiva: λ = λ . (1 − pc ) = λ ⎢1 − ⎥
⎣ 1 − ρc + 1 ⎦
⎧ ρ (c + 1). ρc + 1
⎪ − ρ≠1
⎪ (1 − ρ)
Número promedio de clientes en el sistema: L s = ⎨ 1 − ρc + 1
⎪ k
ρ=1
⎪⎩ 2
⎧ (1 − ρc ) . ρ
⎪L s − 1 − ρc + 1 ρ≠1
⎪
Número promedio de clientes en la cola: L q = L s − (1 − p 0 ) = ⎨
⎪ c . (c − 1) ρ=1
⎪⎩ 2 . (c + 1)
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 105
Number of servers (Numero de servidores): s =
Service rate (Tasa de servicio): μ =
Customer arrival rate (Tasa de llegada de clientes): λ =
Queue capacity (Capacidad de la cola: Por defecto aparece M indicando que es infinita. Cuando la
cola es finita se pone el tamaño máximo de la cola menos el número de servidores (c − 1)
Customer population (Tamaño de la población de clientes): Aparece por defecto M, indicando que
es infinita. En caso de fuente limitada se pone el tamaño de la población.
Busy server cost per hour: Coste del servidor ocupado ≡ c q + c s
Idle server sost per hour: Coste del servidor desocupado ≡ c s .s
Customer waiting cost per hour: Coste de espera de los clientes ≡ c q
Customer being served cost per hour: Coste de los clientes siendo servidos
Cost of customer being balked: Coste por la pérdida de clientes, en el caso que la cola sea finita
Unit queue capacity cost: Coste unitario de capacidad de cada unidad de cola
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 106
En un taller mecánico llegan vehículos para una puesta a punto antes de pasar la ITV, las
llegadas siguen un proceso de Poisson de promedio 18 vehículos/hora.
Las dimensiones del taller sólo permiten que haya 4 vehículos, y las ordenanzas municipales
no permiten esperar en la vía pública. El taller despacha un promedio de 6 vehículos/hora
de acuerdo con una distribución exponencial. Se pide:
a) Probabilidad de que no haya ningún vehículo en el taller. Calcular las probabilidades
b) Promedio de vehículos en el taller
c) ¿Cuánto tiempo pasa por término medio un vehículo en el taller?
d) ¿Cuánto tiempo esperan por término medio los vehículos en la cola?
e) ¿Cuál es la longitud media de la cola?
Solución:
a) Es un modelo de cola M / M / 1 / 4 con c = 4 vehículos
Hay un sola cola, con disciplina FIFO, la capacidad del sistema es limitada, de modo que sólo
puede haber 4 vehículos como máximo en el taller, con lo cual el número máximo de
vehículos en la cola es (4 − 1 = 3) .
Las llegadas siguen un proceso de Poisson de parámetro λ = 18 vehículos/hora , los
tiempos entre llegadas se distribuyen exponencialmente Exp(λ = 18) .
Los tiempos entre servicios se distribuyen exponencialmente Exp(μ = 6) siendo
μ = 6 vehículos/hora el número medio que el taller (servidor) es capaz de atender.
λ 18
El factor de saturación ρ = = = 3 determina como varían las probabilidades pn de
μ 6
que haya n vehículos en el sistema.
⎪ p3 = 3 x 0,008264 = 0,2232
3
⎪p = 34 x 0,008264 = 0,6694
⎩ 4
ρ (c + 1) . ρc + 1 3 5 x 35 3 1215
Ls = − = − = − + = 3,5207 vehículos
1−ρ 1 − ρc + 1 1 − 3 1 − 35 2 242
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 107
Ls
c) Tiempo promedio de un vehículo en el taller: Ws =
λ
Tasa de llegada efectiva:
⎡ (1 − ρ) . ρc ⎤ ⎡ (1 − 3) . 34 ⎤
λ = λ ⎢1 − c+1 ⎥
= 18 ⎢1 − ⎥ = 5,9504 vehículos/hora
⎣ 1 − ρ ⎦ ⎣ 1 − 3 5
⎦
Ls 3,5207
Ws = = = 0,5917 horas
λ 5,9504
d) Tiempo medio de espera de vehículos en la cola:
1 1
Wq = Ws − = 0,5917 − = 0,4250 horas
μ 6
(1 − ρc ) . ρ (1 − 34 ) . 3
o bien, L q = L s − = 3,52 − = 2,5289 vehículos
1 − ρc + 1 1 − 35
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 108
4
Ls = ∑= n . p
n 1
n = 1 x 0,0248 + 2 x 0,0744 + 3 x 0,2231 + 4 x 0,6694 = 3,5207
Ls 3,5207 Lq 2,5289
Ws = = = 0,5917 Wq = = = 0,4250
λ 5,9504 λ 5,9504
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 109
Un asesor fiscal dispone de un local para atender a sus clientes, los cuales se
concentran mayoritariamente entre los meses de mayo y junio. El local tiene una
capacidad máxima de 8 asientos en espera, el cliente se va si no encuentra un asiento
libre, y el tiempo entre llegada de clientes se puede considerar distribuido
exponencialmente según un parámetro λ = 20 clientes por hora en período punta.
El tiempo de una consulta esta distribuido exponencialmente con una media de 12
minutos.
¿Cuántas consultas por hora realizará en promedio?
¿Cuál es el tiempo medio de permanencia en el local?
Solución:
Es un modelo M / M / 1 / 9 c = 8 clientes espera + 1 cliente atendido
60
λ = 20 clientes/hora μ= = 5 clientes/hora
12
λ 20
El factor de saturación ρ = = = 4 determina como varían las probabilidades pn de
μ 5
que haya n clientes en el sistema.
(1 − ρ) (1 − 4)
p0 = c+1
= = 0,0000029
1−ρ 1 − 4 10
(1 − ρ) . ρn (1 − 4) . 4 9
Probabilidades del estado: pn = → p9 = = 0,75
1 − ρc + 1 1 − 4 10
(1 − 4)
O bien, pc = ρc . po = 4 9. = 0,75
1 − 4 9+ 1
Tasa de llegada efectiva: λ = λ .(1 − pc ) = 20 . (1 − 0,75) = 5 clientes/hora
ρ (c + 1) . ρc + 1 4 10 x 4 10
Ls = − = − = 8,6667 clientes
1−ρ 1 − ρc + 1 1 − 4 1 − 4 10
L s 8,6667
Tiempo promedio de estancia en el sistema: Ws = = = 1,7333 horas
λ 5
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 110
9
Promedio de clientes en el sistema: L s = ∑n . p
n=0
n = 8,6667
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 111
9
Promedio de clientes en la cola: L q = ∑ (n − 1) . pn = 7,6667
n=2
8
Tasa efectiva de llegada: λ = λ . ∑p
n=0
n = 20 . 0,2499993 = 5
L s 8,6667
Tiempo promedio de estancia en el sistema: Ws = = = 1,7333 horas
λ 5
Lq 7,6667
Tiempo promedio de estancia en la cola: Wq = = = 1,5333 horas
λ 5
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 112
MODELO DE COLAS M / M / s / c ≡ M / M / s / c / ∞
⎪ s! . sn− s . ⎜ μ ⎟ . p 0 n≥s
⎩ ⎝ ⎠
c s n c n
⎛λ⎞ ⎛λ⎞
∑ ∑ ∑
1 1
pn = . ⎜ ⎟ . p0 + .
n− s ⎜ ⎟
. p0 = 1
n=0 n=0
n! μ
⎝ ⎠ n = s+ 1
s! . s μ
⎝ ⎠
c ⎛ s 1 ⎛λ⎞
n c
⎛λ⎞ ⎞
n
∑ ∑ ∑
1
pn = p 0 ⎜ .⎜ ⎟ + .
n− s ⎜ ⎟
⎟=1
⎜ n = 0 n! ⎝μ⎠ s! . s ⎝ ⎠ ⎟⎠
μ
n=0 ⎝ n = s+ 1
1
p0 = s n c n
⎛λ⎞ ⎛λ⎞
∑ ∑
1 1
.⎜ ⎟ + .
n− s ⎜ ⎟
n=0
n! ⎝μ⎠ n = s+ 1
s! . s ⎝μ⎠
c
λ
Número medio de clientes en el sistema: L s = ∑n=0
n . pn = L q +
μ
c
Número medio de clientes en la cola: L q = ∑
n=s+1
(n − s) . pn
c −1
Tasa efectiva de llegada: λ = λ ∑p
n=0
n
c
λ
Ls = ∑n . p
n=0
n = Lq +
μ
→ λ = μ (L s − L q )
λ
Utilización efectiva del sistema: ρ =
sμ
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 113
Ls
Tiempo medio de clientes en el sistema: Ws =
λ
Lq
Tiempo medio de clientes en la cola: Wq =
λ
s
⎛λ⎞ ρ . p0
Lq = ⎜ ⎟ . . ⎡⎣1 − ρk − s − (k − s) . ρk − s . (1 − ρ) ⎤⎦
⎝ μ ⎠ s! . (1 − ρ)
2
c
Lq
Número promedio clientes en cola para un sistema ocupado: Lb =
pw
pw = ∑p
n≥ s
n
c
Wq
Tiempo promedio clientes en cola para un sistema ocupado: Wb =
pw
pw = ∑p
n≥ s
n
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 114
Number of servers (Numero de servidores): s =
Service rate (Tasa de servicio): μ =
Customer arrival rate (Tasa de llegada de clientes): λ =
Queue capacity (Capacidad de la cola: Por defecto aparece M indicando que es infinita. Cuando la
cola es finita se pone el tamaño máximo de la cola menos el número de servidores (c − s)
Customer population (Tamaño de la población de clientes): Aparece por defecto M, indicando que
es infinita. En caso de fuente limitada se pone el tamaño de la población.
Busy server cost per hour: Coste del servidor ocupado ≡ c q + c s
Idle server sost per hour: Coste del servidor desocupado ≡ c s .s
Customer waiting cost per hour: Coste de espera de los clientes ≡ c q
Customer being served cost per hour: Coste de los clientes siendo servidos
Cost of customer being balked: Coste por la pérdida de clientes, en el caso que la cola sea finita
Unit queue capacity cost: Coste unitario de capacidad de cada unidad de cola
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 115
COLA FINITA: INVESTIGACIÓN
Un grupo de investigadores, formado por seis personas, dispone de dos terminales para realizar
cálculos. El trabajo promedio de cálculo requiere de 20 minutos de tiempo de terminal, y el tiempo
promedio entre solicitudes de servicio es de 30 minutos. Se supone que estas solicitudes están
distribuidas exponencialmente. Se desea saber:
a) Número estimado de investigadores que esperan utilizar una terminal.
b) Tiempo total perdido diariamente si se considera una jornada de 8 horas.
c) Medidas de rendimiento.
Solución:
60 60
Tasa de llegada: λ = = 2 clientes/hora Tasa de servicio: μ = = 3 clientes/hora
30 20
1 1
p0 = s n c n
= = 0,5003
1 ⎛λ⎞ ⎛λ⎞ 1,9988
∑ ∑
1
.⎜ ⎟ + .
n−s ⎜ ⎟
n=0
n! ⎝ μ ⎠ n = s + 1 s! . s ⎝μ⎠
2 n 6 n
1! ⎛ 2 ⎞ ⎛2⎞
∑ ∑
1
.⎜ ⎟ + n−2 ⎜ ⎟
. = 1,8922 + 0,1076 = 1,9998
n=0
n! ⎝ 3 ⎠ n = 3 2! . 2 ⎝3⎠
2 n 2
1! ⎛ 2 ⎞ 2 1 ⎛2⎞
∑
n=0
. ⎜ ⎟ = 1 + + ⎜ ⎟ = 1,8922
n! ⎝ 3 ⎠ 3 2 ⎝3⎠
6 n
⎛2⎞
∑
1
n−2 ⎜ ⎟
. = 0,0741 + 0,0247 + 0,0082 + 0,0003 = 0,1076
n = 3 2! . 2 ⎝ 3 ⎠
⎧ 1 ⎛ λ ⎞n
⎪ . ⎜ ⎟ . p0 n≤2
⎪ n! ⎝ μ ⎠
Probabilidades de cada estado del sistema: pn = ⎨ n
⎪ 1 ⎛λ⎞
⎪ .
n− s ⎜ ⎟
. p0 n≥2
⎩ s! . s ⎝μ⎠
2
1 ⎛2⎞ 1 ⎛λ⎞
p1 = . ⎜ ⎟ . 0,5003 = 0,3336 p2 = . ⎟ . 0,5003 = 0,1112
1! ⎝ 3 ⎠ 2 −2 ⎜
2! . 2 ⎝μ⎠
3 4
1 ⎛2⎞ 1 ⎛2⎞
p3 = . ⎜ ⎟ . 0,5003 = 0,0371 p4 = 2 ⎜ ⎟
. . 0,5003 = 0,0124
2! 2 ⎝ 3 ⎠ 2! . 2 ⎝ 3 ⎠
5 6
1 ⎛2⎞ 1 ⎛2⎞
p5 = 3 ⎜ ⎟
. . 0,5003 = 0,0041 p6 = 4 ⎜ ⎟
. . 0,5003 = 0,0014
2! 2 ⎝ 3 ⎠ 2! 2 ⎝ 3 ⎠
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 116
Tabla estados de probabilidad
n pn n pn (n − 2) pn
0 0,5003
1 0,3336 0,3336
2 0,1112 0,2224
3 0,0371 0,1113 0,0371
4 0,0124 0,0496 0,0248
5 0,0041 0,0205 0,0123
6 0,0014 0,0008 0,0005
0,7454 0,0796
6
Número medio de clientes en el sistema: L s = ∑= n . p
n 0
n = 0,7454
6
Número medio de clientes en cola: L q = ∑ (n − s) . p
n=3
n = 0,0796
5
Tasa efectiva de llegada: λ = 2 ∑p
n=0
n = 2 (1 − 0,0014) = 1,9973
Lq 0,0796
Número promedio clientes en cola para un sistema ocupado: Lb = = = 0,4793 horas
pw 0,1661
Wq 0,0399
Tiempo promedio clientes en cola para un sistema ocupado: Wb = = = 0,2400 horas
pw 0,1661
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 117
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 118
Sstema de pérdidas Erlang B
Modelos de cola sin cola
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 119
MODELO DE COLA M / M / c / c (s = c , Erlang B)
Es un sistema de cola sin cola, estudiado por Agner Krarup Erlang (inventor de la teoría de
colas) y se conoce como sistema de pérdidas Erlang B (Tablas de Erlang B).
Es un caso típico de una red telefónica con "c " líneas, de manera que cuando las líneas se
encuentran ocupadas, el cliente obtiene una señal de ocupado y cuelga.
Al no tener cola, los valores L q = 0 y Wq = 0 , mientras que el tiempo del sistema es igual
1
al tiempo de servicio Ws =
μ
1
Probabilidad de no haber clientes en el sistema: p 0 = c n
1 ⎛λ⎞
n=0
∑
n! ⎜⎝ μ ⎟⎠
n
1 ⎛λ⎞
Probabilidad de n clientes en el sistema: pn = p0
n! ⎜⎝ μ ⎟⎠
λ
Intensidad del tráfico: ρ = Erlangs ‐ Tabla B(c , ρ)
μ
c−1 ⎡ 1 ⎛λ⎞
c ⎤
Tasa efectiva de llegadas: λ = λ ∑p
n=0
n λ = λ ⎢1 − ⎜ ⎟ p 0 ⎥
⎢⎣ c! ⎝ μ ⎠ ⎥⎦
λ
Utilización efectiva del sistema: ρ =
cμ
c
Número promedio de clientes en el sistema: L s = ∑n p
n=1
n λ = μ Ls
Ls
Tiempo promedio de clientes en el sistema: Ws =
λ
Distribución del tiempo en servicio: W(t) = 1 − e − μ t
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 120
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 121
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 122
Dos equipos médicos trabajan en dos salas de investigación, admitiendo ensayos clínicos cuando se
encuentran sin trabajo teniendo las peticiones que abandonar. Las solicitudes de ensayos, indistintamente
de la sala, siguen un proceso de Poisson con tasa de llegada de 8 ensayos/día. Cada sala da salida a los
pedidos con tasa de 4 ensayos/día de acuerdo a una distribución exponencial. Calcular las medidas de
rendimiento.
Solución:
Es una cola M / M / 2 / 2 con lo que el número de ensayos clínicos en cola es c − s = 2 − 2 = 0
λ
Tasa de llegada: λ = 8 ensayos/día Tasa de servicio: μ = 4 ensayos/día =2
μ
1 1
Probabilidad de no haber clientes en el sistema: p0 = = = 0,2
2
1+2+2
∑ n! 2
1 n
n=0
n
1 ⎛λ⎞ 1 1 2
Probabilidad de estado: pn = ⎜ ⎟ p0 p1 = 2 . 0,2 = 0,4 p2 = 2 . 0,2 = 0,4
n! ⎝ μ ⎠ 1! 2!
c −1 1
Tasa efectiva de llegadas: λ= λ ∑p
n=0
n =8 ∑p
n=0
n = 8 . (0,2 + 0,4) = 4,8 ensayos/día
λ 4,8
Utilización efectiva sistema: ρ = = = 0,6
c μ 2.4
2
Número promedio de clientes en el sistema: L s = ∑n p
n=1
n = p1 + 2p2 = 0,4 + 2 . 0,4 = 1,2
L s 1,2
Tiempo promedio de clientes en el sistema: Ws = = = 0,25 días λ = μ L s = 4 . 1,2 = 4,8
λ 4,8
1 1
Denotar que Ws = = = 0,25 días
μ 4
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 123
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 124
En un taller de mantenimiento de aeronaves las llegadas siguen un proceso de Poisson de promedio 3
aeronaves/día. Las dimensiones del recinto permiten albergar a una sola aeronave. El taller despacha un
promedio de 2 aeronaves/día de acuerdo con una distribución exponencial.
Calcular las medidas de rendimiento.
Solución:
Es una cola M / M / 1 / 1 con disciplina FIFO, con una capacidad limitada c = 1 aeronave, con lo que el
número de aeronaves en la cola es (c − 1 = 1 − 1 = 0) . Es una cola sin cola.
Solución:
Tasas de llegada y servicio:
λ 3
λ = 3 aeronaves/día μ = 2 aeronaves/día Factor saturación: ρ =
= = 1,5 Erlangs
μ 2
1 1 1
Probabilidad de que no haya aeronaves en taller: p0 = 1 = 1 = = 0,4
1 ⎛λ⎞
n
1 + 1,5
∑ ∑
1
⎜ ⎟ 1,5n
n=0
n! ⎝ μ ⎠ n=0
n!
L s 0,6
Tiempo promedio de clientes en el sistema: Ws = = = 0,5 días λ = μ L s = 2 . 0,6 = 1,2
λ 1,2
1 1
Denotar que Ws = = = 0,5 días
μ 2
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 125
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 126
Una empresa instala un sistema de comunicación entre sus dos sedes. Las llamadas reciben una señal
de ocupado cuando las líneas están ocupadas. El sistema genera llamadas aleatoriamente, según un
proceso de Poisson con una tasa de 105 llamadas/hora y las llamadas tardan cuatro minutos por término
medio para ser contestadas.
a) La empresa pretende instalar las líneas necesarias para asegurar que la probabilidad de recibir la señal
de ocupado sea inferior a 0,005. ¿Cuántas líneas hacen falta?
b) ¿Cuál será la probabilidad de recibir la señal de congestionado con sólo 10 líneas?
Solución:
a) Tasas de llegada y servicio:
105
λ = 105 llamadas/hora = = 1,75 llamadas/minuto
60
1
μ= llamadas/minuto
4
λ 1,75
Intensidad del tráfico: = = 7 Erlangs
μ 0,25
A partir de la tabla Erlang B (B = 0,5%) hay que encontrar el valor mínimo de "c" que verifique la
desigualdad B(c , 7) ≤ 0,005
Erlang B 0,5%
Servidores
Erlang B 5% x% 10%
7,511 − 6,216 7 − 6,216
Servidores =
0,1 − 0,05 x − 0,05
10 6,216 7 7,511 x = 0,079
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 127
Portal Estadística Aplicada: MODELOS M/M ‐ M/M///k ‐ M/M/1/k ‐ M/M/s /k ‐ M/M/c/c 128
TEORÍA DE COLAS: MODELOS NO EXPONENCIALES
Los modelos de teoría de colas vistos anteriormente se basan en el proceso de nacimiento y
muerte, lo que hace necesario que tanto los tiempos entre llegadas como los de servicio tengan
distribuciones exponenciales. Este tipo de distribuciones de probabilidad proporciona un ajuse
idóneo cuando los tiempos entre llegadas son exponenciales (que implica que las llegadas ocurren
al azar, proceso de entrada de Poisson).
No ocurre lo mismo cuando las llegadas están reguladas o programadas. Por otra parte, las
distribuciones de tiempos de servicio reales con frecuencia se desvían de la distribución
exponencial, en particular cuando los requerimientos de servicio de los clientes son muy parecidos.
Por ello, se hace necesario disponer de modelos de colas que utilizen otras distribuciones de
probabilidad.
λ
Factor de utilización: ρ =
s.μ
⎧ 1 ⎛ λ ⎞n
⎪ ⎜ ⎟ p0 n≤ s
⎪ n! ⎝ μ ⎠
Probabilidad del estado M/M/ s : pn = ⎨ n
⎪ 1 ⎛λ⎞
⎪ s! sn − s ⎜ μ ⎟ p 0 n ≥ s
⎩ ⎝ ⎠
⎛λ⎞ ⎛ 0,04 ⎞
Número promedio de clientes en sistema: L s = L q + ⎜ ⎟ = 0,0019 + ⎜ ⎟ = 0,8019 pruebas
⎝μ⎠ ⎝ 0,05 ⎠
1 1
Tiempo promedio de estancia en el sistema: Ws = Wq + = 0,0465 + = 20,0465 segundos
μ 0,05
1
p1 = 0,8 x 0,4286 = 0,3429
1!
1
p2 = 0,8 2 x 0,4286 = 0,1371
2!
1
p3 = 3−2
x 0,8 x 0,4286 = 0,0549
3
2! 2
1
p8 = 8−2
x 0,8 x 0,4286 = 0,0006
8
2! 2
λ 2 . σ 2 + ρ2
Número promedio de clientes en la cola: L q =
2. (1 − ρ)
λ 2 . σ 2 + ρ2
Lq
Tiempo promedio de espera en la cola: Wq = =
λ 2 . λ . (1 − ρ)
Ls 1
Tiempo promedio de estancia en el sistema: Ws = = Wq +
λ μ
Resulta más beneficioso contratar al primer auxiliar de tierra (Ws1 = 40,20 < Ws2 = 87,50) al ser más
rápido que el segundo.
p 0 = 1 − ρ1 = 1 − 0,6667 = 0,3333
pn = ρn p 0
i
i
i
λ 2 . σ 2 + ρ2 σ2 = 0 ρ2
Número promedio clientes en la cola: L q = ⎯⎯⎯→ L q =
2. (1 − ρ) 2 . (1 − ρ)
ρ2Lq
Tiempo promedio de espera en la cola: Wq = =
λ 2 . λ . (1 − ρ)
Ls 1
Tiempo promedio de estancia en el sistema: Ws = = Wq +
λ μ
λ
Factor de utilización: ρ =
s.μ
⎧ 1 ⎛ λ ⎞n
⎪ ⎜ μ ⎟ p0 n≤ s
⎪ n! ⎝ ⎠
Probabilidad del estado M/M/ s : pn = ⎨ n
⎪ 1 ⎛λ⎞
⎪ s! sn − s ⎜ μ ⎟ p 0 n ≥ s
⎩ ⎝ ⎠
s
1 ⎛λ⎞ ρ
Número promedio clientes en la cola: L q = ⎜ ⎟ p0
s! ⎝ μ ⎠ 2 (1 − ρ)2
⎛λ⎞
Número promedio de clientes en sistema: L s = L q + ⎜ ⎟
⎝μ⎠
Lq
Tiempo promedio de espera en la cola: Wq =
λ
Ls 1
Tiempo promedio de estancia en el sistema: Ws = = Wq +
λ μ
El aeropuerto dispone de un servicio de traslado en el que consiste en llevar a cada empleado que lo
solicite a su casa, hotel o alrededores.
Este servicio puede atender a un empleado cada 7 minutos. El promedio de llegada de empleados es
cada 8 minutos, siguiendo una distribución de Poisson.
a) Encontrar las medidas de eficiencia del servicio.
b) ¿Se podría mejorar el tiempo medio de un empleado en el sistema?
1
Tasa de llegada = = 8 → λ = 0,125 pasajeros / minuto
λ
1
Tasa de servicio = = 7 → μ = 0,1429 pasajeros / minuto
μ
λ 0,1250
Factor de utilización: ρ = = = 0,8750
μ 0,1429
Probabilidad del estado: p0 = 1 − ρ = 1 − 0,8750 = 0,125
ρ2 0,8752
Promedio de empleados en la cola: L q = = = 3,062 empleados
2 . (1 − ρ) 2 x (1 − 0,875 )
1 1
p0 = s−1 n s
= 1 n 2
= 0,3913
1 ⎛λ⎞ 1 ⎛λ⎞ ⎛ 1 ⎞ 1 ⎛ 0,125 ⎞ 1 ⎛ 0,125 ⎞ ⎛ ⎞
∑ n! ⎝⎜ μ ⎠⎟ ∑
1
+ ⎜ ⎟ + ⎜ ⎟ ⎜ ⎟
n=0
s! ⎝⎜ μ ⎠⎟ ⎝⎜ 1 − ρ ⎠⎟ n=0
0! ⎝ 0,1429 ⎠ 2! ⎝ 0,1429 ⎠ ⎝ 1 − 0,4375 ⎠
1
1 ⎛λ⎞ 1
p1 = ⎜ ⎟ p0 = 0,875 x 0,3913 = 0,3424
n! ⎝ μ ⎠ 1!
2
1 ⎛λ⎞ 1
p2 = ⎜ ⎟ p0 = 0,8752 x 0,3913 = 0,1498
n! ⎝ μ ⎠ 2!
n
1 ⎛λ⎞ 1
p3 = n−s ⎜ ⎟ p 0 = 3−2
0,8753 x 0,3913 = 0,0655
s! s ⎝μ⎠ 2! 2
⎛λ⎞ ⎛ 0,125 ⎞
Número promedio de clientes en sistema: L s = L q + ⎜ ⎟ = 0,1036 + ⎜ ⎟ = 0,9786 pasajeros
⎝μ⎠ ⎝ 0,1429 ⎠
Lq 0,1036
Tiempo promedio de espera en cola: Wq = = = 0,8285 minutos
λ 0,125
L s 0,9786
Tiempo promedio de espera en servicio: Ws = = = 7,8285 minutos
λ 0,125
Erlang C
La fórmula supone una infinita población de fuentes, que ofrecen en conjunto, un tráfico
de A Erlangs hacia n servidores.
Si todos los servidores están ocupados cuando una petición llega de una fuente, la
petición es introducida en la cola. Un sinfín de números de peticiones podrían ir a la cola
en este modo simultáneamente.
A ≡ Intensidad total del tráfico ofrecido en unidades de
n An
Erlangs.
(n − A) n!
pw = n − 1 n ≡ Número de servidores [número de troncales].
Ai n An
∑ i! + (n − A) n! Pw ≡ Probabilidad de que un cliente tenga que esperar para
i=0
ser atendido.
Se supone que las llamadas entrantes pueden ser modeladas usando una distribución de
Poisson y que el tiempo de espera de las llamadas son descritas por una distribución
exponencial negativa.
Calcula la probabilidad del tráfico de la cola, suponiendo que las llamadas que son
bloqueadas se queden en el sistema hasta que puedan ser atendidas.
Fórmula de Engset
k! An E ≡ Probabilidad de bloqueo.
(k − n)! n! A ≡ Tráfico de Erlangs generado por cada origen desocupado.
E(n , A, k) = n
k! A i n ≡ Número de servidores.
∑
i = 0 (k − i)! i! k ≡ Número de orígenes.
⎧ E(0 , A, k) = 1
⎪
Se puede calcular por recursión: ⎨ A (k − n + 1) E(n − 1 , A, k)
⎪E(n , A, k) =
⎩ n + A (k − n + 1) E(n − 1 , A, k)
Se supone que las llamadas que llegan pueden ser modeladas por una distribución
Poisson. Sin embargo, debido a que hay un número finito de servidores, la tasa de llegada
de las nuevas llamadas decrece a medida que nuevos orígenes (como abonados
telefónicos) pasan a estar ocupados y, en consecuencia, no pueden originar nuevas
llamadas. Cuando n = k , la fórmula se reduce a una distribución binomial.
Se basa en el método de fases, asume una llegada Poisson sin memoria. El servicio
(servidor) es la concatenación de servicios en serie. Sólo puede haber una llamada en el
sistema completo, por lo que es imposible que haya más de una 'fase' ocupada.
Cada una de las fases tiene un tiempo de servicio distribuido según una exponencial
1
negativa de tasa υ = k. μ , de media , donde:
μ
k 1 1 ⎛1⎞ 1
σ2 = = → σ= .⎜ ⎟ k=
υ k. μ2 ⎝μ⎠ μ x σ2
2 2
k
El tiempo de servicio total es la suma de los tiempos de servicio de cada fase. La suma de k
1 1
variables aleatorias exponenciales negativas con media = da lugar a una variable
υ k. μ
aleatoria con distribución Erlang k.
i Con un único servidor, el modelo M / Ek / 1 de parámetros k y μ , se reduce a un caso
particular del modelo M / G / 1 :
1
σ2 =
k μ2
M / G / 1 ⎯⎯⎯⎯→ M / Ek / 1
1 1
i Cuando σ = → k = 2 2 = 1 , el modelo de cola M / Ek / 1 es un modelo de
μ μ xσ
cola M / M / 1 con tiempo de servicio exponencial.
ρ2 x (1 + 1) ρ2
Promedio de clientes en la cola: Lq = =
2 x 1 x (1 − ρ) (1 − ρ)
MEDIDAS DE RENDIMIENTO:
λ
Factor utilización del sistema: ρ =
μ
λ 2 . (1 + k) (1 + k) . ρ2
Número promedio de clientes en la cola: L q = =
2 . k . μ .(μ − λ ) 2 . k . (1 − ρ)
Lq
Tiempo promedio de espera en la cola: Wq =
λ
Ls 1
Tiempo promedio de estancia en el sistema: Ws = = Wq +
λ μ
Solución:
Rendimiento cola (M / G / 1)
• Con la modificación = Rendimiento cola (M / M / 1)
⎛ 1 + μ2 . σ2 ⎞
⎜ ⎟
⎝ 2 ⎠
se pasa de las medidas de rendimiento de una cola M / G / 1 a medidas de rendimiento de una cola
M/M/1
1 + μ2 . σ2 1 + 0,22. 22
= = 0,58
2 2
Ls 1,3050 Ls 3
+ρ= + 0,75 = 3 = Ls Ws = = = 20 = Ws
0,58 0,58 M/G/1
M/M/1 λ M/G/1 0,15 M/M/1
σ2 = 0
• Con un único servidor, M / G / 1 ⎯⎯⎯→ M / D / 1
λ 2 . σ 2 + ρ2 0,752
Promedio de vehículos en cola: L q = = = 1,125 vehículos
2 . (1 − ρ) 2 x (1 − 0,75)
Promedio de vehículos en sistema: L s = L q + ρ = 1,125 + 0,75 = 1,875 minutos
1,125 Lq
Tiempo promedio de espera en la cola: Wq = =
= 7,5 minutos
λ 0,15
L 1,875
Tiempo promedio de estancia en lavacoches: Ws = s = = 12,5 minutos
λ 0,15
ρ2 . (1 + k) 0,752 x (1 + 6,25)
Lq = = = 1,305 clientes
2 . k . (1 − ρ) 2 x 6,25 x (1 − 0,75)
Promedio de clientes en el sistema: L s = L q + ρ = 1,305 + 0,75 = 2,055 clientes
Lq 1,305
Tiempo promedio en la cola: Wq = = = 8,7 minutos
λ 0,15
Ls 1
Tiempo promedio en el sistema: Ws = = Wq + = 8,7 + 5 = 13,7 minutos
λ μ
1 1
Se observa que cuando σ = → k= 2 = 1 , el modelo de cola M / Ek / 1 es un modelo de
μ μ x σ2
cola M / M / 1 con tiempo de servicio exponencial. En este caso,
ρ2 x (1 + 1) ρ2
Promedio de clientes en la cola: L q = =
2 x 1 x (1 − ρ) (1 − ρ)
9
λ= = 0,15 vehículos/minuto
60
1
Tasa de servicio =
= 5
μ
→ μ = 0,2 pasajeros/minuto
λ 0,15
ρ= = = 0,75
μ 0,2
Modelo Determinista
Sistema de un Aeropuerto: Equilibrio entre
Demanda y Capacidad Operacional
MODELO DETERMINISTA DE COLAS
Cuando la demanda λ es mayor que la capacidad μ , el factor de utilización
λ
ρ = > 1 , entonces para calcular las demoras es necesario utilizar otro
μ
criterio en teoría de colas.
Este fenómeno ocurre en muchos sistemas durante un período corto de tiempo
introduciendo el modelo determinista de colas.
Un modelo determinístico es un modelo matemático donde las mismas entradas
producen invariablemente las mismas salidas, sin contemplar la existencia del azar ni el
principio de incertidumbre.
Estado del sistema ≡ número de clientes en el sistema
L q ≡ Longitud de la cola ≡ número de clientes en cola =
= estado del sistema − número de clientes en servicio
λ ≡ tasa de llegadas
μ ≡ tasa de servicios
λ>μ
⎛ 1⎞
Sí t ∈ ⎜ 0, ⎟ → No ha llegado todavía ningún cliente, en este intervalo de tiempo no
⎝ λ⎠
hay clientes en el sistema.
⎛1 ⎞ ⎛ ⎛ 1 ⎤⎞
Si t ∈ ⎜ , τ ⎟ → N(t) = ⎡⎣llegadas en (0, τ ]⎦⎤ − ⎜ servicios en ⎜ , τ ⎥ ⎟
⎝λ ⎠ ⎝ ⎝ λ ⎦⎠
⎛ 1 ⎤ ⎛ t−1/ λ ⎞ μ
Número de salidas ⎜ , τ ⎥ ≡ ⎜ ⎟ = μt −
⎝ λ ⎦ ⎝ 1/μ ⎠ λ
⎛ μ⎞
En consecuencia, N(t) = λ t − ⎜ μ t − ⎟
⎝ λ⎠
En este caso, siempre que se produce la salida de un cliente hay una llegada
simultánea. Así, no puede suceder que haya simultáneamente una salida y un rechazo,
pues si un cliente sale, se deja una plaza libre en el sistema para el otro cliente que
llega en ese momento.
En consecuencia, el número de clientes es creciente hasta que a partir del
instante τ es constantemente igual a la capacidad (k − 1) del sistema.
Luego, se tiene:
⎧ 1
⎪ 0 t<
λ
⎪
⎪ ⎛ μ⎞ 1
N(t) = ⎨ λ t − ⎜ μ t − ⎟ <t<τ
⎪ ⎝ λ⎠ λ
⎪ k −1 t≥τ
⎪
⎩
Suponiendo que se tiene disciplina FIFO, los valores asociados a cada cliente:
1
Aunque la relación se verifica en general, se observa que S(n) = y si t < τ entonces
μ
Por otra parte, si t ≥ τ un cliente entra al sistema justo cuando sale otro. Por tanto, el
cliente que entra se encuentra con (k − 2) clientes en cola delante de él y un cliente
empezando el servicio. En consecuencia, el tiempo de espera en cola para este nuevo
1
cliente será (k − 2)
μ
⎧ ⎛1 1⎞
⎪(n − 1) ⎜ − ⎟ n = 1, 2, ... , λτ − 1
⎪
En definitiva: Wq(n) = ⎨ ⎝μ λ⎠
⎪ k−2 n ≥ λτ
⎪⎩ μ
1 1
b) ≥
λ μ
El tiempo entre llegadas es mayor o igual que el tiempo de servicio, cuando llega un
cliente es servido y sale del sistema antes de que llegue el siguiente cliente (en la peor
situación, cuando llegue el siguiente cliente si el tiempo entre llegadas es igual al
tiempo de servicio).
Por tanto, considerando que en el instante cero no hay individuos en el sistema, el
número de clientes en el sistema siempre variara entre [0, 1] y el tiempo de espera en
cola de cualquier cliente es cero.
En este caso, una situación particular supone que cuando el sistema se inicia ya hay un
número de clientes C esperando en el mismo.
1 1
• =
λ μ
Si el tiempo entre llegadas coincide con el tiempo de servicio, siempre habrá C clientes
en el sistema y el tiempo de espera en cola del n − ésimo cliente toma la expresión:
⎧ 1
⎪⎪ (n − 1) 1≤n≤ C
μ
Wq = ⎨
(n)
⎪(C − 1) 1 n≥C
⎪⎩ μ
1 1
• >
λ μ
Si el tiempo entre llegadas es mayor que el tiempo de servicio, hay un primer instante
t = 0 no hay clientes N(0) = 0 .
⎡⎛ t ⎞ ⎛ t ⎞ ⎤
N(t) = C − ⎢⎜ ⎟−⎜ ⎟ ⎥ = C − ( μ t − λt )
⎣⎝ 1 / μ ⎠ ⎝ 1 / λ ⎠ ⎦
1 ⎛1⎞ 2
El siguiente cliente llega en el instante τ 2 = τ1 + = λ t ⎜ ⎟ + , después de haber
λ ⎝λ⎠ λ
salido el cliente anterior, después el sistema vuelve a estar vacío en el intervalo
⎡ 1 ⎞ ⎡ 1⎞
⎢ τ1 + μ , τ 2 ⎟ y toma el valor 1 entre ⎢ τ 2 , τ 2 + μ ⎟ .
⎣ ⎠ ⎣ ⎠
⎧
⎪ C − (μ t − λt) 0 ≥ t < τ
⎪
⎪ 1 1 1 k
N(t) = ⎨ 0 τ ≤ t < τ1 τk + ≥ t < τ k + siendo τk = (λτ) +
⎪ μ λ λ λ
⎪ 1
⎪ 1 τk ≥ t < τk +
⎩ μ
Verificando la igualdad Wq(n+1) = Wq(n) + S(n) − T (n) se tiene la siguiente expresión para el
tiempo de espera en cola:
⎧ 1
⎪ (n − 1) 1≤n≤ C
μ
⎪
⎪ 1 k
Wq(n) = ⎨(C + k − 1) − n = C + k 0 ≤ k ≤ λτ
⎪ μ λ
⎪ 0 n = C + k k > λτ
⎪
⎩
λ 45
Utilización promedio del sistema: us = = = 1,125
μ 40
λ 45
Factor de utilización o intensidad tráfico: ρ = = = 0,5625
s μ 2 x 40
λ
Número medio de aviones en sistema: L s = L q + = 0,52 + 1,125 = 1,645
μ
Con una hoja de cálculo se puede elaborar un analisis de sensibilidad del sistema y un
diagrama de la variación de aviones en cola según demanda. En esta línea, para
distintas tasas de llegadas de aviones se obtiene:
λ 20 25 30 35 40 45 50 55
Wq 0,1 0,162 0,245 0,355 0,5 0,694 0, 962 1,344
••••••••
1 1 60 x 60
= = = 80 ⇒ μ = 80 pasajeros/hora seguridad
μ 45 / 60 x 60 45
λ 144
Factor de utilización del sistema de cola: ρ = = = 0,9
s μ 2 x 80
λ
Número medio de pasajeros en sistema: L s = L q += 7,674 + 1,8 = 9, 474
μ
L 7,674
Tiempo promedio espera en la cola: Wq = q = (60 x 60) = 192 segundos
λ 144
L s 9, 474
Tiempo promedio espera en sistema: Ws = = (60 x 60) = 237 segundos
λ 144
En la situación actual con dos aparatos de rayos‐X (servidores) el tiempo medio de
pasajeros en cola es de 192 segundos > 120 segundos, por lo que se necesitaría
incorporar servidores (aparatos rayos‐X)
L q 22,2157
Tiempo promedio espera en la cola: Wq = = (60 x 60) = 347 segundos
λ 230, 4
Siendo Wq = 347 segundos > 120 segundos se necesita incrementar el número de
servidores en el sistema.
L q 1,1846
Tiempo promedio espera en la cola: Wq = = (60 x 60) = 18,51 segundos
λ 230, 4
Wq = 18,51 segundos < 120 segundos , con lo que habría que instalar 4 servidores (4
aparatos de rayos‐X), con número medio de cuatro pasajeros en el sistema:
λ
L s = L q + = 1,1846 + 2,88 = 4,06 ≈ 4
μ
(λ / μ)n (λ / μ)n
Probabilidad del estado n: pn = p0 n ≤ s pn = p0 n > s
n! s! s(n − s)
⎛ 4 (λ / μ)n 8
(λ / μ)n ⎞
P(n > 8) = 1 − P(n ≤ 8) = 1 − ⎜
⎝
∑
⎜ n = 0 n!
p0 +
n=5
s! s∑(n − s)
p 0⎟=
⎟
⎠
= 1 − (0,66937 + 0,24256) = 0, 08807
4
(λ / μ)n
4
2,88n ⎛ 2,882 2,883 2,88 4 ⎞
∑
n=0
n!
p0 = ∑
n=0
n!
X 0,045 = ⎜ 1 + 2,88 +
⎝ 2
+
6
+
24
⎟ X 0,045 = 0,66937
⎠
Con condición de operación FIR, las llegadas aleatorias a la pista son de 24 a la hora,
mientras que el tiempo de servicio definido por separación de llegadas es de 120
segundos con una desviación estándar de 20 segundos.
20 1
σ = 20 segundos = = horas
60 x 60 180
λ 2σ 2 + ρ2 242 (1 / 180)2 + 0, 82
Número de aviones en cola L q = = = 1,644
2(1 − ρ) 2(1 − 0,8)
λ
Número medio de aviones en sistema L s = L q + = 1,644 + 0,8 = 2, 444
μ
L q 1,644
Tiempo medio de aviones en cola Wq = = = 0,0685 horas = 4,11 minutos
λ 24
Wq = = = = 4,11 minutos
λ 2 [1 − (λ / μ)] 2 [1 − (24 / 30)]
L s 2, 444
Tiempo medio aviones en sistema Ws = = x (60) = 6,11 minutos
λ 24
λ
El factor de utilización ρ = > 1 , la demanda es mayor que la capacidad, por lo que es
μ
necesario utilizar un modelo determinístico de teoría de colas.
t
Estado de la cola: L t = ∫ (λ t − μ t ) dt
0
En la práctica, L t = L t − 1 + (λ t − μ t ) Δt
se estima con facilidad en una hoja de
cálculo.
Longitud cola:
L t = distancia vertical (líneas de demanda, capacidad acumuladas)
Tiempo espera:
Wt = distancia horizontal (líneas de demanda, capacidad acumuladas)
Demora total Td
Td
Tiempo promedio de espera (demora) Wq = Nd ≡ número pasajeros afectados
Nd
Td 500
Wq = = = 0,25 horas = 15 minutos
Nd 2000
Td
Promedio de pasajeros en cola: L q = t d ≡ tiempo duración de la cola
td
Td 500
Promedio pasajeros en cola L q = = = 250 pasajeros
td 2
Si los clientes en los nodos 1 y 2 son variables aleatorias independientes se verifica que
rn, m = rn, 0 . r0, m , propiedad que verifica las ecuaciones de equilibrio.
λ λ 2
λ
L red = ∑ (n + m) rn, m = ∑ n rn, 0 + ∑ m r0, m = + = ∑ μi − λ
n, m n m μ 1 − λ μ 2 − λ i =1
2
1
Tiempo medio de clientes sistema en tándem: Wsistema en red = ∑ (μi − λ)
i =1
2
λ2 2
ρi2
Número medio de clientes colas en tándem: L cola en red = ∑ μi (μi − λ) = ∑
i =1 i = 1 (1 − ρi )
2
λ
Tiempo medio de clientes colas en tándem: Wcola en red = ∑ μi (μi − λ)
i =1
L red
Tiempo medio de cliente en la red: Wred = L red = L cola en red + (ρ1 + ρ2 )
λ
⎛ 1 1 ⎞ L cola en red
Tiempo medio de cliente en cola: c Wq = Wred − ⎜ + ⎟ Wq =
⎝ μ1 μ2 ⎠ λ
e
54
a) λ = 54 clientes/hora = = 0,9 clientes/minuto λ = λ1 = λ2 = λ 3
60
λ2 ρi2
Número máximo de clientes en cada nodo: L qi = =
μi (μi − λ ) (1 − ρi )
ρ23
L q3 = = 2 → ρ23 + 2 ρ3 − 2 = 0 ⇒ ρ3 = 0,732 ρ3 = −2,732
(1 − ρ3 )
0,9 0,9
Intensidad de tráfico nodo 3: ρ3 = ⇒ μ3 = = 1,2295 minutos
μ3 0,732
1 1
Ws3 = = = 3,0349 minutos
μi − λ 1,2295 − 0,9
ρ2i 3
b) Número medio clientes en nodos: L cola en red = ∑ = L q1 + L q2 + L q3
i=1 (1 − ρ i )
1
= 1 minuto → μ1 = 1 minuto
μ1
1 30
= 30 segundos = = 0,5 minutos → μ2 = 2 minutos
μ2 60
∑
1 1 1 1
c) Wred = = + + = 13,9441 minutos
i =1
μi − λ 1 − 0,9 2 − 0,9 1,2295 − 0,9
L red 12,5496
Wred = = = 13,9441 minutos
λ 0,9
Solución:
a) La empresa de ITV se puede modelizar como una red de colas en tándem con tres
nodos (subsistemas), el nodo 1 un modelo de cola M / M / 1 y los nodos 2 y 3 un
modelo de cola M / M / 2
L s3
b) Tiempo medio en el sistema de la oficina (cola + servicio): Ws3 =
λ3
1
= 2 minutos/cliente → μ 3 = 0,5 cliente/minuto s 3 = 2 servidores
μ3
λ3 0,95
Factor de utilización o intensidad tráfico oficina: ρ3 = = = 0,95 < 1 con lo
s 3 . μ 3 2 . 0,5
que la oficina (subsistema) no se satura, existe un estado estacionario.
λ 3 0,95
Utilización promedio de la oficina: us3 = = = 1,9
μ3 0,5
Probabilidad de que ningún cliente se encuentre en la cola de la oficina:
1 1
p 03 = 1 = = 0,026
(1,9) (1,9) ⎛
n 2
⎞ 1 + 1,9 + 36,1
∑
2 . 0,5
+ ⎜ ⎟
n=0
n! 2 ⎝ 2 . 0,5 − 0,95 ⎠
Número promedio de clientes en cola de la oficina:
s
1 ⎛ λ3 ⎞ ρ3
3
1 0,95
L q3 = ⎜ ⎟ p 03 = x 1,92 x x 0,026 = 17,833 clientes
s3 ! ⎝ μ 3 ⎠ (1 − ρ3 )2
2 (1 − 0,95)2
c) El tiempo medio de un cliente en la ITV es la suma de los tiempos medios en los tres
3
nodos (subsistemas): Wred = ∑ Ws i
i=1
3
Wred = ∑ Ws i = 20 + 2,223 + 20,771 = 43 minutos
i=1
d)
Atendiendo a la intensidad del tráfico, al ser mayor que 0,85, habría que aumentar la
capacidad de los subsistemas (nodos) 1 y 3, es decir, habría que añadir servidores.
El nodo 1 pasa de ser un modelo de cola M / M / 2 y el nodo 3 a un modelo de cola
M/M/3
Nodo 1:
λ 1 0,95
Utilización promedio de la caseta: us1 = = = 0,95
μ1 1
Nodo 3:
λ 3 0,95
Utilización promedio de la oficina: us3 = = = 1,9
μ3 0,5
λ3 0,95
Intensidad del tráfico de la oficina: ρ3 = = = 0,633
s 3 . μ 3 3 . 0,5
Probabilidad de que ningún cliente se encuentre en la oficina:
1 1
p 03 = s − 1 = 2 =
3
(λ 3 / μ 3 )n (λ 3 / μ 3 )s ⎛ s 3 . μ 3 ⎞ 1,9n 1,93 ⎛ ⎞
∑ ∑
3 x 0,5
+ ⎜ ⎟ + ⎜ ⎟
n=0
n! s 3 ! s .
⎝ 3 3 μ − λ 3 ⎠ n=0
n! 6 ⎝ 3 x 0,5 − 0,95 ⎠
1
= = 0,1278
1,92
1 + 1,9 + + 3,1177
2
Número promedio de clientes en la oficina:
s
1 ⎛ λ3 ⎞ ρ3
3
1 0,633
L q3 = ⎜ ⎟ p 03 = x 1,93 x x 0,1278 = 0,6866
s3 ! ⎝ μ 3 ⎠ (1 − ρ3 )2
6 (1 − 0,633)2
SUPUESTOS CONSIDERADOS
∗ Capacidad infinita en los nodos.
∗ Efecto Bloqueo: Si un cliente ha finalizado su servicio en el nodo i y se dirige a un
nodo j que está al máximo de su capacidad.
El sistema se bloquea con tres posbilidades:
(a) Las llegadas al nodo i se rechazan.
(b) El cliente debe ir inmediatamente a otro nodo en su lugar.
(c) El cliente debe abandonar el sistema.
Las k ecuaciones anteriores forman un sistema lineal con solución única, que se resuelve
para hallar las tasas de llegada a cada nodo Λ i
Las k ecuaciones anteriores forman un sistema lineal con solución única, que se resuelve
para hallar las tasas de llegada a cada nodo Λ i
MEDIDAS DE RENDIMIENTO
Las medidas de rendimiento para cada nodo se calculan según las ecuaciones del modelo
M / M / s , teniendo las siguientes consideraciones:
k
• Distribución del número de clientes en cada nodo: p (n) = ∏ pi (ni ) ∀ n1 , … , nk ≥ 0
i=1
pi (ni ) ≡ probabilidad de que haya ni clientes en el nodo i
k
• λ red = ∑ λ i ≡ Número de llegadas que entran en la red por unidad de tiempo
i=1
• L red ≡ Número medio de clientes en el sistema (cola + servicio), suma del número
k
medio de clientes en cada uno de los nodos: L red = ∑L
i=1
si
• Wred ≡ Tiempo medio en el sistema, tiempo medio que un cliente pasa desde que
L
entra en la red hasta que sale de ella: Wred = red
Λ red
• Vi ≡ Número medio de clientes que visitan el nodo i desde que entran en la red hasta
Λ
que salen: Vi = i ∀ i = 1, 2, ... , k
Λ red
El hecho de que los nodos se comporten como si fueran modelo M / M / s podría
interpretarse que se puede utilizar las distribuciones de los tiempos de espera de estos
modelos. Sin embargo, esto no es necesariamente cierto en las redes de Jackson, donde
se permite la realimentación
L qi Λi ⎛ 1⎞
Tiempo medio espera en cola de nodo i: Wqi = = ⎜ Wsi = Wqi + ⎟
Λi μi (μi − Λ i ) ⎝ μi ⎠
L si 1
Tiempo medio de espera en cada nodo (subsistema): Wsi = =
Λi μi − Λ i
2
Ecuaciones de tráfico o ecuaciones de equilibrio: Λ i = λ i + ∑Λ r
j=1
j ji
⎛ Λ 1 ⎞ ⎛ λ 1 ⎞ ⎛ r11 r21 ⎞ ⎛ Λ 1 ⎞
En forma matricial: ⎜ ⎟ = ⎜ ⎟ + ⎜ ⎟⎜ ⎟
⎝ Λ 2 ⎠ ⎝ λ 2 ⎠ ⎝ r12 r22 ⎠ ⎝ Λ 2 ⎠
⎧ Λ 1 = 20 + 0,2 . Λ 1 + 0,25 . Λ 2
⎨
⎩ Λ 2 = 30 + 0,1 . Λ 1
Λ 1 = 35,484 Λ 2 = 33,548
2
λ red = ∑ λ i = 20 + 30 = 50
i=1
b) Número medio que un proceso visita cada nodo, desde que entra hasta que sale de red:
Λ1 35,484
V1 = = = 0,514 veces/minuto
Λ red 69,032
Λ2 33,548
V2 = = = 0,486 veces/minuto
Λ red 69,032
c) Tiempo medio de un proceso en la red (desde que entra hasta que sale):
L red 2,9907
Wred = = = 0,04332 minutos = 3 segundos
Λ red 69,032
Se pide:
a) Número mínimo de servidores en cada nodo para que la red sea estacionaria.
b) Demoras medias en las colas de la red.
Solución:
a) En cada nodo el flujo de entrada debe ser igual al flujo de salida
⎧ λ 1 = 20 λ 2 = 10 λ 3 = 0
Datos del esquema son: ⎨
⎩r12 = 0,5 r13 = 0,5 r23 = 1 r33 = 0,75
⎧Λ1 = λ1 ⎧ Λ 1 = 20
⎪ ⎪
⎨ Λ 2 = λ 2 + Λ 1 r12 → ⎨ Λ 2 = 10 + 20 x 0,5 = 20
⎪Λ = Λ r + Λ r + Λ r ⎪ Λ = 20 x 0,5 + 20 x 1 + Λ x 0,75 → Λ = 120
⎩ 3 1 13 2 23 33 33 ⎩ 3 3 3
⎧ Λ 1 = 20 Λ 2 = 10 + 0,5 x 20 = 20
⎨ → Λ 1 = 20 , Λ 2 = 20 , Λ 3 = 120
⎩ 0,5 x 20 + 20 + 0,75 Λ 3 ⇒ Λ 3 = 120
Tasas totales de llegada a cada subsistema (desde fuera o de otro nodo):
3
Λ red = ∑ Λ i = 20 + 20 + 120 = 160
i=1
Para que ninguna de las colas del sistema se sature, es preciso que se cumpla para cada
uno de los nodos por separado:
Λ
ρi = i < 1 Λ i ≡ Tasa de procesos que salen del nodo i (Tasas totales llegadas)
si μi
El número mínimo de servidores que verifiquen en cada nodo ρi < 1 :
20 20 20
Nodo 1: ρ1 = < 1 → s1 > → s1 = 2 servidores ⇒ ρ1 = = 0,667
15 s1 15 15 x 2
20 20 20
Nodo 2: ρ2 = < 1 → s2 > → s2 = 2 servidores ⇒ ρ2 = = 0,667
15 s2 15 15 x 2
120 120 120
Nodo 3: ρ3 = < 1 → s3 > → s3 = 9 servidores ⇒ ρ3 = = 0,889
15 s3 15 15 x 9
b) Nodo 1 y Nodo 2: Cola M / M / 2
Λ 1 Λ 2 20
Utilización promedio del nodo 1 o 2: us1 = us2 = = = = 1,333
μ1 μ2 15
Probabilidad que ningún cliente se encuentre en el sistema de cola de cada nodo:
1
p 0i = s − 1 si
(Λ i / μi )n 1 ⎛ Λi ⎞
∑
i
1
+ ⎜ ⎟
n=0
n! s i ! ⎝ μ i ⎠ 1 − ρi
1 1
p 01 = p 02 = 2−1
= 1
=
(20 / 15)n (20 / 15)2 ⎛ 2 x 15 ⎞ (1,333)n
∑ n!
+
2!
⎜
⎝ 2 x 15 − 20
⎟
⎠ n=0 n1 ! ∑
+ 2,665
n=0
1 1
= 1 = = 0,2
(1,333)n 1 + 1,333 + 2,665
∑
n=0
n!
+ 2,665
Nodo 3: Cola M / M / 9
Λ 3 120
Utilización promedio del nodo 3: us3 = = =8
μ3 15
Probabilidad que ningún cliente se encuentre en el sistema de la cola del nodo 3:
1 1
p 03 = 9 − 1 s
= 9−1 =
(Λ 3 / μ3 )n3
⎛ Λ3 ⎞ n3
∑
3
∑
1 1 8 1 1
+ ⎜ ⎟ + 9
x 8 x
n3 ! s3 ! ⎝ μ 3 ⎠ 1 − ρ3 n = 0 n3 ! 9! 1 − 0,889
n3 = 0 3
1
= = 0,0002
1766,33 + 3328,81
8
8 n3
∑
n3 =0
n 3 !
= 1 + 8 + 32 + 85,33 + 170,66 + 273,06 + 364,08 + 416,10 + 416,10 = 1766,33
L q3 5,33
Tiempo medio de espera en la cola de nodo 3: Wq3 = = = 0,045
Λ3 120
Solución:
a) En cada nodo el flujo de entrada debe ser igual al flujo de salida
⎧ λ 1 = 0,3 λ 3 = 0,7
Datos del esquema son: ⎨
⎩ r12 = 0,3 r13 = 0,7 r34 = 0,6 r35 = 0,4 r53 = 0,8
Intuitivamente:
Λ1 = λ1
Λ 2 = Λ 1 r12
Λ 3 = λ 3 + Λ 1 r13 + Λ 5 r53
Λ 4 = Λ 3 r34
Λ 5 = Λ 3 r35
Li 1
Tiempo medio de espera en cada nodo (subsistema): Wsi = =
Λ i μi − Λ i
1 1 1 1
Ws1 = = = 0,3704 Ws2 = = = 0,3436
μ1 − Λ 1 3 − 0,3 μ2 − Λ 2 3 − 0,09
L qi Λi ⎛ 1⎞
Tiempo medio de espera en cola de nodo: Wq i = = ⎜ si
W = Wqi + ⎟
Λi μi (μi − Λ i ) ⎝ μi ⎠
Λ1 0,3 1 1
Wq1 = = = 0,0371 Wq1 = Ws1 − = 0,3704 − = 0,0371
μ1 (μ1 − Λ 1 ) 3(3 − 0,3) μ1 3
Λ2 0,09 1 1
Wq2 = = = 0,0103 Wq2 = Ws2 − = 0,3436 − = 0,0103
μ2 (μ2 − Λ 2 ) 3(3 − 0,09) μ2 3
Λ3 1,338 1 1
Wq3 = = = 0,1257 Wq3 = Ws3 − = 0,3757 − = 0,1257
μ 3 (μ 3 − Λ 3 ) 4 (4 − 1,338) μ3 4
Λ4 0,803 1 1
Wq4 = = = 0,1219 Wq4 = Ws 4 − = 0,4552 − = 0,1219
μ 4 (μ 4 − Λ 4 ) 3(3 − 0,803) μ4 3
Λ5 0,535 1 1
Wq5 = = = 0,0386 Wq5 = Ws5 − = 0,2886 − = 0,0386
μ 5 (μ 5 − Λ 5 ) 4 (4 − 0,535) μ5 4
5 5
⎛ 1⎞
Wred = ∑
i=1
Wsi = ∑
i=1
⎜ qi
⎝
W +
μ i
⎟ = 1,8325
⎠
⎛ Λ 1 ⎞ ⎛ 10 ⎞ ⎛ 0 1 0 0⎞ ⎛ Λ1 ⎞ ⎧ Λ 1 = 10 + Λ 2
⎜ ⎟ ⎜ ⎟ ⎜ 0 0 0 ⎟⎟ ⎜Λ ⎟ ⎪ Λ = 0,2 . Λ
⎜ Λ 2 ⎟ = ⎜ 0 ⎟ + ⎜ 0,2 ⎜ 2⎟ →
⎪ 2
⎨
1
⎜ Λ 3 ⎟ ⎜ 20 ⎟ ⎜ 0,8 . 0,6 0 0 0⎟ ⎜ Λ3 ⎟ ⎪ Λ 3 = 20 + 0,8 . 0,6 . Λ 1
⎜ ⎟ ⎜ ⎟ ⎜ ⎟ ⎜ ⎟
⎝ Λ 4 ⎠ ⎝ 0 ⎠ ⎝ 0,8 . 0,4 0 0 0⎠ ⎝ Λ4 ⎠ ⎩⎪ Λ 4 = 0,8 . 0,4 . Λ 1
Λi
Para que la red no se sature en cada nodo (subsistema): ρi = <1 i = 1, 2, 3 , 4
μi
Λ 1 12,5 Λ 2,5 Λ 26 Λ 4
ρ1 = = = 0,5 ρ2 = 2 = = 0,25 ρ3 = 3 = = 0,65 ρ 4 = 4 = = 0,8
μ1 25 μ2 10 μ 3 40 μ4 5
La red no se satura en ningún nodo, existe una distribución estacionaria.
L
∑ L si
i=1
Tiempo medio de clientes en la red: Wred = red =
Λ red 4
∑ Λi
i=1
ρi Λi
Número medio de trabajos en el sistema (cola + servicio): L si = =
1 − ρi μ i − Λ i
0,5 0,25 0,65 0,8
L s1 = = 1 L s2 = = 0,3333 L s3 = = 1,8571 L s4 = =4
1 − 0,5 1 − 0,25 1 − 0,65 1 − 0,8
∑ n!
+
s!
⎜
s
⎝ 3 3μ − Λ
⎟
3 ⎠ n=0
n!
+∑ 2!
⎜ ⎟
⎝ 2 . 16 − 26 ⎠
n=0
1 1
= 1 = = 0,103
1,625n 1,6252 1 + 1,625 + 7,041
∑
n=0
n!
+
2
. 5,333
L red 10,101
Tiempo medio de un trabajo en red: Wred = = = 0,022 unidades tiempo
Λ red 45
Solución:
5
Las ecuaciones de tráfico o equilibrio son intuitivas: Λ i = λ i + ∑Λ r
j=1
j ji
⎧ Λ 1 = λ 1 = 1,5
⎪ Λ = Λ r = 1,5 . 0,2 = 0,3
⎪ 2 1 12
⎪Λ 3 = Λ 1 r13 = 1,5 . 0,8 = 1,2 6
⎨ Λ red = ∑ Λ i = 5,2
⎪ Λ 4 = Λ 3 r34 = 1,2 . 0,6 = 0,72 i=1
⎪ Λ 6 = λ 6 = 0,5
⎪
⎩ Λ 5 = Λ 3 r35 + Λ 6 r65 = 1,2 . 0,4 + 0,5 . 1 = 0,98
Li 1
Tiempo medio de espera en cada nodo (subsistema): Wsi = =
Λ i μi − Λ i
1 1 1 1
Ws1 = = =2 Ws2 = = = 0,5882
μ1 − Λ 1 2 − 1,5 μ2 − Λ 2 2 − 0,3
1 1 1 1
Ws3 = = = 1,25 Ws4 = = = 0,7812
μ 3 − Λ 3 2 − 1,2 μ 4 − Λ 4 2 − 0,72
1 1 1 1
Ws5 = = = 0,9803 Ws6 = = = 0,6666
μ 5 − Λ 5 2 − 0,98 μ 6 − Λ 6 2 − 0,5
1
Tiempo medio de espera en la cola de cada nodo (subsistema): Wqi = Wsi −
μi
1 1 1
Wq1 = 2 − = 1,5 Wq2 = 0,5882 − = 0,0882 Wq3 = 1,25 − = 0,75
2 2 2
1 1 1
Wq4 = 0,7812 − = 0,2812 Wq5 = 0,9803 − = 0,4803 Wq6 = 0,6666 − = 0,1666
2 2 2
1 k ρni
pn1 , n2 , , nk = ∏
G(N) i = 1 ai (ni )
k
ρni i ⎪⎧ n! n ≤ si
donde, G(N) = ∑ ∏ y ai (n) = ⎨ n−s
n ≥ si
n1 + + nk = N i = 1 ai (ni ) ⎪⎩ si ! si i
El cálculo de G(N) puede resultar costoso cuando N y k son grandes, dado que el
⎛N+ k −1⎞
número de posibles estados es ⎜ ⎟
⎝ N ⎠
JP. Buzen desarrolló un algoritmo recursivo para n = N y m = k donde se observa que
gk (N) = G(N) :
m
ρni i ⎧ n! n ≤ si
g m (n) = ∑ ∏ y ai (n) = ⎨ n − si
n ≥ si
n1 + + nm = n i = 1 ai (ni ) ⎩ si ! si
ρni
denotando f i (n) = para i = 1, ... ,k y n = 0, 1, ... ,N
ai (n)
Se observa que g1 (n) = f1 (n) pudiendo aplicar la ecuación de recurrencia. Por otra
parte, gm (0) = fm (0) = 1 m = 1, ... ,k
Lm
9 Tiempo medio de permanencia de un cliente en un nodo: Wm =
Λm
K
El valor calculado Λ i en las ecuaciones de equilibrio Λ i = ∑Λ r
j=1
j ji es una de las
infinitas soluciones no nulas que relaciona las tasas de entrada, no tienen porqué ser el
valor correcto de las Λ i (son valores proporcionales a los Λ i verdaderos).
Aunque el algoritmo de Buzen hace más cómodo el cálculo de G(N) sigue resultando
costoso.
L i (m)
Λ i (m) = i = 1, 2, ,k m = 1, 2, ,M
Wi (m)
Λ i (m) L (m)
ρi (m) ≡ Utilización del servidor en el nodo i‐ésimo: ρi (m) = = i
μi μi Wi (m)
m Λ i Wi (m)
L i (m) = k
i = 1, 2, ,k m = 1, 2, ,M , L i (0) = 0 i = 1, 2, ,k
∑ Λ i Wi (m)
i=1
Wi (m) ≡ Tiempo medio que cada cliente pasa en el nodo i cuando hay m clientes.
1 + L i (m − 1)
Wi (m) = i = 1, 2, ,k m = 1, 2, ,M
μi
Solución:
Ecuaciones de equilibrio:
k
Λi = ∑Λ r
j=1
j ji
Λ 1 = Λ 2 r21 + Λ 4 r41
Λ 2 = Λ 1 r12
Λ 3 = Λ 1 r13
Λ 4 = Λ 3 r34
⎧ Λ 2 = Λ 1 r12 ⎧ Λ 2 = 0,7
⎪ ⎪
Tomando Λ 1 = 1 : ⎨ Λ 3 = Λ 1 r13 → ⎨ Λ 3 = 0,3
⎪ ⎪ Λ = 0,3
⎩ Λ 4 = Λ 3 r34 ⎩ 4
1 + L i (m − 1)
Tiempo de espera en el nodo: Wi (m) = con μi = 5 , i = 1, 2, 3, 4
μi
m Λ i Wi (m)
Número medio de clientes en el nodo: L i (m) = k i = 1, 2, 3, 4
∑ Λ i Wi (m)
i=1
1+0
L i (0) = 0 i = 1, 2, 3, 4 Wi (1) = = 0,2 i = 1, 2, 3, 4
5
1 . 0,2 0,2
L 1 (1) = = = 0,4348
0,2 + 0,7 . 0,2 + 0,3 . 0,2 + 0,3 . 0,2 0,2 . 2,3
1 + L i (1)
Wi (2) = i = 1, 2, 3, 4
5
1 + L 1 (1) 1,4348 1 + L 3 (1) 1,1304
W1 (2) = = = 0,2870 W3 (2) = = = 0,2261
5 5 5 5
♦ Continúan las iteraciones, con una hoja de cálculo como Excel se obtiene:
Destino
1 2 3 4
Origen
1 0,25 0,15 0,20 0,40
2 0,15 0,35 0,20 0,30
3 0,50 0,25 0,15 0,10
4 0,40 0,30 0,25 0,05
Las dos primeras estaciones (servidores) son biprocesadoras, cada uno con un tiempo
de pocesado exponencial y capacidad de 5 procesos/minuto. Las dos últimas
estaciones son monoprocesadoras y pueden atender por minuto respectivamente a 10
y 15 procesos. Se solicita:
a) Modelizar el proceso
b) Número medio de procesos en la cuarta estación.
c) Tiempo medio que transcurre desde que llega un proceso al servidor cuarto
hasta que finaliza su ejecución.
Solución:
a) Se puede modelizar mediante una red de Jackson, donde los clientes son cada uno
de los tres procesos que recorren el sistema, con:
s1 = 2 s2 = 2 s3 = 1 s4 = 1
N=3 k=4
μ1 = 5 μ2 = 5 μ 3 = 10 μ 4 = 15
♦ Ecuaciones de equilibrio:
K 4
Λi
Λi = ∑
j=1
Λ j rj i = ∑
j=1
Λ j rj i ρi =
μi
→ Λ i = ρi . μ i
1081
Λ 1 = 1,5363 Λ 2 = 1,2716 Λ 3 = 1 Λ 4 = = 1,1537
937
♦ De este modo, la utilización del servidor:
Λ 1 1,5363 Λ 2 1,2716
ρ1 = = = 0,3072 ρ2 = = = 0,2543
μ1 5 μ2 5
Λ 1 Λ 1,1537
ρ3 = 3 = = 0,1 ρ4 = 4 = = 0,0769
μ 3 10 μ4 15
ρni
siendo: fi (n) =
ai (n)
⎧⎪ n! n ≤ si
ai (n) = ⎨ n−s
→ a1 (n) = a2 (n) = 2! 2n − 2 = 2n−1 a3 (n) = a4 (n) = 1
⎪⎩ si ! si i n ≥ si
ρni
fi (n) = n − 1 i = 1, 2 f i (n) = ρni i = 3, 4 y n ≥ 1
2
ρni
• Para m = 1 : g1 (n) = f1 (n) = n − 1 i = 1, 2
2
n = 0 : g1 (0) = f1 (0) = 1
ρ12 0,30722
n = 2 : g1 (2) = f1 (2) = 2 − 1 = = 0,0472
2 2
ρ13 0,30723
n = 3 : g1 (3) = f1 (3) = 3 − 1 = = 0,0072
2 22
3
ρni
• Para m = 2 : g 2 (n) = ∑ f2 (i) . g1 (n − i) f i (n) = n − 1
2
i = 1, 2
i=0
ρ2
f2 (0) = 1 f2 (1) = = ρ2 = 0,2543
21 − 1
ρ22 ρ22 0,25432 ρ23 ρ23 0,25433
f2 (2) = 2 − 1 = = = 0,0323 f2 (3) = 3 − 1 = 2 = = 0,0041
2 2 2 2 2 4
n = 0 : g 2 (0) = f2 (0) = 1
1
n = 1 : g 2 (1) = ∑ f2 (i) . g1 (1 − i) = f2 (0) . g1 (1) + f2 (1) . g1 (0) =
i=0
n
• Para m = 3 : g 3 (n) = ∑ f3 (i) . g2 (n − i) f i (n) = ρni
i=0
n = 0 : g 3 (0) = f3 (0) = 1
1
n = 1 : g 3 (1) = ∑ f3 (i) . g2 (1 − i) = f3 (0) . g2 (1) + f3 (1) . g2 (0) =
i=0
n
• Para m = 4 : g 4 (n) = ∑ f4 (i) . g 3 (n − i) f i (n) = ρni
i=0
n = 0 : g 4 (0) = f4 (0) = 1
1
n = 1 : g 4 (1) = ∑ f4 (i) . g 3 (1 − i) = f4 (0) . g 3 (1) + f4 (1) . g 3 (0) =
i=0
f4 (n k ) . g 3 (N − n k )
pm (n k ) = p• • • n k = n k = 0, 1, 2, 3 m = 1, 2, 3, 4
g 4 (N)
Lm
c) Tiempo medio de permanencia de un proceso en un nodo: Wm =
Λm
4,1985
Λ4 = c . Λ4 → 4,1985 = c . 1,1537 → c= = 3,6391
1,1537
de donde,
Λ 1 = c . Λ 1 → Λ 1 = 3,6391 . 1,5363 = 5,5907
Λ 3 = c . Λ 3 → Λ 3 = 3,6391 . 1 = 3,6391
Las personas a menudo estudian un sistema para medir su ejecución o mejorar su operación, o diseñarlo
si es que no existe. Los gerentes o controladores de un sistema desean tener una ayuda disponible para
las operaciones cotidianas, como decidir qué hacer en una fábrica si una máquina importante se avería.
Muchas veces los analistas de simulación encuentran que el proceso para definir el funcionamiento del
sistema (que debe realizarse antes de que se puede empezar a desarrollar el modelo de simulación)
proporciona una gran perspectiva sobre los cambios que tienen que hacerse.
Muchas veces los analistas de simulación encuentran que el proceso para definir el funcionamiento del
sistema (que debe realizarse antes de que se puede empezar a desarrollar el modelo de simulación)
proporciona una gran perspectiva sobre los cambios que tienen que hacerse.
Los sistemas de colas son modelos de sistemas que proporcionan un servicio. Como modelo, pueden
representar sistema donde clientes llegan buscando un servicio de algún tipo y salen después de que
dicho servicio haya sido atendido.
Entre simuladores, QSIM y a QSS de WinQSB.
QSS es un simulador de WinQSB que permite imitar un sistema de colas con múltiples servidores,
diversos tipos de clientes y varias colas. Cada tipo de clientes puede tener distintos tipos de llegadas y
tiempos de servicio que habrá que especificar.
Cuando un cliente sale de un servidor, éste pasa a la siguiente cola o servidor o acaba el servicio. Si la
cola está llena, el cliente se queda en el servidor que está hasta que haya espacio. Si el cliente se va a
unir a otra cola, se le puede asignar un tiempo fijo de llegada a esa cola (tiempo de transferencia).
QSS siempre supone que la población es infinita, por tanto, no permite simular situaciones donde la
población es finita.
Al abrir un nuevo problema, aparece una pantalla donde hay que dar el número total de
Componentes del sistema de colas (nº de tipos de clientes + nº de colas + nº de servidores) y la
unidad de tiempo en la que se va a trabajar.
A continuación aparece una matriz donde hay que introducir los datos del sistema de colas, por
columnas.
Conviene hacer un gráfico del funcionamiento del sistema antes de comenzar a introducir los datos.
10ª columna (Service time distribution): Para cada servidor hay que dar la distribución del tiempo de
servicio. Puede ocurrir que un servidor tenga tiempos de servicio diferentes para clientes de distinto
tipo. Entonces, se indicará la distribución del tiempo de servicio para cada cliente, separada por
comas: Cliente A/exp/0/0.3, Cliente B/normal/1/5
1 1
Tasa de servicio: = = 0,06667 cliente/minuto cajero
μ 15
σ = 0,01
Capacidad = 15 personas/fila
Se considera que el Banco utiliza un sistema de espera titpo PEPS (FIFO – First In First Out), es decir,
los primeros clientes que entran serán los primeros en ser atendidos.
La columna Distribución de los tamaños de los lotes (Batch Size Distribution), indica si los clientes
llegan de forma agrupada o individual.
En este caso, la columna se rellena con Constant/1, indicando que los Clientes llegan al Banco de uno
en uno.
Para programar la Cola, hay que indicar que los dos Cajeros se alimentarán de ella colocando los
nombres en las casillas correspondientes a la columna Sucesor inmediato (Immediate Follower).
En Disciplina de la cola (Queue Discipline) se marca FIFO y en Capacidad de la Cola (Queue Capacity)
su capacidad (máximo 15 personas en espera).
El Cajero 1 atendió 244 Clientes y el Cajero 2 a 217 Clientes, siendo un total (Customer Processed)
de 461 Clientes atendidos.
La opción Results/Show Queue Analysis facilita información sobre la Cola.
El promedio de Clientes en la
Cola (Average Q. Length) es de
1,4261
El máximo de Clientes en la Cola
(Maximun Q. Length) es de 15
El promedio de automóviles
en la cola (Average Q.
Length) fue de 0.
El número máximo de
automóviles en la cola es de
1 (Maximun Q. Length)