julia - DifferentialEquations.jl 在特定点的输出
问题描述
我正在尝试使用来自 Julia 的 DifferentialEquations.jl。我设法让它工作,但我想知道如何在特定时间点生成输出。文档对此并不清楚,我还没有找到一个这样做的例子。我目前使用的代码来自教程:
using DifferentialEquations
using Plots
function lorenz(du,u,p,t)
du[1] = 10.0*(u[2]-u[1])
du[2] = u[1]*(28.0-u[3]) - u[2]
du[3] = u[1]*u[2] - (8/3)*u[3]
end
u0 = [1.0;0.0;0.0]
tspan = (0.0,100.0)
prob = ODEProblem(lorenz,u0,tspan)
sol = solve(prob)
plot(sol,vars=(1,2,3))
目前它产生 1287 点,我不知道它是如何决定的。我的问题是如果我想在跨度 0 到 100 之间生成 20 个点怎么办?
解决方案
“控制求解器”的 ODE 教程部分演示saveat
了用于此目的的方法。演示是:
sol = solve(prob,reltol=1e-6,saveat=0.1)
在您的示例中,这将保存为 0.0、0.1、...。正下方指出:
更一般地说,saveat 可以是要保存的任何时间点的集合。
因此,例如,我们可以使用 save at 仅在 t=30、60 和 78 时保存,如下所示:
sol = solve(prob,saveat=[30.0,60.0,78.0])
这些例子应该让你朝着正确的方向前进。有关更多详细信息,请参阅文档的输出控件部分。
推荐阅读
- r - 在圆环图中排序图例不起作用ggplot2 R
- python - 有没有类似于 Python 中的“查找”命令的东西?
- sql - PowerQuery 条件最大值
- javascript - javascript 未在浏览器中执行
- oracle - 如何在 oracle 数据库的所有字段/表中查找数据模式匹配项?
- python - 在 Python 中使用 scour 或任何其他优化器优化 svg
- c# - log4net - 记录到数据库表
- python - Django:注释每个对象的重复值数量的计数
- java - ArrayList 的转换
> 列出 - >
- python - 安装 Anaconda 时出错