python - 另外限制 scipy 的根求解器

Question

我正在尝试使用 scipy.optimize.root() 求解非线性方程组

其中一个方程包含类似于 x[5]**epsilon 的幂，其中 epsilon 是一个参数。对于 epsilon=1 我设法找到了根，但对于 1 周围的任何不同 epsilon，我却找不到。我在这个过程中的某个地方得到了 NaN，我认为它来自于求解器尝试负值并为此得到 NaN。

是否有任何方法可以进一步限制求解器，例如告诉他 x[5] >= 0 或类似的？我知道 scipy 的最小化工具可以使用这些边界，但我在 documentation 中找不到类似的东西。对于root。

Answer 1

在目标函数中将其替换为 np.exp(x[5])**epsilon，而不是 x[5]**epsilon。那么x[5]可以是-inf和+inf之间的任何值以及np.exp(x[5]) 永远是正数。

远离 abs(x[5]) 或可能导致导数不连续的类似其他转换。