Estimation

An good estimator for the correlation function5 for deterministic signals could be:

\begin{displaymath}
\hat{r}_{xx}[n] = \frac{1}{L}\sum_{k=0}^{L-\vert n\vert+1} x[k]^* x[\vert n\vert+k]
\end{displaymath}

We define bias of an estimator $\hat{X}$ of a process or random variable $X$ as:

\begin{displaymath}
\hat{\Theta}_{\hat{X}} = \vert E[\hat{X}]-X\vert
\end{displaymath}

To check how good it is:

\begin{displaymath}E\left\lbrace \hat{r}_{xx}[n] \right\rbrace = \frac{1}{L}\sum...
... n\vert+k] \right\rbrace
= \frac{L-\vert n\vert}{L} r_{xx}[n]
\end{displaymath}

Which shows that the estimator of the correlation is biased,If we take the limit as the number of samples $L$ goes to infinity:

\begin{displaymath}
\lim_{L\to\infty} E\left\lbrace \hat{r}_{xx}[n] \right\rbra...
...\lim_{L\to\infty} \frac{L-\vert n\vert}{L} r_{xx}[n] = r_{xx}
\end{displaymath}

So the autocorrelation estimator $\hat{r}_{xx}[n]$ is biased, but asimptotically unbiased. In fact, an unbiased estimator of the correlation would be:

\begin{displaymath}
\check{r}_{xx}[n] = \frac{1}{\frac{L-\vert n\vert}{L}} \hat{...
...n\vert} \sum_{k=0}^{L-\vert n\vert+1} x[k]^* x[\vert n\vert+k]
\end{displaymath}


\begin{displaymath}
E\left\lbrace \check{r}_{xx}[n] \right\rbrace = r_{xx}[n]
\end{displaymath}

But this last one would be seldom used, and we will often use the $\hat{r}_{xx}[n]$ as the estimator for the autocorrelation. $\hat{r}_{xx}[n]$ is called the biased estimator of the correlation function, $\check{r}_{xx}[n]$ is called the unbiased estimator of the correlation function.

The variance of an estimator $\hat{X}$ of a process or random variable $X$ is defined as:

\begin{displaymath}
\sigma_{\hat{X}}^2 = E\left\lbrace (\hat{X} - E\left\lbrace X \right\rbrace)^2 \right\rbrace
\end{displaymath}

Pedro Larroy 2005-04-29