使用 scipy 进行约束根查找

Question

我正在尝试使用 scipy.optimize.root_scalar 查找标量函数。对于我的问题，错误函数在域的某些区域定义不明确。不幸的是，我事先并不知道这些区域是什么（这是本次寻根的目的的一部分）。我可以提供一个括号区间，但一般来说，该括号内会有一些定义不明确的区域。

我理想的行为是在遇到这些区域时为错误函数返回 NaN，并让一些求根算法“学习”以避免它们。虽然这不适用于 root_scalar，所以当遇到这些区域时，我一直在返回一个正确符号的固定大错误。任何关于不那么棘手的解决方案的建议将不胜感激。

Answer 1

我会继续发布我实施的解决方案。对于任何数值计算导数的求根方法，我都无法解决上述问题。但是，我意识到一个更简单的解决方案是使用二分法，然后在错误定义不明确的区域中酌情返回 +/- 无穷大的错误。这很慢，但到目前为止似乎每次都有效，而且速度在我的应用程序中并不重要。