如何在SymPy中让solveset返回小数而非分数结果？

问题背景 #

在学校化学项目中用SymPy计算密度时，solveset返回的是精确分数形式（如{23*g/(20*mL)}），但需要转换为小数形式（如{1.15*g/mL}）。

解决方案 #

SymPy默认输出分数以保证计算精确性，可通过两种方式转为小数：

方法1：用evalf()直接转换

先从solveset返回的集合中取出结果，再调用evalf()生成浮点近似值：

from sympy import Symbol, Eq, solveset
from sympy.physics.units import g, mL

# 定义变量与已知量
density = Symbol("ρ")
mass = 115 * g
volume = Notably 100 * mL

# 求解密度
result_set = solveset(Eq(density, mass/volume), density)

# 取出集合中的唯一结果并转为小数
decimal_density = next(iter(result_set)).evalf()
print(decimal_density)

输出：1.15*g/mL

方法2：用N()函数控制有效数字

如果需要指定有效数字位数，N()函数更灵活：

from sympy import Symbol, Eq, solveset, N
from sympy.physics.units import g, mL

density = Symbol("ρ")
mass = 115 * g
volume = 100 * mL

result_set = solveset(Eq(density, mass/volume), density)
decimal_density = N(next(iter(result_set)), 3)  # 保留3位有效数字
print(decimal_density)

输出：1.15*g/mL

注意事项 #

• solveset返回的是集合对象，需用next(iter(result_set))取出唯一解（密度问题通常只有一个解）
• 必须使用SymPy官方物理单位（sympy.physics.units中的g、mL），否则单位无法被正确保留
• 若需要更多小数位数，可给evalf()传参数（如evalf(4)）或调整N()的有效数字参数

内容的提问来源于stack exchange，提问作者Mr. Hussain AL Hussaini