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 main
focus is to discover more about random processes, a collection
of random signals, we will deal with two random processes in
this discussion, where in this case we will deal with samples
from two different random processes. We
will analyze the expected value of the product of these two
variables and how they correlate to one another, where the
argument to this correlation function will be the time
difference. For the correlation of signals from the same random
process, look at the autocorrelation function.
When dealing with multiple random processes, it is also
important to be able to describe the relationship, if any,
between the processes. For example, this may occur if more
than one random signal is applied to a system. In order to do
this, we use the crosscorrelation function, where
the variables are instances from two different wide sense
stationary random processes.
- Definition 1: Crosscorrelation
if two processes are wide sense stationary, the expected
value of the product of a random variable from one random
process with a time-shifted, random variable from a
different random process
Looking at the generalized formula for the crosscorrelation,
we will represent our two random processes by allowing
U=Ut
U
U
t
and
V=Vt−τ
V
V
t
τ
. We will define the crosscorrelation function as
R
u
v
tt−τ=EUV=∫-∞∞∫-∞∞uvfuvdvdu
R
u
v
t
t
τ
U
V
u
v
u
v
f
u
v
(1)
Just as the case with the autocorrelation function, if our
input and output, denoted as
Ut
U
t
and
Vt
V
t
, are at least jointly wide sense stationary, then the
crosscorrelation does not depend on absolute time; it is just
a function of the time difference. This means we can simplify
our writing of the above function as
R
u
v
τ=EUV
R
u
v
τ
U
V
(2)
or if we deal with two real signal sequences,
xn
x
n
and
yn
y
n
, then we arrive at a more commonly seen formula for
the discrete crosscorrelation function. See the formula below
and notice the similarities between it and the
convolution of two
signals:
R
x
y
nn−m=
R
x
y
m=∑n=-∞∞xnyn−m
R
x
y
n
n
m
R
x
y
m
n
x
n
y
n
m
(3)
Below we will look at several properties of the
crosscorrelation function that hold for two wide
sense stationary (WSS) random processes.
-
Crosscorrelation is not an even
function; however, it does have a unique symmetry
property:
R
x
y
-τ=
R
y
x
τ
R
x
y
τ
R
y
x
τ
(4)
-
The maximum value of the crosscorrelation is not always
when the shift equals zero; however, we can prove the
following property revealing to us what value the
maximum cannot exceed.
|
R
x
y
τ|≤
R
x
x
0
R
y
y
0
R
x
y
τ
R
x
x
0
R
y
y
0
(5)
-
When two random processes are statistically independent
then we have
R
x
y
τ=
R
y
x
τ
R
x
y
τ
R
y
x
τ
(6)
Let us begin by looking at a simple example showing the
relationship between two sequences. Using Equation 3, find the
crosscorrelation of the sequences
xn=
…002-361300…
x
n
…
0
0
2
-3
6
1
3
0
0
…
yn=
…001-241-300…
y
n
…
0
0
1
-2
4
1
-3
0
0
…
for each of the following possible time shifts:
m=
03-1
m
0
3
-1
.
-
For
m=0
m
0
, we should begin by finding the product
sequence
sn=xnyn
s
n
x
n
y
n
. Doing this we get the following sequence:
sn=
…0026241-900…
s
n
…
0
0
2
6
24
1
-9
0
0
…
and so from the sum in our crosscorrelation function
we arrive at the answer of
R
x
y
0=22
R
x
y
0
22
-
For
m=3
m
3
, we will approach it the same was we did
above; however, we will now shift
yn
y
n
to the right. Then we can find the product sequence
sn=xnyn−3
s
n
x
n
y
n
3
, which yields
sn=
…000001-600…
s
n
…
0
0
0
0
0
1
-6
0
0
…
and from the crosscorrelation function we arrive at
the answer of
R
x
y
3=-6
R
x
y
3
-6
-
For
m=-1
m
-1
, we will again take the same approach;
however, we will now shift
yn
y
n
to the left. Then we can find the product sequence
sn=xnyn+1
s
n
x
n
y
n
1
, which yields
sn=
…00-4-126-3000…
s
n
…
0
0
-4
-12
6
-3
0
0
0
…
and from the crosscorrelation function we arrive at
the answer of
R
x
y
-1=-13
R
x
y
-1
-13