使用 np trapz 集成数组

Question

我一直在使用 np.trapz 在数组上进行集成，直到现在还没有遇到任何问题。我得到了一个面积明显小于 1 的分布，因为它的最大值为 0.16，分布的宽度大约为 6，但似乎返回分布下方的面积大于 60。

这是我的代码：

import numpy as np
import matplotlib.pyplot as plt
data = np.load('dist.npy')
thetavals=np.linspace(0,2*np.pi,1000)
plt.xlabel(r'$\theta$')
plt.ylabel(r'$P(\theta)$')
plt.plot(thetavals,data[0:1000])
plt.show()

integralvalue=np.trapz(data)
print('The integral of this distribution results in: ',integralvalue)

Answer 1

使用 numpy trapz，在不选择 x 参数的情况下，假设我们的分布间距是均匀分布的，但是这些应该与最初形成分布的 theta 值相关，使用以下代码：

import numpy as np
import matplotlib.pyplot as plt
data = np.load('dist.npy')
thetavals=np.linspace(0,2*np.pi,1001)
plt.xlabel(r'$\theta$')
plt.ylabel(r'$P(\theta)$')
plt.plot(thetavals,data)
plt.show()

integralvalue=np.trapz(data,thetavals)
print('The integral of this distribution results in: ',integralvalue)

如预期的那样，得到一个小于 1 的数。