two2n.m

多元统计分析是一种应用非常广泛的数据处理方法

function [IndicesN]=Two2N(DimX,Indices2);
%TWO2N Conversion of indices between unfoldings and N-way arrays
%
%
% function [IndicesN]=Two2N(DimX,Indices2);
%
%
% This algorithm requires access to:
% ''
%
% ---------------------------------------------------------
%                Conversion of indices
%         between unfoldings and N-way arrays
% ---------------------------------------------------------
%
% [IndicesN]=Two2N(DimX,Indices2);
%
% DimX     : Dimensions of the N-way array.
% Indices2 : Indices in the unfolded 2-way array.
% indicesN : Indices in the N-way array
%
% This function helps you resolve the correct N-way indices
% to an entry in an unfolded array. If you e.g. want to know
% the N-way indices of the [3 240] element of a 5-way core
% with dimensions [7 6 5 8 7] you would write:
% two2n([7 6 5 8 7],[3 240]) and the answer is
% [3 6 5 8 1].

% Copyright (C) 1995-2006  Rasmus Bro & Claus Andersson
% Copenhagen University, DK-1958 Frederiksberg, Denmark, [e-mail]
%
% This program is free software; you can redistribute it and/or modify it under 
% the terms of the GNU General Public License as published by the Free Software 
% Foundation; either version 2 of the License, or (at your option) any later version.
%
% This program is distributed in the hope that it will be useful, but WITHOUT 
% ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS 
% FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more details.
% You should have received a copy of the GNU General Public License along with 
% this program; if not, write to the Free Software Foundation, Inc., 51 Franklin 
% Street, Fifth Floor, Boston, MA  02110-1301, USA.



% $ Version 0.01 $ Date 11. July 1997 $ Not compiled $

C=size(DimX,2);
a=Indices2(1);
b=Indices2(2);

if a>DimX(1),
   fprintf('two2n.m: The column index cannot be that high!\n');
   return
end;

if b>prod(DimX)/DimX(1),
   fprintf('two2n.m: The row index cannot be that high!\n');
   return
end;

Tb=b;
Par(1)=a;    
for c=C:-1:3,
   factor=prod(DimX(2:c-1));
   Par(c)=floor(Tb/factor)+1;
   if (Tb-(Par(c)-1)*factor)<1,
      Par(c)=Par(c)-1;
   end;
   Tb=Tb-(Par(c)-1)*factor;
end;
Par(2)=Tb;
IndicesN=Par;