硕士论文中SymPy计算太阳能回收期返回ConditionSet，求数值解问题

你的方程 (x = \frac{3}{(1.02)^x}) 属于超越方程，无法通过初等函数得到解析形式的解，这就是solveset返回ConditionSet的原因——它找不到闭式解，但从函数单调性和图像可以判断，实数域内有且仅有一个解。

要得到具体的数值解，改用SymPy的nsolve函数即可，它专门用于求解方程的数值根。修改后的代码如下：

from sympy import symbols, Eq, nsolve

x = symbols("x")
# 构造方程
eq = Eq(x, 3 / (1.02)**x)
# 指定初始猜测值（根据函数趋势选合理值即可，比如x=1）
solution = nsolve(eq, x, 1)
print(f"投资回收期TVT的数值解为：{solution.evalf()}")

运行后会输出近似数值解（约为1.46），代入原方程可验证精度：

• 左边：(1.46)
• 右边：(3/(1.02^{1.46}) \approx 1.46)，误差极小。

如果需要更高精度，可调整evalf的参数，比如solution.evalf(10)会输出10位有效数字的结果。

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