Cross-Correlation and Cross-Covariance

Relating Two Processes

In many engineering systems, we observe the output of a system and want to relate it to the input. The input signal X(t)X(t) and the output signal Y(t)Y(t) are two different processes with a known relationship (e.g., Y(t)=h(t)βˆ—X(t)+N(t)Y(t) = h(t) * X(t) + N(t)). The cross-correlation function captures the statistical dependence between XX and YY at different time instants. It is the key to deriving the input-output relations for LTI systems (Chapter 15) and to designing optimal estimators like the Wiener filter.

Definition:

Cross-Correlation Function

For two second-order processes {X(t)}\{X(t)\} and {Y(t)}\{Y(t)\}, the cross-correlation is rXY(t1,t2)=E[X(t1)Yβˆ—(t2)].r_{XY}(t_1, t_2) = \mathbb{E}[X(t_1)Y^*(t_2)].

For jointly WSS processes, this depends only on Ο„=t1βˆ’t2\tau = t_1 - t_2: rXY(Ο„)=E[X(t)Yβˆ—(tβˆ’Ο„)].r_{XY}(\tau) = \mathbb{E}[X(t)Y^*(t - \tau)].

For discrete-time jointly WSS processes: rxy[k]=E[XnYnβˆ’kβˆ—]r_{xy}[k] = \mathbb{E}[X_n Y_{n-k}^*].

Definition:

Cross-Covariance Function

The cross-covariance of two processes is cXY(t1,t2)=Cov(X(t1),Y(t2))=rXY(t1,t2)βˆ’ΞΌX(t1)ΞΌYβˆ—(t2).c_{XY}(t_1, t_2) = \text{Cov}(X(t_1), Y(t_2)) = r_{XY}(t_1, t_2) - \mu_X(t_1)\mu_Y^*(t_2).

For jointly WSS processes with means ΞΌX\mu_X and ΞΌY\mu_Y: cXY(Ο„)=rXY(Ο„)βˆ’ΞΌXΞΌYβˆ—.c_{XY}(\tau) = r_{XY}(\tau) - \mu_X \mu_Y^*.

Definition:

Jointly Wide-Sense Stationary

Two processes {X(t)}\{X(t)\} and {Y(t)}\{Y(t)\} are jointly WSS if:

  1. Each is individually WSS, and
  2. The cross-correlation rXY(t1,t2)=rXY(t1βˆ’t2)r_{XY}(t_1, t_2) = r_{XY}(t_1 - t_2) depends only on the time difference.

Being individually WSS is necessary but not sufficient for jointly WSS. The cross-correlation must also be shift-invariant.

Theorem: Properties of Cross-Correlation (Jointly WSS)

Let {X(t)}\{X(t)\} and {Y(t)}\{Y(t)\} be jointly WSS. Then:

  1. Conjugate symmetry: rXY(Ο„)=rYXβˆ—(βˆ’Ο„)r_{XY}(\tau) = r_{YX}^*(-\tau). (For real processes: rXY(Ο„)=rYX(βˆ’Ο„)r_{XY}(\tau) = r_{YX}(-\tau).)

  2. Cauchy-Schwarz bound: ∣rXY(Ο„)βˆ£β‰€rXX(0)β‹…rYY(0)|r_{XY}(\tau)| \leq \sqrt{r_{XX}(0) \cdot r_{YY}(0)}.

  3. Generalized bound: ∣rXY(Ο„)∣2≀rXX(0)β‹…rYY(0)|r_{XY}(\tau)|^2 \leq r_{XX}(0) \cdot r_{YY}(0) for all Ο„\tau.

Note: rXY(Ο„)r_{XY}(\tau) is not necessarily even, and ∣rXY(Ο„)∣|r_{XY}(\tau)| need not be maximized at Ο„=0\tau = 0.

Definition:

Orthogonal and Uncorrelated Processes

Two WSS processes {X(t)}\{X(t)\} and {Y(t)}\{Y(t)\} are:

  • Orthogonal if rXY(Ο„)=0r_{XY}(\tau) = 0 for all Ο„\tau.
  • Uncorrelated if cXY(Ο„)=0c_{XY}(\tau) = 0 for all Ο„\tau, i.e., rXY(Ο„)=ΞΌXΞΌYβˆ—r_{XY}(\tau) = \mu_X \mu_Y^* for all Ο„\tau.

If either process has zero mean, orthogonality and uncorrelatedness are equivalent.

Example: Cross-Correlation of LTI Input and Output

Let Yn=βˆ‘k=0Lh[k]Xnβˆ’kY_n = \sum_{k=0}^{L} h[k] X_{n-k} be the output of a causal FIR filter with impulse response h[0],h[1],…,h[L]h[0], h[1], \ldots, h[L] driven by a zero-mean WSS input {Xn}\{X_n\}. Find rYX[m]=E[YnXnβˆ’mβˆ—]r_{YX}[m] = \mathbb{E}[Y_n X_{n-m}^*].

Example: Proper (Circularly Symmetric) Complex Processes

Let Zn=Xn+jYnZ_n = X_n + jY_n be a complex-valued WSS process with E[Zn]=0\mathbb{E}[Z_n] = 0. Define the pseudo-covariance c~ZZ[n,m]=E[ZnZm]\tilde{c}_{ZZ}[n,m] = \mathbb{E}[Z_n Z_m] (without conjugation). Show that if ZnZ_n is proper (c~ZZ[n,m]=0\tilde{c}_{ZZ}[n,m] = 0 for all n,mn,m), then cXX[k]=cYY[k]c_{XX}[k] = c_{YY}[k] and cXY[k]=βˆ’cYX[k]c_{XY}[k] = -c_{YX}[k].

Quick Check

Two zero-mean WSS processes X(t)X(t) and Y(t)Y(t) satisfy rXY(Ο„)=0r_{XY}(\tau) = 0 for all Ο„\tau. Which of the following is true?

X(t)X(t) and Y(t)Y(t) are independent

X(t)X(t) and Y(t)Y(t) are both orthogonal and uncorrelated

rXX(au)=rYY(au)r_{XX}( au) = r_{YY}( au)

Common Mistake: Assuming Cross-Correlation Is Even

Mistake:

Writing rXY(Ο„)=rXY(βˆ’Ο„)r_{XY}(\tau) = r_{XY}(-\tau), by analogy with the autocorrelation.

Correction:

The cross-correlation satisfies rXY(Ο„)=rYXβˆ—(βˆ’Ο„)r_{XY}(\tau) = r_{YX}^*(-\tau), not rXY(Ο„)=rXY(βˆ’Ο„)r_{XY}(\tau) = r_{XY}(-\tau). In general, rXYr_{XY} is neither even nor odd. Only the auto-correlation has the even symmetry property.

Orthogonal vs. Uncorrelated vs. Independent Processes

ConceptConditionZero-mean simplification
OrthogonalrXY(Ο„)=0r_{XY}(\tau) = 0 for all Ο„\tauSame as uncorrelated
UncorrelatedcXY(Ο„)=0c_{XY}(\tau) = 0 for all Ο„\tauSame as orthogonal
IndependentJoint fdds factor for all ordersImplies uncorrelated (not converse)
πŸ”§Engineering Note

System Identification via Cross-Correlation

If a BIBO-stable LTI system with impulse response h[n]h[n] is excited by zero-mean white noise with rxx[k]=Οƒ2Ξ΄[k]r_{xx}[k] = \sigma^2\delta[k], then the input-output cross-correlation is ryx[k]=Οƒ2h[k]r_{yx}[k] = \sigma^2 h[k]. This gives a direct method for estimating h[k]h[k]: inject white noise and cross-correlate the output with the input. In wireless systems, this principle underlies pilot-based channel estimation, where known pilot symbols play the role of the white-noise input.

Cross-Correlation

rXY(Ο„)=E[X(t)Yβˆ—(tβˆ’Ο„)]r_{XY}(\tau) = \mathbb{E}[X(t)Y^*(t-\tau)] for jointly WSS processes. Measures the linear dependence between two processes at different time instants.

Related: Autocorrelation Function

Orthogonal Processes

Two processes with rXY(Ο„)=0r_{XY}(\tau) = 0 for all Ο„\tau. For zero-mean processes, equivalent to being uncorrelated.

Related: Cross-Correlation

Key Takeaway

The cross-correlation rXY(Ο„)r_{XY}(\tau) measures the linear relationship between two processes at different time instants. Unlike the autocorrelation, rXY(Ο„)r_{XY}(\tau) is generally neither even nor maximized at the origin. The key relationship rYX[k]=(hβˆ—rXX)[k]r_{YX}[k] = (h * r_{XX})[k] for an LTI system is the foundation of the Wiener filter and pilot-based channel estimation.