如何优化python中的非线性方程？

Question

我已经写到这里了，但我无法继续。如何定义 x0、x1 和 x2？我该如何解决这个问题？

def nonlin(a,b,c):
    return [3*x0-cos(x1*x2)+ a, x0*x0-81*(x1+0.1)**2+ sin(x2)+b, exp(-x0*x1)+20*x2+c]

Answer 1

您可以尝试使用scipy.optimize.fsolve 解决此问题：

from scipy.optimize import fsolve
from numpy import cos, sin, exp

# Define the non-linear equations
def equations(X, a,b,c):
    x0, x1, x2 = X
    return [3*x0-cos(x1*x2)+ a, x0*x0-81*(x1+0.1)**2+ sin(x2)+b, exp(-x0*x1)+20*x2+c]

a = 1.0
b = 0.5
c = 2.5

# Using (0.0, 0.0, 0.0) as an initial vector
x0, x1, x2 = fsolve(equations, x0=(0.,0.,0.), args=(a,b,c)) 
print(x0, x1, x2)

输出：

-6.826133854836766e-06 -0.03657006638562982 -0.17499998751839313