用Julia求解变参数常微分方程组

用Julia求解变参数常微分方程组,julia,ode,Julia,Ode,想象一下流行病学中的SIR微分方程系统: S'(t)=-b*S(t)I(t)/N I'(t) = b*S(t)I(t)/N - c*I(t) R'(t) = c*I(t) 用b和c常数。我想知道是否有可能用b作为向量而不是常数来求解这样一个系统: S'(t)=-b(t)*S(t)I(t)/N I'(t) = b*S(t)I(t)/N - c*I(t) R'(t) = c*I(t) 假设我们有一个外部函数来计算每个时间步的b值,因此我们最终得到一个向量b=[b(1),b(2),…,b(tspa

想象一下流行病学中的SIR微分方程系统:

S'(t)=-b*S(t)I(t)/N
I'(t) = b*S(t)I(t)/N - c*I(t)
R'(t) = c*I(t)
用b和c常数。我想知道是否有可能用b作为向量而不是常数来求解这样一个系统:

S'(t)=-b(t)*S(t)I(t)/N
I'(t) = b*S(t)I(t)/N - c*I(t)
R'(t) = c*I(t)

假设我们有一个外部函数来计算每个时间步的b值,因此我们最终得到一个向量
b=[b(1),b(2),…,b(tspan)]
,我们是否可以用这个向量作为ODE问题和解算器的输入?

这是对解算器如何工作的严重误解,请参阅中以前的讨论,对于该插值,您可能希望使用这是对解算器如何工作的严重误解,请参见中的前面讨论,对于该插值,您可能希望使用