相关疑难解决方法(0)

有没有办法在Numpy中有效地实现1D阵列的滚动窗口？

例如,我有这个纯Python代码片段来计算1D列表的滚动标准偏差,其中observations是1D值列表,并且n是标准差的窗口长度:

stdev = []
for i, data in enumerate(observations[n-1:]):
    strip = observations[i:i+n]
    mean = sum(strip) / n
    stdev.append(sqrt(250*sum([(s-mean)**2 for s in strip])/(n-1)))

有没有办法在Numpy中完全执行此操作,即没有任何Python循环？标准偏差是微不足道的numpy.std,但滚动窗口部分完全残留我.

我发现这篇关于Numpy滚动窗口的博文,但它似乎不适用于1D阵列.