贝叶斯 ODE 与 Julia

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/