squaregrey.m

该函数是求出N*P矩阵的各行间的1减去平均灰色关联度对称矩阵

%该函数是求出N*P矩阵的各行间的1减去平均灰色关联度对称矩阵，并用矩阵Z输出
function Z=squaregrey(load)
[m,n]=size(load);
p=0.5;
d=zeros(m,m);
f=zeros(m,n);
colu=repmat(load(:,1),1,n);%取矩阵的第一列并形成与矩阵同维的第一列重复矩阵
load1=load./colu; %将矩阵的数据全部除以第一列的值，将数列进行无量纲化并使所有数列有公共交点
for t=1:m
    load_minus=abs(repmat(load1(t,:),m,1)-load1);%第t行与其它行相减的绝对值矩阵
    ma=repmat(maxm(load_minus),m,n);
    mi=repmat(min(min(load_minus)),m,n);
    coe=(mi+p*ma)./(load_minus+p*ma);%形成各点的灰度关联系数矩阵,p值一般取0.5
   for u=1:m
    f(t,u)=mean(coe(u,:));%形成行与行之间的关联度矩阵
    end
end
for t=1:m
    for u=1:m
        d(t,u)=mean([f(t,u),f(u,t)]);%形成均值的对称矩阵
    end
end
Z=1-d;