Skip to content

Energy, Power, and Autocorrelation

the section gives the full continuous- and discrete-time definitions and classification examples. The decisive tests are recalled here because the correct correlation definition depends on the signal class:

A nonzero signal cannot belong to both classes. Some signals belong to neither, and the all-zero signal is normally excluded by the strict positive conditions. For periodic signals whose squared magnitude is integrable or summable over one period,

The starting point is arbitrary because a complete period is averaged.

For the angular-frequency Fourier convention used in this chapter, a CT energy signal with transform X(ω)X(\omega) has energy spectral density (ESD)

For a DT energy signal with DTFT X(ejΩ)X(e^{\mathrm{j}\Omega}),

Ψx(ejΩ)=∣X(ejΩ)∣2,Ex=12π∫−ππΨx(ejΩ) dΩ.\Psi_x(e^{\mathrm{j}\Omega})=|X(e^{\mathrm{j}\Omega})|^2, \qquad E_x=\frac{1}{2\pi}\int_{-\pi}^{\pi} \Psi_x(e^{\mathrm{j}\Omega})\,\mathrm{d}\Omega.

ESD therefore states where a finite total energy lies in frequency.

A power signal generally has infinite total energy, so its ordinary Fourier transform cannot simply be squared as a finite-energy function. Its power spectral density (PSD) is instead defined through autocorrelation. Under the same angular-frequency convention,

with an integral over any 2π2\pi interval for a DT PSD. ESD has units of energy per unit frequency and integrates to energy; PSD has units of mean-square value per unit frequency and integrates to average power.

Periodic power signals have line spectra. If x(t)=∑r=−∞∞Crejrω0tx(t)=\sum_{r=-\infty}^{\infty}C_r e^{\mathrm{j}r\omega_0t}, then

This is the power-signal counterpart of Parseval’s theorem.

Autocorrelation measures similarity between a signal and a delayed copy of itself. This section consistently uses the lag convention x∗(t−τ)x^*(t-\tau) or x∗[n−m]x^*[n-m].

For CT and DT energy signals,

These are inner products between the signal and a shifted version. For real signals the conjugate has no visible effect, but it must be retained for complex baseband signals.

For deterministic power signals, replace total integration or summation by a long-time average:

provided the limits exist. For a CT periodic signal, the equivalent one-period expression is

Rxx(τ)=1T0∫t0t0+T0x(t)x∗(t−τ) dt.R_{xx}(\tau)=\frac{1}{T_0} \int_{t_0}^{t_0+T_0}x(t)x^*(t-\tau)\,\mathrm{d}t.

The DT periodic formula similarly averages over one complete period.

A typical real finite-pulse autocorrelation. It is even, reaches its maximum magnitude at zero lag, and has R_(xx)(0) equal to signal energy. Other signals need not have this triangular shape.

A typical real finite-pulse autocorrelation. It is even, reaches its maximum magnitude at zero lag, and has Rxx(0)R_{xx}(0) equal to signal energy. Other signals need not have this triangular shape.

For a complex random process X(t)X(t), the general second-order correlation is

RX(t1,t2)=E ⁣{X(t1)X∗(t2)}.R_X(t_1,t_2)=\mathbb{E}\!\left\{X(t_1)X^*(t_2)\right\}.

If the process is wide-sense stationary (WSS), its mean is constant and its correlation depends only on the lag τ=t1−t2\tau=t_1-t_2:

Expectation is an ensemble average. A time average from one realization equals it only under an appropriate ergodicity assumption; stationarity by itself does not guarantee that equality.

Autocorrelation must also be distinguished from autocovariance. If the WSS mean is μX\mu_X, then

CX(τ)=RX(τ)−∣μX∣2.C_X(\tau)=R_X(\tau)-|\mu_X|^2.

Thus RX(0)=E{∣X(t)∣2}R_X(0)=\mathbb{E}\{|X(t)|^2\} is mean square, not variance unless the mean is zero. A nonzero mean creates a constant floor in autocorrelation and a DC spectral line in the PSD; mean removal is often necessary before estimating periodicity or delay.

The following properties hold when the defining energy, time-average power, or WSS second moment exists.

  1. Value at zero: Zero lag equals the relevant signal measure:
Rxx(0)=Exfor an energy signal,Rxx(0)=Pxfor a power signal.R_{xx}(0)=E_x\quad\text{for an energy signal}, \qquad R_{xx}(0)=P_x\quad\text{for a power signal}.

The DT statements use Rxx[0]R_{xx}[0]; for a WSS process it is the mean-square value.

  1. Hermitian symmetry: Autocorrelation obeys
Rxx(−τ)=Rxx∗(τ),Rxx[−m]=Rxx∗[m].R_{xx}(-\tau)=R_{xx}^*(\tau), \qquad R_{xx}[-m]=R_{xx}^*[m].

It is real and even for a real signal. For a complex signal it need not be real or even; its real part is even and its imaginary part is odd.

  1. Maximum magnitude at zero: Cauchy–Schwarz gives
∣Rxx(τ)∣≤Rxx(0),∣Rxx[m]∣≤Rxx[0].|R_{xx}(\tau)|\leq R_{xx}(0), \qquad |R_{xx}[m]|\leq R_{xx}[0].

This is a magnitude bound, not an ordering of complex values. Equality can also occur at a nonzero lag when a shifted signal is proportional to the original, as for periodic signals; zero lag need not be the unique peak.

  1. Periodicity: If x(t)x(t) has period T0T_0, then Rxx(τ)R_{xx}(\tau) is also periodic with period T0T_0. The analogous statement holds for a DT period N0N_0.

  2. Nonnegative definiteness: For arbitrary complex constants cic_i and lags τi\tau_i,

∑i∑jcicj∗Rxx(τi−τj)≥0.\sum_i\sum_j c_i c_j^*R_{xx}(\tau_i-\tau_j)\geq0.

An autocorrelation may still take negative pointwise values; this condition is what guarantees a nonnegative PSD.

  1. Ideal white noise: A zero-mean ideal white-noise process has
Rww(τ)=N02δ(τ),Rww[m]=σw2δ[m].R_{ww}(\tau)=\frac{N_0}{2}\delta(\tau), \qquad R_{ww}[m]=\sigma_w^2\delta[m].

Distinct lags are uncorrelated. The impulse correlation and flat infinite-bandwidth PSD describe an idealization, not every random process; physical white noise is bandwidth-limited.

The Fourier transform of autocorrelation gives spectral density. For a CT energy signal this is the ESD; for a deterministic power signal or WSS random process it is the PSD:

For an energy signal, direct substitution gives

Sxx(ω)=∣X(ω)∣2.S_{xx}(\omega)=|X(\omega)|^2.

The DT counterpart is the Fourier-series pair

The PSD is real and nonnegative. Setting lag to zero recovers total energy or average power by integrating the applicable density.

For a finite NN-sample record, define circular autocorrelation by

rxx(c)[m]=∑n=0N−1x[n]x∗[(n−m)N].r_{xx}^{(c)}[m]=\sum_{n=0}^{N-1}x[n]x^*[(n-m)_N].

It obeys the exact DFT pair

rxx(c)[m]⟷∣X[k]∣2.r_{xx}^{(c)}[m]\longleftrightarrow |X[k]|^2.

An FFT and IFFT therefore compute circular correlation efficiently. Linear correlation requires sufficient zero padding to prevent wraparound, exactly as for linear convolution.

Cross-Correlation and Cross-Spectral Density

Section titled “Cross-Correlation and Cross-Spectral Density”

Cross-correlation compares two different signals. For CT and DT energy signals,

Power-signal versions use time averages. Jointly WSS random processes use the ensemble correlation

RXY(τ)=E ⁣{X(t)Y∗(t−τ)}.R_{XY}(\tau)=\mathbb{E}\!\left\{X(t)Y^*(t-\tau)\right\}.

Unlike autocorrelation, cross-correlation is generally neither even nor real:

Its Fourier transform is the cross-spectral density Sxy(ω)=F{Rxy(τ)}S_{xy}(\omega)=\mathcal{F}\{R_{xy}(\tau)\}; for energy signals, Sxy(ω)=X(ω)Y∗(ω)S_{xy}(\omega)=X(\omega)Y^*(\omega).

The correlation peak locates relative delay, but its sign depends on order and lag convention. With the convention used here, if y(t)=x(t−τ0)y(t)=x(t-\tau_0), then Rxy(τ)R_{xy}(\tau) peaks at −τ0-\tau_0, while Ryx(τ)R_{yx}(\tau) peaks at +τ0+\tau_0.

Correlation is convolution with a conjugated, time-reversed second signal:

Thus the same direct or FFT-based convolution routine can compute correlation after conjugate reversal of one input.

For finite, nonzero signal energies, a dimensionless coefficient is

ρxy(τ)=Rxy(τ)ExEy,∣ρxy(τ)∣≤1.\rho_{xy}(\tau)=\frac{R_{xy}(\tau)}{\sqrt{E_xE_y}}, \qquad |\rho_{xy}(\tau)|\leq1.

For power signals or WSS processes, replace Ex,EyE_x,E_y by the corresponding zero-lag mean-square values. Several practical qualifications matter:

  • subtract sample means when the aim is waveform similarity rather than a shared DC level;

  • for finite records, use lag-dependent norms over the samples that actually overlap, or explicitly accept edge bias from a fixed denominator;

  • the peak may repeat for periodic data and may broaden or move because of filtering, noise, or multipath;

  • normalization is undefined when either denominator is zero, and high correlation alone does not establish causation.

Autocorrelation detects hidden periodicity, estimates pitch and symbol or clock timing, and converts naturally to PSD. Cross-correlation estimates relative delay and similarity in radar and sonar ranging, communication synchronization, matched filtering, channel sounding, and pattern matching.