Lomb normalized periodogram

버전 1.2.0.0 (6.57 KB) 작성자: Christos Saragiotis
Both functions calculate the Lomb-Scargle periodogram (aka Gauss-Vanicek/Least-squares spectrum)
다운로드 수: 3.7K
업데이트 날짜: 2008/12/4

라이선스 보기

Both functions caculate the Lomb normalized periodogram (aka Lomb-Scargle, Gauss-Vanicek or Least-Squares spectrum) of a vector x with coordinates in t, which is essentially a generalization of the DFT for unevenly sampled data.

The codes are transcriptions from Fortran of the subroutines found in Section 13.8 (pp. 569-577) of "Numerical recipes in Fortran 77: the art of scientific computing", 2nd ed., vol. 1, Cambridge University Press, NY, USA, 2001 by WH Press, SA Teukolsky, WT Vetterling and BP Flannery,

However, Matlab's characteristics have been taken into account in order to make it fast for Matlab.

FASTLOMB is much faster than LOMB (especially when the length of the input increases) but even LOMB is faster than any other implementation I found in FileExchange. Also they both do not suffer from memory problems (I tested them both for inputs of 100,000 samples).

I'd also like to acknowledge file ID: 20004 (for some reason I can't get two file IDs in the acknowledgements)

인용 양식

Christos Saragiotis (2024). Lomb normalized periodogram (https://www.mathworks.com/matlabcentral/fileexchange/22215-lomb-normalized-periodogram), MATLAB Central File Exchange. 검색됨 .

MATLAB 릴리스 호환 정보
개발 환경: R2008a
모든 릴리스와 호환
플랫폼 호환성
Windows macOS Linux
카테고리
Help CenterMATLAB Answers에서 Fortran with MATLAB에 대해 자세히 알아보기
도움

받음: lombscargle.m

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!
버전 게시됨 릴리스 정보
1.2.0.0

I corrected a formula where a sqrt was supposed to exist and it didn't. The results don't differ much though.

Sorry about this...

1.1.0.0

After running the Matlab profiler, I realized that the waitbar I had in the lomb.m function was the main factor of slowness of the lomb.m function. So I removed the waitbar and updated the info in the "Other requirements" section.

1.0.0.0