[Pyhton] SimPy 離散事件模擬框架詳解 —— 以一個簡單的汽車充電排隊模擬為例

一、背景知識
二、SimPy 講解
三、後續
參考連結
附錄二

一、背景知識

人們在生產活動和社會活動中，經常遇到一類複雜的系統，這類系統中有許多事件時而出現，時而消失，時而動作，時而停止，而啟動和停止都發生在一些離散的時刻，並帶有一定的隨機性。例如，港口中船舶的停靠碼頭、生產線上機床的啟停、電話的接通和斷開、計算機系統中某項作業的進行和退出，凡此種種，都帶有上述特點，這類系統叫做離散事件動態系統(DEDS)。隨著生產和科技的日益發展，以及人類社會交往的日趨頻繁，這類系統的數量和種類也越來越多。

設計此類系統時，往往需要模擬來評估演演算法或方案的效能，常見的軟體有：

商業軟體：

名字	簡介	圖示
AnyLogic	通用多方法建模工具。結合了基於代理、系統動力學和離散事件建模。
Arena	一種離散事件模擬程式，也允許對連續過程進行建模。
Care pathway simulator	專門為服務行業（如醫療保健）設計的離散事件模擬程式。
Enterprise Dynamics	一個模擬軟體平臺，用於模擬和分析幾乎任何製造、材料處理和物流挑戰。
ExtendSim	通用模擬軟體包
DELMIA	3DEXPERIENCE 平臺的一部分
FlexSim	拖拖拉拉做離散事件模擬，3D
GoldSim	將動態離散事件模擬嵌入到 Monte Carlo 框架
GPSS	離散事件模擬語言。供應商可以提供不同的實現
Micro Saint Sharp	通用離散事件建模工具，使用拖放介面和C#程式語言
MS4 Modeling Environment	基於離散事件和混合模型的通用DEVS方法的軟體環境
Plant Simulation	能夠模擬和優化生產系統和流程的軟體
ProModel AutoCAD® Edition	在數位孿生Autodesk®工具集（AutoCAD®和Inventor®）中快速構建空間精確的佈局和過程模擬模型。
Simcad Pro	實時變化，零程式碼有介面，支援 VR
SimEvents	向MATLAB / Simulink環境新增離散事件模擬。
SIMUL8	基於物件的模擬軟體
VisualSim	基於時序、功耗和功能的電子、嵌入式軟體和半導體的基於模型的系統架構探索
WITNESS	可在桌面和雲端使用 VR 進行離散事件模擬

開源軟體：

名字	語言	型別	License	簡介
JaamSim	Java	App	Apasche 2.0	JaamSim是一款自由開源的離散事件模擬軟體，包括拖放式使用者介面、互動式3D圖形、輸入和輸出處理以及模型開發工具和編輯器。
CPN Tools	BETA	App	GPLv2	用於分析所有型別應用中的物流/排隊模型的工具。
DESMO-J	Java	Lib	Apasche 2.0	Java離散事件模擬框架，支援混合事件/過程模型，並提供2D和3D動畫。
Facsimile	Scala	Lib	LGPLv3	離散事件模擬/模擬庫。
PowerDEVS	C++	App	AFL, GPLv2	基於DEVS形式的混合系統建模和模擬整合工具。
Ptolemy II	Java	Lib	BSD	支援面向角色設計實驗的軟體框架。
SIM.JS	JavaScript	Lib	LGPL	JS是一個完全用JavaScript編寫的通用離散事件模擬庫。在瀏覽器中執行，支援基於GUI的建模工具。
SimPy	Python	Lib	MIT	SimPy是基於標準Python的基於過程的離散事件模擬框架。
Simula	Simula	Language		一種專門為模擬而設計的程式語言。
SystemC	C++	Lib	Apache 2.0	提供事件驅動模擬核心。

注：軟體示意圖見《附錄二》

二、SimPy 講解

2.1 SimPy 概述

1）SimPy 是 python 的離散事件模擬框架。
2）SimPy 中的流程由 Python 生成器函數定義。例如，可以用於為客戶、車輛或代理等活動元件建模。SimPy 還提供各種型別的共用資源來模擬容量有限的擁塞點（如伺服器、結賬櫃檯和隧道）。
3）模擬可以要多快有多快、實時、或者手動步進三種方式進行。
3）儘管理論上 SimPy 可以用於連續模擬，但它沒有任何功能來幫助您實現這一點。另一方面，SimPy 對於具有固定步長的模擬來說是過火的，在這種情況下，您的流程不會相互互動或共用資源。

模擬兩個時鐘在不同時間間隔滴答作響的簡短範例如下：

>>> import simpy
>>>
>>> def clock(env, name, tick):
...     while True:
...         print(name, env.now)
...         yield env.timeout(tick)
...
>>> env = simpy.Environment()
>>> env.process(clock(env, 'fast', 0.5))
<Process(clock) object at 0x...>
>>> env.process(clock(env, 'slow', 1))
<Process(clock) object at 0x...>
>>> env.run(until=2)
fast 0
slow 0
fast 0.5
slow 1
fast 1.0
fast 1.5

2.2 基本概念

1）SimPy 是一個離散事件模擬庫。活動元件（如車輛、客戶或訊息）的行為是用流程建模的。所有程序都存在於一個環境中。它們通過事件與環境和彼此互動。
2）流程由簡單的 Python generator 描述。你可以將他們稱為過程函數或過程方法，取決於它是函數還是類的方法。在其整個生命週期內，他們產生事件等待被觸發。
3）當一個過程產生一個事件時，該程序就會被掛起。當事件發生時（我們說事件被觸發），SimPy 恢復該過程。多個程序可以等待同一個事件。SimPy 以它們產生該事件的相同順序恢復它們。
4）一個最重要的事件型別就是 Timeout 類事件。它允許在程序給定的時間內休眠（或保持其他狀態）。

2.3 一個汽車開開停停的例子

下面是一個簡單的汽車走走停停的例子，列印其走停的時間戳：

>>> def car(env):
...     while True:
...         print('Start parking at %d' % env.now)
...         parking_duration = 5
...         yield env.timeout(parking_duration)
...
...         print('Start driving at %d' % env.now)
...         trip_duration = 2
...         yield env.timeout(trip_duration)

>>> import simpy
>>> env = simpy.Environment()
>>> env.process(car(env))
<Process(car) object at 0x...>
>>> env.run(until=15)
Start parking at 0
Start driving at 5
Start parking at 7
Start driving at 12
Start parking at 14

2.4 在走走停停過程中增加充電過程（過程互動）

我們在上面汽車例子基礎上引入充電的過程：車走一段時間，停下來充電，電充好了，才能繼續走。
這裡引入了 charge_duration 過程，在該過程中簡單寫了一個超過的掛起事件：

>>> class Car(object):
...     def __init__(self, env):
...         self.env = env
...         # Start the run process everytime an instance is created.
...         self.action = env.process(self.run())
...
...     def run(self):
...         while True:
...             print('Start parking and charging at %d' % self.env.now)
...             charge_duration = 5
...             # We yield the process that process() returns
...             # to wait for it to finish
...             yield self.env.process(self.charge(charge_duration))
...
...             # The charge process has finished and
...             # we can start driving again.
...             print('Start driving at %d' % self.env.now)
...             trip_duration = 2
...             yield self.env.timeout(trip_duration)
...
...     def charge(self, duration):
...         yield self.env.timeout(duration)

>>> import simpy
>>> env = simpy.Environment()
>>> car = Car(env)
>>> env.run(until=15)
Start parking and charging at 0
Start driving at 5
Start parking and charging at 7
Start driving at 12
Start parking and charging at 14

如果我們不想等充電結束，而是想中斷充電過程並開始駕駛，可以使用 SimPy 的 interrupt() 方法來中斷正在執行的程序：

>>> def driver(env, car):
...     yield env.timeout(3)
...     car.action.interrupt()

由於原來的充電過程被中斷會報異常，因此我們要對例外處理下：

...     		try:
...                 yield self.env.process(self.charge(charge_duration))
...             except simpy.Interrupt:
...                 # When we received an interrupt, we stop charging and
...                 # switch to the "driving" state
...                 print('Was interrupted. Hope, the battery is full enough ...')

再次執行：

>>> env = simpy.Environment()
>>> car = Car(env)
>>> env.process(driver(env, car))
<Process(driver) object at 0x...>
>>> env.run(until=15)
Start parking and charging at 0
Was interrupted. Hope, the battery is full enough ...
Start driving at 3
Start parking and charging at 5
Start driving at 10
Start parking and charging at 12

2.5 共用資源

SimPy 提供三種型別的資源，用於解決建模中多個進行希望使用有限資源的問題（例如：加油站汽車場景中的燃油泵）或典型的生產者-消費者問題。

我們還用汽車充電的例子：汽車開到充電樁旁 a battery charging station (BCS)，向兩個充電樁申請使用其一進行充電，如果兩個樁都在被使用，它將會等待直到可用，然後開始充電，然後開走。

>>> def car(env, name, bcs, driving_time, charge_duration):
...     # Simulate driving to the BCS
...     yield env.timeout(driving_time)
...
...     # Request one of its charging spots
...     print('%s arriving at %d' % (name, env.now))
...     with bcs.request() as req:
...         yield req
...
...         # Charge the battery
...         print('%s starting to charge at %s' % (name, env.now))
...         yield env.timeout(charge_duration)
...         print('%s leaving the bcs at %s' % (name, env.now))

備註： bcs.request() 將會產生一個事件，該事件會阻塞直到資源可用，一般情況下使用資源後需要呼叫 release 對資源進行釋放，這裡的 with 語句意味著自動釋放。

我們建立有兩個充電樁的資源：

>>> import simpy
>>> env = simpy.Environment()
>>> bcs = simpy.Resource(env, capacity=2)

然後我們建立 4 輛車：

>>> for i in range(4):
...     env.process(car(env, 'Car %d' % i, bcs, i*2, 5))

最後，我們可以開始模擬了。由於汽車程序在此模擬中都自行終止，因此我們無需指定直到時間——當沒有更多事件時，模擬將自動停止：

>>> env.run()
Car 0 arriving at 0
Car 0 starting to charge at 0
Car 1 arriving at 2
Car 1 starting to charge at 2
Car 2 arriving at 4
Car 0 leaving the bcs at 5
Car 2 starting to charge at 5
Car 3 arriving at 6
Car 1 leaving the bcs at 7
Car 3 starting to charge at 7
Car 2 leaving the bcs at 10
Car 3 leaving the bcs at 12

注意到前兩輛車到達BCS後可以立即開始充電，而2號車和3號車需要等待，符合預期。