From: <cd...@us...> - 2012-08-07 06:08:51
|
Revision: 10834 http://octave.svn.sourceforge.net/octave/?rev=10834&view=rev Author: cdf Date: 2012-08-07 06:08:45 +0000 (Tue, 07 Aug 2012) Log Message: ----------- beautify code syntax in lombnormcoeff Modified Paths: -------------- trunk/octave-forge/extra/lssa/inst/lombnormcoeff.m Modified: trunk/octave-forge/extra/lssa/inst/lombnormcoeff.m =================================================================== --- trunk/octave-forge/extra/lssa/inst/lombnormcoeff.m 2012-08-07 06:03:45 UTC (rev 10833) +++ trunk/octave-forge/extra/lssa/inst/lombnormcoeff.m 2012-08-07 06:08:45 UTC (rev 10834) @@ -23,11 +23,13 @@ ## @end deftypefn -function coeff = lombnormcoeff(T,X,omega) -tau = atan2( sum( sin( 2.*omega.*T)), sum(cos(2.*omega.*T))) / 2; -coeff = ( ( sum ( X .* cos( omega .* T - tau ) ) .^ 2 ./ sum ( cos ( omega .* T - tau ) .^ 2 ) - + sum ( X .* sin ( omega .* T - tau ) ) .^ 2 / sum ( sin ( omega .* T - tau ) .^ 2 ) ) - / ( 2 * var(X) ) ); +function coeff = lombnormcoeff (T, X, omega) + +tau = atan2 (sum (sin (2 .* omega .*T)), sum (cos (2 .* omega .* T))) / 2; +coeff = ((sum (X .* cos (omega .* T - tau)) .^ 2 / sum (cos (omega .* T - tau) .^ 2) + + sum (X .* sin (omega .* T - tau)) .^ 2 / sum (sin (omega .* T - tau) .^ 2 )) / + (2 * var(X))); + endfunction %!shared t, x, o, maxfreq This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site. |