SymPy中同次多项式LT与LC运算结果差异的技术问询

可复现性确认 #

该现象可以稳定复现。在SymPy环境中运行提供的代码，会得到完全一致的浮点数差异结果：

import sympy
x = sympy.var('X')
A = sympy.Poly(1 / 9 * x, x, domain='RR')
B = sympy.Poly(1 / 7 * x, x, domain='RR')
print(sympy.LT(A) / sympy.LT(B) == sympy.LC(A) / sympy.LC(B))  # 输出False
print((sympy.LT(A) / sympy.LT(B)).evalf(60))
print((sympy.LC(A) / sympy.LC(B)).evalf(60))

问题本质：浮点数精度特性与内部实现差异 #

这并非SymPy的bug，而是双精度浮点数运算的固有特性结合SymPy内部函数实现差异导致的：

• 首先，RR域对应双精度浮点数域，1/9和1/7无法用二进制浮点数精确表示，存储时会带有微小的舍入误差。
• LC(A)直接提取多项式的首项系数，即1/9的双精度近似值；LC(B)同理是1/7的双精度近似值，二者相除是两个近似系数直接做除法运算。
• LT(A)返回的是完整的首项多项式（(1/9)*X），LT(A)/LT(B)的计算逻辑是先处理多项式的除法：由于两个多项式次数相同，结果为系数的比值，但SymPy在处理多项式除法时的内部计算路径，和直接调用LC取系数再除法的路径存在细微差别，导致最终浮点数的舍入方向不同——一个是7/9的向下取整近似，一个是就近取整近似，最终表现为末尾的微小数值差异。

这种差异是浮点数运算的正常现象：数学上等价的表达式，在浮点数近似计算中，不同的计算路径会引入不同的舍入误差，导致结果出现极小的偏差。

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