Index of /~gari/CODE/Lomb

      Name                    Last modified       Size  Description


[DIR] Parent Directory        19-May-2009 15:42      -  
[   ] LSP.m                   10-Aug-2008 12:29     2k  
[   ] Liv_LOMB_demo_1.m       10-Aug-2008 12:29     6k  
[   ] fLSPw.m                 10-Aug-2008 12:29     4k  
[   ] lomb.zip                10-Aug-2008 12:37   108k  
[   ] testLSP.mat             10-Aug-2008 12:29   100k  
[   ] vLSP.m                  10-Aug-2008 12:29     3k  

The code in this directory, contributed by Tony Fisher of Liverpool 
University in the UK, provides various versions of the Lomb-Scargle
Periodogram and described below by Tony in his own words ....

--- Files provided:
README                - This file
LSP.m                 - Original Lomb-Scargle Periogoram
vLSP.m                - Vectorized LSP
Liv_LOMB_demo_1.m     - Demo Script - Run this
fLSPw.m               - Numerical Recipes in C fasper.c
testLSP.mat           - Test data
lomb.zip              - Zip file of all the above code & data
----

vLSP.m (vectorised + stats) and fLSPw.m (Press & Rybiki + stats)

I have done a function script Liv_LOMB_demo.m which runs these over 4 HRV datasets (in testLSP.mat).
Note:
 
Just vectorising (even very carefully as in vLSP.m) is not much good other than in function loops (with 'persistent' working variables) and for small data sets ... the intermediate matrices grow very quickly and soon overwhelm the computer. On the other hand, the Press & Rybiki (fLSPw.m) is terrific: it can do v big sets and is v quick.
 
The output plots and the command line returns show the story. I have not returned all the stats as they are rather comprehensive  .. for details just use help vLSP and help fLSPw from the command line and you get to see all the comment code (as normal).
 
Summary  ... run Liv_LOMB_demo from the command line: note figures are modal and wait for you to kill them to move on.