R语言常微分方程数值解海强作业

最新推荐文章于 2022-08-04 23:40:10 发布

原创

最新推荐文章于 2022-08-04 23:40:10 发布 · 7.3k 阅读

24 ·

CC 4.0 BY-SA版权

文章标签：

#R语言

本文介绍了如何使用R语言的deSolve包求解常微分方程，包括一阶微分方程的人口增长模型和一阶多元线性微分方程组的Lorenz模型。此外，还探讨了欧拉方法及其改进的Adams-Bashforth方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

调包解法

**
一阶微分方程：一元以人口增长为例
在这里插入图片描述
我们使用R package deSolve ODEs 函数。
A simplified form of the syntax for solving ODEs is: ode(y, times, func, parms, …)
where times holds the times at which output is wanted, y holds the initial conditions, func is the name of the R function that describes the differential equations,andparmscontainsthe parametervalues(oris NULL)
dy/dx.
Ode()的y相当于y(0) times相当与dx，function 相当与ry(1-y/k)

r <- 1
K <- 10
yini <- 2
derivs <- function(t, y, parms) list(r * y * (1-y/K))
library(deSolve)
times <- seq(from = 0, to = 20, by = 0.2)
out <- ode(y = yini, times = times, func = derivs, parms = NULL)
我们换个y（0）=12，输出out2.
r <- 1
K <- 10
yini <- 12
derivs <- function(t, y, parms) list(r * y * (1-y/K))
librar