python - DiracDelta 没有给出正确的结果

Question

我必须在一个复杂的积分中使用狄拉克三角洲，并希望看到它如何处理一个简单的情况，但它返回了错误的答案。任何线索我在下面做错了什么？

from sympy import DiracDelta
from scipy import integrate

def f(x):
     return x*DiracDelta(x-1)

b, err = integrate.quad(f, 0, 5)    
print b

这会返回 0.0 而它不应该。

Answer 1

似乎 sympy 函数与 scipy 集成不兼容。需要使用 sympy 集成。以下给出正确答案

from sympy import *
x = Symbol('x')
print integrate(x*DiracDelta(x-1), (x, 0, 5.0))    

不过，我不确定 sympy.integrate 是否像 scipy.integrate 一样强大和通用。