gcm.asv

灰色管理法的实现

%%------将数据均一化-------%%
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:在原始矩阵中，首行为总的值，
                       %均值的计算是列均值，要计算同一指标在不同情况下的指标值
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

%%--------计算关联系数----------%%
[bR,bL]=size(b);%计算矩阵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)));%计算关联度
corre=mean(kesai');%计算出关联系数
disp(corre');