如何用Python数值求解非线性方程？Sympy求函数交点遇报错求助

问题原因 #

你遇到的Cannot create mpf from x报错，本质是两个问题：

• nsolve作为数值求解器必须提供初始猜测值，否则无法启动迭代计算
• 代码中使用了浮点数（如0.5、0.1），SymPy在混合浮点数与符号变量时容易出现类型兼容问题

修正后的代码 #

import sympy
from sympy import S

# 限定x为正数，匹配函数定义域（x**-0.5要求x>0）
x = sympy.symbols('x', positive=True)

# 用SymPy有理数替代浮点数，避免符号计算精度问题
f_x = x**(-S(1)/2)
g_x = -sympy.log(S(1)/10 + x**(-S(1)/2), 10)

# 使用nsolve时传入方程、变量、初始猜测值（这里选5，对应函数交点的大致区间）
sol = sympy.nsolve(f_x - g_x, x, 5)
num = float(sol)
print(num)

关键改动说明 #

• 限定变量定义域：通过positive=True确保x在有效范围内，避免求解过程中出现无意义的计算
• 替换浮点数为有理数：用S(1)/2、S(1)/10替代0.5、0.1，让SymPy的符号计算更稳定，避免浮点数精度干扰
• 添加初始猜测值：nsolve需要一个迭代起点，这里选择5是因为通过函数趋势分析，交点在x=1到x=10之间，你也可以根据实际情况调整初始值

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