使用SymPy解方程时遇TypeError: cannot create mpf from x错误求助

错误原因 #

触发TypeError: cannot create mpf from x的核心原因有两个：

1. 使用了math.e而非SymPy提供的符号化自然常数sympy.E。math.e是浮点数，和SymPy符号变量x混合运算后，会导致表达式无法被nsolve正确解析为符号方程。
2. SymPy的Sum默认仅支持整数步长的求和操作，代码中使用0.5作为求和起始值或上下限，不符合SymPy求和规则，会导致符号表达式解析异常。

修正方案 #

1. 替换自然常数 #

导入sympy.E，替换代码中所有的math.e。

2. 修正求和逻辑 #

原代码的求和索引逻辑混淆了整数索引与数值，应直接基于amount_of_days的元素构建求和项，而非通过非整数的i范围遍历。

修正后的代码 #

from sympy import symbols, Eq, nsolve, E

# 测试变量
cr_user = 0.25 
amount_of_days = [0.108219, 0.6082191781]
md_user = 0.6082191781
price_user = 98.37

x = symbols("x")
n = len(amount_of_days)

if n % 2 == 0:
    # 前半部分元素求和
    lhs_terms = [cr_user * E**(-x*(md_user - day)) for day in amount_of_days[:n//2]]
    lhs = sum(lhs_terms)
    # 后半部分元素求和
    rhs_terms = [-cr_user * E**(-x*(md_user - day)) - (cr_user+100)*E**(-x*md_user) + price_user for day in amount_of_days[n//2:]]
    rhs = sum(rhs_terms)
else:
    # 前半部分元素求和
    lhs_terms = [cr_user * E**(-x*(md_user - day)) for day in amount_of_days[:(n-1)//2]]
    lhs = sum(lhs_terms)
    # 后半部分元素求和
    rhs_terms = [-cr_user * E**(-x*(md_user - day)) - (cr_user+100)*E**(-x*md_user) + price_user for day in amount_of_days[(n-1)//2:]]
    rhs = sum(rhs_terms)

# 提供初始猜测值帮助求解器收敛
print(nsolve(Eq(lhs, rhs), x, 0.1))

关键说明 #

• 用sum()直接计算列表项的和，比Sum()更适合这种基于已有元素的求和场景，避免了索引步长的问题。
• nsolve需要提供初始猜测值（示例中用0.1），帮助数值求解器找到收敛方向。

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