5
例如，当处理大型雅可比矩阵时，如何处理数字代码中的复杂性？
我正在求解一个耦合方程组的非线性系统，并且已经计算出离散系统的雅可比行列式。结果是非常复杂，下面是（仅！）的第一一个的3列3 × 93×93\times 9矩阵， （之所以会出现复杂性，部分原因是数值方案需要指数拟合才能保持稳定性。） 关于使用Jacobian实现数字代码，我有一个普遍的问题。 我可以继续在代码中实现此矩阵。但是我的直觉告诉我，由于纯粹的复杂性和不可避免的引入错误，期望几天（也许几周！）乏味的调试。似乎不可避免地如何应对这样的复杂性？您是否使用从符号包中自动生成代码（然后手动调整代码）？ 首先，我计划用有限差分近似调试解析雅可比行列式，我是否应该注意任何陷阱？您如何处理代码中的类似问题？ 更新资料 我正在用Python进行编码，并使用sympy生成了Jacobian。也许我可以使用代码生成功能？