Periodogram

Periodogram

The periodogram is an estimate of the spectral density of a signal. The term was coined by Arthur Schuster in 1898Schuster, A., "On the investigation of hidden periodicities with application to a supposed 26 day period of meteorological phenomena," "Terrestrial Magnetism and Atmospheric Electricity", 3, 13-41, 1898.] as in the following quote: [cite web | url = http://books.google.com/books?id=v-EOAAAAIAAJ&q=schuster+%22convenient+to+have+a+word+for+some+representation+of+a+variable+quantity%22&dq=schuster+%22convenient+to+have+a+word+for+some+representation+of+a+variable+quantity%22&ei=4OwaR7WSOI7wpwKd5Zm4CQ&pgis=1 | title = The Annals of Mathematical Statistics | year = 1972 ]

cquote|THE PERIODOGRAM. 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:Let:::frac{T}{2}a = int_{t_1}^{t_1+T}f(t)cos(kt)dt:::frac{T}{2}b = int_{t_1}^{t_1+T}f(t)sin(kt)dtwhere "T" may for convenience be chosen to be equal to some integer multiple of :::frac{2pi}{k},and plot a curve with {2pi}/k as abscissæ and :::r = sqrt{a^2+b^2}as ordinates; this curve, or, better, the space between this curve and the axis of abscissæ, represents the periodogram of f(t).

Note that the term "periodogram" may also be used to describe the quantity r^2,Box, George and Jenkins, Gwilym (1970) "Time series analysis: Forecasting and control", San Francisco: Holden-Day.] which is its common meaning in astronomy (as in "the modulus-squared of the discrete Fourier transform of the time series (with the appropriate normalisation)"Vaughan and Uttley (2006), "Advances in Space Research", Volume 38, Issue 7, p. 1405-1408] ). See Scargle (1982) for a detailed discussion in this context.Scargle (1982), "Astrophysical Journal", Part 1, vol. 263, Dec. 15, 1982, p. 835-853]

A spectral plot refers to a smoothed version of the periodogram. [ [http://www.itl.nist.gov/div898/handbook/eda/section3/eda33r.htm Spectral Plot] , from the NIST Engineering Statistics Handbook.] [ [http://www.itl.nist.gov/div898/software/dataplot/refman1/ch2/spectral.pdf Short explanation] of the relation between the spectral plot and the periodogram.] Smoothing is performed to reduce the effect of measurement noise.

In practice, the periodogram is often computed from a finite-length digital sequence using the fast Fourier transform (FFT). The raw periodogram is not a good spectral estimate because of spectral bias and the fact that the variance at a given frequency does not decrease as the number of samples used in the computation increases.

The spectral bias problem arises from a sharp truncation of the sequence, and can be reduced by first multiplying the finite sequence by a window function which truncates the sequence gracefully rather than abruptly.

The variance problem can be reduced by smoothing the periodogram. Various techniques to reduce spectral bias and variance are the subject of spectral estimation.

One such technique to solve the variance problems is also known as the method of averaged periodogramsEngelberg, S. (2008), "Digital Signal Processing: An Experimental Approach", Springer, Chap. 7 p. 56] . The idea behind it is, to divide the set of samples into N sets of M samples, compute the DFT of each set, square it to get the power spectral density and compute the average of all of them. This leads to a decrease in the standard deviation as frac{1}{sqrt{N .

See also

*Fourier transform
*Discrete Fourier transform
*Least-squares spectral analysis
*Spectral density
*Spherical harmonics
*Welch method

References


Wikimedia Foundation. 2010.

Игры ⚽ Поможем написать курсовую

Look at other dictionaries:

  • periodogram — noun A graph used in the harmonic analysis of spectra, time series etc …   Wiktionary

  • periodogram — pe·ri·od·ogram …   English syllables

  • periodogram — ˌpirēˈädəˌgram, ˈpirēədə noun Etymology: International Scientific Vocabulary period (I) + o + gram : a curve exhibiting graphically the periodicity of any natural or physical phenomenon …   Useful english dictionary

  • Bayesian Kepler periodogram — is an automatic instrument installed in Markov chain Monte Carlo observatory in British Columbia, Canada. It is designed to find more planets in single or multiple planetary systems. P.C. Gregory operates this instrument to find evidence for any… …   Wikipedia

  • Least-squares spectral analysis — (LSSA) is a method of estimating a frequency spectrum, based on a least squares fit of sinusoids to data samples, similar to Fourier analysis. [cite book | title = Variable Stars As Essential Astrophysical Tools | author = Cafer Ibanoglu |… …   Wikipedia

  • Welch's method — In physics, engineering, and applied mathematics, Welch s method, named after P.D. Welch, is used for estimating the power of a signal vs. frequency, reducing noise compared to the methods it is based upon. Welch s method is based on the concept… …   Wikipedia

  • Bartlett's method — In physics, engineering, and applied mathematics, Bartlett s method, named after M.S. Bartlett, is used for estimating power spectra. Bartlett s method provides a way to reduce the variance of the periodogram in exchange for a reduction of… …   Wikipedia

  • Doppler spectroscopy — Diagram showing how a smaller object (such as an extrasolar planet) orbiting a larger object (such as a star) could produce changes in position and velocity of the latter as they orbit their common center of mass (red cross). Doppler spectroscopy …   Wikipedia

  • HD 11964 d — Planetbox begin name = HD 11964dPlanetbox star star = HD 11964 constell = Cetus ra = RA|1|57|09.6064 dec = DEC| 10|14|32.739 dist ly = 110.83 dist pc = 33.98 class = G5IVPlanetbox orbit semimajor = 1.13±0.04 eccentricity = 0.63+0.34 0.17 period …   Wikipedia

  • HD 208487 c — Planetbox begin name = HD 208487 c Planetbox star star = HD 208487 constell = Grus RA=RA|21|57|19.8477 DEC=DEC| 37|45|49.037 dist ly=144.4 dist pc=44.30 class=G2V: Planetbox orbit semimajor = 1.87±|0.13|0.14cite journal |… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”