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

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

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

%初值象
x = x ./ x(1);
y = y ./ y(1);

%始点零化象
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(len));
yy = abs(ys + 1/2 * sy(len));
zz = abs(xx - yy);

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