SymPy调用solve方法求解自定义超越函数实根时无法得到结果的原因

• sp.solve()的适用场景不匹配。你构造的函数是带tanh嵌套平方根的超越隐式方程，这类方程本身不存在可以用初等函数写出的解析闭式解，通用符号求解器无法输出有效结果，运行时会因为遍历海量求解分支出现卡顿、返回空列表的问题。
• 符号变量缺少物理场景的定义域约束。代码中定义的符号变量x没有提前限定为正实数，求解器会默认遍历所有复数域、负实数域的无意义解分支，进一步冗余计算流程，甚至触发求解逻辑异常。
• 求解策略和问题类型错配。提前硬编码所有物理常数数值的写法，没有给求解器提供匹配物理意义的变量边界提示，进一步拉低了求解器的运行效率。

因为目标方程无闭式解析解，需要将求解策略从“符号解析求解”调整为“数值迭代求解”，具体调整点：

• 给符号变量增加明确的物理属性约束，将其定义为正实数，从根源上排除无意义的求解分支
• 放弃通用解析求解接口sp.solve()，替换为SymPy内置的专门针对非线性方程的数值求解接口sp.nsolve()，匹配这类超越方程的求解场景
• 保留原本设计的初始猜测值输入逻辑，作为数值迭代的起始点，提升求解收敛速度和稳定性

修正后可直接运行的代码 #

import sympy as sp
import matplotlib.pyplot as plt

# 物理参数配置
Cd = 0.25
g = 9.81
pf = 10**(-6)  # 原有冗余变量不影响运行，无需删除
t = 4
v = 36
xr = [int(input('initial guess : '))]
i = 0
Ea = 1
Es = 0.01

# 定义带正实数约束的质量符号变量，匹配物理场景
x = sp.Symbol('x', positive=True, real=True)

# 目标求解函数
def f(m):
    return sp.sqrt(g * m / Cd) * sp.tanh(sp.sqrt(g * Cd / m) * t) - v

# 调用数值求解器计算实根
ans = sp.nsolve(f(x), xr[0])
print(f"目标函数的实根计算结果：{ans.evalf():.4f}")

代码运行后，结合给定的物理参数，最终求得的实根约为74.21，符合带空气阻力的自由下落终端速度模型的参数计算逻辑。

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