Main Content

Spectral Analysis

Background Information

The goal of spectral estimation is to describe the distribution (over frequency) of the power contained in a signal, based on a finite set of data. Estimation of power spectra is useful in a variety of applications, including the detection of signals buried in wideband noise.

The power spectral density (PSD) of a stationary random process x(n) is mathematically related to the autocorrelation sequence by the discrete-time Fourier transform. In terms of normalized frequency, this is given by

Pxx(ω)=12πm=Rxx(m)ejωm.

This can be written as a function of physical frequency f (for example, in hertz) by using the relation ω = 2πf / fs, where fs is the sampling frequency:

Pxx(f)=1fsm=Rxx(m)ej2πmf/fs.

The correlation sequence can be derived from the PSD by use of the inverse discrete-time Fourier transform:

Rxx(m)=ππPxx(ω)ejωmdω=fs/2fs/2Pxx(f)ej2πmf/fsdf.

The average power of the sequence x(n) over the entire Nyquist interval is represented by

Rxx(0)=ππPxx(ω)dω=fs/2fs/2Pxx(f)df.

The average power of a signal over a particular frequency band [ω1ω2], 0 ≤ ω1 ≤ ω2 ≤ π, can be found by integrating the PSD over that band:

P¯[ω1,ω2]=ω1ω2Pxx(ω)dω=ω2ω1Pxx(ω)dω.

You can see from the above expression that Pxx(ω) represents the power content of a signal in an infinitesimal frequency band, which is why it is called the power spectral density.

The units of the PSD are power (e.g., watts) per unit of frequency. In the case of Pxx(ω), this is watts/radian/sample or simply watts/radian. In the case of Pxx(f), the units are watts/hertz. Integration of the PSD with respect to frequency yields units of watts, as expected for the average power.

For real–valued signals, the PSD is symmetric about DC, and thus Pxx(ω) for 0 ≤ ω ≤ π is sufficient to completely characterize the PSD. However, to obtain the average power over the entire Nyquist interval, it is necessary to introduce the concept of the one-sided PSD.

The one-sided PSD is given by

Pone-sided(ω)={0,πω<0,2Pxx(ω),0ωπ.

The average power of a signal over the frequency band, [ω1,ω2] with 0 ≤ ω1 ≤ ω2 ≤ π, can be computed using the one-sided PSD as

P¯[ω1,ω2]=ω1ω2Pone-sided(ω)dω.

Spectral Estimation Method

The various methods of spectrum estimation available in the toolbox are categorized as follows:

  • Nonparametric methods

  • Parametric methods

  • Subspace methods

Nonparametric methods are those in which the PSD is estimated directly from the signal itself. The simplest such method is the periodogram. Other nonparametric techniques such as Welch's method [8], the multitaper method (MTM) reduce the variance of the periodogram.

Parametric methods are those in which the PSD is estimated from a signal that is assumed to be output of a linear system driven by white noise. Examples are the Yule-Walker autoregressive (AR) method and the Burg method. These methods estimate the PSD by first estimating the parameters (coefficients) of the linear system that hypothetically generates the signal. They tend to produce better results than classical nonparametric methods when the data length of the available signal is relatively short. Parametric methods also produce smoother estimates of the PSD than nonparametric methods, but are subject to error from model misspecification.

Subspace methods, also known as high-resolution methods or super-resolution methods, generate frequency component estimates for a signal based on an eigenanalysis or eigendecomposition of the autocorrelation matrix. Examples are the multiple signal classification (MUSIC) method or the eigenvector (EV) method. These methods are best suited for line spectra — that is, spectra of sinusoidal signals — and are effective in the detection of sinusoids buried in noise, especially when the signal to noise ratios are low. The subspace methods do not yield true PSD estimates: they do not preserve process power between the time and frequency domains, and the autocorrelation sequence cannot be recovered by taking the inverse Fourier transform of the frequency estimate.

All three categories of methods are listed in the table below with the corresponding toolbox function names. More information about each function is on the corresponding function reference page. See Parametric Modeling for details about lpc and other parametric estimation functions.

Spectral Estimation Methods/Functions

MethodDescriptionFunctions

Periodogram

Power spectral density estimate

periodogram

Welch

Averaged periodograms of overlapped, windowed signal sections

pwelch, cpsd, tfestimate, mscohere

Multitaper

Spectral estimate from combination of multiple orthogonal windows (or “tapers”)

pmtm

Yule-Walker AR

Autoregressive (AR) spectral estimate of a time-series from its estimated autocorrelation function

pyulear

Burg

Autoregressive (AR) spectral estimation of a time-series by minimization of linear prediction errors

pburg

Covariance

Autoregressive (AR) spectral estimation of a time-series by minimization of the forward prediction errors

pcov

Modified Covariance

Autoregressive (AR) spectral estimation of a time-series by minimization of the forward and backward prediction errors

pmcov

MUSIC

Multiple signal classification

pmusic

Eigenvector

Pseudospectrum estimate

peig