x(1,:)=[83 0.191 12.9 7.2 89.4 0.432 6.33];
x(2,:)=[75 0.189 11.6 9.1 82.3 0.453 5.87];
x(3,:)=[64 0.165 11.9 10.3 69.3 0.512 6.31]; %列出各数值，可修改
x(4,:)=[63 0.165 12.8 9.7 68.2 0.455 6.6];
x(5,:)=[56 0.211 13.2 12.6 77.5 0.317 7.12];
m=5;n=7;
x0=[83 0.211 13.2 7.2 89.4 0.317 5.87]; %参考序列
for i=1:n
avg(i)=0; %均值初始化
end
for i=1:m
for j=1:n
avg(j)=avg(j)+x(i,j);
end
end %求均值序列
for i=1:n
avg(i)=avg(i)/m;
end
for j=1:m
for i=1:n
x(j,i)=x(j,i)/avg(i); %均值化
end
end
for i=1:n
x0(i)=x0(i)/avg(i); %参考序列均值化end
for j=1:m
for i=1:n
delta(j,i)=abs(x(j,i)-x0(i)); %求序列差
end
end
max=delta(1,1);
for j=1:m
for i=1:n
if delta(j,i)>max
max=delta(j,i);
end
end
end
%求两极差
min=0;
for j=1:m
xgd(j)=0;
for i=1:n
glxs(j,i)=0.5*max/(0.5*max+delta(j,i)); %计算关联系数及相关度
xgd(j)=xgd(j)+glxs(j,i);
end
xgd(j)=xgd(j)/n;
end
xgd
%因此，A—E区与参考序列（最佳指标）的相关度分别为0.8489 0.6983 0.5588 0.5858 0.7105
