[英]R - deSolve,parameters
我是这个网站的新手。 我正在R中实现脚本,因此我需要利用deSolve,为了做到这一点,我必须定义在编写方程式的“变化率”时有用的参数。
这是我的问题,对于使用软件包我必须定义这些参数,但实际上我并不需要它们,我的意思是方程式的变化率我没有任何参数,所以我要避免定义的方法。
我部分解决了以这种方式定义参数的问题(不是很优雅):
parameters <- c(1)
但我不喜欢这种解决方案。
预先感谢您的帮助,如果我的信息不太清楚,请向我报告,对不起,但英语不是我的母语!
那是代码:
parameters <- c(1)
e=0.2056
ini <- c(q1 =1-e, q2 = 0, p1 = 0, p2 = sqrt ((1+e)/(1-e)) )
rhs <- function(t, ini,parameters)
{
with(as.list(c(ini,parameters)),{
# rate of change
dq1 <- p1
dq2 <- p2
dp1 <- -q1/((q1^2+q2^2)^(3/2))
dp2 <- -q2/((q1^2+q2^2)^(3/2))
# return rate of change
list(c(dq1,dq2,dp1,dp2))
} )
}
###### EULER ######
library(deSolve)
times <- seq(0,40, by = 0.0005)
out<- ode(y = ini, times = times, func = rhs, parms = parameters,method="euler")
head(out)
'ode'需要这些参数,如果我不放置参数,它会提示错误
parms
传递给实际执行求解的方法,然后传递给提供的函数,因此必须提供。 但是它不需要一个值,您的函数也不需要使用它。
# Accept and ignore third argument
rhs <- function(t, ini,...)
{
with(as.list(c(ini)),{
# rate of change
dq1 <- p1
dq2 <- p2
dp1 <- -q1/((q1^2+q2^2)^(3/2))
dp2 <- -q2/((q1^2+q2^2)^(3/2))
# return rate of change
list(c(dq1,dq2,dp1,dp2))
} )
}
# Pass NULL to parms
out2 <- ode(y = ini, times = times, func = rhs, parms = NULL,method="euler")
identical(out, out2)
## [1] TRUE
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.