Mercurial > hg > aimmat

diff aim-mat/tools/@frame/getweightedintervalsum.m @ 0:74dedb26614d
Initial checkin of AIM-MAT version 1.5 (6.4.2011).
author: tomwalters
date: Fri, 20 May 2011 12:32:31 +0100
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/aim-mat/tools/@frame/getweightedintervalsum.m	Fri May 20 12:32:31 2011 +0100
@@ -0,0 +1,39 @@
+% method of class @frame
+% 
+%   INPUT VALUES:
+%  
+%   RETURN VALUE:
+%
+% 
+% (c) 2003, University of Cambridge, Medical Research Council 
+% Stefan Bleeck (stefan@bleeck.de)
+% http://www.mrc-cbu.cam.ac.uk/cnbh/aimmanual
+% $Date: 2003/01/17 16:57:46 $
+% $Revision: 1.3 $
+
+function sumsig=getweightedintervalsum(cframe,gauss)
+% usage: spectal_activity=getweightedspektralsum(cframe,gauss)
+% returns the sum over frequency as function of interval
+% the frame is weighted by "gauss"
+
+
+number=length(cframe);
+if number==1
+    val=cframe.values;
+    sval=val';
+    gava=getvalues(gauss);
+    for i=1:getnrpoints(cframe)
+        sval(i,:)=sval(i,:).*gava';
+    end
+    intervals=sum(sval');
+
+    sumsig=signal(intervals);
+    sumsig=setsr(sumsig,getsr(cframe));
+    sumsig=setname(sumsig,sprintf('Weighted sum of Frame: %s',getname(cframe)));
+    sumsig=setstarttime(sumsig,getminimumtime(cframe));
+    
+    
+else
+        not implemented yet
+end
+