julia - Julia 的贝叶斯 ODE

Question

我一直在尝试实现贝叶斯 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) 

Answer 1

您没有将任何参数传递到您的 ODE。 prob1 = ODEProblem(Arps!,u0,tspan,p)。

关于贝叶斯部分，看教程:

https://turing.ml/dev/tutorials/10-bayesiandiffeq/