tomwalters@0: % method of class @signal
tomwalters@0: % 
tomwalters@0: %   INPUT VALUES:
tomwalters@0: %  
tomwalters@0: %   RETURN VALUE:
tomwalters@0: %
tomwalters@0: % 
tomwalters@0: % (c) 2003-2008, University of Cambridge, Medical Research Council 
tomwalters@0: % Maintained by Tom Walters (tcw24@cam.ac.uk), written by Stefan Bleeck (stefan@bleeck.de)
tomwalters@0: % http://www.pdn.cam.ac.uk/cnbh/aim2006
tomwalters@0: % $Date: 2008-06-10 18:00:16 +0100 (Tue, 10 Jun 2008) $
tomwalters@0: % $Revision: 585 $
tomwalters@0: 
tomwalters@0: function w=getcenterofmass(sig)
tomwalters@0: 
tomwalters@0: vals=getvalues(sig);
tomwalters@0: nr=getnrpoints(sig);
tomwalters@0: 
tomwalters@0: su1=0;
tomwalters@0: su2=0;
tomwalters@0: 
tomwalters@0: for i=1:nr
tomwalters@0:    su1=su1+vals(i)*i;
tomwalters@0:    su2=su2+vals(i);
tomwalters@0: end
tomwalters@0: 
tomwalters@0: if su2~=0
tomwalters@0: 	w=su1/su2;
tomwalters@0: else
tomwalters@0: 	w=nr/2;
tomwalters@0: end
tomwalters@0: