您可以使用root:
import numpy as np
from scipy.optimize import root
def your_funcs(X):
x, y, a, b = X
f = [np.cos(x) / x - 0.2 * a - 0.1,
np.cos(y) / y - 0.2 * b - 0.1,
a + b - 1,
1.03 * np.sinc(x) - 1.02 * np.sinc(y)]
return f
sol2 = root(your_funcs, [0.1, 0.1, 0.3, 0.1])
print(sol2.x)
将打印
[ 1.30301572 1.30987969 0.51530547 0.48469453]
您的函数必须以评估为0的方式进行定义,例如a b-1而不是a b = 1.
快速检查:
print(your_funcs(sol2.x))
给
[-1.9356960478944529e-11, 1.8931356482454476e-11, 0.0, -4.1039033282785908e-11]
因此,解决方案应该可以(请注意e-11基本上为0).
或者,您也可以使用fsolve:
from scipy.optimize import fsolve
sol3 = fsolve(your_funcs, [0.1, 0.1, 0.3, 0.1])
这会给您相同的结果:
[ 1.30301572 1.30987969 0.51530547 0.48469453]
您可以使用args参数传递其他参数&#