Guía de Laboratorio: 06
Simulación de Procesos
FACULTAD DE INGENIERÍA INDUSTRIAL
Tema: Simulación con R
Fecha Fecha Inicio del Laboratorio: 13/09/2021
Actualización: Fecha Terminación del Laboratorio: 13/09/2021
13/09/2021
Docente
Maryory Patricia Villamizar León
Estudiante
DAIRY YAMILE CHAPARRO VILLAMIZAR
OBJETIVO DEL LABORATORIO
Practica de modelos de línea de espera con R
Ejercicios D/D/1
El departamento de exámenes médicos de un hospital cita a los pacientes cada 5 min comenzando a las 9:00 a.m. para practicarles una revisión general, que
tarda 7 min. Durante un lapso de 2h de funcionamiento:
a) ¿Cuál será el número promedio de pacientes en espera y siendo atendidos?
b) ¿De qué tamaño será la línea de espera a las 2h?
Hora Paciente en espera Paciente en revisión
9:00 - A
9:05 B A
9:07 - B
9:10 C B
9:14 - C
9:15 - C
9:20 D, E C
9:21 E D
9:25 E, F D
9:28 F E
9:30 F, G E
9:35 G, H F
9:40 G, H, I F
9:42 H, I G
9:45 H, I, J G
9:49 I, J H
9:50 I, J, K H
9:55 I, J, K, L H
9:56 J, K, L I
10:00 J, K, L, M I
10:03 K, L, M J
10:05 K, L, M, N J
10:10 L, M, N, O K
10:15 L, M, N, O, P K
10:17 M, N, O, P L
10:20 M, N, O, P, Q L
10:24 N, O, P, Q M
10:25 N, O, P, Q, R M
10:30 N, O, P, Q, R, S M
10:31 O, P, Q, R, S N
10:35 O, P, Q, R, S, T N
10:38 P, Q, R, S, T O
10:40 P, Q, R, S, T, U O
10:45 Q, R, S, T, U, V P
10:50 Q, R, S, T, U, V, W P
10:52 R, S, T, U, V, W Q
10:55 R, S, T, U, V, W, X Q
10:59 S, T, U, V, W, X R
11:00 S, T, U, W, X, Y R
a) En la tabla se observa que durante el lapso de las dos horas hay en todo momento un paciente en revisión.
Sin embargo, en la línea de espera hay minutos sin ningún paciente, incluso algunos periodos de hasta siete
pacientes, por tanto, para obtener el promedio se deben cuantificar los minutos de cero hasta siete
pacientes de las 9:00 am a las 11:00 horas. Al efectuar esto se tendrá:
Cero pacientes: 5+3+1= 9 [minutos]
Un paciente: 2+4+5+4+2= 17 [minutos]
Dos pacientes: 1+3+5+5+3+1= 18 [minutos]
Tres pacientes: 2+4+5+4+2= 17 [minutos]
Cuatro Pacientes: 1+3+5+5+3+1= 18 [minutos]
Cinco Pacientes: 4+5+2+4+2= 17 [minutos]
Seis Pacientes: 1+3+5+5+3+1= 18 [minutos]
Siete Pacientes: 2+4= 6 [minutos]
Aplicando la ecuación:
0∗9+1∗17+2∗18+ 3∗17+ 4∗18+5∗17+6∗18+7∗6
X M= =¿
9+ 17+18+17+18+17+ 18+6
El número promedio de pacientes en espera es: 3.
b) De la tabla inicial se ve que a las 11:00h, la línea de espera cuenta con: 7, incluyendo al Y, que se incorpora
justo a esa hora.
11:00 S, T, U, W, X, Y R
Ejercicio D/D/K
Al departamento de control de calidad de la empresa Envases Plásticos llega cada hora 4 lotes de envases para su
inspección, que se realiza por parte de dos encargados, cada uno de los cuales tarda 34 minutos en la revisión.
En un turno de 8 horas:
a) ¿Cuántos lotes de envases en espera de ser revisados habrá al final de este?
b) ¿Cuál será el número promedio de lotes en la línea de espera a lo largo del turno?
Solución: para comprobar ambas preguntas lo primero que se hará será simular la operación del sistema durante un
turno, cuya situación se presenta a continuación:
Tiempo de turno en Lotes en espera de Lote en revisión Lote en revisión
horas revisión inspector 1 inspector 2
0:00 C, D A B
0:34 - C D
1:00 E, F, G, H C D
1:08 G, H E F
1:42 - G H
2:00 I, J, K, L G H
2:16 K, L I J
2:50 - K L
3:00 M, N, O, P K L
3:24 O, P M N
3:58 - O P
4:00 Q, R, S, T O P
4:32 S, T Q R
5:00 S, T, U, V, W, X Q R
5:06 U, V, W, X S T
5:40 W, X U V
6:00 W, X, Y, Z, A1, B1 U V
6:14 Y, Z, A1, B1 W X
6:48 A1, B1 Y Z
7:00 A1, B1, C1, D1, E1, F1 Y Z
7:22 C1, D1, E1, F1 A1 B1
7:56 E1, F1 C1 D1
8:00 E1, F1, G1, H1,I1, J1 C1 D1
Tiempo con cero lotes: 26+18+10+2= 56 [minutos]
Tiempo con dos lotes: 34+ 34+ 34+ 34+ 28+ 20+ 12+ 4= 200 [minutos]
Tiempo con cuatro lotes: 8+16+24+32+34+34+34= 182 [minutos]
Tiempo con seis lotes: 6+14+22= 42 [minutos]
0∗56+2∗200+4∗182+6∗42
X M= =2.875≈ 3
56+200+182+ 42
a) El número promedio de lotes en la línea de espera a lo largo del turno es de: 3
b) Lotes de envases en espera de ser revisados habrá al final de este turno: 6
8:00 E1, F1, G1, H1,I1, J1 C1 D1
Modelo de línea de espera de un solo canal con llegadas de Poisson y
tiempos de servicio exponencial. Modelo M/M/1: Ejemplo
Recuerden que para el problema del Corral tuvimos una tasa media de llegadas de λ=0.75 clientes por
minuto y una tasa media de servicio de μ=1 cliente por minuto. Por tanto, con μ > λ, pueden usarse las
ecuaciones para brindar las características operativas para la línea de espera de un solo canal:
rm(list=ls())
lambda<-45/60
u<-60/60
n<-4
Características Sigl Código Respuesta
Operativas a
Probabilidad de que no P0 #Probabilidad de que no haya unidades en el sistema 0.25
haya unidades o clientes P0_MM1<-round(1-(lambda/u),4)
en el sistema P0_MM1
Cantidad promedio de Lq #Cantidad promedio de unidades en línea de espera 2
unidades en la línea de Lq_MM1<-round((lambda^2)/(u*(u-lambda)),0) [clientes]
espera Lq_MM1
Cantidad promedio de L #Cantidad promedio de unidades en el sistema 3
unidades en el sistema L_MM1<-round(Lq_MM1+(lambda/u), 0) [clientes]
L_MM1
Tiempo promedio que Wq #Tiempo promedio que pasa una unidad en la línea de 2.6667
pasa una unidad en la espera [minutos]
línea de espera Wq_MM1<-round(Lq_MM1/lambda, 4)
Wq_MM1
Tiempo promedio que W #Tiempo promedio que pasa una unidad en el sistema 3.6667
pasa una unidad en el W_MM1<-round(Wq_MM1+(1/u),4) [minutos]
sistema W_MM1
Probabilidad que tiene Pw #Probabilidad de que una unidad que 0.75
una unidad que llega #llega tenga que esperar por el servicio
tenga que esperar por el Pw_MM1<-round(lambda/u,4)
servicio Pw_MM1
Probabilidad de n Pn #Probabilidad de n unidades en el sistema
unidades en el sistema Pn_MM1<-c() 0.18750000
para: n= 1,2,3,4 for (n in 0:n)
{Pn_MM1[n]<-((lambda/u)^n)*P0_MM1
} 0.14062500
Pn_MM1
0.10546875
0.07910156
Modelo de línea de espera con canales múltiples, con llegadas de Poisson y
tiempos de servicio exponencial. Modelo M/M/k. Ejemplo
Para ilustrar el modelo de línea de espera con canales múltiples, suponga que la administración desea
evaluar la conveniencia de abrir una segunda edición de procesamiento de pedidos ya que hay cambio de
personal en el área de servicio de modo que puede servirse simultáneamente a dos clientes. Suponga una
línea de espera única con el primer cliente en la línea pasando al primer empleado disponible. Evaluamos a
continuación el sistema de dos canales, k = 2, para la tasa media de llegadas de λ=0.75 clientes por minuto y
una tasa media de servicio de 40 clientes en 60 minutos μ=0.666 para cada canal, obtenemos las
características operativas:
rm(list=ls())
lambda<-45/60
u<-40/60
k<-2
Características Sigl Código Respuesta
Operativas a
Probabilidad de P0 #Probabilidad de que no haya unidades en el sistema
que no haya 0.28
unidades o k1<-k-1
clientes en el suP0<-c()
sistema for (i in 0:k1)
{suP0[i]<-((lambda/u)^i)/factorial(i)
}
P0_MMK<-round(1/(1+sum(suP0)+
(((lambda/u)^k)/factorial(k)*(k*u)/(k*u-lambda))),4)
P0_MMK
Cantidad Lq #Cantidad promedio de unidades en línea de espera 1
promedio de Lq_MMK<-round((((lambda/u)^k)*(lambda*u))/((factorial(k- [clientes]
unidades en la 1)*(k*u-lambda)^2))*P0_MMK,0)
línea de espera Lq_MMK
Cantidad L #Cantidad promedio de unidades en el sistema 2
promedio de L_MMK<-round(Lq_MMK+(lambda/u), 0) [clientes]
unidades en el L_MMK
sistema
Tiempo promedio Wq #Tiempo promedio que pasa una unidad en la línea de espera 1.3333
que pasa una Wq_MMK<-round(Lq_MMK/lambda, 4) [minutos]
unidad en la línea Wq_MMK
de espera
Tiempo promedio W #Tiempo promedio que pasa una unidad en el sistema 2.8333
que pasa una W_MMK<-round(Wq_MMK+(1/u),4) [minutos]
unidad en el W_MMK
sistema
Probabilidad que Pw #Probabilidad de que una unidad que
tiene una unidad #llega tenga que esperar por el servicio 0.405
que llega tenga Pw_MMK<-round((1/factorial(k))*((lambda/u)^k)*((k*u)/(k*u-
que esperar por el lambda))*P0_MMK,4)
servicio Pw_MMK
Probabilidad de n Pn #Probabilidad de n unidades en el sistema
unidades en el n<-4 0.3150
sistema para: n= Pn_MMK<-c()
1,2,3,4 for (n in 1:n)
if(n<k){ 0.1772
Pn_MMK[n]<-round(((lambda/u)^n)/(factorial(n))*P0_MMK,4)
}else{ 0.0997
Pn_MMK[n]<-round(((lambda/u)^n)/(factorial(k)*k^(n-
k))*P0_MMK,4) 0.0561
}
Pn_MMK
