灰色关联分析好的

1、 灰色关联分析用途：考虑到影响*因素的指标个数之多，并且彼此之间存在着一定的相关性，因此上海市就业是一个多因素复杂的系统，我们采用灰色关联理论对各因素与城镇就业人数之间的关系进行分析研究。灰色关联分析反映了曲线间的关联程度，反映了各相关因素对体统特征行为的接近次序，其中关联度最大的为最优因素，因此灰色关联分析对于一个系统发展变化态势提供了量化的度量，非常适合动态历程分析。1）建立原始数列的因变量参考数列和自变量比较数列因变量参考数列又称母序列，记作 x ，= ，自变, x(k)xx , x , x(k )(k)(1)0(2)0(3)000量参考数列又称母序列，记作 x ，x= = 。, x ,

2、i 1,2, nx , x , x(2)(k )(k)(1)(3)(k)iiiii2）将原始序列进行初始化法、均值化法的无法量纲处理，目的是消除数量级大小不同的影响，以便于进行计算和比较分析，我们采用了这两种方法对数据进行了处理。3）计算每个时刻点上母序列与各子序列差的绝对值，并从中取得最大差和最小差序列：d(k) x=-x (i 1,2, ,n)=(k)，则差序列为：(k)i0i()i 1,2, ,nd = d (1),d (2),d (3), d (k) ， =iiiii（21）其 中 ， 最 大 差 ： d max = max max x - x； 最 小 差 ：(k)(k)0iiid

3、min = min min x - x 。(k)(k)0iii4）计算灰色关联度系数ld min+ d max利用公式 l =计算灰色关联度系数，其中 l 是第 k 个点的(k)(k )ld + d max子因素与母因素的相对差值， 为分辨系数，一般在 0 与 1 之间选取，通常取为0i0il0.5。5）计算灰色关联度为求总的关联度，需要考虑不同的观测点在总体观测中的重要性程度，则需要确定各点的权重，我们采用算数平均的方法计算灰色关联度 r ，公式为0i1nr =l(k)0in0ik=16）关联度排序根据 r 的大小安排关联序的先后顺序，关联度越接近于 1，说明关联程度0i越大，根据经验，当

4、= 0.5时，两因素的关联度大于 0.6，便认为其关联性显著l。我们利用 matlab 软件编程（代码详见附录 3）计算可得，各指标的关1314联度如小表所示，从结果来看，所有的因素关联度均大于 0.6，则说明选取的指标均对上海的就业有影响。 表 6 各指标灰色关联度对照表xxxxxxxxx123456789关联度 0.81 0.79 0.93 0.75 0.71 0.75 0.71 0.89 0.73难点：样本数据和方阵数据的含义及其使用方法（未解决）定义：对于两个系统之间的因素，其随时间或不同对象而变化的关联性大小的量度，称为关联度。在系统发展过程中，若两个因素变化的趋势具有一致性，即同步

5、变化程度较高，即可谓二者关联程度较高；反之，则较低。因此，灰色关联分析方法，是根据因素之间发展趋势的相似或相异程度，亦即“灰色关联度”，作为衡量因素间关联程度的一种方法。就业人数进出口 城市差额 化（万亿） 水平27.60 0.5934.72 0.6140.26 0.6134.64 0.6133.17 0.6231.56 0.6227.74 0.6317.30 0.6419.60 0.6523.24 0.6619.65 0.6622.16 0.6732.11 0.6734.36 0.6833.53 0.6819.32 0.691992 806.91 0.01 1114.32 0.03 5.39

6、 4273 55.90 357.381993 853.11 0.01 1519.23 0.02 5.54 5650 53.10 653.911994 850.04 0.01 1990.86 0.02 5.66 7401 53.50 1123.29 22.87 0.701995 855.72 0.01 2499.43 0.02 5.81 9279 53.40 1601.79 41.29 0.711996 851.21 0.01 2957.55 0.02 5.90 10663 50.70 1952.05 42.13 0.711997 847.25 0.01 3438.79 0.02 5.97 11

7、425 51.70 1977.59 46.84 0.72 1998 836.21 0.01 3801.09 0.02 6.03 12059 50.60 1964.83 66.10 0.731999 812.09 0.01 4188.73 0.02 6.10 14147 45.20 1856.72 -10.34 0.742000 828.35 0.01 4771.17 0.02 6.18 15420 44.50 1869.67 -40.02 0.752001 752.26 0.02 5210.12 0.01 6.26 17764 43.40 1994.73 -56.42 0.752002 792

8、.04 0.02 5741.03 0.01 6.37 19473 39.40 2187.06 -85.54 0.762003 813.05 0.02 6694.23 0.01 6.42 22160 37.20 2452.11 -154.33 0.782004 836.87 0.02 8072.83 0.01 6.54 24398 36.37 3084.66 -129.86 0.812005 863.32 0.02 9247.66 0.01 6.62 26823 35.90 3542.55 -48.81 0.842006 885.51 0.03 10572.24 0.01 6.74 29569

9、35.60 3925.09 -3.43 0.862007 1024.33 0.03 12494.01 0.01 6.86 34707 35.50 4458.61 48.83 0.872008 1053.24 0.03 14069.87 0.01 6.95 39502 36.60 4829.45 165.62 0.872009 1064.42 0.03 15046.45 0.01 7.02 42789 35.00 5273.33 60.97 0.88clear;clc;yangben=load(yangben.txt);fangzhen=load(fangzhen.txt); %待判数据rows

10、,cols=size(fangzhen);p=0.5;%分辨系数m,n=size(yangben);yangben=(yangben./max(yangben);for i=1:cols% 因素3/9越小越好if (i=9 | i=3)clear value;clear minvalue;value = fangzhen(:,i);minvalue=min( value );for j=1:rowsfangzhen(j,i) = minvalue / fangzhen(j,i);endelseif(i=6)clear value;clear minvalue;clear maxvalue;va

11、lue = fangzhen(:,i);minvalue=min( value );maxvalue=max( value );for j=1:rowsfangzhen(j,i) = ( fangzhen(j,i) - minvalue ) / ( maxvalue - minvalue );endelseclear value;clear minvalue; value = fangzhen(:,i);maxvalue=max( value );for j=1:rowsfangzhen(j,i) = fangzhen(j,i) / maxvalue;endendendfangzhen=fan

12、gzhen;% 方阵的列行数据发生值发生交换for i=1:colsfor j=1:rowsfangzhen(i,j)=abs( yangben(1,j)-fangzhen(i,j) );endendmatrixmax=max(max(fangzhen);matrixmin=min(min(fangzhen);totalvalue=matrixmin + matrixmax;for i=1:colsfor j=1:rowsfangzhen(i,j)=totalvalue/(fangzhen(i,j)+matrixmax);endendr = sum(fangzhen,2)/rows或灰色关联度

13、 matlab 源程序最近几天一直在写算法，其实网上可以下到这些算法的源程序的，但是为了搞懂，搞清楚，还是自己一个一个的看了，写了，作为自身的积累，而且自己的的矩阵计算类库也迅速得到补充，以后关于算法方面，基本的矩阵运算不用再重复写了，挺好的，是种积累，下面把灰关联的 matlab 程序与大家分享。灰色关联度分析法是将研究对象及影响因素的因子值视为一条线上的点,与待识别对象及影响因素的因子值所绘制的曲线进行比较,比较它们之间的贴近度,并分别量化,计算出研究对象与待识别对象各影响因素之间的贴近程度的关联度,通过比较各关联度的大小来判断待识别对象对研究对象的影响程度。关联度计算的预处理，一般初值化

14、或者均值化，根据我的实际需要，本程序中使用的是比较序列与参考序列组成的矩阵除以参考序列的列均值等到的，当然也可以是其他方法。 %注意：由于需要，均值化方法采用各组值除以样本的各列平均值clear;clc;yangben=47.924375 25.168125 827.4105438 330.08875 1045.164375 261.37437516.3372 6.62 940.2824 709.2752 962.1284 84.87455.69666667 30.80333333 885.21 275.8066667 1052.42 435.81; %样本数据fangzhen=36.27 1

15、4.59 836.15 420.41 1011.83 189.5464.73 35.63 755.45 331.32 978.5 257.8742.44 23.07 846 348.05 1025.4 296.6959.34 39.7 794.31 334.63 1016.4 317.2752.91 17.14 821.79 306.92 1141.94 122.044.21 4.86 1815.52 2584.68 963.61 0.006.01 2.43 1791.61 2338.17 1278.08 30.873.01 1.58 1220.54 956.14 1244.75 3.9125.65 7.42 790.17 328.88 1026.01 92.82115.80 27 926.5 350.93 1079.49 544.3812.63 8.75 1055.50 1379.00 875.10 1.65; %待判数据fpxer mro)=n 0or,;d w.n s5i; =,;r cso wio =zl 1srdydlfe (=:yaion yt=jdr osfmar ;nw=a_oi% szrgxnw=分beiyg2=s1 e(jyzhum:辨nh;a1imfer )a;