邓氏关联度及其代码

%邓氏关联度,灰色系统及其应用第七版73页
%读取数据,第一行为参考序列，其余行为比较序列
data=xlsread("2.xls");
r = size(data,1);
c = size(data,2);
%第一步：求解各序列的初值像
data1=data(:,:)./data(:,1);
%第二步，求X0与Xi的初值像对应分量之差的绝对值序列
deta=abs(data1(1,:)-data1(2:end,:));
%第三步，最大值与最小值
m = min(min(deta)); %全局最小值
M = max(max(deta)); %全局最大值
%计算关联系数
r0=0.5;%r0的取值范围属于0到1
gamma=(m+r0*M)./(deta+r0*M);
%第五步，最后求出关联系数的平均值既是所要的关联度
degree=sum(gamma,2)./c;