gusucode.com > ​常见的三种灰关联分析matlab源码程序 > GRA2.m

    %计算两向量绝对灰关联度 方式2
%完成时间:2014.5.27

function [gabd] = GRA2(x,y)
[r,c] = size(x);
len = max(r,c);
if(size(x) ~= size(y))
    error('错误:序列长度不等');
    quit;
end

%始点零化象
sx = x - x(1);
sy = y - y(1);

%具体公式
xs = sum(sx) - sx(1) - sx(len);
ys = sum(sy) - sy(1) - sy(len);

xx = abs(xs + 1/2 * sx(1));
yy = abs(ys + 1/2 * sy(1));
zz = abs(xx - yy);

gabd = (1 + xx + yy) / (1 + xx + yy + zz);