关于灰色关联度的程序与算例 doc

最近在写论文，当利用到灰色关联度的时候发现网上其实很多，但是新手们往往猴急，急于求成，所以估计也没心思静下心看程序。我自己收集了一个M文件，但测试的时候发现不是我想要的那样子，因此在利用Excel手动验算一遍后，发现改动后的程序与Excel计算的一样，因此想将一点体会写出来，以供后来的师弟师妹们参考。不过，我得首先感谢那些在论坛里之前出现的零碎的代码供给着，他们的思路也许是对后来者的一种有价值的学习机会。

例子：比如，计算北部湾与南宁、北海、防城港、钦州、玉林、崇左的关联度大小。

一般我推荐用综合关联度分析，因为它比相对关联度分析更全面。

First,利用Excel计算一遍，相关过程参考中国知网或百度文库里面关于相对关联度的计算方法，这里暂不列出。

对各序列进行始点零化像

绿色部分是北部湾与其他城市之间的关联度大小（即上述原始数列矩阵），其中0代表北部湾，字母i依次代表南宁、北海、防城港、钦州、玉林、崇左城市名称。



用Matlab计算也是一样的结果，只是上述Excel算的其实计算量也不是很大，但如果出现巨大数量的矩阵的时候，估计你只好想到编程求解了，以下提供综合关联度的计算程序。将以下文本拷贝成M文件，即可在软件中调用，

function com_xy = greyrelation(x,y) %调用的时候不用function英文前缀！！！我几次都弄错了，看来是Matlab基础薄弱啊~~~

y_row = size(y,1);%;%计算矩阵y的行数

y_col =size(y,2);%;%计算矩阵y的列数

x_col = size(x,2);%;%计算x的列数

if y_col ~= x_col

error(message('MATLAB:greyrelation:wrong in input data'));

end

temp_y = y;%绝对关联度中比较序列中的数据处理后的矩阵

temp_x = x;%x数据处理后的矩阵

for i =1:x_col

temp = x(i)-x(1);

temp_x(i)=temp;

end

for i =1:y_row

for j=1:y_col

temp = y(i,j) - y(i,1);

temp_y(i,j)=temp;

end

end%处理过程

%temp_x;

%temp_y;

s0 = abs(sum(temp_x)-0.5*temp_x(x_col));

abs_xy =[];

for i=1:y_row

si = abs(sum(temp_y(i,:))-0.5*temp_y(i,y_col));

si_s0 = abs(si-s0);

abs_xy(i,1) =(1+s0+si)/(1+s0+si+si_s0);

end%计算绝对关联度%该步骤在编程过程中存在的问题是公式理解上存在偏差，公式没理解好

%下面开始计算相关关联度

temp_y2 = y;

temp_x2 = x;

for i =1:x_col

temp = x(i)/x(1);

temp_x2(i)=temp-1;

end

for i =1:y_row

for j=1:y_col

temp = y(i,j) / y(i,1);

temp_y2(i,j)=temp-1;

end

end

s02 = abs(sum(temp_x2)-0.5*temp_x2(x_col));

rela_xy=[];

for i=1:y_row

si2 = abs(sum(temp_y2(i,:))-0.5*temp_y2(i,y_col));

si2_s02 = abs(si2-s02);

rela_xy(i,1) =(1+s02+si2)/(1+s02+si2+si2_s02);

end

%上述步骤是为了计算相关关联度，

%下面计算综合关联度

com_xy = 0.5*abs_xy +(1-0.5)*rela_xy;%返回的是综合关联度

%关于如何返回相对关联度的方法是：将函数名改为function rela_xy = greyrelation(x,y),

%然后在最后一行写上rela_xy，如此返回的才是相对关联度



以上是M程序，通过程序运行，得到的结果是（相对关联度）

rela_xy=

是不是和Excel算的一样呢（当然，保留4位小数即可）

com_xy= PS:这个是综合关联度，至于做论文或分析时，自己看情况选取吧。

不知道你们看懂没哦，给出的程序主要求解com_xy，但如果有同学想求解rela_xy，即相对关联度，将程序第一行改为function rela_xy = greyrelation(x,y) ,同时最后一行的com_xy = 0.5*abs_xy +(1-0.5)*rela_xy删除，替换为rela_xy即可。就是说首尾呼应。



以上是程序和例子说明，我只提供灰色关联度分析的程序和例子，自己今天总结的，仅供参考！不足之处，还请程序高手批评指正！

5/18/2020 1:44:22 AM

本文来源：https://www.2haoxitong.net/k/doc/5cd77463ef630b1c59eef8c75fbfc77da26997da.html