特定点的DifferentialEquations.jl输出

Question

我试着用朱莉娅的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分怎么办？

Answer 1

关于“控制解决方案”的ODE教程部分演示了为此目的使用saveat。示威活动是：

sol = solve(prob,reltol=1e-6,saveat=0.1)

这将节省0.0，0.1，.在你的例子中。在下面，人们注意到：

更普遍地说，稀释剂可以是任何时间点的集合来保存。

例如，我们可以使用save来保存t=30、60和78，如下所示：

sol = solve(prob,saveat=[30.0,60.0,78.0])

这些例子应该会让你走上正确的方向。有关更多细节，请参见文档的输出控件部分。