【数学建模】熵值法与多指标评价系统9 x$ U. u2 Q7 J) h9 ?6 M
一.熵的基本概念' q5 v5 Z; L" w' b; a' D/ B, L
熵最初是热力学中的一个概念,后来延伸到信息领域。熵是对系统无序度的度量,系统越无序,熵越大。熵的统计学定义是:4 K# a. K: G1 ^
& {, [4 z# x7 g
某个宏观状态的微观状态数,取对数lnWlnW! l9 g0 A, g. B8 I
举个例子,抛10次硬币,正面朝上的次数是3是宏观状态,对应的微观状态数,也就是使得正面朝上的排列组合数为3。可以看出熵越大,该熵对应的状态出现的概率也越大,弄清楚定义之后也就不难理解为什么熵越大系统会越稳定了。9 G" Z" m- [( ^# G# `
" o' W; q% p( H6 {, |7 ?1 R, v0 f0 g二.熵值法在多指标评价系统中的应用# D1 [6 X R# W, E7 m7 P
对于多指标系统评价的赋权问题除了采用层次分析法(AHP)等主观分析方法外,还可以通过熵值法来确定权重。 ( |4 x8 `/ O* W5 y/ L: O1 b/ Y, M
熵值法的核心思想是用信息的无序度来衡量信息的效用值。信息的无序度越低(越不稳定),该信息的效用值就越大。换句话说,越稳定的信息越无用。 + y, X0 H! q; q5 X7 ?" z* N在多指标系统中应用熵值法是基于以下假定:7 N" ?, X F# `& [
5 H ^+ H0 ]* S6 }. p* O
那些很稳定基本不怎么变化的指标对最终评价造成的影响也很小。信息的离散程度越大越重要。 9 S' _: _5 ^- Y& U8 h, V* o) n# ]- N* o7 N5 Y
国内外关于熵值法的资料不是很多,我认为主要原因还是这个假定不是很有说服力,因为实际中存在大量比较稳定但也同样重要的信息,这种信息即使是很小一点的变动也会引起结果的很大变化,在非线性动力系统中有大量这样的例子。 5 p ?1 m) m9 y我认为,信息的离散程度只是信息的某一方面的反映,信息的重要性也是信息某方面的反映,单单从信息的离散程度或是重要性都无法还原信息本身,而通过信息的某个侧面去推测信息的另一个侧面这种方法确实不太靠谱。 3 ]1 ]% ^0 D7 [0 ` D' V% K( E1 @8 r
综上,这种方法虽然是一种客观的评分方法,但由于理论上的硬伤,还是不能在多指标评价系统中出于优势地位。 % s- c6 z% {" u" {$ i
不过由于评价这种东西本来就很主观,也没有客观的标准去评价某个方法的有效性,所以对多指标评价的方法也很难分出孰优孰劣。如果是有专业知识作为指导还是最好选择主观赋权的方法。- }% \+ d. n$ O- B; a& i, H4 r# {* n& L
( ]5 N5 \+ c ^6 ?7 N三.熵值法算法步骤! K( m& \ I: I5 J( p W
选取n个样本,m个指标,则 & z& f9 ?( J3 c( z1 |8 }( |, t7 r- [8 H9 r& ?
xijxij' b) a9 h" ?' |2 L$ c
为第i个样本的第j个指标的数值(i=1, 2…, n; j=1,2,…, m); / {3 @" |" P5 i+ W/ [/ u( [. w L% k2 m$ i0 w( j% F. p
指标的归一化处理:异质指标同质化- y0 V/ C% n) H, p$ p# ]
3 d! e% _3 P& ?" V- i3 S" I( r
由于各项指标的计量单位并不统一,因此在用它们计算综合指标前,先要对它们进行标准化处理,即把指标的绝对值转化为相对值,并令 ; a# C4 b) i2 L9 u5 `# S; {8 L3 u7 ?5 l3 \% P, |2 f: ~
xij=|xij|xij=|xij|" k; l. s" W( F1 H! q
从而解决各项不同质指标值的同质化问题。而且,由于正向指标和负向指标数值代表的含义不同(正向指标数值越高越好,负向指标数值越低越好),因此,对于高低指标我们用不同的算法进行数据标准化处理。其具体方法如下: 5 P; `4 y) T# {$ V, o- y) A! r+ }. D7 K7 {( Z9 g+ g
对于正向指标:6 ^1 \+ U4 O9 k
2 v7 a1 t* y( ]+ u2 X, g
zij=xij−min(xj)max(xj)−min(xj)zij=xij−min(xj)max(xj)−min(xj)0 T$ s' K/ v: [
对于负向指标:4 `/ T1 l& }( j8 C& p) ]
! B' j) u2 D6 z* X$ b& q6 O5 B1 Wzij=max(xj)−xijmax(xj)−min(xj)zij=max(xj)−xijmax(xj)−min(xj)0 r' R: N9 S% w) {$ w4 _" W! |
则zijzij为第i个样本的第j个指标的数值(i=1, 2…, n; j=1, 2,…, m)。 6 D" I# K- t* d% a* v& L! Z$ p' h ]3 Q, V h' V
计算第j项指标下第i个城市占该指标的比重:. J" ~9 e4 w8 L( N V
5 W: x# E1 R" m2 ]pij=zij∑ni=1zij,i=1,...,n,j=1,...,mpij=zij∑i=1nzij,i=1,...,n,j=1,...,m + D/ X: E* K+ w# n+ [ Q& Z计算第j项指标的熵值: 8 s6 b4 r2 P/ L# J' p. y% q6 G% C. u1 q: z2 K* m( U
ej=−k∑ni=1pijln(pij)ej=−k∑i=1npijln(pij)# `9 ]# K1 l9 [6 ]9 F j
其中k=1/ln(n)>0k=1/ln(n)>05 b* \! L# C7 Y$ l8 ]/ p
满足ej≥0ej≥0 " B( v( h) o: w! Q/ |计算信息效用值: - r% i/ f/ D0 H" i. z# e+ q* V : P6 Q. O- p* Y5 `dj=1−ejdj=1−ej+ z4 v8 J/ d4 `+ q7 l% t2 r9 k! B
计算各项指标的权值:" a, u+ R6 J7 ~& E# e
/ X9 D3 a8 c1 {0 X: ?. Wwj=dj∑mj=1djwj=dj∑j=1mdj3 z. ?2 q8 l+ a
多层评价系统的评价' Z7 o3 |& ?3 \5 H% {( ^' V/ k
/ }) Z0 e4 }' z" O0 E根据熵的可加性,可以利用下层结构的指标信息效用值按比例确定对应于上层结构的权重WjWj的数值。对下层结构的每类指标的效用值求和,得到上层指标的效用值和,记作Dk(k=1,2,...,K)Dk(k=1,2,...,K)。进而得到全部指标效用值的总和0 \3 v; K2 U C6 A$ H+ x- {; A
7 p ^2 x. \# u7 {( iD=∑Kk=1DkD=∑k=1KDk3 K/ e% Z, H- k6 q4 t
则相应主题的权重为:/ H; U% o: {+ |4 S6 Q3 |. Y
! o2 ]. x8 H/ d( ?; C# qWk=DkD" @0 I& o" E' d0 n& a
2 Z0 L8 i1 Y4 C V, n; A. \* `
0 j# u9 t1 E6 @$ [0 l& I ~ B
% G& ^* A- t y: G" j' Y