MATLAB File Help: lomb | View code for lomb | WFDB Contents |

lomb

function varargout=lomb(varargin) [Pxx,F]=lomb(x,dcOffset,smooth) Wrapper to WFDB LOMB: http://www.physionet.org/physiotools/wag/lomb-1.htm Transforms a real-valued time series 'x' into a power spectrum 'X', using a technique known as the Lomb periodogram. The input is a Nx2 matrix containing a sampled time series, presented as two columns of numbers (the sample times in the first column and the sample values in the second). The intervals between consecutive samples need not be uniform. Input Parameters: x Nx2 vector of doubles. First column is sample time index (in seconds), and second column is the sample value of the signal at that time. dcOffset (Optional) Booelan. If present add constant to input samples ( x(:,2) ), such that the mean values of the time series is zero (default=1). smooth (Optional) Boolean String specifying the if the output should be smoothed (default =1). Output Parameters: Pxx Mx1 Double. Estimated power spectrum. F Mx1 Double. Frequency of the estimated power spectrum (Hz). CITING CREDIT: To credit this function, please cite the following paper at your work: 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 . Additional References: 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 . MATLAB wrapper written by Ikaro Silva, 2013 Last Modified: - Version 1.0 Since 0.9.0 %Example: Heart Rate Spectral Analysis: [signal,Fs,tm]=rdsamp('mitdb/100',1); [ann]=rdann('mitdb/100','atr'); [Pxx,F]=lomb([tm(ann) signal(ann)]); plot(F,Pxx);grid on;hold on See also RDANN, TACH, SQRS, WQRS