sympy.lcm未正确扩展术语

Question

我试图用交感找到两个项的LCM。我使用的代码是：

import sympy
s=sympy.Symbol("s")
A=(s+0.0004)*(s+0.02)
B=(s+0.02)
sympy.lcm(A,B) 

回传

1.0*s**2 + 0.0204*s + 7.99999999999992e-6

然而，通过检查lcm应该是

1.0*s**2 + 0.0204*s + 8.0e-6

我如何才能得到sympy.lcm来计算精确的解？

谢谢。

Answer 1

这里的问题是，0.02和0.004没有一个有限的精确表示作为浮动，所以你有问题的四舍五入甚至在你开始之前。您可以使用sympy.Rationals保持这些数字的精确性：

A = (s + sympy.Rational(4, 1000))*(s + sympy.Rational(2, 100))
B = (s + sympy.Rational(2, 100))
sympy.lcm(A, B).expand()

这现在给出了一个确切的结果：

s**2 + 3*s/125 + 1/12500