如何解决SymPy:难以简化包含因子一的表达式
对不起,对于新手问题,但到目前为止我还没有找到答案。
假设我们有一个符号,然后将其乘以浮点数1:
import sympy
x = sympy.Symbol('x')
x = 1.0 * x
无论我做什么,所有包含此x
的表达式也都将包含因子1。
from pprint import pprint
pprint(x)
pprint(sympy.simplify(x))
输出:
1.0⋅x
1.0⋅x
在上面的示例中,我想输出x
而不是1.0⋅x
。是否可以用浮点数1求解所有乘法来简化表达式?我应该将此浮点数转换为其他数据类型吗?
解决方法
这似乎是由于浮点数与整数1的关系-如果您可以避免使用1,请不要使用1.0?
您可以检查它是否与if f.as_integer_ratio() == (1,1)
完全相同,如果不是,则不要在它们前面加上不必要的'1.0':
>>> xf = 1.0 * x
>>> xi = 1 * x
>>> xf
1.0*x
>>> xi
x
>>> pprint(xf)
1.0*x
>>> pprint(xi)
x
这可能是由于浮点数不精确(我不太了解sympy,不能说更多)。
编辑:针对您的问题,float.as_integer_ratio()
的运行时很小:
>>> python -m timeit -n 1000000 '(1.0).as_integer_ratio() == (1,1)'
1000000 loops,best of 5: 172 nsec per loop
>>> python -m timeit -n 1000000 '(1.1).as_integer_ratio() == (1,best of 5: 327 nsec per loop
>>> python -m timeit -n 1000000 '(1.01).as_integer_ratio() == (1,best of 5: 331 nsec per loop
>>> python -m timeit -n 1000000 '(0.999999999999).as_integer_ratio() == (1,best of 5: 329 nsec per loop
,
我找到了解决方案。有一个特殊的功能nsimplify()
可以完成工作。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。