出生死亡过程未完成的代码

时间:2016-07-12 18:57:47

标签: r queue simulation markov

我有一个出生&死亡过程代码问题。我有4个状态S = {0,1,2,3}

在州0,没有客户。 在州1,有1名顾客正在接受治疗。 在状态2中,有1个客户在队列中被处理+1。 在状态3中,有1个客户在队列中被处理+2。

如果有更多的客户来自状态3,他们就会离开并稍后回来。

Lamba是每小时到达的客户数量。 穆是每小时接受治疗的客户的数量。

最初我们从州0开始。

目前我只想了解下的代码。在我用Caps Lock编写的代码中,我需要帮助。

如果您需要更多信息来理解这个问题,请写下评论。

bd_process <- function(lambda, mu, initial_state = 0, steps = 100) {
time_now <- 0
state_now <- initial_state
time <- 0
state <- initial_state

for (i in 1:steps) {


if (state_now == 3) {
lambda_now <- 0
}else {
  lambda_now <- lambda
}


if (state_now == 0) {
mu_now <- 0
}else {
mu_now <- mu
}


#? WHAT DOES TIME_TO_TRANSISTION DO? WHAT SHOULD BE INSTEAD OF "..."?
time_to_transition <- ...


 #? WHAT DOES THIS CODE DO? WHAT SHOULD BE INSTEAD OF "..."?
if (...) {
  state_now <- state_now - 1
} else {
  state_now <- state_now + 1
}


time_now <- time_now + time_to_transition #WHAT IS TIME_NOW?
time <- c(time, time_now) #WHAT DOES THIS VECTOR CONSIST OF?
state <- c(state, state_now) #WHAT DOE THIS VECTOR CONSIST OF?
}
list(tid = tid, state = state) }

0 个答案:

没有答案