python - Python 中加速度样本数组的速度

Question

我有一个相当大的加速样本列表(大约 300k)，我用 matplotlib 绘制了它们。我的目标是绘制速度并获得其最大值。

有没有办法整合加速度列表以获得速度列表？

Answer 1

我认为您正在寻找scipy.integrate.cumtrapz。它会比上面的 for 循环实现快得多。请参阅此答案 here 。基本思想是取普通梯形公式的累积和(而不是简单的和):

def cumtrapz_example(y, x):
     return np.cumsum(0.5*((x[1:]-x[:-1])*(y[1:]+y[:-1])))

而不是

def trapz_example(y, x):
     return np.sum(0.5*((x[1:]-x[:-1])*(y[1:]+y[:-1])))

这个lesson关于 Python 中的数值积分也很有帮助。