TOPICS
Search

Autocorrelation


Autocorrelation of a real-valued stochastic process {X_t} measures the correlation of values separated by a given lag. For a weakly stationary process with mean mu, the autocovariance and autocorrelation functions are

gamma(h)=E[(X_t-mu)(X_(t+h)-mu)]
(1)
rho(h)=(gamma(h))/(gamma(0)),
(2)

where E[Y] denotes the expectation value of Y and gamma(0)>0 is assumed when defining rho(h). They depend only on the lag h, not on the time t. Their sample counterpart is the sample autocorrelation.

For a finite sequence whose indices are taken modulo its length, the corresponding notion is periodic autocorrelation. The non-wrapping counterpart is aperiodic autocorrelation.

For a square integrable complex function f, the autocorrelation is defined by

f*f=int_(-infty)^inftyf(tau)f^_(tau-t)dtau
(3)
=int_(-infty)^inftyf^_(tau)f(tau+t)dtau,
(4)

where * denotes cross-correlation and f^_ is the complex conjugate (Bracewell 1965, pp. 40-41).

Note that the notation rho_f(t) is sometimes used for f*f and that the quantity

 R_f(t)=lim_(T->infty)1/(2T)int_(-T)^Tf(tau)f(t+tau)dtau
(5)

is sometimes also known as the autocorrelation of a continuous real function f(t) when the limit exists (Papoulis 1962, p. 241).

The autocorrelation discards phase information, returning only the power, and is therefore an irreversible operation.

There is also a somewhat surprising and extremely important relationship between the autocorrelation and the Fourier transform known as the Wiener-Khinchin theorem. Let F_t[f(t)](omega)=F(omega), and F^_ denote the complex conjugate of F, then the inverse Fourier transform of the absolute square of F(omega) is given by

 F_omega^(-1)[|F(omega)|^2](t)=int_(-infty)^inftyf^_(tau)f(tau+t)dtau.
(6)

For the square integrable complex function f, Schwarz's inequality gives

 |[f*f](t)|<=int_(-infty)^infty|f(tau)|^2dtau=[f*f](0).
(7)

Thus its complex modulus is maximum at the origin. For a square integrable real function f, f*f itself is maximum at the origin; in other words,

 int_(-infty)^inftyf(tau)f(tau+x)dtau<=int_(-infty)^inftyf^2(tau)dtau.
(8)

To see this, let epsilon be a real number. Then

 int_(-infty)^infty[f(tau)+epsilonf(tau+x)]^2dtau>=0
(9)
 int_(-infty)^inftyf^2(tau)dtau+2epsilonint_(-infty)^inftyf(tau)f(tau+x)dtau+epsilon^2int_(-infty)^inftyf^2(tau+x)dtau>=0
(10)
 int_(-infty)^inftyf^2(tau)dtau+2epsilonint_(-infty)^inftyf(tau)f(tau+x)dtau+epsilon^2int_(-infty)^inftyf^2(tau)dtau>=0.
(11)

Define

a=int_(-infty)^inftyf^2(tau)dtau
(12)
b=2int_(-infty)^inftyf(tau)f(tau+x)dtau.
(13)

Then the preceding inequality is aepsilon^2+bepsilon+a>=0 for every real number epsilon. If a=0, the result is immediate. If a>0, the discriminant satisfies b^2-4a^2<=0, so |b|<=2a and hence b/2<=a. It follows that

 int_(-infty)^inftyf(tau)f(tau+x)dtau<=int_(-infty)^inftyf^2(tau)dtau,
(14)

with the equality at x=0. This proves that f*f is maximum at the origin.


See also

Aperiodic Autocorrelation, Autocovariance, Average Power, Convolution, Correlation, Cross-Correlation, Expectation Value, Periodic Autocorrelation, Quantization Efficiency, Recurrence Plot, Sample Autocorrelation, Stochastic Process, Weak Stationarity, Wiener-Khinchin Theorem

Explore with Wolfram|Alpha

References

Bracewell, R. "The Autocorrelation Function." The Fourier Transform and Its Applications. New York: McGraw-Hill, pp. 40-45, 1965.Papoulis, A. The Fourier Integral and Its Applications. New York: McGraw-Hill, 1962.Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. "Correlation and Autocorrelation Using the FFT." §13.2 in Numerical Recipes in FORTRAN: The Art of Scientific Computing, 2nd ed. Cambridge, England: Cambridge University Press, pp. 538-539, 1992.Zwillinger, D. (Ed.). CRC Standard Mathematical Tables and Formulae. Boca Raton, FL: CRC Press, p. 223, 1995.

Referenced on Wolfram|Alpha

Autocorrelation

Cite this as:

Weisstein, Eric W. "Autocorrelation." From MathWorld--A Wolfram Resource. https://mathworld.wolfram.com/Autocorrelation.html

Subject classifications