求大神解答 代码如下
m=6; %子序列
k=6;%序列长度
X0=zeros(6,6);%差序列矩阵
R0=zeros(6,6);%关联系数矩阵
a0=[4788,5383,6421,7166.19,7583.36,7987]; %母序列
a1=[715,676,677,679,682,690];
a2=[20464,22368,24998,28567,32145,35227];
a3=[3805652,4594829,6024387,7827193,9301099,10485229];
a4=[672,814,1020,925,1212,957];
a5=[1398062,1858234,2563024,2438708,3185557,2889915]; %子序列
x0=a0./a0(1);
x1=a1./a1(1);
x2=a2./a2(1);
x3=a3./a3(1);
x4=a4./a4(1);
x5=a5./a5(1);
x6=a6./a6(1);
X=[x1',x2',x3',x4',x5',x6'];
for i=1:m
for ii=1:k
X0(ii,i)=abs(x0(ii)-X(ii,i));
end
end
Max=max(max(X0));
Min=min(min(X0));
for i=1:m
for ii=1:k
R0(ii,i)=(Min+0.5*Max)/(0.5*Max+X0(ii,i));
end
end
r=sum(R0)/k;
disp('关联度依次为:')
disp(r)
disp(R0)
disp(Max)
disp(Min)
m=6; %子序列
k=6;%序列长度
X0=zeros(6,6);%差序列矩阵
R0=zeros(6,6);%关联系数矩阵
a0=[4788,5383,6421,7166.19,7583.36,7987]; %母序列
a1=[715,676,677,679,682,690];
a2=[20464,22368,24998,28567,32145,35227];
a3=[3805652,4594829,6024387,7827193,9301099,10485229];
a4=[672,814,1020,925,1212,957];
a5=[1398062,1858234,2563024,2438708,3185557,2889915]; %子序列
x0=a0./a0(1);
x1=a1./a1(1);
x2=a2./a2(1);
x3=a3./a3(1);
x4=a4./a4(1);
x5=a5./a5(1);
x6=a6./a6(1);
X=[x1',x2',x3',x4',x5',x6'];
for i=1:m
for ii=1:k
X0(ii,i)=abs(x0(ii)-X(ii,i));
end
end
Max=max(max(X0));
Min=min(min(X0));
for i=1:m
for ii=1:k
R0(ii,i)=(Min+0.5*Max)/(0.5*Max+X0(ii,i));
end
end
r=sum(R0)/k;
disp('关联度依次为:')
disp(r)
disp(R0)
disp(Max)
disp(Min)