我正在用simmer做一个DES模型。而不是指数或均匀分布的到达过程,我想要一个具有两个高峰小时的到达过程。现在我只有这两个函数:
add_generator("traffic" , trajectory = weekday , function() rnorm(100, 15, 5))%>%
add_generator("traffic2" , trajectory = weekday , function() rexp(100, 1/5))%>% 我想要的是10辆车/小时和50辆车在上午8-10点的持续输入
发布于 2021-04-14 06:43:27
下面的代码每小时在hstart和hstop之间生成cph汽车
library(simmer)
weekday <- trajectory()
cph <- 10
hstart <- 8
hstop <- 10
simmer() %>%
add_generator(paste0(cph,"carsperhour"), weekday , at(seq(from=hstart,to=hstop, by = 1/cph))) %>%
run() %>%
get_mon_arrivals()
#> name start_time end_time activity_time finished replication
#> 1 10carsperhour0 8.0 8.0 0 TRUE 1
#> 2 10carsperhour1 8.1 8.1 0 TRUE 1
#> 3 10carsperhour2 8.2 8.2 0 TRUE 1
#> 4 10carsperhour3 8.3 8.3 0 TRUE 1
#> 5 10carsperhour4 8.4 8.4 0 TRUE 1
#> 6 10carsperhour5 8.5 8.5 0 TRUE 1
#> 7 10carsperhour6 8.6 8.6 0 TRUE 1
#> 8 10carsperhour7 8.7 8.7 0 TRUE 1
#> 9 10carsperhour8 8.8 8.8 0 TRUE 1
#> 10 10carsperhour9 8.9 8.9 0 TRUE 1
#> 11 10carsperhour10 9.0 9.0 0 TRUE 1
#> 12 10carsperhour11 9.1 9.1 0 TRUE 1
#> 13 10carsperhour12 9.2 9.2 0 TRUE 1
#> 14 10carsperhour13 9.3 9.3 0 TRUE 1
#> 15 10carsperhour14 9.4 9.4 0 TRUE 1
#> 16 10carsperhour15 9.5 9.5 0 TRUE 1
#> 17 10carsperhour16 9.6 9.6 0 TRUE 1
#> 18 10carsperhour17 9.7 9.7 0 TRUE 1
#> 19 10carsperhour18 9.8 9.8 0 TRUE 1
#> 20 10carsperhour19 9.9 9.9 0 TRUE 1
#> 21 10carsperhour20 10.0 10.0 0 TRUE 1https://stackoverflow.com/questions/67078197
复制相似问题