Before diving into a more complex statistical analysis of random signals and
processes, let us quickly review the idea of correlation. Recall that
the correlation of two signals or variables is the expected
value of the product of those two variables. Since our focus
will be to discover more about a random process, a collection of
random signals, then imagine us dealing with two samples of a
random process, where each sample is taken at a different point
in time. Also recall that the key property of these random
processes is that they are now functions of time; imagine them
as a collection of signals. The expected value of the product of these two
variables (or samples) will now depend on how quickly they
change in regards to time. For example, if
the two variables are taken from almost the same time period,
then we should expect them to have a high correlation. We will
now look at a correlation function that relates a pair of random
variables from the same process to the time separations between
them, where the argument to this correlation function will be
the time difference. For the correlation of signals from two
different random process, look at the crosscorrelation function.
The first of these correlation functions we will discuss is
the autocorrelation, where each of the random
variables we will deal with come from the same random process.
- Definition 1:
Autocorrelation
the expected value of the product of a random variable or
signal realization with a time-shifted version of itself
With a simple calculation and analysis of the autocorrelation
function, we can discover a few important characteristics
about our random process. These include:
-
How quickly our random signal or processes changes with
respect to the time function
-
Whether our process has a periodic component and what the
expected frequency might be
As was mentioned above, the autocorrelation function is simply
the expected value of a product. Assume we have a pair of
random variables from the same process,
X
1
=X
t
1
X
1
X
t
1
and
X
2
=X
t
2
X
2
X
t
2
, then the autocorrelation is often written as
R
xx
t
1
t
2
=E
X
1
X
2
=∫-∞∞∫-∞∞
x
1
x
2
f
x
1
x
2
d
x
2
d
x
1
R
xx
t
1
t
2
X
1
X
2
x
1
x
2
x
1
x
2
f
x
1
x
2
(1)
The above equation is valid for stationary and nonstationary
random processes. For
stationary processes, we can generalize
this expression a little further. Given a wide-sense
stationary processes, it can be proven that the expected
values from our random process will be independent of the
origin of our time function. Therefore, we can say that our
autocorrelation function will depend on the time difference
and not some absolute time. For this discussion, we will let
τ=
t
2
-
t
1
τ
t
2
t
1
, and thus we generalize our autocorrelation
expression as
R
xx
tt+τ=
R
xx
τ=EXtXt+τ
R
xx
t
t
τ
R
xx
τ
X
t
X
t
τ
(2)
for the continuous-time case. In most DSP course we will be
more interested in dealing with real signal sequences, and thus
we will want to look at the discrete-time case of the
autocorrelation function. The formula below will prove to be
more common and useful than
Equation 1:
R
xx
nn+m=∑n=-∞∞xnxn+m
R
xx
n
n
m
n
x
n
x
n
m
(3)
And again we can generalize the notation for our
autocorrelation function as
R
xx
nn+m=
R
xx
m=EXnXn+m
R
xx
n
n
m
R
xx
m
X
n
X
n
m
(4)
Below we will look at several properties of the
autocorrelation function that hold for
stationary random processes.
-
Autocorrelation is an even function for τ
τ
R
xx
τ=
R
xx
-τ
R
xx
τ
R
xx
τ
-
The mean-square value can be found by evaluating the
autocorrelation where
τ=0
τ
0
, which gives us
R
xx
0=X2¯
R
xx
0
X
2
-
The autocorrelation function will have its largest value
when
τ=0
τ
0
. This value can appear again, for example in
a periodic function at the values of the equivalent
periodic points, but will never be exceeded.
R
xx
0≥|
R
xx
τ|
R
xx
0
R
xx
τ
-
If we take the autocorrelation of a period function,
then
R
xx
τ
R
xx
τ
will also be periodic with the same frequency.
Sometimes the whole random process is not available to us.
In these cases, we would still like to be able to find out
some of the characteristics of the stationary random
process, even if we just have part of one sample function.
In order to do this we can estimate the
autocorrelation from a given interval,
0 0 to T
T seconds, of the sample function.
Ř
xx
τ=1T-τ∫0T-τxtxt+τdt
Ř
xx
τ
1
T
τ
t
T
τ
0
x
t
x
t
τ
(5)
However, a lot of times we will not have sufficient
information to build a complete continuous-time function of
one of our random signals for the above analysis. If this
is the case, we can treat the information we do know about
the function as a discrete signal and use the discrete-time
formula for estimating the autocorrelation.
Ř
xx
m=1N-m∑n=0N-m-1xnxn+m
Ř
xx
m
1
N
m
n
N
m
1
0
x
n
x
n
m
(6)
Below we will look at a variety of examples that use the
autocorrelation function. We will begin with a simple example
dealing with Gaussian White Noise (GWN) and a few basic
statistical properties that will prove very useful in these
and future calculations.
We will let
xn
x
n
represent our GWN. For this problem, it is
important to remember the following fact about the mean of a
GWN function:
Exn=0
x
n
0
Along with being zero-mean, recall that
GWN is always independent. With these
two facts, we are now ready to do the short calculations
required to find the autocorrelation.
R
xx
nn+m=Exnxn+m
R
xx
n
n
m
x
n
x
n
m
Since the function,
xn
x
n
, is independent, then we can take the product of
the individual expected values of both functions.
R
xx
nn+m=ExnExn+m
R
xx
n
n
m
x
n
x
n
m
Now, looking at the above equation we see that we can break
it up further into two conditions: one when
m m and
nn are equal and one when they
are not equal. When they are equal we can combine the
expected values. We are left with the following piecewise
function to solve:
R
xx
nn+m=
ExnExn+mifm≠0Ex2nifm=0
R
xx
n
n
m
x
n
x
n
m
m
0
x
n
2
m
0
We can now solve the two parts of the above equation. The
first equation is easy to solve as we have already stated
that the expected value of
xn
x
n
will be zero. For the second part, you should
recall from statistics that the expected value of the square
of a function is equal to the variance. Thus we get the
following results for the autocorrelation:
R
xx
nn+m=
0ifm≠0σ2ifm=0
R
xx
n
n
m
0
m
0
σ
2
m
0
Or in a more concise way, we can represent the results as
R
xx
nn+m=σ2δm
R
xx
n
n
m
σ
2
δ
m