用Python解决非线性方程组（scipy.optimize.fsolve）

Question

我正在尝试解决以下简单的非线性方程组（Source（第二个示例））：

(I) y - x^2 = 7 - 5x
(II) 4y - 8x = -21

它应该只有一个解（x = 3.5，y = 1.75）。

我当前使用scipy堆栈的方法如下：

from scipy.optimize import fsolve

def equations(p):
    x, y = p
    return (y - x**2 -7 + 5*x, 4*y - 8*x + 21)

x, y =  fsolve(equations, (5, 5))

print(equations((x, y)))

并产生以下结果（不是结果）：

(0.0, 0.0)

我已经尝试过不同的初始估算，但是它没有提供正确的解决方案。

我的方法有什么问题？我想念什么吗？

提前致谢！

Answer 1

这工作得很好：

In [1]: %paste
from scipy.optimize import fsolve

def equations(p):
    x, y = p
    return (y - x**2 -7 + 5*x, 4*y - 8*x + 21)

x, y =  fsolve(equations, (5, 5))

print(equations((x, y)))

## -- End pasted text --
(0.0, 0.0)

In [2]: x
Out[2]: 3.5000000414181831

In [3]: y
Out[3]: 1.7500000828363667

equations(x, y)是(0, 0)表示两个 y - x**2 -7 + 5*x 和4*y - 8*x + 21均为0。

也许您感到困惑并打算print(x, y)代替print(equations(x, y))？