Function[result,sum0]=zjjc(s,opt)
[m,n]=size(s)
p=max(s(:,1)) *运行时 p 为小数
q=n-1
sum0=zeros(p,q) *要求 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
if(opt= =1)
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
else
minsum0=min(sum0)
for kk=1:q bn
modmin=mod(find(sum0= =minsum0(kk)),p)
if modmin= =0
modmin=p
end
result(2,kk)=(modmin) *左右数组维度不一样
end
end
运行不了
arr=reshape(A,prod(size(A)),1);
pd=max(arr)-min(arr);%极差
sk=skewness(arr);%偏度
ku=kurtosis(arr);%峰度
