朱莉娅和常微分方程组

时间:2018-07-28 07:50:09

标签: julia ode differential-equations

我想尝试求解一个常微分方程组,也许是并行化的并且遇到了Julia和DifferentialEquations.jl。系统看起来像

x'(t) = f(t)*z(t)
y'(t) = g(t)*z(t)
z'(t) = f(t)*(1-2*x(t))/(2) -g(t)*y(t)

超过10 ^ 2

x(10^14) == 0
y(10^14) == 0
z(10^14) == 0

有人可以向我解释如何在茱莉亚设置这个问题吗?我检查了文档,仅能找到u0作为参数,但是它没有提供有关选择右手边界条件集的详细信息,非常感谢!

1 个答案:

答案 0 :(得分:1)

您正在寻找解决边界值问题(BVP)的方法。尽管此领域目前不如DifferentialEquations.jl的其他领域发达,但the tutorial on solving BVPs中显示了一些针对此领域的方法。 SPECIES方法可以尝试。

但是我会注意到您的时间尺度很大,可能会导致数字错误。可能需要使用更高精度的数字(BigFloat,ArbFloat,DoubleFloat)来处理该范围,或者您可能希望重新调整方程式中的时间,以便使其更适合于标准的双精度浮点数(MIRK4