%該函數(shù)是求出向量j和矩陣subload的間的灰色關(guān)聯(lián)度，并用Grey輸出
function Grey=grey_correlation(j,subload)
jsubload=[j;subload];%合并兩矩陣
[m,n]=size(jsubload);
p=0.5;

colu=repmat(jsubload(:,1),1,n);%取矩陣的第一列并形成與矩陣同維的第一列重復(fù)矩陣
jsubload1=jsubload./colu; %將矩陣的數(shù)據(jù)全部除以第一列的值，將數(shù)列進(jìn)行無量綱化并使所有數(shù)列有公共交點(diǎn)

    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);%形成各點(diǎn)的灰度關(guān)聯(lián)系數(shù)矩陣,p值一般取0.5
   Grey=mean(mean(coe,2));%第一行與其他行的平均灰色關(guān)聯(lián)度輸出