灰色關聯度Matlab代碼

load x.txt %把原始數據存放在純文本文件x.txt中，其中把數據的"替換替換成.
for i=1:40
x(i,:)=x(i,:)/x(i,1); %標准化數據
end
data=x;
n=size(data,2); %求矩陣的列數，即觀測時刻的個數
ck=data(1,:); %提出參考數列
bj=data(2:end,:); %提出比較數列
m2=size(bj,1); %求比較數列的個數
    for j=1:m2
    t(j,:)=bj(j,:)-ck;
    end
mn=min(min(abs(t'))); %求最小差
mx=max(max(abs(t'))); %求最大差
rho=0.5; %分辨系數設置
ksi=(mn+rho*mx)./(abs(t)+rho*mx); %求關聯系數
r=sum(ksi')/n %求關聯度
[rs,rind]=sort(r,'descend') %對關聯度進行排序