function [y, factor] = normalize(x)
%
% Normalize to make the sum of magnitude in a burst 
% to be one in spatial domain.
%   [m n] = size(x);   leng=length(x);   if m ~= leng      temp=m;      m = n;      n = temp ;      x = x.';   end   factor=sum(x);   sumx=ones(m,1)*factor;   y=x./sumx;