grey_correlation.m

该函数是求出向量j和矩阵subload的间的灰色关联度

%该函数是求出向量j和矩阵subload的间的灰色关联度，并用Grey输出
function Grey=grey_correlation(j,subload)
jsubload=[j;subload];%合并两矩阵
[m,n]=size(jsubload);
p=0.5;

colu=repmat(jsubload(:,1),1,n);%取矩阵的第一列并形成与矩阵同维的第一列重复矩阵
jsubload1=jsubload./colu; %将矩阵的数据全部除以第一列的值，将数列进行无量纲化并使所有数列有公共交点

    load_minus=abs(repmat(jsubload1(1,:),m-1,1)-jsubload1(2:m,:));%第1行与其它行相减的绝对值矩阵
    ma=repmat(maxm(load_minus),m-1,n);%两级最大矩阵
    mi=repmat(min(min(load_minus)),m-1,n);%两级最小矩阵
    coe=(mi+p*ma)./(load_minus+p*ma);%形成各点的灰度关联系数矩阵,p值一般取0.5
   Grey=mean(mean(coe,2));%第一行与其他行的平均灰色关联度输出