3

与此问题类似，我正在尝试使用与时间相关的输入参数来解决此 ODE 。它由一系列离散的回调组成。在某些时候，参数会改变（不是状态！）。时间和值存储在nx2 Array. 但是我无法获取affect在指定时间找到对应参数值的函数。在给定的示例中，分配给的值u[1]通常是恒定的。考虑这个 MWE（使用非常类似于 Matlab 的方法），它可以在没有回调的情况下正常工作：

using DifferentialEquations
using Plots

function odm2prod(dx, x, params, t)
    k_1, f_1, V_liq, X_in, Y_in, q_in = params

    rho_1 = k_1*x[1]
    q_prod = 0.52*f_1*x[1]
    # Differential Equations
    dx[1] = q_in/V_liq*(X_in - x[1]) - rho_1
    dx[2] = q_in/V_liq*(Y_in - x[2])
end

x0      = [3.15, 1.5]
tspan   = (0.0, 7.0)
params  = [0.22, 43, 155, 249, 58, 0]
prob    = ODEProblem(odm2prod, x0, tspan, params)

input   = [1.0 60; 1.1 0; 2.0 60; 2.3 0; 4.0 430; 4.05 0]
dosetimes = input[:,1]
function affect!(integrator)
    ind_t = findall(integrator.t == dosetimes)
    integrator.p[6] = input[ind_t, 2]
end
cb = PresetTimeCallback(dosetimes, affect!)
sol = solve(prob, Tsit5(), callback=cb, saveat=1/12)

plot(sol, vars=[1, 2])

这没用。错误出现在第 22 行，因为在 Julia 中似乎没有定义将向量与标量进行比较，或者存在我不知道的特殊语法。

我知道可以在 Julia 中使用与时间相关的参数，但我想这仅适用于连续函数，而不适用于离散变化！？我已经查看了 的帮助interpolate，但我不确定如何将它用于我的具体情况。

有人可以告诉我如何让它工作吗？应该可能只需要几行代码。另外，我不一定想dosetimes成为 的一部分sol.t，除非它们重合。

4

1 回答 1

3

你用findall错了，文档说

findall(f::Function, A)

返回whereI的索引或键的向量。Af(A[I])true

然后你必须考虑到搜索“all”的结果是一个列表。如您所料，它只有一个元素，请仅使用第一个元素

function affect!(integrator)
    ind_t = findall(t -> t==integrator.t, dosetimes)
    integrator.p[6] = input[ind_t[1], 2]
end

你得到了情节

ode 解的图

于 2020-01-25T15:01:25.713 回答