对于僵硬的系统，boost :: numeric :: odeint，如果我们无法获得分析的Jacobean矩阵怎么办？

Question

我正在使用boost :: numeric :: odeint rosenbrock4积分器。问题是我无法获得分析的Jacobean矩阵，有没有解决方案来解决这个问题？

Answer 1

您可以使用数值方案来计算近似雅可比。但这必须由您自己或其他一些图书馆实施。 odeint不提供此功能。

出于好奇：你的系统究竟是什么，你找不到雅可比人？如果rhs不可微分，则可能会遇到数字问题的其他问题。