coredian.m

来自「多维数据处理：MATLAB源程序用于处理多维数据」· M 代码 · 共 23 行

function [DIA_Rel,DIA_Abs]=CoreDian(C,Fac)
%[DIA_Rel,DIA_Abs]=CoreDian(C,Fac)
%
%DIA_Rel : Diagonality in percent
%DIA_Abs : Absolute sum of squares of
%          the diagonal elements 
%
% Claus Andersson

c=size(Fac,2);
d=ones(1,c);
l=min(Fac);

C=C.^2;

DIA_Abs=0;
for j=1:l,
  [i,j]=GetIndxn(Fac,j*d);
  DIA_Abs = DIA_Abs + C(i,j);
end;

DIA_Rel=100*DIA_Abs/sum(sum(C));