Я хотел бы выполнить автокорреляцию по сигналу, показанному ниже. Время между двумя последовательными точками составляет 2,5 мс (или частота повторения 400 Гц).

Это уравнение для оценки автоакреляции, которое я хотел бы использовать (взято изhttp://en.wikipedia.org/wiki/Autocorrelation, раздел Оценка):

Какой самый простой метод определения предполагаемой автокорреляции моих данных в python? Есть ли что-то похожее наnumpy.correlate что я могу использовать?

Или я должен просто рассчитать среднее значение и дисперсию?

Редактировать:

С помощьюunutbu, Я написал:

from numpy import *
import numpy as N
import pylab as P

fn = 'data.txt'
x = loadtxt(fn,unpack=True,usecols=[1])
time = loadtxt(fn,unpack=True,usecols=[0]) 

def estimated_autocorrelation(x):
    n = len(x)
    variance = x.var()
    x = x-x.mean()
    r = N.correlate(x, x, mode = 'full')[-n:]
    #assert N.allclose(r, N.array([(x[:n-k]*x[-(n-k):]).sum() for k in range(n)]))
    result = r/(variance*(N.arange(n, 0, -1)))
    return result

P.plot(time,estimated_autocorrelation(x))
P.xlabel('time (s)')
P.ylabel('autocorrelation')
P.show()