数学建模社区-数学中国

标题: 1989年度西山矿务局五个生产矿井经济效益综合评价问题 [打印本页]
作者: 小丫头妮妮呵呵    时间: 2010-9-4 16:41
标题: 1989年度西山矿务局五个生产矿井经济效益综合评价问题
1989年度西山矿务局五个生产矿井经济效益综合评价问题
1 问题的提出
1989年度西山矿务局五个生产矿井实际资料如表1所示，对西山矿务局五个生产矿井1989年的企业经济效益综合评价(至少用两种方法)。
表1 1989年度西山矿务局五个生产矿井技术经济指标实现值
指  标
白家庄矿
杜尔坪矿
西铭矿
官地矿
西曲矿
原煤成本
99.89
103.69
97.42
101.11
97.21
原煤利润
96.91
124.78
66.44
143.96
88.36
原煤产量
102.63
101.85
104.39
100.94
100.64
原煤销售量
98.47
103.16
109.17
104.39
91.90
商品煤灰分
87.51
90.27
93.77
94.33
85.21
全员效率
108.35
106.39
142.35
121.91
158.61
流动资金周转天数
71.67
137.16
97.65
171.31
204.52
资源回收率
103.25
100
100
99.13
100.22
百万吨死亡率
171.2
51.35
15.90
53.72
20.78
2 问题的假设
1 假设题中提供的实际资料都是可靠准确的；
2 假设除了题中提供的指标外，经济效益不再受其它因素影响。
3 符号的约定
评价对象的个数；
评价指标的个数；
参考数列；
比较数列， ；
表示时刻；
关联系数；
分辨系数， ；
第 个评价对象对理想对象的灰色关联度；
指标的个数；
矩阵；
相关系数矩阵；
单位特征向量， ；
原煤成本；
原煤利润；
原煤产量；
原煤销售量；
商品煤灰分；
全员效率；
流动资金周转天数；
资源回收率；
百万吨死亡率；
4 问题的分析
这是一个综合评价问题，评价方法有很多种。层次分析法、主成分分析法、模糊评价法、灰色关联分析法等。我们可以利用其中的方法对题目进行求解。
5 模型的建立与求解
5.1灰色关联分析模型
基本定义
设评价对象为 个，评价指标为 个，参考数列为
，其中 表示时刻
比较数列为
则称
，其中 为分辨系数
8 E' ?5 @1 \. |/ P3 G0 |
* N* g/ m2 J0 `. c! c+ x& ~! a

( u4 {: \. f- \, e( F
为比较数列 对参考数列 在 时刻的关联系数。
称 、 分别为两级最小差及两级最大差。
一般来讲，分辨系数 越大，分辨率越大； 越小，分辨率越小。
计算灰色关联度
灰色关联度的计算公式为
                                 
这里 为第 个评价对象对理想对象的灰色关联度。
评价分析
根据灰色关联度的大小，对各评价对象进行排序，可建立评价对象的关联序，关联序越大其评价结果越好。
模型的求解
原煤利润、原煤产量、原煤销售量、商品煤灰分、全员效率、流动资金周转天数、资源回收率指标属于效益型指标，原煤成本、百万吨死亡率属于成本型指标。现分别对上述指标进行规范化处理。见表2
表2 比较数列和参考数列值
评价指标
待选生产矿井
最优生产矿井
白家庄矿
杜尔坪矿
西铭矿
官地矿
西曲矿
指标1
0.5864
0
0.9676
0.3981
1.0000
1
指标2
0.3931
0.7526
0
1.0000
0.2828
1
指标3
0.5307
0.3227
1.0000
0.0800
0
1
指标4
0.6196
0.3480
0
0.2768
1.0000
1
指标5
0.7478
0.4452
0.0614
0
1.0000
1
指标6
0.9625
1.0000
0.3114
0.7028
0
1
指标7
1.0000
0.5070
0.8044
0.2500
0
1
指标8
0
0.7888
0.7888
1.0000
0.7354
1
指标9
0
0.7717
1.0000
0.7565
0.9686
1
取 ，计算 及 的值。见表3
表3 关联系数和关联度值
评价指标
待选生产矿井
白家庄矿
杜尔坪矿
西铭矿
官地矿
西曲矿
指标1
0.4602
1.0000
0.3407
0.5567
0.3333
指标2
0.5599
0.3992
1.0000
0.3333
0.6388
指标3
0.4851
0.6078
0.3333
0.8621
1.0000
指标4
0.4466
0.5896
1.0000
0.6437
0.3333
指标5
0.4007
0.5290
0.8906
1.0000
0.3333
指标6
0.3419
0.3333
0.6162
0.4157
1.0000
指标7
0.3333
0.4965
0.3833
0.6667
1.0000
指标8
1.0000
0.3879
0.3879
0.3333
0.4047
指标9
1.0000
0.3932
0.3333
0.3979
0.3405
0.5586
0.5263
0.5873
0.5788
0.5982
由表3知，按灰色关联度排序可看出， ，由于生产矿井5(即西曲矿)的关联度与最优生产矿井的关联度最大，亦即生产矿井5(即西曲矿)的经济效益优于其他生产矿井的经济效益。
5.2 主成分分析模型
基本定义与基本思想
主成分分析也称为主分量分析，是一种通过降维来简化数据结构的方法：如何把多个变量(指标)化为少数几个综合变量(综合指标)，而这几个综合变量可以反映原来多个变量的大部分信息。为了使这些综合变量所含的信息互不重叠，应要求它们之间互不相关。
基本步骤
设有 个样品，每个样品观测 个指标，将原始数据写成矩阵
1)将原始数据标准化，即对同一变量减去其均值再除以标准差，以消除量纲影响，得到标准化的矩阵；
2)根据标准化后的矩阵计算相关系数矩阵 ；
3)求相关系数矩阵 的特征值及对应的单位特征向量 ，
其中 ， ；
4)写出主成分。
其中
模型的求解
运用SPSS软件求出相关系数矩阵的特征值、方差贡献率以及累计贡献率。见表4
表4 相关系数矩阵的特征值和方差、累计贡献率
特征值
方差贡献率( )%
累计贡献率 %
1
3.222
35.800
35.800
2
2.816
31.291
67.091
3
2.651
29.457
96.548
4
.311
3.452
100.000
5
2.127E-16
2.364E-15
100.000
6
8.023E-17
8.915E-16
100.000
7
-1.612E-16
-1.791E-15
100.000
8
-2.341E-16
-2.601E-15
100.000
9
-3.161E-16
-3.513E-15
100.000
由表4可以看出，前3个特征值累计贡献率已达96.548%，说明前3个主成分基本包含了全部指标具有的信息，我们取前3个特征值，并计算出相应的特征向量：
表5 前3个特征值的特征向量
第一特征向量
第二特征向量
第三特征向量
0.25
0.42
-0.28
0.04
0.48
-0.35
0.32
-0.17
0.47
0.27
0.31
0.44
0.15
0.41
0.38
-0.46
-0.23
0.25
-0.51
0.18
-0.16
0.33
-0.43
-0.22
0.41
-0.20
-0.33
因而前三个主成分分别为
第一主成分
第二主成分
第三主成分
在第一主成分的表达式中，除第2项和第五项指标之外，其余变量的系数均在0.25~0.5之间，说明第一主成分是其余七个变量的综合。在第二主成分的表达式中，第1、2、4、5和第8项指标的影响大。尤其是第1和第2项指标影响尤大。在第三主成分的表达式中，第3和第4项指标影响尤大。
计算综合评价函数
通过SPSS软件得出前三个主成分的主成分函数值、综合函数值以及企业排名。见表6
表6 前三个主成分的主成分函数值、综合函数值以及企业排名
排名
白家庄矿
119.28
40.06
28.20
63.54
4
杜尔坪矿
46.77
102.07
66.01
68.12
5
西铭矿
-.66
20.56
74.92
28.26
2
官地矿
16.96
114.87
14.05
46.15
3
西曲矿
-38.04
82.64
40.00
24.02
1
附录
clc, clear
a=[99.89
0 p* n- q1 a& M9 v4 g3 B103.69! \/ @5 f1 s5 ?( L  x- |
97.423 K/ @; E5 {) N- {
101.11
, I* y5 X9 h8 w( i! o2 _97.21
96.91
, a; k6 N& b+ \( c5 m0 Q124.78
: r: C: `2 a0 f66.44
0 O" D- u. f) ?' J3 w8 b143.962 L% h8 _5 P- V
88.36
102.63
8 Z$ t# x/ Q$ v* i  Q101.85
& e+ ^$ n+ Z) v# q: a* G104.39
3 Q* `: @& n$ a& a: Y100.94
8 p; ?3 A) x3 x7 B$ b# H$ ~9 y, o100.64
98.47) E& x9 J/ h, M. C* N
103.16! x$ ~' e8 i1 M7 s
109.172 v, V5 W4 X3 v1 V: S
104.39
3 q7 n' x: E4 d- D9 k! J, O91.90
87.51
1 o8 ~1 I: G# j( i. B4 G+ K, F) u90.27  y4 G: L/ {% w; R: S. H
93.77
" l& ]7 l+ l' r" Q94.33# P4 C4 I1 i  ^8 g: U8 ?4 d
85.21
108.35
% K7 Y) j% c) F  v$ ^1 B9 V+ ^- k106.399 T! M8 Q/ d0 p# N0 H; X, F
142.356 l! `0 b- |8 Z* ?0 m0 _' q2 v
121.918 d, b. j( y4 y& ~; Y
158.61
71.67
7 c" `9 g' {; }; _/ G9 i137.16) x. ^3 {6 c7 ?4 {$ {' L
97.65
" d' |4 n* }: V9 z171.312 W% `" y8 G* i) M3 J
204.52
103.25  I  w4 Z# X" O( ~7 s% z
100. r3 C0 K$ E- t) m" s7 n
100
, Y4 c" L2 Y/ W/ v; w6 W99.13
# {" w% t- B: K4 e# V, x" x  Y" d100.22
171.2
8 S) |  `3 G- i# A% ~4 S1 v9 \51.35( h/ v# D! P% K! R
15.90
4 r/ I- A  R8 x# X53.72) e6 T- [: u2 K6 R
20.78
];
for i=[1 4:9]
a(i,=(a(i,-min(a(i,))/(max(a(i,)-min(a(i,));
end
for i=2:3
a(i,=(max(a(i,)-a(i,)/(max(a(i,)-min(a(i,));
end
[m,n]=size(a);
cankao=max(a')';
t=repmat(cankao,[1,n])-a;
mmin=min(min(t));
mmax=max(max(t));
rho=0.5;
xishu=(mmin+rho*mmax)./(t+rho*mmax);
guanliandu=mean(xishu)
[gsort,ind]=sort(guanliandu,'descend')
作者: 小丫头妮妮呵呵    时间: 2010-9-4 16:45
呵呵
作者: 794150845l    时间: 2010-9-4 17:23
顶一下！！！！
作者: 小丫头妮妮呵呵    时间: 2010-9-4 17:24
回复 794150845l 的帖子4 E# N5 S  c9 r* W1 _
- Q9 [  F' Q; q$ ]! d# W, M2 g
谢啦！呵呵6 a. o* p* g5 X: g
   
作者: liubin6p6p    时间: 2010-9-4 20:00
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  
作者: xunavy    时间: 2010-9-5 08:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: 戈朗    时间: 2010-9-5 12:00
来报道!!!!!!!!!!!
作者: 李问问    时间: 2010-9-5 15:00
我来了~~~~~~~~~ 闪人~~~~~~~~~~~~~~~~  
作者: yeppy    时间: 2010-9-5 15:19
喔的天( ⊙ o ⊙ )啊！这也行
作者: 小丫头妮妮呵呵    时间: 2010-9-5 17:43
回复 liubin6p6p 的帖子" N; S6 S4 T: `; m! w5 a* u

- ?: P9 F3 H/ G+ d0 j! F. z4 u0 hThank you!呵呵
3 A5 f5 T* P; ^   
作者: 小丫头妮妮呵呵    时间: 2010-9-5 17:43
回复 xunavy 的帖子. H' d. @0 \) I5 P( n, H
谢谢！
3 M$ i. Z& ?# a$ L3 t7 E' s% |, [  I7 }7 v9 X7 K
   
作者: web_star    时间: 2010-9-5 20:00
留个脚印```````
作者: mixin    时间: 2010-9-6 08:00
来报道!!!!!!!!!!!
作者: qiaoxlong    时间: 2010-9-6 12:00
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  
作者: mengxue11180611    时间: 2010-9-6 15:00
顶顶更健康，越顶吃的越香。
作者: loverose    时间: 2010-9-6 20:00
强人，佩服死了。呵呵，不错啊
作者: normanlai    时间: 2010-9-7 08:00
看起来好~~像啊~~~~~
作者: huanglonghaoliu    时间: 2010-9-7 12:00
来报道!!!!!!!!!!!
作者: zhnagbei_2004    时间: 2010-9-7 15:00
鉴定完毕！  
作者: ydnzd    时间: 2010-9-7 20:00
试试运气啦~~~~~~~~~~~
作者: abcd9876    时间: 2010-9-8 08:00
强烈支持。楼主万岁
作者: jfwiang    时间: 2010-9-8 12:00
强人，佩服死了。呵呵，不错啊
作者: lgfree    时间: 2010-9-8 15:00
我基本上是采用看英语文章的办法，先泛读，再精读，再一句一句看，最后再提纲挈领，总算是明白一点了，当然，也可能还是领悟错了。最后要说的一句话是：楼主，你很牛叉，希望你不是真的有病。   
作者: dajun539    时间: 2010-9-8 20:00
强人，佩服死了。呵呵，不错啊
作者: yefeng    时间: 2010-9-9 08:00
呵呵 大家好奇嘛 来观看下~~~~  
作者: cow990191    时间: 2010-9-9 08:00
不错不错,我喜欢看  
作者: xiaogui_zhou    时间: 2010-9-9 12:00
楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   
作者: lvyongzhi    时间: 2010-9-9 15:00
鉴定完毕！  
作者: byxh111    时间: 2010-9-9 20:00
鉴定完毕！  
作者: zjnuzhc    时间: 2010-9-9 20:00
强烈支持。楼主万岁
作者: 362514    时间: 2010-9-10 08:00
我回不抢呢 考虑再三 还是不抢了吧 ^_^
作者: caobin    时间: 2010-9-10 11:39
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  
作者: mxyo3111    时间: 2010-9-10 11:51
强人，佩服死了。呵呵，不错啊
作者: rchjia    时间: 2010-9-10 12:00
强人，佩服死了。呵呵，不错啊
作者: cdly1    时间: 2010-9-10 15:00
强烈支持。楼主万岁
作者: xuefei_310    时间: 2010-9-10 20:00
哦～～
作者: 星期九    时间: 2010-9-11 08:00
呵呵 大家好奇嘛 来观看下~~~~  
作者: ken_he    时间: 2010-9-11 12:00
强人，佩服死了。呵呵，不错啊
作者: algebra    时间: 2010-9-11 15:00
楼主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  
作者: Ready    时间: 2010-9-11 20:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: 橘子先生    时间: 2010-9-12 08:00
声明一下：本人看贴和回贴的规则，好贴必看，精华贴必回。
作者: brainyrain    时间: 2010-9-12 12:00
楼主那种裂纸欲出的大手笔，竟使我忍不住一次次的翻开楼主的帖子……   
作者: Beatrice    时间: 2010-9-12 15:01
哦～～
作者: kxhu    时间: 2010-9-12 20:00
哦～～
作者: 星竹石    时间: 2010-9-13 12:00
我回不抢呢 考虑再三 还是不抢了吧 ^_^
作者: talent212121    时间: 2010-9-13 15:00
鉴定完毕！  
作者: hyxie616    时间: 2010-9-13 20:00
楼主，你写得实在是太好了。我惟一能做的，就只有把这个帖子顶上去这件事了
作者: ljk_770    时间: 2010-9-14 12:00
试试运气啦~~~~~~~~~~~
作者: dajun539    时间: 2010-9-14 15:00
我要把这个帖子一直往上顶，往上顶！
作者: futures    时间: 2010-9-14 20:00
我要把这个帖子一直往上顶，往上顶！
作者: 流氓兔    时间: 2010-9-15 15:00
哦～～
作者: sjs914915    时间: 2010-9-15 20:00
试试运气啦~~~~~~~~~~~
作者: chxq666    时间: 2010-9-15 20:00
哦～～
作者: lylhx    时间: 2010-9-16 12:00
顶顶更健康，越顶吃的越香。
作者: xzhg    时间: 2010-9-16 15:00
鉴定完毕！  
作者: wwzwu    时间: 2010-9-16 20:00
顶顶更健康，越顶吃的越香。
作者: gnemgnayfight    时间: 2011-5-1 01:18
老看不见符号和公式！！！2 V8 @0 L- H1 D$ e, E" m# c3 E
作者: gnemgnayfight    时间: 2011-5-1 01:24
楼主给个原版吧？？？？谢谢啦
作者: 李贤    时间: 2011-5-27 21:39
主的帖子实在是写得太好了。可是我立刻想到，这么好的帖子，倘若别人看不到，那么不是浪费楼主的心血吗？经过痛苦的思想斗争，我终于下定决心，牺牲小我，奉献大我。我要拿出这帖子奉献给世人赏阅，我要把这个帖子一直往上顶，往上顶！顶到所有人都看到为止！  
作者: super-wem    时间: 2012-7-12 14:24
看看。。。。



欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5