关于地区经济发展水平的综合评价方法

关于地区经济发展水平的综合评价方法

摘要

本文以多元统计分析中主成分分析的方法为基础，对我国16个省份经济发展水平的7个主要指标进行分析评价，为了处理不确定环境下的综合评价问题，进一步明确各项指标的权重系数，提出了一种基于蒙特卡罗仿真思想的随机模拟综合评价方法, 将点赋值的序关系法推广至区间赋值的情形。基于随机模拟的思想, 依据经典的综合评价过程, 模拟反映被评价对象之间优劣关系的优胜度矩阵，在此基础上, 给出各被评价对象之间的最佳排序。通过实证分析，得到该种方法比主成分分析更为有效。

关键词：主成分分析，随机模拟，综合评价

Abstract

The paper is based on the principal component analysis which is a kind of the multivariate statistics analysis methods. It aims at studying and comparing with the 16 provinces’economic development level with the seven major indicators. But for solving the uncertain problem ,we want to make sure The weight of each indicator coefficient further, so a comprehensive evaluation method based on Monte Carlo simulation which is therefore presented, where the dot assignment can be extended to the interval assignment. Based on the stochastic simulation and according to the classical comprehensive evaluation process, the superiority matrix reflecting the relation between the superiority and inferiority of the object s being evaluated is simulated. Then, the optimal ranking of all the objects to be evaluated is given. Finally, using the case analysis, we draw a conclusion that stochastic simulation is more efficient.

Key Words：principal component analysis, stochastic simulation, comprehensive evaluation

目录

摘要................................................... I Abstract.................................................. I 第1节问题背景 .. (1)

第2节主成分分析 (1)

第3节综合评价的相关理论及随机模拟的实现 (2)

3.1基本概念与方法原理 (2)

3.2随机模拟的软件实现步骤 (4)

3.3综合评价值的确定 (5)

第4节实证分析与研究 (6)

4.1数据来源 (6)

4.2指标选择 (6)

4.3综合评价排序 (6)

第5节程序输入及输出结果 (8)

略 (8)

第6节总结 (8)

参考文献 (9)

附录一 (9)

第1节 问题背景

综合评价即是对被评价对象完成时态的状态进行客观、公正、合理的全面评价, 其理论与方法在社会科学与自然科学的多个领域有着广泛的应用。指标权重系数的确定一直是综合评价理论研究中的一个核心问题。文献3提出了一种简捷有效的赋权方法—序关系分析法。它的主要特点有: 不用构造判断矩阵,更无需一致性检验; 方法简便、直观, 便于应用; 对元素或指标的个数没有限制; 具有保序性。序关系分析法一经提出, 就在综合评价理论中得到了广泛的应用。

随着综合评价问题研究的深入和推广, 不确定环境下的综合评价问题引起了人们的关注, 并已取得了丰硕的成果。由于客观事物的复杂性、不确定性以及人类思维的模糊性, 在使用序关系分析法时, 专家在给出任意两个指标间的重要性程度之比时, 往往很难给出精准的数值。从实际应用出发, 专家会更偏爱给出一个相对宽松的比值区间, 这样操作的优点是通过降低评价的精确度以提高评价的效率, 并能从整体上提高评价的效果。但同时,不确定序关系分析过程也使得问题求解的难度大幅增加针对这种情况, 本文给出了一种基于随机模拟的综合评价方法,为了说明这种新方法的优劣，本文考虑利用多元统计知识中的主成分分析法对其进行验证。其理论是将多个指标化为少数综合性指标,这些综合指标是原指标的线性组合，且彼此之间互不相关。再根据综合指标进行评价。它可以在力保原始数据信息丢失最少的情况下，对高维变量空间进行降维处理，是一种较新的综合评价方法。

第2节 主成分分析

主成分分析的基本思想]1[

对于较多变量(设为p )的多元统计问题，由于彼此间存在多重共线性，因而使得观测数据有一定成都的信息重叠。希望用较少的几个（设为p m <个）综合变量在（主成分）来代替原来较多的变量，使得几个综合变量彼此互不相关，且尽可能多的反映原变量的信息。p 个变量的n 次观测可以简化为m 个主成分的n 次观测，即主成分得分。主成分分析可归纳为如下几步： （1）设原始变量为n x x x ,,,21 ，主成分分析后得到的综合变量为

)(,,,21n m z z z m < ，为了排除数量级和量纲不同带来的影响，首先对原始数据进行标准化处理．p j n i x x x i i ij ij ,,2,1;,,2,1,/)(* ==-=σ。其中ij x 为第i 个指标第j 个分区的原始数据；i x 和i σ分别为样本均值与标准差． （2）由标准化数据表，计算相关系数矩阵n n ij r R ?=)(，其中

j i j kj n

k i ki ij x x x x n r σσ/)()(11

--=∑=。

（3）计算R 的特征值和特征向量．将特征值从大到小排列n λλλ≥≥≥ 21，同时可得对应的单位特征向量n u u u ,,,21 。

（4）计算贡献率∑==n

i i i i e 1

/λλ和累积贡献率∑∑===n

i i j m j j E 1

1

/λλ。

（5）计算主成分∑∑===

p j n

i ij ij j x u z 11

*。 （6）综合分析．一般取%85>j E 的最小)(n m m <，得m 个主成分，对m 个主

成分进行综合分析．

第3节 综合评价的相关理论及随机模拟的实现

3.1 基本概念与方法原理]2[

设由n 个被评价对象(或方案) n o o o ,,21，m 个指标(属性)m x x x ,,,21 组成的多指标评价系统：),,2,1;,,2,1)((m j n i o x x i j ij ===为被评价对象

i o 关于指标j x 的观测值。评价数据矩阵(决策矩阵)可表示为：

?

?

???

?

??????==?nm n n m m n

n ij x x x x x x x x x x A 2

1

2222111211][

记{}{})3,(,,2,1,,,2,1≥==n m n N m M 不妨设m 个指标均为极大型(效益型)的,A 中数据为预处理后的数据。

综合评价的过程可描述为:取一正变换函数f ,使),,,,(21im i i i x x x f y =i y 为被评价对象i o 的综合评价值,对n o o o ,,21按照n y y y ,,,21 的大小(做降序)排序,进而完成n o o o ,,21的优劣比较。通常,将f 取为线性函数,

即有 ： ∑==m

j ij j i x w y 1 （1）

式中：j w 为指标j x 对应的权重系数,∑==≥m

j j j w w 1

1,0。

确定权重的方法很多, 序关系分析法]3[是一种简便优良的主观赋权方法，其基本操作步骤如下:

3.1.1 确定指标重要性的序关系

定义1.若评价指标i x 相对于某评价准则(或目标)的重要性程度大于(或不小于)j x 时,则记为j i x x >.

定义2.若评价指标m x x x ,,,21 相对于某评价准则(或目标)具有关系式：

*

*2*1m

x x x >>> 时,则称评价指标m x x x ,,,21 之间按"">确立了序关系。这里*i x 表示}{i x 按序关系"">排定顺序后的第i 个评价指标),,2,1(m i =,为书写方便

且不失一般性,记为： m x x x >>> 21 （2）

3.1.2 给出相邻两个指标的相对重要程度的比较判断

设专家关于评价指标1-k x 与k x 的重要性程度之比k k w w /1-的理性判断为:2,,1,,/1 -==-m m k r w w k k k 。

k r 的赋值:指标1-k x 与指标k x 具有同样重要性,1=k r ；指标1-k x 与k x 稍微重要, 2.1=k r ；指标1-k x 比指标k x 明显重要, 4.1=k r ；指标1-k x 比指标

k x 强烈重要6.1=k r ,指标1-k x 比指标k x 极端重要, 8.1=k r 。给出k r 的点值形式的判断,则指标相应的权重m x x x ,,,21 可由下述定理得到。

3.1.3 算权重系数

定理1 若专家(或决策者)给出了指标两两之间重要性的比值判断k r ,则m w 为：12)1(-==∑∏+=m

k m

k i i m r w ,而2,,1,,1 -==-m m k w r w k k k 。

在实际评价中,由于问题的复杂性,人类认识的局限性及思维的不确定性,专家通常不能给出),1,,2(m m k r k -= 的精确点值形式,而更有把握给出k r 的区间值,此时综合评价问题的求解显得比较复杂,本文将给出一种仿真模拟的算法。

定义3 若专家给出的关于评价指标1-k x 与k x 的重要性程度之比k r 为区间数,则称k

r 为不确定比值判断,记为),1,,2(~m m k r k

-= 。

设u k l k u k l k u k l k k r r r r r r r ,,],,[~

≤可参照表4在区间]8.1,1[上连续选取,即]8.1,1[,∈u k l k r r 。

定义4 设).,2(*m k r k =为不确定比值，判断k r ~区间上服从均匀分布的随

机变量,称),,,(*

1*2*1*m w w w w =为序关系随机模拟权,这里

,)1(12*

-==∑∏+=m

k m

k i i m r w 2,,1,,*

**1 -==-m m k w r w k k k 。 定义5 对任意2个被评价对象)(,**i i o o i i ≠'',称)(*i i o o s >'为i o '对*i o 的优胜度,有： )(5.0)()(1

1

**1

1

**

***∑∑∑∑=='==''=+>=>m

j m

j j j i j

j i m

j m

j j

j i j

j i i i w x w x P w x w x P o o s （3）

式中集合函数P 为事件概率。

对优胜度)(*i i o o s >'有如下的一些基本性质: 性质1:有界性.即1)(0*≤>≤'i i o o s 性质2:归一性.即1)()(**=<+>''i i i i o o s o o s

性质3:等价性.即5.0)()()(***=<=>?-'''i i i i i i o o s o o s o o s 性质4:自反性.即5.0)(*=>'i i o o s

3.2 随机模拟的软件实现步骤

本文运用蒙特卡罗仿真的方法求解被评价对象两两之间的优胜度及随机模拟权,步骤如下:

步骤1. 选取被评价对象集中任意两个被评价对象*,i i o o '(若有n 个被评价对象, 共需进行n n 2-次比较).

步骤2. 设置仿真次数监控变量count(初始化时count 置1). 步骤3. 运用数学工具包（SAS 或Matlab ）中的随机数发生器产生给定区间),,2(~m k r k =上均匀分布的不确定比值),,2(*m k r k =,进而确定出序关系随机模拟权*w .

步骤4. 设置计数变量t e s r r r ,,( 初始化时t e s r r r ,,均置0) 分别表示

"","",""***i i i i i i o o o o o o <=>'''的次数.

若)

,(),(*1*

1*

i i w i i w j m

j j i j

m j j i '''>'''∑∑=='λλ,则1+=s s r r ; 若)

,(),(*1*

1*

i i w i i w j m

j j i j

m

j j i '''='''∑∑=='λλ,则1+=e e r r ; 若)

,(),(*1

*

1*

i i w i i w j m

j j i j

m

j j i '''<'''∑∑=='λλ,则1+=t t r r ; 步骤5. 1+=count count 若sum count =(sum 为决策者给定的仿真总次数,一般指标数量越多,sum 值越大),则转步骤6,否则转步骤2;

步骤6. 统计优势度)(*i i o o s >'的仿真值,(sum r r o o s e s i i /)5.0()(*

+=>',保存

数值,若已完成被评价对象之间的两两对比,则退出程序,否则转步骤1.

通过模拟仿真,得到n 个被评价对象的优胜度矩阵S ,记),(j i ij o o s s =,有

????????????==?nn n n n n n n ij s s s s s s s s s s S 2

12222111211][ 式中: S 对角线上的元素均为5.0,理论上,应满足1=+ji ij s s , 故只需得

到对角线以上(或以下)，2/n)(n 2-个元素的值便可确定出S 中的所有元素,但是从精确性角度考虑,可独立地模拟出ij s 的值,此时,若1-+ji ij s s 足够小,表明

仿真次数足够,结果是精确可信的;否则,需要增加仿真次数，即提高sum 值.

模拟出的S 在满足精度要求后,需要对其中的元素进行调整,调整公式为2/)]1([ji ij ij s s s -+=,调整后的元素严格满足1=+ji ij s s .仍记调整后的优胜度矩阵为S .

3.3 综合评价值的确定

定义6 i z 为被评价对象i o 的相对综合评价值(从优胜度矩阵S 中导出

的),T

n z z z z ）

（ ,,21=为相对综合评价值向量,满足c c z z z n ,21=+++ 为一常数,则称 )(ji ij j i ij s s z z ---=σ （4）

为被评价对象i o 与j o 的比较函数.

定义7 关于n 个被评价对象之间的整体比较函数为：

∑∑∑∑====---==n

i n

j ji ij j i n

i n

j ij s s z z z F 11

211

2

)]([)(σ

（5） 从偏差最小的角度出发,)(z F 越小越好,最佳的T n z z z z ）（ ,,21=将使)(z F 取最小值。

规划模型为: ??

?

??=+++---==∑∑∑∑====c z z z t s s s z z z F n n i n

j ji ij j i n i n j ij 2111

2112

..)]([)(min σ （6） 定理2 模型(6)达到最优时, ()n

c

s s n z n j ji ij i +-=∑=11 （7）

证明：构造拉格朗日函数：)()(),(L 1∑=-+=n

i i c z z F z λλ

令,0=??i z L 得到：0)]([21

=+---∑

=n

j ji ij j i s s z z λ即

0)(][21

=+---∑=λn

j ji ij i s s c nz

（8） 式( 8) 两边对i 求和,,0])([21

11

=+---∑∑∑===n i n i n

j ji ij i n s s nc z λ推得：

∑∑==-=n i n

j ji ij s s n 11

)(2λ （9）

将式(9)代入式（8)，可得

()n

c

s s n s s n z n j ji ij n i n j ji ij i +-+-=∑∑∑===11121)(1

（10） 因为，0)(1112=-∑∑==n i n j ji ij s s n 所以()n

c s s n z n j ji ij i +-=∑=11，证毕

由式(7)可知,c 的选取对T

n z z z z ）（ ,,21=中元素的大小顺序没有影响,

因而也不会影响被评价对象最终的排序。

第4节 实证分析与研究

4.1 数据来源

本数据来源于2010中国统计年鉴，详见附录一。

4.2 指标选择

为客观评价2010年16个省份的经济发展水平，本着评价的科学性、可操作

性、全面性和可比性原则，选取相关指标。选取了反映各地区经济发展状况的7个评价指标，以较为全面地反映各地区发展的实际情况，其中衡量地区经济实力的地区生产总值GDP,在这里选用GDP 增长率为一指标，它同时也反映了该地区的经济增长速度；及地方财政收入，和衡量地区的质量和效益的指标全社会固定资产投资和衡量地区富裕程度的指标人均GDP 、社会消费品零售总额、人均可支配收入和衡量该地区开发程度的指标进出口总值]4[。记号如下：

4.3 综合评价排序

4.3.1 应用主成分分析对于对评价结果验证

1、利用SAS 软件（见附录）求相关阵R 的特征值λ和特征向量u ，并计算出贡献率e 和累积贡献率E ，结果见表2．

2%37.87，并

以该精度将变量由7维降到2维，可依2个主成分进行综合评价．

根据这2个主成分对应的方差贡献率为权数计算综合得分值，结果如表3．

22

12

1211z e e e z e e e Z +++=

（11）

表格 3 16个省、市经济发展水平综合评价及排序

4.3.1 应用随机模拟对16省经济发展水平的评价

1、评价指标无量纲化处理：对于观测数据X ,为了消除量纲的影响，决策之

前首先应将评价指标无量纲化处理。)min()max()min(*

j j j ij ij x x x x x --=

，其中)min(j x 为第j 个指标中最小的观测值，)max(j x 为最大值，处理后，便得到16个省市的经济发展水平归一化数据表，见表4。

2、由定义1和定义2,确定评价指标的序关系，将各指标排列成：m x x x >>> 21。

3、 由定义3,经过模拟，假定专家根据问题给出的不确定比值判断k r ~

在]8.1,1[上连续取值为],5001.1,3103.1[],417.1,2103.1[],3558.1,1103.1[],2956.1,0103.1[ ]7706.1,5103.1[],5467.1,4103.1[。

4、由定义4和定义5,求得随机模拟权向量,,124,0.18200.2718,0.2(*=w ,0.0371)971,0.06340.1362,0.0，求解出被评价对象两两之间的优胜度矩阵S

( 调整后) ：???

?

???

?

???

??

???=5.03333.05

.016667.05.05.015.05.05

.010

005.0

S 5、由定义6, 定义7 及定理2 得到被评价对象）（16

,,21=n o o o n 的相对综合评价值向量,在这里，n c =，则()n

s s n z n j ji ij i 1

11+-=∑=，那么：

,0.9281,

,1.4428994,1.21415,1.35,0.9,0.5,1.187875,1.02080.3125,0.6(=z

）

88915,0.849,0.3916,0.746,1.27641.2764,1.6。 6、根据相对综合评价值向量, 得到被评价对象的综合排序:

4.3.3 对表3和表5 结果的分析说明

比较两种方法的排序结果，大多数省份的排名还是一致的，因此，随机模

拟的方法在一定程度上可以说明该地区的发展状况。

第5节 程序输入及输出结果

略

第6节 总 结

本文研究了16个省份的经济发展水平，经过查询资料，确定了七个影响指标：GDP 增长率，地方财政收入，全社会固定资产投资，人均GDP ，社会消费品零售总额，人均可支配收入，进出口总值；在主成分分析中，选取累积贡献率大于%85的主成分对其进行评价，结合我们所了解的省份发展的具体情况，广东，江苏，浙江，山东，福建都属于沿海省份，进出口贸易频繁，劳动力密集，经济发展水平较高，两种方法得出的结论较为一致；而湖北作为我国南方重工业，高新科技地区，又有着较为有利的地理条件，发展水平应该在河北，河南这些农业省份之上，山西的煤矿业所带来的经济效益也应该刚高于内陆省份河北，江西。根据各方面的经济实力综合，随机模拟得出的结论更符合我国的实际情况。从而得到如下结论：主成分法是一种客观的评价方法,而随机模拟的

方法本质上属于一种主观评价方法, 这是二者的根本区别，它虽然降低了信息的精确度, 但是却提高了信息准确度。从方法结构上看,该方法可拓展性很强，仿真模拟部分是独立的,函数可任意选取,因而可将基本思想及方法推广用于解决类似的很多问题。

参考文献

[1] 张晓冉.统计分析及其SAS实现.清华大学出版社.2011,5.

[2]郭亚军，易平涛，李玲玉.基于随机模拟的综合评价方法及应用.东北大学学报(自然科学版).2010,31(10):1499-1503.

[3] 郭亚军,潘德惠.一类决策问题的新算法[J].决策与决策支持系统,1992,2( 3):56- 62.

[4] 张吉献,张艳玲.中原经济区各城市经济发展水平评价及发展对策.经济师.2009,1:277-278.

[5]国家统计局网站.

https://www.360docs.net/doc/4512242601.html,/tjsj/ndsj/2011/indexch.htm.2012-06-08.

附录一

1x :GDP 增长率; 2x :地方财政收入; 3x :全社会固定资产投资; 4x :人均GDP;

5x :社会消费品零售总额; 6x :人均可支配收入; 7x :进出口总值