Simulation d'événements discrets en Python
Diogo Costa (PhD, MSc)
Adjunct Professor, University of Saskatchewan, Canada & CEO of ImpactBLUE-Scientific
Description d'un système naturel ou humain à l'aide de concepts et du langage mathématiques.
Les modèles se classent en :
Les modèles mathématiques peuvent être simples ou complexes.
George Box : « All models are wrong, but some are useful. »

Modèle d'écoulement d'une rivière

Ce modèle prédit le débit à la station hydrologique de Goaya, en Chine.
La prédiction repose sur les précipitations utilisées comme intrants du modèle.
Autres exemples
Prédire l'inflation

Autres exemples
# Define model parameters
processes = {"process_1": 5,
"process_2": 2,
"process_3": 3}
# Simulation period
simulation_time = 365
# Run model
discrete_model(processes, simulation_time)
def discrete_model(processes, simulation_time): # 1) Run end-condition while (time < simulation_time): process_names = list(processes.keys())# 2) Loop over all processes for p in range(len(process_names)): process_name_p = process_names[p]# 3) Account for effect of each process time += processes[process_name_p]
Exemple de sortie d'un modèle à événements discrets pour une activité manufacturière
=> START OF SIMULATION (Time = 0 days)
Time = 6.00 days | Process Complete: Transport of raw material
Time = 9.00 days | Process Complete: Building components
Time = 11.00 days | Process Complete: Assembling parts
Time = 14.00 days | Process Complete: Selling product
=> COMPLETED: Supply-Chain cycle #1 | Time = 15.5 days
Time = 21.50 days | Process Complete: Transport raw material
Time = 24.50 days | Process Complete: Building components
Time = 26.50 days | Process Complete: Assembling parts
Time = 29.50 days | Process Complete: Selling product
=> COMPLETED: Supply-Chain cycle #2 | Time = 31.0 days
Visualisation graphique des résultats : utile pour repérer des motifs et des seuils de basculement du système.
La visualisation doit être adaptée aux objectifs de simulation.
De nombreuses bibliothèques sont utiles : matplotlib, seaborn, plotly.

Exemples
y) selon le temps correspondant (x)plt.plot(x, y, color='green', marker='o',
markersize=12, linestyle='dashed',
linewidth=2)
x et compter le nombre de valeurs par classeplt.hist(x, 50, density=True,
facecolor='g', alpha=0.75)
Simulation d'événements discrets en Python