灰色系统理论

一.关联分析

1.消除量纲差异

包括初值化变换（除x1），均值化变换，百分比变换（除xmax），倍数变换（除xmin），归一化变换

注意，成本型指标是效益型指标的倒数

 2.代码

clc,clear

load x.txt %把原始数据存放在纯文本文件 x.txt 中

for i=1:N1

x(i,:)=x(i,:)/x(i,1); %标准化数据

end

for i=N1+1:N2

x(i,:)=x(i,1)./x(i,:); %标准化数据

end

data=x;

n=size(data,1);

ck=data(1,:);m1=size(ck,1);  %取出第1行数据，对第一行（成绩）进行因素分析

bj=data(2:n,:);m2=size(bj,1);

for i=1:m1

for j=1:m2

t(j,:)=bj(j,:)-ck(i,:);

end

jc1=min(min(abs(t')));jc2=max(max(abs(t')));

rho=0.5;

ksi=(jc1+rho*jc2)./(abs(t)+rho*jc2);

rt=sum(ksi')/size(ksi,2);  %各个时刻的关联度取一个平均值

r(i,:)=rt;

end

r

[rs,rind]=sort(r,'descend')

说明：得到的r越大，则r对于总成绩的贡献度越大

            size(A,1)该语句返回的时矩阵A的行数， c=size(A,2) 该语句返回的时矩阵A的列数