Misplaced Pages

Autocorrelation technique

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
The topic of this article may not meet Misplaced Pages's general notability guideline. Please help to demonstrate the notability of the topic by citing reliable secondary sources that are independent of the topic and provide significant coverage of it beyond a mere trivial mention. If notability cannot be shown, the article is likely to be merged, redirected, or deleted.
Find sources: "Autocorrelation technique" – news · newspapers · books · scholar · JSTOR (February 2012) (Learn how and when to remove this message)
This article includes a list of references, related reading, or external links, but its sources remain unclear because it lacks inline citations. Please help improve this article by introducing more precise citations. (February 2012) (Learn how and when to remove this message)

The autocorrelation technique is a method for estimating the dominating frequency in a complex signal, as well as its variance. Specifically, it calculates the first two moments of the power spectrum, namely the mean and variance. It is also known as the pulse-pair algorithm in radar theory.

The algorithm is both computationally faster and significantly more accurate compared to the Fourier transform, since the resolution is not limited by the number of samples used.

Derivation

The autocorrelation of lag 1 can be expressed using the inverse Fourier transform of the power spectrum S ( ω ) {\displaystyle S(\omega )} :

R ( 1 ) = 1 2 π ∫ − π π S ( ω ) e i ω d ω . {\displaystyle R(1)={\frac {1}{2\pi }}\int _{-\pi }^{\pi }S(\omega )e^{i\,\omega }d\omega .}

If we model the power spectrum as a single frequency S ( ω )   = d e f   δ ( ω − ω 0 ) {\displaystyle S(\omega )\ {\stackrel {\mathrm {def} }{=}}\ \delta (\omega -\omega _{0})} , this becomes:

R ( 1 ) = 1 2 π ∫ − π π δ ( ω − ω 0 ) e i ω d ω {\displaystyle R(1)={\frac {1}{2\pi }}\int _{-\pi }^{\pi }\delta (\omega -\omega _{0})e^{i\,\omega }d\omega }
R ( 1 ) = 1 2 π e i ω 0 {\displaystyle R(1)={\frac {1}{2\pi }}e^{i\,\omega _{0}}}

where it is apparent that the phase of R ( 1 ) {\displaystyle R(1)} equals the signal frequency.

Implementation

The mean frequency is calculated based on the autocorrelation with lag one, evaluated over a signal consisting of N samples:

ω = ∠ R N ( 1 ) = tan − 1 ⁡ im { R N ( 1 ) } re { R N ( 1 ) } . {\displaystyle \omega =\angle R_{N}(1)=\tan ^{-1}{\frac {{\text{im}}\{R_{N}(1)\}}{{\text{re}}\{R_{N}(1)\}}}.}

The spectral variance is calculated as follows:

var { ω } = 2 N ( 1 − | R N ( 1 ) | R N ( 0 ) ) . {\displaystyle {\text{var}}\{\omega \}={\frac {2}{N}}\left(1-{\frac {|R_{N}(1)|}{R_{N}(0)}}\right).}

Applications

External links

Categories: