G( k' Y3 h+ b' i2 l. `9 z9 \% 数据预处理,每一个元素除所在列的均值 {) O1 m y! Q9 o1 F
Mean = mean(gdp); y8 a( _) M& w0 R3 F; {2 l
gdp = gdp ./ repmat(Mean,size(gdp,1),1);/ |5 j- t3 Y# t2 F! U& z% q
, x- t, a$ Q- t$ F' p4 k3 z
% 确定子母序列 " W( v2 H Y$ m; aY = gdp(:,1);) m: I( C* L e1 f: l
X = gdp(:,2:end);/ M& y8 n8 W5 Y! M. m, G0 k
) \. P6 b7 K% u; N* e! |1 u% 计算两级最小差,和最大差3 }) W& o8 F. V$ g x
val_min = 100;3 E) q+ J9 n, V0 z- W6 p
val_max = 0;- Z- r# J: q5 }8 r) B0 u
for i = 1:31 p0 L: u) M" `/ x. _6 `6 G
val_min = min(val_min,min(abs(X(:,i)-Y)));6 ~! E& `( ^# Z4 A/ ?* |$ L
val_max = max(val_max,max(abs(X(:,i)-Y))); ! W* C8 q5 _' p5 Pend 1 `7 V+ a/ ]; R% Q Z9 M8 D6 [8 ?$ K* w8 B. o; W
% 计算灰度关联系数 5 U M3 z0 n4 ]% m& c9 Z) jp = 0.5; % 分辨系数取0.5 1 d' l3 Z9 z% J, n1 i& wabsX0_Xi = abs(X - repmat(Y,1,size(X,2))); % 计算|X0-Xi|矩阵& \- F$ C2 Z8 L) r
g = (val_min+p*val_max) ./ (absX0_Xi+p*val_max); % 各指标与母序列关联系数5 m2 m y3 T2 h. j- x
1 K2 W4 b9 s% H: \! |2 N) g) U2 ?4 u% 计算灰色关联度. }- [. p1 Z2 O I4 F0 Q
res = mean(g) 5 s+ t1 v& c8 t" T/ z该模型也可以用来确定权重 % L% f+ u4 r$ k! l) Z4 A, M S2 |7 } ]& l" Q4 K! N* v