python求解具有解> 0的联立方程

Question

我有联立方程（有很多方程），这是没有解的，我需要找到一个最佳解。解的每个元素必须>0（不为0或负数）。

np.linalg.lstsq可以得到最好的解决方案，但也可能出现负面因素。

scipy.optimize.nnls可以在没有负数的情况下得到最佳解，但也可能出现0。

如何限制每个元素的最佳解决方案>0 或者是否有任何其他功能可以用来解决问题？

Answer 1

对于每个参数p，将p_ = p - eps 定义为非常小的eps。

求解{p_}，然后赋值p = p_ + eps，你会得到所有大于零的参数。

在实践中，只需将 epsilon 添加到解决方案中即可。

如果这还不够，请解释原因。