宜居城市建模分析

《宜居城市建模分析》由会员分享，可在线阅读，更多相关《宜居城市建模分析（28页珍藏版）》请在金锄头文库上搜索。

1、五一数学建模竞赛承诺书我们仔细阅读了五一数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与本队以外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其它公开的资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们愿意承担由此引起的一切后果。我们授权五一数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒

2、体进行正式或非正式发表等）。参赛题号：C参赛队号：03357参赛组别：本科所属学校：西南交通大学参赛队员：1. 戴潘宁2. 何逸茗3. 卢晓孟日期：2017 年 05 月 01 日获奖证书邮寄地址：西南交通大学犀浦校区邮政编码：611756收件人姓名：何逸茗联系电话：五一数学建模竞赛题题 目目淮海经济区核心区宜居城市评价模型关键词：熵权法 层次分析法 主成分分析法 灵敏度 灰色分析法 粗糙集摘要摘要本文主要探讨了淮海经济区内的 8 个城市的宜居评价问题， 建立了包括综合评价模型与灵敏度分析的数学模型，并设计算法进行了求解。针对问题一，本文首先选取了城市宜居指数作为评价城市宜居情况的一级指标体系

3、，逻辑上肯定了客观评价指标和主观评价指标的总要性。结合国内外相关文献资料，通过统计分析，从性质上选取了安全性、健康性、便捷性、舒适性、可持续性共五个指标作为二级指标体系，以人均 GDP、城市交通网络运输能力、年空气达标率等共 45 个指标作为三级指标体系，建立了针对评价城市宜居的三级评价体系；其次，针对 45 个指标，在国内相关统计年鉴上收集相应数据， 并进行了必要的数据预处理，利用熵权法确定了 45 个指标的权重系数，得到了关于评价宜居城市的综合评价模型。针对问题二， 本文利用问题一中建立的城市宜居评价模型， 代入 2013 年至 2015 年共三年各项指标的标准化期望数据进行计算求解，得到

4、了8 个城市宜居得分（满分分），按得分大小顺序依次排列为：济宁、徐州、连云港、淮北、枣庄、商丘、宿迁、宿州。针对问题三，本文选取主成分分析法，先从问题二模型中筛选出部分有效指标，再对数据标准化预处理，用 SPSS 计算相关系数矩阵、特征值和特征向量。通过求解主成分载荷矩阵，找出各指标对主成分的贡献率，判断出对宜居城市排名产生显着性影响的评价指标， 按灵敏度大小顺序依次排列为： 人均医院拥有数、 人均卫生机构人员拥有数、人均卫生机构床位数、年空气达标率、地表水质量、城镇居民人均可支配收入。针对问题四，本文运用灰色预测法和熵权法相结合，提出了基于不确定性理论的宜居城市评价模型。由于宏观政策调整的自

5、由度具有强烈随机性，所以对宜居城市的影响程度不确定性大，无法进行定量建模分析， 所以本文引入了两个不确定因素：突发自然灾害和房价大幅波动。然后，本文选取了淮海经济区常发生的自然灾害：洪涝、台风、雪灾等，按照自然灾害发生的年数频次，对其进行科学风险赋值。赋值后，考虑各种突发自然灾害的发生概率和灾后对社会的经济与环境造成的损失，与原评价指标多元组合，运用粗糙集简化指标，重新计算出不确定因素的对宜居城市的影响。针对问题五，基于以上计算结果，将徐州市的各项指标与淮海经济区八个城市指标的最优值、 最差值以及平均值进行了比较。 本文建议徐州市政府加快转变城市发展方式，改善城市环境质量， 继续完善城市基础设

6、施的建设， 为提高徐州市宜居水平再接再厉， 。一、问题提出一、问题提出城市宜居性是当前城市科学研究领域的热点议题之一， 也是政府和城市居民密切关注的焦点。 建设宜居城市已成为现阶段我国城市发展的重要目标,对提升城市居民生活质量、完善城市功能和提高城市运行效率具有重要意义。我国宜居城市的排名每年都是热门话题，不同机构对宜居城市的排名结果也不尽相同。中科院在2016年发布了中国宜居城市研究报告，在被调查的40个城市中，排名前十的城市分别为：青岛、昆明、三亚、大连、威海、苏州、珠海、厦门、深圳、重庆。而美世人力资源咨询公司（ WilliamMercer）公布的2016年全球宜居城市排行中大陆前十名分

7、别为上海、北京、广州、成都、南京和深圳（并列）、西安、重庆、青岛、沈阳、吉林。由此可见，宜居城市评价指标体系不同，宜居城市排名结果也会发生变化。问题1. 通过查阅资料，筛选评价宜居城市的主要指标，并阐述这些指标的合理性。根据所筛选的主要指标，建立评价宜居城市的数学模型。问题2. 利用你构建的评价宜居城市的数学模型， 对淮海经济区内的8个城市 （宿迁、连云港、宿州、商丘、济宁、枣庄、徐州、淮北）进行合理性研究，给出宜居城市排名。问题3. 以问题2为例，定量分析你所建立的模型中，哪些评价指标的变化会对宜居城市排名产生显着的影响。问题4. 一些不确定性的因素（如突发自然灾害、房价大幅波动、宏观政策的

8、重大调整等）会对宜居城市的某些指标产生重大影响。建立基于某些不确定性因素的评价宜居城市的数学模型，并重新讨论问题2。问题5. 根据上述定量分析的结果， 请有针对性地给出进一步提高徐州市宜居水平的政策建议。二、问题分析问题一的分析问题一要求我们构建出一个科学合理的评价宜居城市的数学模型。 由于宜居城市的评价与社会发展的诸多因素有关， 本文希望能够得到一个综合的、 较为客观的评价模型。首先，本文考虑宜居城市主要与城市安全性、便捷性、舒适性、健康性、可持续性共五个方面相关，分别由客观评价指标和主观评价指标反映。但由于主观评价反映如居民满意度不易调查，因此本文只从客观评价指标对宜居城市的五个方面进行评

9、价。其次，本文针对这五个方面选定适当的分指标对宜居城市情况进行充分评价， 收集各个分指标所对应的数据，并进行必要的数据预处理。考虑需要评价的是城市宜居情况，本文采用层次分析法建立数学模型。在将城市宜居指数作为目标的前提下，我们查阅了国内外研究宜居城市问题的相关文献，对它们选择的二三级指标进行了统计，并绘制了表格（附录1）。根据统计结果，我们筛选出了五个二级指标，四十五个三级指标。所以，本文将以城市宜居指数为一级指标，以安全性、健康性、便捷性、舒适性、可持续性为二级指标，以人均GDP、年空气达标天数等四十五个指标为三级指标，构建了整个模型框架。然后，需要选取较为客观的方法，判断各个指标对宜居城市

10、影响的大小，即确定各个指标的权值，方可对城市宜居进行综合评价。本文认为，我们所做的评价宜居城市指标不仅应切实地注意到城市的物质丰足性，生活便利性，同时也应注意在城市中生活的人们自我提高与未来发展问题，所以，本文中所给出的模型指标需要是相对科学、合理的。问题二的分析问题二要求我们用问题一所构造的数学模型，对淮海经济区内的 8个城市（宿迁，连云港，宿州，商丘，济宁，枣庄，徐州，淮北）进行城市宜居排名。我们在问题一的基础上广泛搜集资料，查询了淮海经济区八个城市 2013-2015年的大量指标数据，并使用插补法，填补了少数缺失的数据。由此，本文得到了较为详尽的数据表（附录 2）。在已得数据的基础上，先

11、对数据预处理，再运用相对客观的熵权法，对已确定的指标进行权值地划分。将预处理后的标准化数据与对应指标的权值相乘再求和，就得到了淮海经济区八个城市的宜居得分与排名。问题三的分析问题三要求我们以问题二的结果为基础， 定量分析哪些指标对于城市宜居性产生明显的影响。首先，本文选取主成分分析法，按照数理统计原则，筛选出部分有效指标。用SPSS先对有效指标的原始数据进行正态标准化处理，计算其相关系数矩阵，再计算出特征值和特征向量。通过求解主成分载荷矩阵，找出各指标对主成分的贡献率，也就可以判断出对宜居城市排名产生显着性影响的评价指标，并且可以相应地分析。问题四的分析问题四要求我们在已有模型的基础上对如自然

12、灾害，房价大幅波动，宏观政策调整等突发事件进行分析， 讨论不确定因素对模型的影响。 本文运用粗糙集和熵权法相结合，提出了基于不确定性理论的宜居城市评价模型，对淮海经济区内的 8个城市宜居水平进行了合理评价，最终确定了在不确定因素影响下宜居城市得分与排名，同时考虑了评价所设计的随机性和模糊性因素，使其评价结果更客观、有效。问题五的分析问题五要求我们在前四个问题结果的基础上，对徐州市的城市宜居进行具体的分析，分析其优势点与劣势点，再指出徐州市该如何进一步挺高自身的城市宜居水平。三、基本假设三、基本假设1、假设：各个统计网站上提供的相关数据真实有效；2、假设：本文选取的四十五个三级指标符合实际，客观

13、、科学的反映出城市宜居情况；3、假设：本文选取的四十五个三级指标简单明了，均可以量化。4、假设：本文采用的灰度预测法对城市未来不确定因素发生概率基本正确四、符号说明及名词定义四、符号说明及名词定义符号说明符号意义wj各个单因素指标经过熵权法计算后得到的权重ei第i个城市第j项指标标准化后的得分xij第i个城市第j项指标yi第i个城市城市宜居性得分Xij第i个城市第j项指标正态标准化后的得分M相关系数矩阵Xj第j项指标正态标准化后的权重名词解释城市宜居性： 宜居城市是指对城市适宜居住程度的综合评价。 其特征是： 环境优美，社会安全，文明进步，生活舒适，经济和谐，美誉度高。宜居城市建设是城市发展到

14、的产物，是指宜居性比较强的城市，是具有良好的居住和空间环境、人文社会环境、生态与自然环境和清洁高效的生产环境的居住地。1996 年联合国第二次人居大会提出了城市应当是适宜居住的人类居住地的概念。城市基础设施：城市基础设施（urban infrastructure）是生存和发展所必须具备的工程性基础设施和社会性基础设施的总称， 是城市中为顺利进行各种和其他社会活动而建设的各类设备的总称。城市公共服务：城市公共服务就是指城市公共部门面向城市公众提供的公共产品和服务，包括城市基础设施的投资和维护，提供和加强就业岗位，社会保障服务，兴办和支持教育、科技、文化、医疗卫生、体育等公共事业，及时发布有个社会

15、信息，为社会公众生活质量的提高和参与公共事务提供有力的保障和创造相关的条件。城市生态环境： 是城市人类与周围生物和相互作用而形成的一类具有一定功能的网络结构，也是在改造和适应的基础上建立起来的特殊的。五、模型准备五、模型准备数据获取根据本文在后面建立的针对城市宜居性的评价指标体系，本文需要查找的有以下数据：空气质量达标率，地表水、地下水、饮用水合格率，城市噪声环境质量，城市辐射环境质量，城市绿化覆盖率，日平均气温在15-25 度之间的天数等。本文需对 8个城市宜居性情况进行评价，因此需要收集 2013 年至 2015 年以上 44 项指标的数据。本文的数据主要来源于网络公开资料整理及以下网站：

16、1.中国城市统计年鉴山东省统计年鉴：江苏省统计年鉴：安徽省统计年鉴： 河南省统计： 数据预处理计算指标增长比例本文建立的是关于淮海经济区八个城市宜居性的评价模型，单独评价各项指标的数量并不能说明淮海经济区八个城市宜居性的发展状况。基于此，我们进行了部分指标的增长比例的计算，增长比例能更直接地反映淮海经济区八个城市宜居性的发展状况。增长比例计算公式如下：XijUijUi(j 1)Ui(j 1)其中：Uij为第i项指标第j年的数据，Xij为第i项指标第j年相较于第j 1年的增长。数据标准化各个评价指标之间由于各自的度量单位及数量级的差别，存在不可公度性。为了消除不同变量量纲的影响，大部分统计学模型

17、均须要求数据进行无量纲化处理。本文中均以极大极小值法对数据无量纲化至0,1的区间内：Fi Fi minFiFi max Fi min其中：Fi为各指标无量纲化后的数据，Fi max，Fi min分别为该指标中的极小、极大值。六、模型的建立与求解六、模型的建立与求解问题一模型建立与求解问题一要求我们建立评价宜居城市的数学模型，本文通过以下的步骤建模并求解：步骤一：城市宜居性指标体系的选取与建立。步骤二：用熵选法确定各个指标的权重步骤三：建立基于熵权法的综合评价模型。问题一模型的建立城市宜居性评价指标的确定问题一要求对对目前淮海经济区八个城市宜居性情况进行评价， 本文考虑从客观指标进行科学评价。由

18、于宜居指标涉及范围很广，想要完整而客观的选取充分的指标并不容易，因此本文考虑利用统计学院里， 选取典型并能够涵盖足够信息的指标，指标的选取过程如下：宜居城市评价二级指标在各体系中出现次数及相应情况指标所在体系指标名称出现次数/重复出现比率国内/国年研究学者或机构外国外国外国外国内生态环境8/%国内国内国内国内国外国外国内交通及生活便捷度7/%国内国内国内国内国外国外居住条件6/50国外国内国内国内国外安全性6/50份197大卫.史密斯4200英国经济学家智囊团4201美国中老年人协会5199宁月敏9200周志田4200中国城市科学研究会6200建设部6200李嘉菲8199Knox5201美国中

19、老年人协会5200周志田4200中国城市科学研究会6200建设部6200张文忠6200董晓峰91995Knox1998Gideon2015美国中老年人协会1999宁越敏2006建设部2006张文忠199Knox国内国内国内国内国内国外社会环境4/国内国内国内国外基础设施及公共服务4/国外国内国内国外经济环境4/国内国内国内国外教育情况3/25%国外国内就业情况2/%国外国外5200620062006200820092004200420062008199820041999200619982004200620081998200420061998201中国城市科学研究会建设部张文忠李嘉菲董晓峰英国经

20、济学家智囊团周志田中国城市科学研究会李嘉菲Gideon英国经济学家智囊团宁越敏建设部Gideon周志田中国城市科学研究会李嘉菲Gideon . Omuta英国经济学家智囊团建设部Gideon . Omuta美国中老年人协会5200国外英国经济学家智囊团4健康情况2/%200国内张文忠6197历史1/%国外大卫.史密斯4199负面事件1/%国外Gideon . Omuta8基于以上统计结果，我们将其按照性能分类，依次为健康性、方便性、舒适性、经济发展性安全性五大类。查阅资料，每类也可继续细分。因此，结合资料，本文建立了五个三级 45 个四级指标的评价指标体系如下：一级二级三级年空气达标天数率地表

21、水地下水健康性饮用水达标率城市区域昼间噪声（A）城市交通干线昼间噪声（A）垃圾无害化处理率?城市交通网络运输能力（万人）?公路密度两个指标反映?卫生机构数（所/万人）卫生机构人员（人/万人）方便性卫生机构床位（个/万人）学校数（所/万人）各级各类学校数（所）每平方公里拥有超市数目燃气普及率有线电视普及率自来水普及率日平均气温在 15-25舒适性城市绿化覆盖率人均公园绿地面积（平方米）人均淡水量（立方米）人均住房面积人均 GDP（元）城镇居民人均可支配收入城镇失业率经济可持续第三产业在所有产业中所占的比重度第三产业增长率居民消费价格总指数（以上年为100）城市宜居度粮食总产量增长率油料总产量增长

22、率进出口总额增长率人口自然增长率（）人均平均房价平均工资与房价比人均工业产值（万元）对公共安全事件处理成功安全性安全教育表宜居城市评价指标结合上表可知，对淮海经济区内 8 个城市宜居度的评价指标可分为三层，描述城市的健康性指标E1、 方便性指标E2、 舒适性指标E3、 经济发展性指标E4、 安全性指标E5。E E1,E2,E3,E4,E5因此， 一级指标体系， 其中E1包括年空气达标天数率e1、 地表水e2、地下水e3、饮用水达标率e4、城市区域昼间噪声（A）e5、城市交通干线昼间噪声（A）e6、垃圾无害化处理率e7；E2包括城市交通网络运输能力（万人）e8、公路密度e9、卫生机构数（所/万人

23、）e10、卫生机构人员（人/万人）e11、卫生机构床位（个/万人）e12、学校数（所/万人）e13、各级各类学校数（所）e14、每平方公里拥有超市数目e15；E3包括燃气普及率e16、有线电视普及率e17、自来水普及率e18、日平均气温在 15-25e19、城市绿化覆盖率e20、人均公园绿地面积（平方米）e21、人均淡水量（立方米）e22、人均住房面积e23；E4包括人均 GDP（元）e24、城镇居民人均可支配收入e25、城镇失业率e26、第三产业在所有产业中所占的比重e27、第三产业增长率e28、居民消费价格总指数（以上年为 100）e29、粮食总产量增长率e30、油料总产量增长率e31、进

24、出口总额增长率e32、人口自然增长率（）e33、人均平均房价e34、平均工资与房价比e35、人均工业产值（万元）e36；E5包括对公共安全事件处理成功率e37、安全教育e38，因此二级指标集为E1e1,e2,e3,e4,e5,e6,e7、E2e8,e9,e10,e11,e12,e13,e14,e15、E3e16,e17,e18,e19,e20,e21,e22,e23、E4e24,e25,e26,e27,e28,e29,e30,e31,e32,e33,e34,e35,e36E5e37,e38。熵权法确定各个指标的权重评价体系中各指标权重的确定一般包括主观赋权法和客观赋权法两种。 主观赋权法主观随

25、意性大，可靠性不足。因此，本文选择了主观性相对较小，能充分利用数据特征的熵值赋权法。熵权法可根据各项指标的变异度，利用信息熵计算出各指标权重的客观赋权法，信息熵越小，则表明指标的变异程度越大、重要程度越大。计算步骤如下：Step 1:计算标准化后的指标Xij的比重Rij，即：RijXijXi1nijStep 2:计算第j项指标的熵值ejk，即： 1 nejk RijlnRijlnni1Step 3:计算第j项指标差异性系数gjk，即：gjk1ejkStep 4:计算指标Xij的权属wjk，计算公式为：gjkwjkgi1n, j 1,2,3,.; k 1,2,3,.jkStep 5:利用 Mat

26、lab 编程计算权重集合：w |j 1,2,3,.; k 1,2,3,.jk基于熵权法对城市宜居度综合评价模型的建立基于熵选法本文可以分别得到每个评价指标的权重， 因此可以建立城市宜居性情况的多指标综合评价模型如下：yiwjk Xji,i 1,2,3,.; j 1,2,3,.j1n其中，ej为第j项指标标准化后的得分，wj为第j项指标通过熵权法所得到的权重。至此，基于熵权法的综合评价模型求解完毕。问题一模型的求解问题一要求建立评价宜居城市的数学模型，由于评价指标具有多样化，因此运用数理统计的方法，调查国内外城市宜居性评价指标， 选取具有代表性的一部分指标，并统计出各个评价指标在不同评价体系里的

27、频数频率， 最终根据频率得出本文评价宜居城市的评价指标，统计结果见附录 1。问题一结果的分析及验证从表和表可知，本文所建立的评价宜居城市的数学模型中，评价指标的选取符合统计学中统计指标应具有反映现象总体数量特征的原则。为验证指标的有效性，本文从以下四个方面说明：（1）目的性：解决了问题一；（2）科学性：达到了正确、科学的要求；（3）度量性：满足可测定和可计量原则；（4）可比性：保证了指标的相对稳定性。问题二模型建立与求解问题二要求我们基于问题一所建立的评价模型， 去计算淮海经就去 8 个城市的宜居得分，本文通过以下步骤建模并求解：步骤一：计算 2013 年至 2015 年各项标准化的指标数据的

28、期望步骤二：将各项指标处理后的数据代入综合评价模型，计算淮海经济区八个城市宜居情况的得分，并对计算结果进行分析。问题二模型的建立基于熵权法对城市宜居度综合评价模型的建立基于熵权法， 本文可以分别得到 45 个指标的权重， 因此可以建立针对 2013 至 2015年淮海经济区内城市宜居情况的多指标综合模型如下：yiwjej,i 1,2,.,8j145其中，ej为第j项指标标准化后的得分，wj为第j项指标通过熵权法所得到的权重问题二模型的求解本文根据问题一的模型对城市宜居度进行评价。 首先从公开网络资料和权威统计年鉴收集到八个城市 2013 年至 2015 年 45 个指标数据，对八个城市这三年的

29、各项指标数据求期望，得到了 45*8 的原始数据矩阵，经插补法获取少量缺失数据后，得到的原始数据原始数据见附表：对上述所得数据依次进行标准化和归一化，得到了标准化数据。利用熵权法，得到各指标权重系数为wj，见附录。综上所述，得出基于熵权法的淮海经济区城市宜居性的多指标综合评价模型，评价得分表如下表：城宿迁连云港宿州商丘济宁枣庄徐州淮北市三年期望值表 6. 2 三年综合评价得分模型注：以上得分满分为。问题二结果的分析及验证近三年淮海地经济区城市平均宜居性得分为分，发展状况较好。按照得分高低，排名依次为济宁、徐州、商丘、连云港、淮北、枣庄、宿迁、宿州。个体上，济宁在第三产业上表现突出，徐州在交通网

30、络运输能力表现突出，商丘教育设施普及度的值最高，枣庄人均拥有卫生机构值最高，宿迁在声污染方面存在问题，宿州在人均公路面积方面存在问题。 总体上， 相较于 2013 年和 2014 年的淮海地经济区城市平均宜居性得分， 2015年得分有所增长。随着经济的不断发展以及人们对生活品质追求的提高，淮海经济区内8 个城市平均宜居性水平是在不断提高。2015 年得分最高，其中不乏客观因素， 本文首先通过分析 2015 年 45 个三级评价指标，发现可持续性里的多数指标值都是较高的，且这些基本评价指标的权重系数占比很大。其次，在 2015 年，淮海经济区核心区一体化建设推动九项重点工作这则总体规划纲要对拓展

31、区域产业投资合作、拓展区域商贸物流等经贸合作、拓展区域医疗合作覆盖范围、推动区域基础设施互联互通、 健全区域生态环境保护机制、强化区域警务合作等九个方面做出了重要的指示。 其中，基础设施互通、产业协同发展、 公共服务共享、生态环境共治四个方面第一时间得到了全民的关注。这四个个因素均间接导致了 2015 年淮海经济地区 8 个城市宜居情况的提高。问题三模型建立与求解问题三要求我们基于问题二，定量分析对宜居城市排名产生显着性影响的评价指标，本文通过以下的步骤建模并求解：步骤一：对 2013 年至 2015 年各项指标数据的期望做正态标准化处理。步骤二：基于问题二熵权法所求得的结果， 对指标进行有效

32、筛选，建立主成分分析法评价模型。步骤三：对评价指标显着性影响模型进行求解，并计算出结果。问题三模型的建立问题三要求找出对宜居城市排名有显着影响的评价指标， 其关键是将部分相关性很高的指标转化为彼此相互独立或不相关指标。本文在问题一中给出了 45 个描述城市宜居性的三级指标，现在需要确定哪些指标对排名有显着影响。本文基于在力保数据信息丢失最少的原则下，对多变量的截面数据表进行最佳综合简化。对原始数据进行标准化处理的模型建立由问题一可知，通过数据筛选，将基于熵权法模型下权重极小的评价指标舍弃，保留剩下的 29 个。下文中需要进行主成分分析的的指标变量有 29 个，分别为x1,x2,.xj, j1,

33、2,.,29，共有 8 个评价对象，第i个城市的第j个评价指标的值为xij。将各指标值xij转换成标准化指标值Xij。计算步骤如下：Step 1:计算第j个指标的样本均值j，即：18jxijni1Step 2:计算第j个指标的样本标准差Sj，即：18Sjxijj, j 1,2,.,45n1i1Step 3:将各指标xij转换成标准化指标值Xij，即：Xijxijj,i 1,2,.,8; j 1,2,.,45SjStep 4:计算标准化指标变量，即：xjxjj, j 1,2,.,mSj建立求解相关系数矩阵M的特征值及特征向量的模型先利用 SPSS 软件分析已标准化数据，用软件计算得到相关系数矩阵

34、M mijm m，mijXk18ki Xkj,i j 1,2,.,mSj其中：mii1,mij mji，mij是第i个指标和第j个指标的相关系数。再通过 SPSS 软件计算得到主成分载荷矩阵，在其基础上通过软件计算得到主成分特征向量矩阵，相关系数矩阵M的特征值mm1,2,., m及其对应的特征向量mm1,2,., m。 由特征向量组成个新的指标变量：y1 u11x1u21x2.um1xm,y2 u12x1u22x2.um2xm,ym u1mx1u2mx2.ummxm,最终得到各指标对各成分的贡献率以及各成分对于总体评价指标的信息贡献率。选择主成分根据各主成分的信息贡献率，我们得到各个指标在主成

35、分分析法下的新权重。根据各个指标占权重之比，比较得出人均医院拥有数、人均卫生机构人员拥有数、 人均卫生机构床位数、年空气达标率、 地表水质量、城镇居民人均可支配收入这八个指标对宜居城市排名有显着影响。问题三模型的求解利用 SPSS 软件，计算得到的解释总方差表如下：解释的总方差解释的总方差成份合计12345初始特征值方差的 %累积 %合计提取平方和载入方差的 %累积 %67表 解释的总方差通过上表我们可以看出前 5 个特征根的累计贡献率就达到了以上， 主成分分析的效果很好，故我们选取前 5 个主成分进行综合评价，得到的主成分综合评价模型：Z=+前 5 个特征根对应的主成分载荷矩阵如下：主成分X

36、2X3X4X7X8X12X13X141主成分载荷矩阵2345（部分表格，详情请见附录 1）利用上述主成分载荷矩阵， 将各指标除以各自特征根的算术平方根就可以得出各主成分的特征向量，计算结果如下：主成分X2X3X4X7X8X12X13X141主成分特征向量2345（部分表格，详情请见附录 1）最后，运用软件计算出每个指标在主成分模型下的权重，计算结果如下：主成分X2X3X4X7X8X12权重X13X14X16X18X20X21X22X23X24X27X28X29X31X33X34X35X38X39X40X41X42X43X44结合实际，对评价指标进行灵敏度分析。将其按影响度由大到小进行排序， 得

37、出结果如下：人均医院拥有数、人均卫生机构人员拥有数、人均卫生机构床位数、年空气达标率、地表水质量、城镇居民人均可支配收入。问题三结果的分析及验证本文通过主成分分析法，通过比较各指标灵敏度的大小，得出了六个对排名影响显着的三级指标，包括医疗、空气、水质、可支配收入。其中，医疗条件的三个条件位居其中， 说明评价一座城市是否宜居， 医疗占很大因素。 近年来， 由于空气污染日渐严重，因此空气质量对城市宜居也有巨大的影响。总体来说，影响最大的六个因素囊括了居民的衣食住行。问题四模型建立与求解问题四模型的建立本文在已有模型的基础上，引入了两个不确定因素：突发自然灾害，房价大幅波动率。突发自然灾害首先基于淮

38、海经济区过去十年的灾害情况， 大致拟合估计出旱灾、 洪涝、 山体滑坡、泥石流、台风、风雹灾害和低温冷冻和雪灾等常发的自然灾害的发生周期，从而得到各个突发自然灾害的发生概率。淮海经济区主要突发自然灾害评价结果如下表：灾害种类发生周期 T人员受灾直接经济生态环境受灾风险评价（天）数 R损失 Z面积 H（公顷） 值 L（万人）（万元）旱灾360洪涝、山体滑295坡、 泥石流和台风风雹灾害375低温冷冻和雪450灾R Z HT然后，在灾害风险程度的比较过程中，首先将上述每类灾害按照相同的指标赋予风险评价值。灾害的风险评价值的确定主要考虑如下几个指标，分别是灾害对人的影响程度即人员受灾，其次是灾害对社会

39、经济以及环境的影响。根据近几年来淮海经济区的自然灾害数据，对每类自然灾害进行赋权，评价方法主要是通过先通过熵权法算评价指标的权向量。下一步，对经常发生的灾害赋值为 1，几十年甚至上百年才发生一次的灾害赋值为10。对这些灾害风险评价指标赋值后，考虑各种突发自然灾害的发生概率和灾后对社会经济以及环境造成的损失，与原评价指标组合，按照粗集中简约的知识简化指标，重新计算各突发自然灾害的对城市的灾害风险权重。最后，根据风险权重和近几年的灾害数据期望值得出城市的自然灾害受损率。房价大幅波动首先， 基于淮海经济区八个城市城市过去十年的房价波动情况对城市房价波动情况进行灰色预测， 得到未来城市的房价大幅波动可

40、能性 （本文以波动高于 20%作为房价大幅波动的标准）。本文中得到的可能性值为 Pij（i,j 分别为表格的行与列），具体值如下表：计算公式：L=城宿迁市房%价大幅上涨房价大%连云港宿州商丘济宁枣庄徐州淮北%幅下降然后，以近几年的数据为基准，选取淮海经济区八个城市，对房价与其他二级指标进行相关性分析，得到房价与其他二级指标的相关系数， 取平均数，并选取产生重大影响的指标作为问题四模型建立指标 （本文选取了： 居民消费指数， 城镇居民可支配收入，人口自然增长率，失业率）。本文中得到的相关系数为 Gm，具体数据表格如下：指标居民消费城镇居民可支人口自然增长失业率配收入率%影响度 （房价%上升每个百

41、分点对该指数的影响）下一步，将 Step1,Step2 得到的数据进行运算并替换入问题 1 得到的矩阵中对每个城市每一项计算公式如下：X4 X1(1 20Gm)(P1j P2 j)(1 P1j P2 j)最后，对得到的新的数据进行标准化运算，求熵值，求权重。问题四模型的求解添加不确定因素后各城市的评分：城市宿迁评分连云港宿州商丘济宁枣庄徐州淮北问题五模型建立与求解问题五模型的建立Step1：本文对淮海经济区八个城市的数据进行归一化处理，得到一个每一个数据值都处于0,1的表格Step2: 基于 Step1 得到的表格，对 45 个二级指标进行筛选，选出每一组指标中，最大的值，即 Mi，并且保留徐

42、州市每一个二级指标的数据值为 Xi，得到一个新的表格Step3: 基于 Step2 得到的表格，求出徐州市单项指标与淮海经济区八个城市指标最大值的差异值，计算公式为：Ji Xi/Mi问题五模型的求解根据在中建立的模型， 本文选取了徐州市单项指标与淮海经济区八个城市单项最大指标差异值较大的一些量，进行分析、归纳、总结，得出了一些合理的、可执行行强的政策建议以提高徐州市的城市宜居度：(1)不要满足于自给自足，而是要增加对外销售量，出台一些城市出口促进政策，以增强与外界城市的商业联系，以期提高城市的出口率、人均消费水平，带动经济发展，提高城市宜居度；(2)加快转变城市发展方式，创造现代化的城市发展方

43、式，大力整治城市生态环境，着力于改善城市空气质量、城市水环境质量、城市绿化覆盖情况以及人均饮用水数量，以期提高城市的环境水平，提高居民健康程度，提高城市宜居性；(3)加快推进现代化农业发展模式，不以牺牲第一产业为代价而发展第二第三产业，提高农业的机械化水平，提高城市粮食自给率，提高市民生活水平，以提高城市宜居性；(4)提高城市基础设施建设水平，提高城市基础设施建设的拨款，大力进行教育事业、卫生与交通设施的建设， 提高城市每万人学校数量、 医院人均床位数以及城市道路密度，以期提高市民的生活便利度，提高城市宜居性水平；(5)出台强力政策控制房价，改变城市房价居高不下、工资与房价比低的情况，努力提高

44、城市人均住房面积，提高市民生活幸福度，提高城市宜居性水平。七、模型的评价与推广七、模型的评价与推广模型的评价问题一模型的评价1. 模型一的优点:（1）本文在正确、清晰地分析了题意的基础上，建立了综合的、较为客观的评价模型对淮海经济区八个城市宜居性状况进行评价。（2）分别从安全性、健康性、便捷性、舒适性四个方面出发对我国养老设施进行评价，并从各大权威网站获取数据，同时对数据进行必要的预处理。（3）本文将现在的评价指标得分与前两年年对应指标进行对比，较为客观的得出各个指标淮海经济区八个城市宜居性的影响大小。（4）本文在建立城市宜居性评价模型时考虑到了创新与发展项，考虑到了市民的未来发展，较为客观的

45、得到各个城市宜居性的情况。2. 模型一的缺点:（1）数据查找不完善，存在数据缺失的情况，虽采用插补法得到缺失的数据，但依然缺少一定的准确性。（2）模型的评价指标数较少，覆盖面不够广，不能够体现本模型的普遍适用性。问题三模型的评价模型三优点：由于评价指标过多， 采用主成分分析法，能够在保证保留大多数信息情况下，用少数几个指标代替。 主成分分析法的计算较为规范，克服了某些评价方法中认为确定权数的缺陷。模型三缺点：主成分的解释含义在一定程度上具有模糊性，不像原始变量的含义那么清楚、确切，这是变量降维过程中不得不付出的代价。问题四模型的评价模型四的优点:本文为了得到突发自然灾害、房价大幅波动等不确定因

46、素的概率，从各大权威统计网站查取了数据，得到了未来房价大幅波动、 自然灾害发生概率的表格，并且对自然灾害进行了分级化处理，得到的数据较为客观，可靠。本文为了较为客观的表示不确定性因素对淮海经济区八个城市城市宜居性各个的影响，分情况选取了典型年份对不确定因素与各个指标标准化处理后的相关系数进行了分析，制表，得到了灰色不确定因素对城市宜居性的影响，较为客观。本文将自然灾害分为了：旱灾、洪涝、山体滑坡、泥石流和台风、风雹灾害四类，并对各种自然灾害的影响数据进行了标准化处理后分析，得到的结果较为全面。模型四的缺点:（1）数据查找不完善，存在数据缺失的情况，虽采用插补法得到缺失的数据，但依然缺少一定的准

47、确性。（2）模型的评价指标数较少，覆盖面不够广，不能够体现本模型的普遍适用性。（3）在本文的自然灾害分类中，难免带有一些主观因素，分类不够客观。问题五模型的评价模型五的优点:（1）本文将徐州市的各项指标与淮海经济区八个城市指标的最优值进行了对比，并求出了差异度，分析过程较为直观，易懂。模型五的缺点:（1）本文对徐州市与其他城市客观差异性了解不是非常全面，提出的政策建议难免带有较大主观性。模型的推广针对研究过程中的研究与发现， 我们可以发现淮海经济区评价宜居城市的数学模型同样适用于评价宜居村镇。七、参考文献七、参考文献1姜启源等, 数学模型（第三版），北京：高等教育出版社，2003 年，50-5

48、2.2崔杰，吴昊, 误差处理的方法研究，误差理论分析，23(1):50-52,2003 年.3胡伏湘，胡希军，城市宜居性评价指标体系构建J.生态经济，2014（8）4张文忠.宜居城市的内涵即评价指标体系探讨J.城市规划学科，2007（3）5张磊.新农村建设评价体系研究J.经济纵横，2009（7）6赵之枫.乡村宜居性建设的构想J.生态经济，2001（5）7江苏、山东、安徽、河南四省 2013-2015 的统计年鉴八、附录八、附录附录清单附录 1：问题一数据指标表格附录 2：求解问题二，问题四的 matlab 程序及命令附录 3：求解问题一的中间数据附录 4：问题二的完整数据结果附录 5：问题三的

49、 SPSS 过程截图及结果表格附录 6：问题五的分析结果表格附录正文附录 1：问题一数据指标表格经济居民消费价格总粮食总产油料总产量增长率 指数（以上年为量增长率100）%第三产业在所有产业中所占城镇失业率第三产业增长率的比重%环境水环境水进出口总额增长 城率%122222122NA122232222222222城市噪声环境质量达标率 地表城市绿化覆盖率日平均气温在 15-25饮用水地下水城市区域昼间噪声（A） 城市交通干线昼间噪声（A）达标率100%100%100%100%100%100%160151133128136132类类类类类类类类类类类类类类类类类类类类类类类类类类类类类类类类类类

50、类类类类类类类类类类类类类类100%100%100%100%100%100%100%100%100%100%100%100%100%100%100%100%100%100%122132117132128152137124118142129128150139123150126125生活便利度城市宿迁宿迁宿迁连云港连云港连云港宿州宿州宿州商丘商丘商丘济宁济宁济宁枣庄枣庄年份20132014201520132014201520132014201520132014201520132014201520132014城市交通网络运输能力（万人）66236734596753815433470072107719

51、51711631488188262103825471371577593260公路密度枣庄徐州徐州徐州淮北淮北淮北20152013201420152013201420152431148141506313347923534711841城市宿迁宿迁宿迁连云港连云港连云港宿州宿州宿州商丘商丘商丘济宁济宁济宁枣庄枣庄枣庄徐州徐州徐州淮北淮北淮北年份201320142015201320142015201320142015201320142015201320142015201320142015201320142015201320142015卫生机构人卫生机构床位 人口（万人）各级各类学校学校数（所/万人）卫

52、生机构数员数数（所）65017411776266036781415342713741381185512644904642345246224261774189518806631686922082375451046196817057182108023862261792305617400199002141636563464825125351827167962005821023397004700051700100431044210796190852024023290190351660017800194713276646723754224436714732171861915035500462004790

53、0102791171311852824000000001.000000卫生机2.000000城市基础设施附录 2：求解问题二，问题四的 matlab 程序及命令程序function W,S=shang(A,f)n,m=size(A);k=1/log(n);X=zeros(n,m);for j=1:m temp=sort(A(:,j); if(f(j) for i=1:nX(i,j)=(A(i,j)-temp(1)/(temp(n)-temp(1);%(A(i,j)-min(A(:,j)/(max(A(:,j)-min(A(:,j);+1 if(X(i,j)=1) X(i,j)=; elseif

54、(X(i,j)=0) X(i,j)=; end end else for i=1:nX(i,j)=(temp(n)-A(i,j)/(temp(n)-temp(1);%(max(A(:,j)-A(i,j)/(max(A(:,j)-min(A(:,j);%+1 if(X(i,j)=1) X(i,j)=; elseif(X(i,j)=0) X(i,j)=; end end endendfor i=1:n%j=1:m for j=1:m%i=1:n R(i,j)=X(i,j)/sum(X(:,j); endendfor j=1:m W(j)=g(j)/sum(g);endS=W*X;问题二命令：x=2

55、24100 462710.133581 2426261792329053. 47.22133596722233161100 5832400 45132650 2616230561903551. 42.221324700257281813953423214161947132.12236302713382280365633276640.1323572294100 461711.558518274436762. 53.3311717761880117517115236782660152826211834272580118371527877113100 405013.400210231915054.

56、49.225792283100 649510.394517004790059.226219101100 5045398 421185249. 54.1425690s,w=shang(x)问题四命令：x=附录附录 3 3：求解问题一的中间数据：求解问题一的中间数据对数据归一化求标准值结果：归一化后宿迁连云0港.宿州0.商丘0.济宁0.枣庄0.徐州2128212312518410.0 0. .1281381196024311231264478013347125464 718 107960.0.0 0. .0.0.0.淮北附录 4：问题二的完整数据结果问题二各项指标权重：对公生人人共医活平均均 城有

57、安安院垃均自淡人工 市人燃线居全全、圾平工来水均业 噪口气电水民指 事教卫无均资水量住产 声（第普视环城消标 件育生害房与普（房值 环万人人三及普日卫卫境 日市费处人院化价房及立面（ 境人卫均均产城率及平各生生平交价人理均数处价率方积万 质）公学生公业镇率均级机机均通人粮油格进 口成量理比米元 量 城园卫校机园在 第居气各卫构构气网均食料总出 自功率）市绿卫生数构空绿所城 三民温类生人床温络公G总总指口 然率绿地生机（数气地有镇 产人在学机员位在运路D产产数总 增达面产失 业均城化面机构所（1校构（1输密P量量（额 长标积业业 增可城市覆积构床/ 所5数人人个5能度（增增以增 率率（中率 长支

58、市交盖（数位万/-（员/ /-力元长长上长 （数人万平所 率配区通率平2所万万2（）率率年率 饮方）人方占收域干5）人人5万为）用米）的入昼线地地米）人1水）比间昼表下）0达）重噪间水水0标声噪）率（声A （）A）权0重得分计算结果：城市宿迁宿州商丘济宁枣庄徐州淮北连云港最终分数0.附录 5：问题三的 SPSS 过程截图及结果表格主成分1X20X28X21X3X43X22X35X18X13X2X8X12X33X23X38X29X34X24X14X42X31X16X41X40X4X39X7X44X272主成分载荷矩阵345主成分1X20X28X21X3X43X22X35X18X13X2X8X12

59、X33X23X38X29X34X24X14X42X31X16X41X40X4X39X7X44X27-0.0.-0.-0.2-0.30.0.主成分特征向提取方法 :主成份。a 已提取了 7个成份。附录 6：问题五的分析结果表格徐州市各项指标与最好指标差值（已降序排序）：居民第城消城人人三卫卫 市费对卫市均均产城人日 日生城价公 生 人生 交人各学生交公公业镇均平平 平机医市格共 活 均机 通进粮口级油人校机通 城 园卫 园在第居有淡均 人均 均安料人 构院 自区总卫安 垃 工构 饮网出食自各均数空城构燃干 市 绿生 绿所三民线卫公水工 均气 气全总口 人、 来域指生全 圾 平 业床 用络口总然类

60、G（地 地 气镇数气线 绿 地机 地有产人电生城路量资 住温 温教产（ 员卫 水昼数机事 无 均 产位 水运总产增学D所下 表 达失（视普昼 化 面构业均机面产市量长校密（与 房在 在育量万 （生 普间（构件 害 房 值（ 达输额P/ 水 水 标业所及间 覆 积床 积业增可普构度立房 面15 15人增人 人院 及噪以人处 化 价 （个 标能增增率数（万率率/率噪 盖 （位长支及数（中方价 积-2 -2均长） / 数 率声上员理 处 万/率力元长长（人万声 率 平数 平所率配率米比5 5万量（年成 理 元万 （率率所率）人（方方占收） 人A）率 ）人 万）为功）A 米米的入） 人1率）比）0重0）差0.00000000000000 0值