numpy数组元素的累积积分

Question

我想进行以下类型的集成：

假设我有 2 个数组

a = np.array[1,2,3,4]
b = np.array[2,4,6,8]

我知道如何使用以下方法集成这些：

c = scipy.integrate.simps(b, a)

c = 15 用于上述数据集。

我想做的是将每个数组的第一个元素相乘并添加到名为 d 的新数组中，即 a[0]*b[0] 然后将前 2 个元素整合到数组中，然后整合前 3 个元素，依此类推。所以最终这个数据集，我会得到

d = [2 3 8 15]

我尝试了一些东西，但没有运气；我对编写代码很陌生。

Answer 1

如果我正确理解了您的需求，您可以执行以下操作：

import numpy as np
from scipy import integrate

a = np.array([2,4,6,8])
b = np.array([1,2,3,4])
d = np.empty_like(b)
d[0] = a[0] * b[0]
for i in range(2, len(a) + 1):
    d[i-1] = integrate.simps(b[0:i], a[0:i])
print(d)