我正在用simmer
制作DES模型。我想要一个有两个高峰时间的到达过程,而不是指数或均匀分布的到达过程。现在我只有这两个功能:
add_generator("traffic" , trajectory = weekday , function() rnorm(100, 15, 5))%>%
add_generator("traffic2" , trajectory = weekday , function() rexp(100, 1/5))%>%
我想要的是10辆车/小时和上午8-10点50辆车的恒定输入。
以下代码在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 1