zj.m

来自「MATLAB科学计算与工程分析源代码源程序3」· M 代码 · 共 34 行

s=[  1     1  1    1  857;
     1     2   2   2  951;
     1     3   3   3  909;
     2     1   2   3  878;
     2     2   3   1  973;
     2     3   1   2  899;
     3     1   3   2  803;
     3     2   1   3  1030;
     3     3   2   1   927];
 [m,n]=size(s);
 p=max(s(:,1));
 q=n-1;
 sum0=zeros(p,q);
 for i=1:q
    for k=1:m 
         for j=1:p
          if(s(k,i)==j)
              sum0(j,i)=sum0(j,i)+s(k,n);
          end
        end
     end
 end

maxdiff=max(sum0)-min(sum0);
result(1,:)=maxdiff;
maxsum0=max(sum0);
for kk=1:q
    modmax=mod(find(sum0==maxsum0(kk)),p);
    if modmax==0
        modmax=p;
    end
    result(2,kk)=(modmax);
end
result