Periodogram

From HandWiki
Revision as of 16:01, 30 June 2023 by NBrush (talk | contribs) (fix)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Short description: Estimate of the spectral density of a signal

In signal processing, a periodogram is an estimate of the spectral density of a signal. The term was coined by Arthur Schuster in 1898.[1] Today, the periodogram is a component of more sophisticated methods (see spectral estimation). It is the most common tool for examining the amplitude vs frequency characteristics of FIR filters and window functions. FFT spectrum analyzers are also implemented as a time-sequence of periodograms.

Definition

There are at least two different definitions in use today.[2] One of them involves time-averaging,[3] and one does not.[4] Time-averaging is also the purview of other articles (Bartlett's method and Welch's method). This article is not about time-averaging. The definition of interest here is that the power spectral density of a continuous function, [math]\displaystyle{ x(t), }[/math]  is the Fourier transform of its auto-correlation function (see Cross-correlation theorem, Spectral density, and Wiener–Khinchin theorem):

[math]\displaystyle{ \mathcal{F}\{x(t)\circledast x^*(-t)\} = X(f)\cdot X^*(f) = \left| X(f) \right|^2. }[/math]

Computation

A power spectrum (magnitude-squared) of two sinusoidal basis functions, calculated by the periodogram method.
Two power spectra (magnitude-squared) (rectangular and Hamming window functions plus background noise), calculated by the periodogram method.

For sufficiently small values of parameter T, an arbitrarily-accurate approximation for X(f) can be observed in the region  [math]\displaystyle{ -\tfrac{1}{2T} \lt f \lt \tfrac{1}{2T} }[/math]  of the function:

[math]\displaystyle{ X_{1/T}(f)\ \triangleq \sum_{k=-\infty}^{\infty} X\left(f - k/T\right), }[/math]

which is precisely determined by the samples x(nT) that span the non-zero duration of x(t)  (see Discrete-time Fourier transform).

And for sufficiently large values of parameter N,  [math]\displaystyle{ X_{1/T}(f) }[/math] can be evaluated at an arbitrarily close frequency by a summation of the form:

[math]\displaystyle{ X_{1/T}\left(\tfrac{k}{NT}\right) = \sum_{n=-\infty}^\infty \underbrace{T\cdot x(nT)}_{x[n]}\cdot e^{-i 2\pi \frac{kn}{N}}, }[/math]

where k is an integer. The periodicity of  [math]\displaystyle{ e^{-i 2\pi \frac{kn}{N}} }[/math]  allows this to be written very simply in terms of a Discrete Fourier transform:

[math]\displaystyle{ X_{1/T}\left(\tfrac{k}{NT}\right) = \underbrace{\sum_{n} x_{_N}[n]\cdot e^{-i 2\pi \frac{kn}{N}},}_{DFT}\quad \scriptstyle{\text{(sum over any }n\text{-sequence of length }N)}, }[/math]

where [math]\displaystyle{ x_{_N} }[/math] is a periodic summation:  [math]\displaystyle{ x_{_N}[n]\ \triangleq \sum_{m=-\infty}^{\infty} x[n-mN]. }[/math]

When evaluated for all integers, k, between 0 and N-1, the array:

[math]\displaystyle{ S\left(\tfrac{k}{NT}\right) = \left| \sum_{n} x_{_N}[n]\cdot e^{-i 2\pi \frac{kn}{N}} \right|^2 }[/math]

is a periodogram.[4][5][6]

Applications

Periodogram for Proxima Centauri b is shown at the bottom.[7]

When a periodogram is used to examine the detailed characteristics of an FIR filter or window function, the parameter N is chosen to be several multiples of the non-zero duration of the x[n] sequence, which is called zero-padding (see § Sampling the DTFT).[upper-alpha 1]  When it is used to implement a filter bank, N is several sub-multiples of the non-zero duration of the x[n] sequence (see § Sampling the DTFT).

One of the periodogram's deficiencies is that the variance at a given frequency does not decrease as the number of samples used in the computation increases. It does not provide the averaging needed to analyze noiselike signals or even sinusoids at low signal-to-noise ratios. Window functions and filter impulse responses are noiseless, but many other signals require more sophisticated methods of spectral estimation. Two of the alternatives use periodograms as part of the process:

  • The method of averaged periodograms,[8]  more commonly known as Welch's method,[9][10]  divides a long x[n] sequence into multiple shorter, and possibly overlapping, subsequences. It computes a windowed periodogram of each one, and computes an array average, i.e. an array where each element is an average of the corresponding elements of all the periodograms. For stationary processes, this reduces the noise variance of each element by approximately a factor equal to the reciprocal of the number of periodograms.
  • Smoothing is an averaging technique in frequency, instead of time. The smoothed periodogram is sometimes referred to as a spectral plot.[11][12]

Periodogram-based techniques introduce small biases that are unacceptable in some applications. Other techniques that do not rely on periodograms are presented in the spectral density estimation article.

See also

Notes

  1. N is designated NFFT in the Matlab and Octave applications.

References

  1. Schuster, Arthur (January 1898). "On the investigation of hidden periodicities with application to a supposed 26 day period of meteorological phenomena". Terrestrial Magnetism 3 (1): 13–41. doi:10.1029/TM003i001p00013. Bibcode1898TeMag...3...13S. http://iranithenticate.ir/download/On%20the%20investigation%20of%20hidden%20periodicities%20with%20application%20to%20a%20supposed%2026%20day%20period%20of%20meteorological%20phenomena/iranithenticate-ir.pdf. "It is convenient to have a word for some representation of a variable quantity which shall correspond to the ‘spectrum’ of a luminous radiation. I propose the word periodogram, and define it more particularly in the following way.". 
  2. McSweeney, Laura A. (2004-05-14). "Comparison of periodogram tests". Journal of Statistical Computation and Simulation (online ($50)) 76 (4): 357–369. doi:10.1080/10629360500107618. 
  3. "Periodogram—Wolfram Language Documentation". https://reference.wolfram.com/language/ref/Periodogram.html. 
  4. 4.0 4.1 "Periodogram power spectral density estimate - MATLAB periodogram". https://www.mathworks.com/help/signal/ref/periodogram.html. 
  5. Oppenheim, Alan V.; Schafer, Ronald W.; Buck, John R. (1999). Discrete-time signal processing (2nd ed.). Upper Saddle River, N.J.: Prentice Hall. p. 732 (10.55). ISBN 0-13-754920-2. https://archive.org/details/discretetimesign00alan. 
  6. Rabiner, Lawrence R.; Gold, Bernard (1975). "6.18". Theory and application of digital signal processing. Englewood Cliffs, N.J.: Prentice-Hall. pp. 415. ISBN 0-13-914101-4. https://archive.org/details/theoryapplicatio00rabi/page/415. 
  7. "Do-it-yourself Science — is Proxima c hiding in this graph?". https://www.eso.org/public/images/potw1737a/. 
  8. Engelberg, S. (2008), Digital Signal Processing: An Experimental Approach, Springer, Chap. 7 p. 56
  9. Welch, Peter D. (June 1967). "The use of Fast Fourier Transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms". IEEE Transactions on Audio and Electroacoustics AU-15 (2): 70–73. doi:10.1109/TAU.1967.1161901. Bibcode1967ITAE...15...70W. 
  10. "Welch's power spectral density estimate - MATLAB pwelch". https://www.mathworks.com/help/signal/ref/pwelch.html. 
  11. Spectral Plot, from the NIST Engineering Statistics Handbook.
  12. "DATAPLOT Reference Manual". National Institute of Standards and Technology (NIST). 1997-03-11. http://www.itl.nist.gov/div898/software/dataplot/refman1/ch2/spectral.pdf. "The spectral plot is essentially a “smoothed” periodogram where the smoothing is done in the frequency domain." 

Further reading

  • Box, George E. P.; Jenkins, Gwilym M. (1976). Time series analysis: Forecasting and control. San Francisco: Holden-Day. 
  • Scargle, J.D. (December 15, 1982). "Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data". Astrophysical Journal, Part 1 263: 835–853. doi:10.1086/160554. Bibcode1982ApJ...263..835S. 
  • Vaughan, Simon; Uttley, Philip (2006). "Detecting X-ray QPOs in active galaxies". Advances in Space Research 38 (7): 1405–1408. doi:10.1016/j.asr.2005.02.064. Bibcode2006AdSpR..38.1405V.