Next: LWCAT(1) Up: WFDB Applications Guide Previous: LOG10(1)On This Page

Name

lomb - estimate power spectrum using the Lomb periodogram method

Synopsis

lomb [ options ... ] input-file

Description

lomb transforms a real-valued time series (from the specified input-file, or from the standard input if input-file is specified as "-"; input-file must be in text form) into a power spectrum (on the standard output), using a technique known as the Lomb periodogram.

The input is a text file containing a sampled time series, presented as two columns of numbers (the sample times and the sample values). The intervals between consecutive samples need not be uniform (in fact, this is the most significant advantage of the Lomb periodogram over other methods of power spectral density estimation). lomb writes the Lomb periodogram (the power spectral density estimate derived from the input time series) on the standard output, in two columns (frequency and power). If the units of the sample times in the input file are seconds, the units of the frequencies in the output are Hz.

Options are:

-h
Print a usage summary.
-P
Generate a power spectrum (print squared magnitudes).
-s
Smooth the output.
-z
Add a constant to each input sample, chosen such that the mean value of the entire series is zero.

Among many other applications, this program can be used to estimate heart rate power spectra, in combination with ihr(1) . The Lomb method is ideal for analysis of any time series with missing or noisy data (the noisy data may be removed from the time series and need not be replaced, as would be necessary if conventional PSD estimation algorithms were employed).

See Also

fft(1) , hrfft(1) , memse(1)
Lomb, N.R.
Least-squares frequency analysis of unequally spaced data. Astrophysics and Space Science 39:447-462 (1976).
Press, W.H, and Rybicki, G.B.
Fast algorithm for spectral analysis of unevenly sampled data. Astrophysical J. 338:277-280 (1989).
Press, W.H. Teukolsky, S.A., Vetterling, W.T., and Flannery, B.P.
Numerical Recipes in C: the Art of Scientific Computing, pp. 575-584 (Cambridge Univ. Press, 1992).
Moody, G.B.
Spectral analysis of heart rate without resampling. Computers in Cardiology 1993, pp. 715-718 (IEEE Computer Society Press, 1993). http://www.physionet.org/physiotools/lomb/lomb.html .

Author

George B. Moody (george@mit.edu)

Source

http://www.physionet.org/physiotools/wfdb/psd/lomb.c


Table of Contents

Up: WFDB Applications Guide


Please e-mail your comments and suggestions to webmaster@physionet.org, or post them to:

PhysioNet
MIT Room E25-505A
77 Massachusetts Avenue
Cambridge, MA 02139 USA

Updated 28 May 2015