我对具有不同初始条件的 ODE 的多个解进行平均，因此所有解同时具有值很重要；例如，增量为 0.01。
我一直在使用数值配方 3 (nr3) 中的 ODE 例程。他们进行自适应尺寸步长并使用计算值进行相同的插值顺序。我不能使用它们，因为它们与 boost 冲突。还有其他类似的套路吗？
我查看了 GSL，它非常好，但它没有内置插值。我可以做到的一种方法是使用自适应大小求解 ODE，而不是运行 Akima 插值。但似乎 nr3 解决方案会更快、更准确。
问问题
417 次
2 回答
1
我使用了来自http://www.unige.ch/~hairer/software.html的 DOPRI5，输出密集 = 插值。我发现它可靠。我使用了原始版本（在 Fortran 中）；在同一个网页上还有一个 C 版本，我自己没有使用过，但我似乎记得人们对此很满意。
于 2011-10-05T08:39:11.817 回答
1
您可以使用odeint。它具有 Dopri5、Rosenbrock4 和 Burlish-Stoer 用于密集输出。
于 2011-10-31T21:27:47.527 回答