python-3.x - Scipy quad 与 cdf 作为差分集成

Question

我现在很困惑我写的一个小函数。我想将以下函数嵌入到我的程序中：

其中是泊松分布的累积分布函数的卷积。

如果我做对了，我会像这样计算卷积，其中 k 是间隔中的事件数，Mu 只是平均值，例如 5，j 是卷积数：

from scipy.stats import poisson
poisson.cdf(k, (Mu * (1 + i)))

函数 D 如下所示：

这是在 python 中完成的，如下所示：

def D(k, y, i):
    y = y - k
    return (2 * max(0, y) - 10 min(0, y)) * poisson.cdf(k, (Mu * (1 + i)))

现在我遇到的棘手部分是，如何让积分与 cdf 一起工作。目前我正在尝试这样做，但我不知道如何继续：

from scipy.integrate import quad
def Integrate(i, y):
    return quad(D, 0, np.inf, args=(y, i))[0]

如果有人知道如何做到这一点，我将不胜感激:)

编辑： 为了检查 cdf 的卷积是否正确，我这样做了：

print(poisson.pmf(0,(Mu * (1 + i))) + poisson.pmf(1,(Mu * (1 + i))))
>>> 0.445679641365
print(poisson.cdf(1,(Mu * (1 + i))))
>>> 0.445679641365