如何解决评估Sympy表达式不会产生正确的结果
为什么在以下代码中fdeno
和deno2
之间的差异如此之大?
from sympy import symbols
sign2 = symbols('sign2')
deno = sign2*(sign2 - 1.6)*(sign2 - 1.5)*(sign2 - 1.4)*(sign2 - 1.3)*(sign2 - 1.2)*(sign2 - 1.1)*(sign2 - 1.0)*(sign2 - 0.9)*(sign2 - 0.8)*(sign2 - 0.7)*(sign2 - 0.6)*(sign2 - 0.5)*(sign2 - 0.4)*(sign2 - 0.3)*(sign2 - 0.2)
"""not include the term (sign2 -0.1)"""
fdeno = float(deno.evalf(subs={sign2:round(0.1,1)},n=15))
deno2 = 0.1*(0.1-1.6)*(0.1-1.5)*(0.1-1.4)*(0.1 - 1.3)*(0.1 - 1.2)*(0.1 - 1.1)*(0.1 - 1.0)*(0.1 - 0.9)*(0.1 - 0.8)*(0.1 - 0.7)*(0.1 - 0.6)*(0.1 - 0.5)*(0.1 - 0.4)*(0.1 - 0.3)*(0.1 - 0.2)
print(fdeno,"//",deno2)
# Output:
0.0 // -0.00013076743679999998
解决方法
似乎您的浮动精度已用尽。
要求解,请将sign2定义为带有30个小数位的0.1:
fdeno = float(deno.evalf(subs={sign2: Float('0.1',30)}))
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。