如何在Python中求解非線性方程組？

Question

我有一個方程組：

for i [1, N]:

        |A_i x (X - B_i)|
y_i = ------------------------
           |A_i|

the goal: find X such that it minimizes the target function:
sum_{i in [1, N]} (y_i)^2 -> min

其中A_i, X, B_i是3x1向量， *是標量乘法， |v| 是v歐式范數，並且x是交叉乘法。

如何使用Python（scipy.optimize？）求解此方程組？ 我以前僅使用numpy.linalg.solve解決了Ax = b ，所以我有點困惑。

我認為我應該使用Nelder-Mead單純形算法 ，聽起來正確嗎？

Answer 1

基本上，我最終使用了此SciPy文檔中的代碼：

import numpy as np
from scipy.optimize import minimize

# Target function
def rosen(x):
    """The Rosenbrock function"""
    return sum(100.0*(x[1:]-x[:-1]**2.0)**2.0 + (1-x[:-1])**2.0)

x0 = np.array([1.3, 0.7, 0.8, 1.9, 1.2])
res = minimize(rosen, x0, method='nelder-mead',
               options={'xtol': 1e-8, 'disp': True})

print(res.x)