灰色关联分析法与TOPSIS评价法

灰色关联度评价法一、灰色关联分析(GRA)方法

灰色关联分析是一种多因素统计分析方法,它是以各因素的样本数据为依据用灰色联度来描述因素间关系的强弱、大小和次序的.如果样本数据列反映出两因素变化的态势(方向、大小、速度等)基本一致,则它们之间的关联度较大；反之关联度较小.

与其他传统的多因素分析方法相比,灰色关联分析对数据要求较低且计算量小,便于广泛应用.GRA分析的核心是计算关联度,下面通过一个例子来说明计算关联度的思路和方法.

下表是某地区2000-2005年国内生产总值的统计资料.现在提出这样的问题：该地区产业中,哪一产业的变化与该地区国内生产总值(GDP)的变化态势更一致？也就是哪一产业与GDP的关联度最大呢？年份国内生产总值第一产业第二产业第三产业200019883868397632001206140884680820022335422960953200327504821258101020043356511157712682005380656118931352表

某地区国内生产总值统计资料(百万元)

一个自然的想法就是分别将三次产业产值的时间序列与GDP的时间序列进行比较,为了能够比较,先对各序列进行无量纲化,这里采用均值化法.各序列的均值分别为：2716,461.5,1228.83,1025.67,上表中每列数据除以其均值可得均值化序列(如表所示)年份tGDP

x0(t)一产业x1(t)二产业x2(t)三产业x3(t)20000.73200.83640.68280.744020010.75880.88190.68850.787820020.85970.91440.78120.929120031.01251.04441.02370.984720041.23561.10731.28331.236320051.40131.21561.54051.3182两序列变化的态势是表现在其对应点的间距上.如果各对应点间距均较小,则两序列变化态势的一致性强,否则,一致性弱.分别计算各产业产值与GDP在对应期的间距(绝对差值),结果见表所示.年份t20000.10440.04920.011920010.12310.07040.028920020.05470.07850.069420030.03190.01120.027820040.12840.04770.000620050.18570.13920.0832接下来应该是三个绝对值序列分别求平均再进行比较,就可以解决问题了.但仔细观察表中的数据会发现绝对差值数据序列的数据间存在着较大的数量级差异(最大为0.1857,最小的为0.0006,相差300多倍),不能直接进行综合,还需要对其进行一次规范化.设和分别表示表中绝对值的最大数和最小数，则因而由于在一般情况下,可能为零(即某个为零)故将上式改进为显然越大,说明两序列(xi和x0)的变化态势一致性弱,反之,一致性强,因此可考虑将取倒反向,为了规范化后数据在[0,1]内,可考虑在0和1之间取值.上式可变形为称为序列xi和序列x0在第t期的灰色关联系数(或简称为关联系数).由(6.1)式可以看出,取值的大小可以控制对数据转化的影响,取较小的值,可以提高关联系数间差异的显著性,因而称为分辨系数.利用(6.1)对表6-3中绝对差值进行规范化,取

结果见表6-4,以

计算为例：同样可计算出表6-4中其余关联系数.年份t20000.41910.60670.868720010.37960.51780.725720020.58080.49030.521320030.70550.87610.733820040.36960.61411.00020050.28810.35100.4758表6-4最后分别对各产业与GDP的关联系数序列求算术平均可得r0i称为序列x0和xi(i=1,2,3)的灰色关联度.由于因而第三产业产值与GDP的关联度最大,其次是第二产业、第一产业.可以看出,灰色关联分析需要经过以下几个步骤：1.确定分析序列在对研究问题定性分析的基础上,确定一个因变量因素和多个自变量因素.设因变量数据构成参考序列,各自变量数据构成比较序列n+1个数据序列成成如下矩阵：其中N为变量序列的长度.一般情况下,原始变量序列具有不同的量纲或数量级,为了保证分析结果的可靠性,需要对变量序列进行无量纲化.无量纲化后各因素序列形成如下矩阵：2.对变量序列进行无量纲化计算(6.3)中第一列(参考序列)与其余各列(比较序列)对应期的绝对差值,形成如下绝对差值矩阵：3.求差序列、最大差和最小差其中绝对差值阵中最大数和最小数即为最大差和最小差：对绝对差值阵中数据作如下变换：4.计算关联系数得到关联系数矩阵：式中分辨系数在(0,1)内取值,一般情况下依据(6.10)中数据情况多在0.1至0.5取值,越小越能提高关联系数间的差异.关联系数是不超过1的正数,越小,越大,它反映第i个比较序列Xi与参考序列X0在第k个期关联程度.比较序列Xi与参考序列X0的关联程度是通过N个关联系数(即(6.10)中第i列)来反映的,求平均就可得到Xi与X0的关联度5.计算关联度6.依关联度排序对各比较序列与参考序列的关联度从大到小排序，关联度越大，说明比较序列与参考序列变化的态势越一致．灰色关联度综合评价法

利用灰色关联分析进行综合评价的步骤是：1．根据评价目的确定评价指标体系，收集评价数据。2．确定参考数据列参考数据列应该是一个理想的比较标准，可以以各指标的最优值（或最劣值）构成参考数据列，也可根据评价目的选择其它参照值．记作3．对指标数据进行无量纲化

无量纲化后的数据序列形成如下矩阵：

4．逐个计算每个被评价对象指标序列（比较序列）与参考序列对应元素的绝对差值

即（,为被评价对象的个数）．5．确定

与6．计算关联系数由（12－5）式，分别计算每个比较序列与参考序列对应元素的关联系数．7．计算关联序对各评价对象（比较序列）分别计算其个指标与参考序列对应元素的关联系数的均值，以反映各评价对象与参考序列的关联关系，并称其为关联序，记为：8．如果各指标在综合评价中所起的作用不同，可对关联系数求加权平均值即

9．依据各观察对象的关联序，得出综合评价结果．2.灰色关联分析的应用举例

例1：利用灰色关联分析对6位教师工作状况进行综合评价1．评价指标包括：专业素质、外语水平、教学工作量、科研成果、论文、著作与出勤．2．对原始数据经处理后得到以下数值，见下表

编号专业外语教学量科研论文著作出勤1898752927875738397966474688843658669838689576483．确定参考数据列：

4．计算，见下表编号专业外语教学量科研论文著作出勤1101237022124161302032524311146351330061610422515．求最值6．依据（12－5）式，取计算，得 同理得出其它各值，见下表编号10.7781.0000.7780.6360.4670.3331.00020.636

0.778

0.636

0.467

0.636

0.368

0.778

31.000

0.636

1.000

0.538

0.538

0.412

0.636

40.538

0.778

0.778

0.778

0.412

0.368

0.538

50.778

0.538

0.538

1.000

0.778

0.368

0.778

60.778

1.000

0.467

0.636

0.538

0.412

0.778

7．分别计算每个人各指标关联系数的均值（关联序）：

8．如果不考虑各指标权重（认为各指标同等重要），六个被评价对象由好到劣依次为1号，5号，3号，6号，2号，4号．即

Topsis法(对多个对象作评价)TOPSIS（Techniquefororderpreferencebysimilaritytoidealsolution）法，即逼近理想解排序法，意为与理想方案相似性的顺序选优技术，是系统工程中有限方案多目标决策分析的一种常用方法。它是基于归一化后的原始数据矩阵，找出有限方案中最优方案和最劣方案（分别用最优向量和最劣向量表示），然后分别计算诸评价对象与最优方案和最劣方案的距离，获得各评价对象与最优方案的相对接近程度，以此作为评价优劣的依据。步骤：1.设有n个评价对象、m个评价指标，原始数据可写为矩阵X＝(Xij)n×m

2.对极大、极小指标分别进行一致化、归一化变换

3.归一化得到矩阵Z＝(Zij)n×m，其各列最大、最小值构成的最优、最劣向量分别记为Z＋＝(Zmax1

Zmax2…

Zmaxm)

Z－＝(Zmin1

Zmin2…

Zminm)

4.第i个评价对象与最优、最劣方案的距离分别为5.第i个评价对象与最优方案的接近程度Ci为例某儿童医院1994～1998年7项指标的实际值，用Topsis法比较该医院这5年的医疗质量年份出院人数病床使用率平均住院日病死率抢救成功率治愈好转率院内感染率19942158476.77.31.0178.397.52.019952437286.37.40.8091.198.02.019962204181.87.30.6291.197.33.219972111584.56.90.6090.297.72.919982463390.36.90.2595.597.93.6变换后，得到矩阵

平均住院日、病死率、院内感染率为极小指标，其余为极大指标，一致化、归一化变换计算各列最大、最小值构成的最优、最劣向量分别为Z＋＝(0.48330.48050.46340.81780.47760.44870.5612)Z－＝(0.41420.40810.43210.20240.39160.44550.3118)计算各年与最优、最劣向量的距离（以94年为例）C