1 function [AEF,SHIFTMAT3d] = order(A,AE,N,L);
  2 
  3 SHIFTMAT3d = zeros(N,L,L);
  4 SHIFTMAT = zeros(N,L);
  5 
  6 for k=0:L-1,
  7         for l=0:L-1,
  8                 for shift=0:N-1,
  9                         SHIFTMAT(shift+1,l+1) = sum(abs(A(:,k+1) - circshift(AE(:,l+1),shift)).^2);
 10                 end;
 11         end;
 12         SHIFTMAT3d(:,:,k+1) = SHIFTMAT;
 13 end;
 14 
 15 AEF = zeros(N,L);
 16 
 17 for k=0:L-1,
 18         [val,shiftcol] = min(SHIFTMAT3d(:,:,k+1));
 31         AEF(:,k+1) = circshift(AE(:,shiftcol(2)),shiftcol(1));
 36 end;