Revision: 10832
http://octave.svn.sourceforge.net/octave/?rev=10832&view=rev
Author: cdf
Date: 2012-08-07 05:55:09 +0000 (Tue, 07 Aug 2012)
Log Message:
-----------
fix test syntax
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-05 23:05:36 UTC (rev 10831)
+++ trunk/octave-forge/extra/lssa/inst/lombnormcoeff.m 2012-08-07 05:55:09 UTC (rev 10832)
@@ -22,20 +22,19 @@
##
## @end deftypefn
-%!test
-%!shared t, x, o, maxfreq;
-%! maxfreq = 4 / ( 2 * pi );
-%! t = linspace(0,8);
-%! x = ( 2.*sin(maxfreq.*t) + 3.*sin((3/4)*maxfreq.*t) - 0.5 .* sin((1/4)*maxfreq.*t) - 0.2 .* cos(maxfreq .* t) + cos((1/4)*maxfreq.*t));
-%! o = [ maxfreq , 3 / 4 * maxfreq , 1 / 4 * maxfreq ];
-%!assert( lombnormcoeff(t,x,o(1)),63.3294946603949,5e-10 );
-%!assert( lombnormcoeff(t,x,o(2)),73.2601360674868,5e-10 );
-%!assert( lombnormcoeff(t,x,o(3)),53.0799752083903,5e-10 );
-
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
+%! t = linspace(0,8);
+%! maxfreq = 4 / ( 2 * pi );
+%! o = [ maxfreq , 3 / 4 * maxfreq , 1 / 4 * maxfreq ];
+%! x = ( 2.*sin(maxfreq.*t) + 3.*sin((3/4)*maxfreq.*t) - 0.5 .* sin((1/4)*maxfreq.*t) - 0.2 .* cos(maxfreq .* t) + cos((1/4)*maxfreq.*t));
+%!assert( lombnormcoeff(t,x,o(1)),63.3294946603949,5e-10 );
+%!assert( lombnormcoeff(t,x,o(2)),73.2601360674868,5e-10 );
+%!assert( lombnormcoeff(t,x,o(3)),53.0799752083903,5e-10 );
This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.
|
