Question

我想弄清楚为什么sympy.diff无法按预期区分sympy多项式。通常，如果定义了符号变量并且未使用sympy.diff定义多项式，sympy.Poly就可以正常工作。但是，如果使用sympy.Poly定义函数，则sympy.diff似乎不会实际计算导数。下面的代码示例显示了我的意思：

import sympy as sy

# define symbolic variables
x = sy.Symbol('x')
y = sy.Symbol('y')

# define function WITHOUT using sy.Poly
f1 = x + 1
# define function WITH using sy.Poly
f2 = sy.Poly(x + 1, x, domain='QQ')

# compute derivatives and return results
df1 = sy.diff(f1,x)
df2 = sy.diff(f2,x)
print('f1:  ',f1)
print('f2:  ',f2)
print('df1:  ',df1)
print('df2:  ',df2)

这将打印以下结果：

f1:   x + 1
f2:   Poly(x + 1, x, domain='QQ')
df1:   1
df2:   Derivative(Poly(x + 1, x, domain='QQ'), x)

为什么sympy.diff不知道如何区分多项式的sympy.Poly版本？有没有办法区分sympy多项式，或者将sympy多项式转换为允许区分它的形式？

注意：我尝试使用不同的域名（即domain='RR'代替domain='QQ'），输出不会改变。

Answer 1

这似乎是一个错误。您可以直接在diff实例上调用Poly来绕过它。理想情况下，从顶级sympy模块调用函数diff应该产生与调用方法diff相同的结果。

In [1]: from sympy import *

In [2]: from sympy.abc import x

In [3]: p = Poly(x+1, x, domain='QQ')

In [4]: p.diff(x)
Out[4]: Poly(1, x, domain='QQ')

In [5]: diff(p, x)
Out[5]: Derivative(Poly(x + 1, x, domain='QQ'), x)

In [6]: diff(p, x).doit()
Out[6]: Derivative(Poly(x + 1, x, domain='ZZ'), x)

为什么sympy.diff没有像预期的那样区分交感多项式？

1 个答案: