问在deSolve中更改时间步长的参数值

EN

这是对您的函数的一个非常简单的修改。如果你有兴趣知道你做错了什么，你可以看看下面。

mod <- function(times, yini, param) { 

  dt = times[2] - times[1]
  with(as.list(c(yini, param)), { 

    gamma0 <- ifelse(times <= 10*dt, 5, 0) 

    ## print(gamma0) 

    dalpha0 <- - a*alpha0 + gamma0 
    dbeta0 <- a*alpha0 - b*beta0 
    return(list(c(dalpha0, dbeta0))) 

  })} 

编辑

与G5W的答案一样，问题在于您将times与什么进行比较。

当你在写作时

times %in% seq(0,10,1)

你指的不是时间步骤。您只需引用times的值。

因此，如果您想在最初的10个时间步骤中使用它，只需使用我的代码或任何考虑dt的内容即可。

但有个问题要问你：

如果您不需要gamma0根据times进行更改，并且希望在最初的11 (10)个时间步骤中为5，那么为什么要将它与times进行比较呢？为什么不简单地将这些时间步骤设置为5呢？

我得到的结果和你略有不同。如果我取消注释print(gamma0)，它会打印5次，然后打印出513个零。虽然你可能比我在这里提供的更多，但用一种表面的方式来寻找原因并不难。

如果您有(注释掉的)语句print(gamma0)，则放置行：

cat("g:", gamma0, "  t:", times, "\n")

并运行代码。您将看到它显示的前两次是0。因为这些都在您的列表中，所以seq(0,10,1) gamma0是5。在那之后，时间值显示了变化。注意，打印出来的没有一个是从原来的时间列表( seq(from = 0, to = 10, by = 1/24) )中打印出来的，也没有一个是整数，所以没有一个符合将gamma0设置为5的条件。ode正在用时间(内插吗？)但这并不是简单地使用您提供的值。实际上，它没有打印出241个gamma0值和时间。它打印出515个这样的值。我注意到结果out确实有241个值。

我认为，从您的问题中，您认为ode实际上将在您的times上评估该函数。事实并非如此。它将times视为一个连续变量。但你的情况

gamma0 <- ifelse(times %in% seq(0,10,1), 5, 0) 

仅对11个特定值进行测试，而不是对值的范围进行测试。连续变量不太可能精确地命中这些值。

M-的答案对我没用，如果你只是尝试一下呢？

mod <- function(times, yini, param) {
  with(as.list(c(yini, param)), {

    if (times < 10) {
      gamma0 = 5
    } else if (times >= 10) {
      gamma0 = 0
    }

    dalpha0 <- -a * alpha0 + gamma0
    dbeta0 <- a * alpha0 - b * beta0
    return(list(c(dalpha0, dbeta0)))

  })
}

这是作为works类型函数工作的。