如何解决与 Julia 的贝叶斯 ODE
我一直在尝试实现贝叶斯 ODE。在石油行业,我们使用以下等式来拟合生产数据然后进行预测:
ODE 方程描述为:
其中 0 我的初始代码: 不知道如何解决这个问题:using DiffEqFlux,OrdinaryDiffEq,Flux,Optim,Plots,AdvancedHMC
function Arps!(du,u,p,t)
y = u
K,n = p
du = (y * K * y^n)
end
tspan=(1.0,200.0)
tsteps = range(1,200,length = 200)
u0 = [5505.99]
p=[0.17,0.87]
prob1 = ODEProblem(Arps!,u0,tspan)
sol_ode = solve(prob1,Vern7(),saveat = tsteps)
MethodError: no method matching iterate(::DiffEqBase.NullParameters)
解决方法
您没有将任何参数传递到 ODE 中。 prob1 = ODEProblem(Arps!,u0,tspan,p)
。
关于贝叶斯部分,看教程:
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。