使用nsimplify而不是simplify可以修复此问题。

我之所以使用SymPy，是因为我需要精确的值
症状是针对symbolic mathematics的。如果你真的需要"精确值"到可用的精度，也许mpmath会更合适。
也就是说，Symy有工具来测试黄金分割率的近似值，在本例中，它是Expr. equals（）函数：

from sympy import simplify
from sympy.geometry import RegularPolygon, Point
from sympy import GoldenRatio

sides = 5
polygon = RegularPolygon(Point(0, 0), 1, sides)
diagonal = polygon.vertices[0].distance(polygon.vertices[2])
print("Ratio:", simplify(diagonal/polygon.length))
print(simplify(diagonal/polygon.length).equals(GoldenRatio))

Ratio: sqrt(sqrt(5) + 5)/sqrt(5 - sqrt(5))
True

再说一次，这是符号数学，没有人有计算机能力，如果有的话，Python也没有精确比较无理数的能力。
例如：

>>> 1/3
0.3333333333333333
>>> 1/3 * 3
1.0
>>> 0.3333333333333333 * 3
1.0

CPython很难定义0.9~ == 1。
所以这些比较是行不通的：

>>> from sympy import GoldenRatio, simplify
>>> simplify(diagonal/polygon.length) == GoldenRatio
False
>>> import math
>>> GoldenRatio.equals(math.sqrt(math.sqrt(5)+5)/math.sqrt(5-math.sqrt(5)))
False

另一个例子：

>>> from sympy import GoldenRatio, simplify
>>> simplify(diagonal/polygon.length) - GoldenRatio
-GoldenRatio + sqrt(sqrt(5) + 5)/sqrt(5 - sqrt(5))
>>> (simplify(diagonal/polygon.length) - GoldenRatio) == 0
False
>>> (simplify(diagonal/polygon.length) - GoldenRatio).equals(0)
True

我很抱歉，我无法避免这一点：simplify, man