Simulación de eventos discretos

Mecanismo de avance del tiempo

Mecanismo siguiente evento

(eventos futuros)

Prof. Luis Ulfe

 Evento Discreto y Continuo
Evento Discreto.
• Acción instantánea que ocurre en un punto en el tiempo. Un
cliente arribando a una agencia, una pieza que será procesada,
un auto en un peaje, etc.
• La ocurrencia de estos eventos puede causar que los estados
del sistema cambien.
• Se define la variable “reloj de la simulación”, el cual avanza con
la ocurrencia de un evento.

Evento Continuo.
• Es una acción sin finalización. Este continúa ininterrumpido con
respecto al tiempo, generalmente representado por ecuaciones
diferenciales.
• Tenemos como ejemplo el flujo de jarabe en una fábrica de
bebidas gaseosas, petróleo en un oleoducto
• La mayoría de procesos continuos pueden ser representados a
manera de sumatorias como procesos discretos.
 Simulación de Eventos Discretos
• Modelo en el cual la el estado de las variables
cambian instantáneamente en un punto del tiempo.
• Los eventos definidos como una ocurrencia
instantánea que pueden cambiar el estado del
sistema.
Ejem.
• Colas.
• Variable del estado.
• Estadísticos.
• Arribos - servicios - cambian el estado del sistema.
 Simulación de Eventos Discretos
Mecanismos de avance del tiempo
• Necesitamos rastrear el tiempo simulado y
determinar el siguiente tiempo
• Reloj de la simulación: Variable en el modelo de
simulación que nos da el valor actual del tiempo
simulado
• Especificar las unidades de tiempo

Dos enfoques principales:

• Mecanismo de avance del siguiente evento
• Mecanismo de avance de incremento fijo
 Simulación de Eventos Discretos
Mecanismo del siguiente evento:
• El reloj de la simulación es inicializado en cero.
• Se determinan los tiempos de ocurrencia de eventos
futuros.
• El reloj de la simulación se avanza al tiempo de ocurrencia
del evento más inminente (primero de los eventos futuros).
• En este punto el estado del sistema es actualizado para
contabilizar la ocurrencia del evento.
• Los tiempos de ocurrencia de eventos futuros es
actualizada.
• Luego, el reloj de la simulación avanza al tiempo del (nuevo)
evento más inminente.
• El estado del sistema es actualizado y se determinan los
eventos futuros etc.
• Este proceso continua hasta satisfacer alguna condición de
parada.
Simulación de Eventos Discretos
Obs:
• Los periodos de inactividad son omitidos, saltando el
reloj de un tiempo del evento a otro
• El mecanismo del tiempo de incremento fijo no omite
estos periodos de inactividad, consumiendo más
tiempo de computador.
 Simulación de Eventos Discretos
Ejem. Para un mecanismo de cola simple
• Ti = tiempo de arribo del i - ésimo cliente (t0 = 0)
• Ai = ti – t(i – 1) = tiempo entre arribos entre el (i - 1)
Cliente y el i cliente
• Si = tiempo requerido por el servidor para el cliente i
(no considera demora en cola )
• Ci = ti +Di + Si = tiempo en que el cliente i completa
su servicio y se retira
• Di = demora en la cola del cliente i
• ei =tiempo de ocurrencia del evento i de cualquier tipo
(el valor i que el reloj de la simulación toma
excluyendo el valor e0 = 0)
Simulación de Eventos Discretos
• Generar variables aleatorias de tiempos entre arribos
y de servicios.
FA → 1 ,  2....
Fs → S1, S 2.....
e0 e1 e2 e3 e4 e5

t0 t1 t2 C1 t3 C2

A A A
1 2 3

S1 S1
Simulación de Eventos Discretos
Simulación de un sistema de colas con un servidor
• Sean las Variables Aleatorias (VA) del tiempo entre
arribos y servicio respectivamente, dados abajo.
• Independientes e idénticamente distribuidos (IID).
• Política de cola FIFO.
• La simulación empezara en el estado "vacío o
desocupado"
• La simulación termina cuando se tengan n
observaciones de demora en cola.
 Simulación de un sistema de colas con un servidor
Para n = 6 observaciones de demoras en cola e54312 =
8==ttt2t464t379865==
===1.6
3.8
2.1
Tiempo entre arribos: A1=0.4, A2=1.2, A3=0.5, A4=1.7,
(e 873 =
(e3211210
11 5.6)+
==2.1)
0.4)
1.6)
2.1)
3.8)
4)
5.8)
7.2) ++(A(A
(A243=1.2)
987654=1.7)
=0.5)
=1.7) ee2e13
=1.6)
=0.2)
=1.4)
=1.9) 976=
=
10
11 3.8
5.6
5.8
4
5.6
7.2
9.1
A5=0.2, A6=1.6, A7=0.2, A8=1.4, A9=1.9 (e
(e 1 = 0.4)
(e41975 =
= 2.4)
0.4) +
+ (S
(S21543=
(S 1==0.7)
2)
2) eeee213
0.2)
1.1)
3.7) 986==
1154310 ==ccc1c52143=5=
==2.4
2.4
3.1
8.6
3.3
4.9
8.6
Q(t) 3.1)
3.8)
4.9) + 12
Tiempo de servicios: S1= 2, S2=0.7, S3=0.2,
S4=1.1, S5=3.7, S6=0.6
3

0 t
Arribos

A1 A2 A3 A4 A5 A6 A7 A8
0.4 1.2 0.5 1.7 0.2 1.6 0.2 1.4
e1 e2 e3 e7 e8 e10 e11 e12
t1 = 0.4 t2 = 1.6 t3 = 2.1 t4 = 3.8 t5 = 4 t6 =5.6 t7 =5.8 t8 =7.2
Salidas

S1 S2 S3 S4 S5
2 0.7 0.2 1.1 3.7
e4 e5 e6 e9 e13
c1 = 2.4 c2=3.1 c3=3.3 c4=4.9 c5=8.6
B(t)
1
t
0
PERFORMANCE
DEL SISTEMA
 D1 = 0
• Demora promedio de espera en la cola d(n): D2 = C1 – t2 = 2.4 - 1.6 = 0.8
n
D3 = C2 – t3 = 3.1 – 2.1 = 1
 Di d ( 6) =
5 .7
= 0.95 D4 = 0
Q(t)
d ( n) = i =1
6 D5 = C4 – t5 = 4.9 - 4 = 0.9
n D6 = C5 – t6 = 8.6 - 5.6 = 3
3

0 t
Arribos

A1 A2 A3 A4 A5 A6 A7 A8
0.4 1.2 0.5 1.7 0.2 1.6 0.2 1.4
e1 e2 e3 e7 e8 e10 e11 e12
t1 = 0.4 t2 = 1.6 t3 = 2.1 t4 = 3.8 t5 = 4 t6 =5.6 t7 =5.8 t8 =7.2
Salidas

S1 S2 S3 S4 S5
2 0.7 0.2 1.1 3.7
e4 e5 e6 e9 e13
c1 = 2.4 c2=3.1 c3=3.3 c4=4.9 c5=8.6
B(t)
1
t
0
 T(n): tiempo requerido para observar
Promedio clientes en cola q(n):
 “n” demoras en cola.
  iT i

q ( 6) =
9. 9
= 1.15
Ti : tiempo total durante la simulación
q ( n) = i =1
T (n )
8. 6 en que la cola es de longitud i.

T = (1.6 − 0.0) + (4.0 − 3.1) + (5.6 − 4.9) = 3.2
 iTi = (0 * 3.2) + (1* 2.3) + (2 *1.7) + (3 *1.4) = 9.9
0

i =0 T = (2.1 − 1.6) + (3.1 − 2.4) + (4.9 − 4.0) + (5.8 − 5.6) = 2.3

1

Q(t) T = (2.4 − 2.1) + (7.2 − 5.8) = 1.7

2

3 T = (8.9 − 7.2) = 1.4

3

(T = 0, i  4)
i
2

0 t
Arribos

A1 A2 A3 A4 A5 A6 A7 A8
0.4 1.2 0.5 1.7 0.2 1.6 0.2 1.4
e1 e2 e3 e7 e8 e10 e11 e12
t1 = 0.4 t2 = 1.6 t3 = 2.1 t4 = 3.8 t5 = 4 t6 =5.6 t7 =5.8 t8 =7.2
Salidas

S1 S2 S3 S4 S5
2 0.7 0.2 1.1 3.7
e4 e5 e6 e9 e13
c1 = 2.4 c2=3.1 c3=3.3 c4=4.9 c5=8.6
•Utilización del Servidor: U(n)
• Proporción de tiempo durante la simulación (desde el
tiempo 0 hasta T(n) que el servidor esta ocupado. (no
ocioso) 0  U (n)  1

• Un estimador de U(n) es U (n)

U (n) = la proporción de tiempo observado durante la
B(t)
simulación que el servidor esta ocupado

e1 = 0.4 e6 = 3.3 e7 = 3.8 e13 = 8.6 = T(6)

 (3.3 − 0.4) + (8.6 − 3.8) 7.7

U ( n) = = = 0. 9
8.6 8.6