蒙特卡罗模拟之排队上厕所问题

2024-08-29 07:38

本文主要是介绍蒙特卡罗模拟之排队上厕所问题,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

蒙特卡罗模拟之排队上厕所问题


'''
电影结束后会有20人上厕所
20个人会在0-10分钟内全部到达厕所
每个人上厕所时间在1-3分钟
模拟只有一个厕所到达时间,等待时间,开始上厕所时间,结束时间
'''
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from matplotlib.patches import Circlearrivingtime = np.random.uniform(0,20,size=20)
arrivingtime.sort()
workingtime = np.random.uniform(1,3,size=20)startingtime = [0 for i in range(20)]
finishtime = [0 for i in range(20)]
waitingtime = [0 for i in range(20)]
emptytime = [0 for i in range(20)]startingtime[0] = arrivingtime[0]
finishtime[0] = startingtime[0]+workingtime[0]
waitingtime[0] = startingtime[0]-arrivingtime[0]print(startingtime[0],workingtime[0],finishtime[0],waitingtime[0])for i in range(1,len(arrivingtime)):if finishtime[i-1] > arrivingtime[i]:startingtime[i] = finishtime[i-1]else:startingtime[i] = arrivingtime[i]emptytime[i]=arrivingtime[i]-finishtime[i-1]finishtime[i] = startingtime[i] + workingtime[i]waitingtime[i] = startingtime[i] - arrivingtime[i]print('第%d个人,到达时间 开始时间 工作时间 结束时间 等待时间\n' %i,arrivingtime[i],startingtime[i],workingtime[i],finishtime[i],waitingtime[i], '\n')#添加个排队人数计算
#plt.plot(arrivingtime,'--k')
#plt.plot(startingtime,'*')
#plt.fill_between(arrivingtime,startingtime,0)
#plt.grid()end=[]
y=[]
for i in np.arange(0,40.0,0.5):num=0for j in range(0,20):if i>arrivingtime[j] and i<startingtime[j]:num=num+1dic={}dic['source']=idic['targe']=numend.append(dic)y.append(num)x = [i for i in np.arange(0,40.0,0.5)]
y = y
plt.plot(x, y)
plt.xlabel('x')
plt.ylabel('y')
plt.show()#等待时间作图
fig = plt.figure(figsize = (6,4))
plt.plot(waitingtime,'go')
plt.grid(True,linestyle='--',color='gray',linewidth='0.8')
plt.title('排队上厕所等待时间')
plt.show()

 

这篇关于蒙特卡罗模拟之排队上厕所问题的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1117254

相关文章

Kotlin Map映射转换问题小结

《KotlinMap映射转换问题小结》文章介绍了Kotlin集合转换的多种方法,包括map(一对一转换)、mapIndexed(带索引)、mapNotNull(过滤null)、mapKeys/map... 目录Kotlin 集合转换:map、mapIndexed、mapNotNull、mapKeys、map

nginx中端口无权限的问题解决

《nginx中端口无权限的问题解决》当Nginx日志报错bind()to80failed(13:Permissiondenied)时,这通常是由于权限不足导致Nginx无法绑定到80端口,下面就来... 目录一、问题原因分析二、解决方案1. 以 root 权限运行 Nginx(不推荐)2. 为 Nginx

解决1093 - You can‘t specify target table报错问题及原因分析

《解决1093-Youcan‘tspecifytargettable报错问题及原因分析》MySQL1093错误因UPDATE/DELETE语句的FROM子句直接引用目标表或嵌套子查询导致,... 目录报js错原因分析具体原因解决办法方法一:使用临时表方法二:使用JOIN方法三:使用EXISTS示例总结报错原

Windows环境下解决Matplotlib中文字体显示问题的详细教程

《Windows环境下解决Matplotlib中文字体显示问题的详细教程》本文详细介绍了在Windows下解决Matplotlib中文显示问题的方法,包括安装字体、更新缓存、配置文件设置及编码調整,并... 目录引言问题分析解决方案详解1. 检查系统已安装字体2. 手动添加中文字体(以SimHei为例)步骤

SpringSecurity整合redission序列化问题小结(最新整理)

《SpringSecurity整合redission序列化问题小结(最新整理)》文章详解SpringSecurity整合Redisson时的序列化问题,指出需排除官方Jackson依赖,通过自定义反序... 目录1. 前言2. Redission配置2.1 RedissonProperties2.2 Red

nginx 负载均衡配置及如何解决重复登录问题

《nginx负载均衡配置及如何解决重复登录问题》文章详解Nginx源码安装与Docker部署,介绍四层/七层代理区别及负载均衡策略,通过ip_hash解决重复登录问题,对nginx负载均衡配置及如何... 目录一:源码安装:1.配置编译参数2.编译3.编译安装 二,四层代理和七层代理区别1.二者混合使用举例

怎样通过分析GC日志来定位Java进程的内存问题

《怎样通过分析GC日志来定位Java进程的内存问题》:本文主要介绍怎样通过分析GC日志来定位Java进程的内存问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、GC 日志基础配置1. 启用详细 GC 日志2. 不同收集器的日志格式二、关键指标与分析维度1.

Java 线程安全与 volatile与单例模式问题及解决方案

《Java线程安全与volatile与单例模式问题及解决方案》文章主要讲解线程安全问题的五个成因(调度随机、变量修改、非原子操作、内存可见性、指令重排序)及解决方案,强调使用volatile关键字... 目录什么是线程安全线程安全问题的产生与解决方案线程的调度是随机的多个线程对同一个变量进行修改线程的修改操

Redis出现中文乱码的问题及解决

《Redis出现中文乱码的问题及解决》:本文主要介绍Redis出现中文乱码的问题及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1. 问题的产生2China编程. 问题的解决redihttp://www.chinasem.cns数据进制问题的解决中文乱码问题解决总结

全面解析MySQL索引长度限制问题与解决方案

《全面解析MySQL索引长度限制问题与解决方案》MySQL对索引长度设限是为了保持高效的数据检索性能,这个限制不是MySQL的缺陷,而是数据库设计中的权衡结果,下面我们就来看看如何解决这一问题吧... 目录引言:为什么会有索引键长度问题?一、问题根源深度解析mysql索引长度限制原理实际场景示例二、五大解决