%%------將數(shù)據(jù)均一化-------%%
a=[2.6	2.3	2	2.4	2.8
2.5	2	3	2.5	2.6
2	2.1	1.6	1.5	2.1
3.2	3	3.2	3.1	3.5
3.6	4	3.2	3.8	4
4.5	4.3	4	4.3	4.1
18.4	17.7	17	17.6	19.1

];%輸入原始矩陣
                       %notes:在原始矩陣中，首行為總的值，
                       %均值的計(jì)算是列均值，要計(jì)算同一指標(biāo)在不同情況下的指標(biāo)值
a=a';
[aR,aL]=size(a);%原始矩陣的行列
aM=mean(a);%原始矩陣的列均值
b=ones(aR,aL);
for j=1:aL
    b(:,j)=a(:,j)/aM(j);
            %for k=1:aR
             %   b(k,i)=a(k,i)/aM(j);
            end
       
%for  i=1:aL
   % for j=1:aR
%disp('均一化后的矩陣:');
b=b';
%b

%%--------計(jì)算關(guān)聯(lián)系數(shù)----------%%
[bR,bL]=size(b);%計(jì)算矩陣b的大小
c=b(bR,:);%首行為c矩陣
d=b(1:bR-1,:);%b除去c后剩下的就是d矩陣
cL=length(c);%c的大小
[dR,dL]=size(d);%d的大小
rou=0.5;
for i=1:cL
    for j=1:dL
        if(i==j)
            for k=1:dR
               temp(k,j)=abs(c(i)-d(k,j));
            end
        end
    end
end
kesai=(min(min(temp))+rou*max(max(temp)))./(temp+max(max(temp)));%計(jì)算關(guān)聯(lián)度
corre=mean(kesai');%計(jì)算出關(guān)聯(lián)系數(shù)
disp(corre');