1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
"""bank17: Plotting a Histogram of Monitor results"""
from SimPy.Simulation import *
from SimPy.SimPlot import *
from random import expovariate, seed
## Model components ------------------------
class Source(Process):
""" Source generates customers randomly"""
def generate(self, number, rate):
for i in range(number):
c = Customer(name="Customer%02d" % (i))
activate(c, c.visit(timeInBank=12.0))
yield hold, self, expovariate(rate)
class Customer(Process):
""" Customer arrives, is served and leaves """
def visit(self, timeInBank):
arrive = now()
#print("%8.4f %s: Arrived " % (now(), self.name))
yield request, self, counter
#print("%8.4f %s: Got counter " % (now(), self.name))
tib = expovariate(1.0 / timeInBank)
yield hold, self, tib
yield release, self, counter
#print("%8.4f %s: Finished " % (now(), self.name))
t = now() - arrive
Mon.observe(t)
## Experiment data -------------------------
maxTime = 400.0 # minutes
counter = Resource(1, name="Clerk")
Mon = Monitor('Time in the Bank')
N = 0
## Model ----------------------------------
def model(SEED=393939):
seed(SEED)
initialize()
source = Source()
activate(source,
source.generate(number=20, rate=0.1), at=0.0)
simulate(until=maxTime)
## Experiment ----------------------------------
model()
Histo = Mon.histogram(low=0.0, high=200.0, nbins=20)
plt = SimPlot()
plt.plotHistogram(Histo, xlab='Time (min)',
title="Time in the Bank",
color="red", width=2)
plt.mainloop()
|
