如前所述,“numpy”中的“cumsum”效果很好。确保您的数据是 正确的PDF(即总和为1),否则CDF将不会以统一结束应该. 在这里 是一个最小的工作示例:
import numpy as np
from pylab import *
# Create some test data
dx = 0.01
X = np.arange(-2, 2, dx)
Y = exp(-X ** 2)
# Normalize the data to a proper PDF
Y /= (dx * Y).sum()
# Compute the CDF
CY = np.cumsum(Y * dx)
# Plot both
plot(X, Y)
plot(X, CY, 'r--')
show()