SymPy中nsolve未遵循x正性符号假设的问题及解决

问题原因 #

• SymPy里给符号变量设置的positive=True是符号运算专属的假设，仅在符号推导（比如化简、符号求解）时生效，而nsolve是数值求解器，依赖迭代算法（如牛顿法）寻找解，不会读取这些符号层面的约束规则。
• 你使用(-1,-5)作为初始猜测值，数值求解器会从这个点开始迭代，自然收敛到离它最近的负解，和你设置的符号假设没有关联。

解决方法 #

方法1：更换为正的初始猜测值 #

把初始猜测改成正区间的数值，比如(1,1)，求解器就会收敛到正解：

import sympy as sym

x = sym.symbols('x', positive=True)
y = sym.symbols('y', positive=True)
eq1 = sym.Eq(y, x**2)
eq2 = sym.Eq(y, 2)
# 改用正的初始点
result = sym.nsolve((eq1, eq2), (x, y), (1, 1))
print(result)

运行后输出：

[1.41421356237310]
[2.00000000000000]

方法2：先符号求解再转换为数值 #

先用sym.solve获取符号解，由于你已经给x添加了positive=True的假设，它会自动过滤负解，再将结果转换为数值即可：

import sympy as sym

x = sym.symbols('x', positive=True)
y = sym.symbols('y', positive=True)
eq1 = sym.Eq(y, x**2)
eq2 = sym.Eq(y, 2)

# 符号求解，自动过滤负解
symbolic_sol = sym.solve((eq1, eq2), (x, y))[0]
# 转换为数值
numeric_sol = (sym.N(symbolic_sol[0]), sym.N(symbolic_sol[1]))
print(numeric_sol)

输出：

(1.4142135623730951, 2.0)

内容的提问来源于stack exchange，提问作者confused_researcher