python - 如何使用 SimPy 获得随机时间设置最小阈值?
问题描述
我正在尝试从 SimPy 的模拟程序中获取车辆的随机生成器。我使用的代码是从http://phillipmfeldman.org/Python/discrete_event_simulation/traffic_sim.py中提取并改编的。问题是我希望车辆以不少于特定的车距(例如 4 秒)到达。下面的代码显示车辆随机到达,但有时车头时距太小(例如,车辆#2 到达 2.03,车辆#3 到达 2.45,此差异小于 1 秒)。我想知道是否有办法为模拟设置特定阈值,谢谢...
from collections import deque # double-ended queue
from numpy import random
import simpy
from simpy.util import start_delayed
class Struct(object):
def __init__(self, **kwargs):
self.__dict__.update(kwargs)
random.seed([1, 2, 3])
# Total number of seconds to be simulated:
end_time= 3600.0
# Cars cars arrive at the traffic light according to a Poisson process with an
# average rate of (0.00028-0.5) per second:
demand_per_hour = 1800 #veh/h
sh = 4 # sat headway 4 seconds,
arrival_rate= demand_per_hour/(3600*sh)
t_interarrival_mean= 1.0 / arrival_rate
queue= deque()
arrival_count = departure_count= 0
def arrival():
global arrival_count, env, queue
while True:
arrival_count+= 1
print("Vehicle #%d arrived at time "
"%.3f." % (arrival_count, env.now))
# Schedule next arrival:
yield env.timeout( random.exponential(t_interarrival_mean))
print("\nSimulation of Cars Arriving at Intersection Controlled by a Traffic "
"Light\n\n")
env= simpy.Environment()
t_first_arrival= random.exponential(t_interarrival_mean)
start_delayed(env, arrival(), delay=t_first_arrival)
env.run(until=end_time)
解决方案
有几种方法可以做到这一点。
- 截断您的指数到达间隔时间分布:
t=max(threshold_value, random.exponential(t_interarrival_mean))
- 改变指数 IAT 分布(也称为 2 参数指数分布):
t = threshold_value + random.exponential(t_interarrival_mean)
- 使用不同的分布。你确定指数是最好的吗?三角分布怎么样?
推荐阅读
- android - IONIC 4 上的命令失败,退出代码为 enoent
- jira - QSH 在 Atlassian 中声明 Zephyr API 不匹配
- jquery - 添加新行时如何自动增加输入ID?
- mongodb - 如何在 Scala 的 MongoDB 聚合查询中使用“$let”?
- spring-boot - 在 Spring Boot 属性文件中外部化密码
- mysql - MySQL:根据if条件更新行
- reactjs - 将 ErrorBoundary 应用于没有根组件的反应项目
- powershell - 电源外壳 For Loop 不循环
- python - 为什么 Sublime Text 3 给我一个“python:找不到命令”错误
- scala - 通过将 Unit 作为类型传递给泛型特征的可选方法参数