为什么SymPy无法计算Dirac delta（狄拉克δ函数）的极限？

Dirac δ 是广义函数（也叫分布），不属于SymPy的limit函数默认支持的普通实值函数范畴。SymPy的limit实现依赖于级数展开、极限路径判定等针对普通函数设计的逻辑，遇到δ函数时会进入无意义的符号化简循环，因此会卡住无返回。
你预期的x趋近于非0点时δ(x)极限为0是完全正确的，但目前SymPy没有针对分布类函数的极限做特殊适配，不会自动识别δ的支集仅在x=0的特性直接返回结果。另外δ函数在x→0时本身不存在普通实值极限，只有分布意义下的运算结果，这也是SymPy没有为其适配普通limit逻辑的原因之一。

• 方法1：优先判断极限点是否在δ函数的支集范围内。Diracδ的非零点取值恒为0，当极限点不等于0时，不需要调用limit函数，直接返回0即可。
• 方法2：使用subs方法替代limit查询非支集点的取值，DiracDelta(x).subs(x, -1)会直接返回0，不会出现卡顿问题。
• 方法3：如果需要处理通用极限场景，可以先调用DiracDelta的eval方法预判断：

from sympy import DiracDelta, Symbol
x = Symbol('x', real=True)
a = -1
# 预判断非支集点直接返回0
res = DiracDelta(x).eval(x, a)
if res is not None:
    print(res)
else:
    # 仅当极限点为0时再处理分布意义下的极限逻辑
    print("极限点为δ支集点，按需处理对应运算")

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