sympy 中求解结果的数值

Posted

技术标签:

【中文标题】sympy 中求解结果的数值【英文标题】:Numeric value of solve result in sympy 【发布时间】:2021-02-19 10:44:52 【问题描述】:

我试图计算solve 函数结果的数值。在 WolframAlpha 中,它被称为“近似形式”。

这是我的方程式:

>>> solve([Eq(u, r * i), Eq(p, u * i), Eq(p, 20), Eq(r, 4)])
[i: -√5, p: 20, r: 4, u: -4⋅√5, i: √5, p: 20, r: 4, u: 4⋅√5]

我试过N,但似乎不起作用:

>>> N(solve([Eq(u, r * i), Eq(p, u * i), Eq(p, 20), Eq(r, 4)]), 2)
Traceback (most recent call last):
  File "<string>", line 1, in <module>
  File "/base/data/home/apps/s~sympy-live-hrd/66.426491309333028408/sympy/sympy/core/evalf.py", line 1544, in N
    return sympify(x, rational=True).evalf(n, **options)
AttributeError: 'list' object has no attribute 'evalf'

执行此操作的正确语法是什么?

【问题讨论】:

【参考方案1】:

N 函数只能用于 sympy 表达式,不能用于列表或字典。您可以使用列表/字典推导来重建解决方案数据结构,如下所示:

In [22]: [k: N(e) for k, e in s.items() for s in sol]
Out[22]: [i: -2.23606797749979, p: 20.0, r: 4.0, u: -8.94427190999916, i: 2.23606797749979, p: 20.0, r: 4.0, u: 8.94427190999916]

【讨论】:

以上是关于sympy 中求解结果的数值的主要内容,如果未能解决你的问题,请参考以下文章

Sympy:在有限域中求解矩阵

求解Haversine公式的分量时,sympy挂起[重复]

如何传递用于scipy的sympy表达式?

如何用python进行符号求解?

求解可分微分方程

数值计算 --求解连续微分系统和混沌系统