python - 将时间增量数组转换为瞬时数组

Question

如果我有一个时间增量数组，例如:

intervals = np.random.normal(loc=1,scale=0.1,size=100)

创建相应的时刻数组的一种可能方法是创建一个列表并手动求和:

Sum=0.
instants=[]
for k in range(len(intervals)):
    Sum+=intervals[k]
    instants.append(Sum)

instants=np.array(instants)

所以，我刚刚从 dt(i) 数组切换到 t(i) 数组。

但通常 python 提供了使用 for 循环的优雅替代方案。有更好的方法吗？

Answer 1

你这里描述的是累计和。您可以使用 .cumsum() [numpy-doc] 来计算此值:

intervals<b>.cumsum()</b>

例如:

>>> intervals
array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])
>>> intervals.cumsum()
array([ 0,  1,  3,  6, 10, 15, 21, 28, 36, 45])