svdize

0001 function svdwindows = svdize(data,meanSubtract)
0002 
0003 [r,c] = size(data);
0004 svdwindows = cell(r,c);
0005 for ri = 1:r
0006   for ci = 1:c
0007     windows = data{ri,ci};
0008     nwin = length(windows);
0009     svdwindowsTrial = cell(1,nwin);
0010     for w = 1:nwin
0011       if meanSubtract
0012     windows{w} = windows{w} - ...
0013         repmat(mean(windows{w},2),1,size(windows{w},2));
0014       end
0015       [u,s] = svd(windows{w});
0016       %[u,s]=sorteig(windows{w} * windows{w}');
0017 
0018       for col = 1:size(u,2)
0019     if u(1,col) < 0
0020       u(:,col) = u(:,col) * -1;
0021     end
0022       end
0023       svdwindowsTrial{w} = u;
0024     end
0025     svdwindows{ri,ci} = svdwindowsTrial;
0026   end
0027 end