[Octave-cvsupdate] SF.net SVN: octave:[10836] trunk/octave-forge/extra/lssa/inst/lombcoeff.m

[<cd...@us...>]

Revision: 10836
          http://octave.svn.sourceforge.net/octave/?rev=10836&view=rev
Author:   cdf
Date:     2012-08-07 06:39:02 +0000 (Tue, 07 Aug 2012)
Log Message:
-----------
more syntax beautification

Modified Paths:
--------------
    trunk/octave-forge/extra/lssa/inst/lombcoeff.m

Modified: trunk/octave-forge/extra/lssa/inst/lombcoeff.m
===================================================================
--- trunk/octave-forge/extra/lssa/inst/lombcoeff.m	2012-08-07 06:27:57 UTC (rev 10835)
+++ trunk/octave-forge/extra/lssa/inst/lombcoeff.m	2012-08-07 06:39:02 UTC (rev 10836)
@@ -25,14 +25,16 @@
 
 function coeff = lombcoeff (T, X, o)
 
-  theta = atan2 (sum (sin (2 .* o .* T )), 
-                 sum (cos (2 .* o .* T))) / (2 * o);
+  oT = o .* T;
 
-  coeff = (sum (X .* cos (o .* T - theta)) ^2 / 
-           sum (cos (o .* T - theta) .^2) + 
-           sum (X .* sin (o .* T - theta)) ^2 / 
-           sum (sin (o .* T - theta) .^2));
+  theta = atan2 (sum (sin (2 * oT)), 
+                 sum (cos (2 * oT))) ./ (2 * o);
 
+  coeff = (sum (X .* cos (oT - theta)) ^2 / 
+           sum (cos (oT - theta) .^2) + 
+           sum (X .* sin (oT - theta)) ^2 / 
+           sum (sin (oT - theta) .^2));
+
 endfunction
 
 

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.