五種灰色關(guān)聯(lián)度分析matlab代碼

1、精選優(yōu)質(zhì)文檔-傾情為你奉上灰色鄧氏關(guān)聯(lián)度分析% P12 - The Study on the Grey Relational Degree and Its Application function r1 = gld_deng(x)s = size(x);len = s(2);num = s(1);ro = 0.5;for i = 1: num x(i,:) = x(i,:)./x(i,1);enddx(num,len) = 0;for i = 2 : num for k = 1 : len dx(i,k) = abs(x(1,k) - x(i,k); endendmax_dx = max(ma

2、x(dx);min_dx = min(min(dx);r(1,1:len-1) = 1;for i = 2 : num for k = 1 : len r(i,k) = (min_dx + ro*max_dx)/(dx(i,k) + ro*max_dx); endendr1 = sum(r(2:num,:),2)/(len);改進(jìn)灰色絕對(duì)關(guān)聯(lián)度分析% P11 - The Study on the Grey Relational Degree and Its Application function r1 = gld_gjjd(x)s = size(x);len = s(2);num = s(1

3、);for i = 1: num x(i,:) = x(i,:)./x(i,1);enddx(num,len-1) = 0;for i = 1 : num for j = 1 : len - 1 dx(i,j) = x(i,j+1) - x(i,j); endendc = 1;beta(1,1:len-1) = 0;w(1,1:len-1) = 0;for i = 2 : num temp = sum(abs(x(i,:) - x(1,:),2); for k = 1 : len - 1 beta(i,k) = atan(dx(i,k) - dx(1,k)/(1 + dx(i,k)*dx(1,

4、k); if beta(i,k) < 0 beta(i,k) = pi + beta(i,k); end w(i,k) = 1 - abs(x(i,k) - x(1,k)/temp; endendr = c./(c + tan(beta./2);wr = w.*r;r1 = sum(wr(2:num,:),2)/(len - 1);灰色絕對(duì)關(guān)聯(lián)度分析% P18 - The Study on the Grey Relational Degree and Its Application function r1 = gld_jd(x)s = size(x);len = s(2);num = s

5、(1);for i = 1: num x(i,:) = x(i,:)./x(i,1);enddx(num,len-1) = 0;for i = 1 : num for j = 1 : len - 1 dx(i,j) = x(i,j+1) - x(i,j); endendr(1,1:len-1) = 1;for i = 2 : num for k = 1 : len - 1 r(i,k) = 1/(1 + abs(dx(1,k) - dx(i,k); endendr1 = sum(r(2:num,:),2)/(len - 1);灰色T型關(guān)聯(lián)度分析% P19 - The Study on the

6、Grey Relational Degree and Its Application function r1 = gld_t(x)s = size(x);len = s(2);num = s(1);dx(num,len-1) = 0;for i = 1 : num for j = 1 : len - 1 dx(i,j) = abs(x(i,j+1) - x(i,j); end d_x = sum(dx(i,:),2)/(len - 1); x(i,:) = x(i,:)./d_x;enddx(num,len-1) = 0;for i = 1 : num for j = 1 : len - 1

7、dx(i,j) = x(i,j+1) - x(i,j); endendr(1,1:len-1) = 1;for i = 2 : num for k = 1 : len - 1 if dx(1,k)*dx(i,k) = 0 r(i,k) = sign(dx(1,k)*dx(i,k); else r(i,k) = sign(dx(1,k)*dx(i,k)*min(abs(dx(1,k),abs(dx(i,k) / max(abs(dx(1,k),abs(dx(i,k); end endendr1 = sum(r(2:num,:),2)/(len - 1);灰色斜率關(guān)聯(lián)度分析% P20 - The Study on the Grey Relational Degree and Its Application function r1 = gld_xl(x)s = size(x);len = s(2);num = s(1);for i = 1: num x(i,:) = x(i,:)./x(i,1);enddx(num,len-1) = 0;for i = 1 : num for j = 1 : len - 1 dx(i,j) = x(i,j+1) - x(i,j); endendr(1,1:len-1) = 1;for i = 2 : num