美国数学建模竞赛论文C题论文Word格式文档下载.docx

1、标准化后的健康值相关系数矩阵CiKiB，Ya聚类系数节点i的邻接点数据矩阵健康指数的临界值注：部分未说明的符号在模型会给与介绍。四、问题分析根据对题目的分析，我们知道题目要求我们考虑地球上各系统间的交叉影响的复杂性及地区条件和全球系统之间的影响，建立一个动态的全球网络模型，明确节点的实体和节点间的链接特性，根据地区条件变化引起全球变化提供预警，并能预测未来地球健康状况的变化，为政策制定者提供有效的政策。针对问题，以国家为节点，根据生态学选取地球健康因子,在世界银行数据库中查得所需的数据。具体分析思路如下：（1）根据所查的数据求出各个节点及世界的健康指数。可以通过主成分分析法对所取20年的数据进

2、行综合评价；（2）接着，求解各节点间的影响及内部的相互作用。在主成分分析的基础上，因为数据量较大，在21个国家中挑选出8个国家来构成“小世界”作为真实世界的模拟。对于选出来的8个国家，将他们构成有向复杂网络来考虑各节点间的影响及内部的相互作用。（3）对未来健康状况的预测。根据全球1991年到2010年的健康指数，以20年的数据为基础，采用GM（1,1）进行预测，并提出预警情况。五、模型的建立与求解5.1 节点和指标的选取（1）节点的选取考虑地理位置和国家发展状况，选取了21个典型的国家，这21个国家可以代表全球。选取的21个国家：Argentina，Australia，Brazil，Canad

3、a，China，Egypt，Arab Rep.，France，Germany，Iran，Islamic Rep.，Israel，Italy，Japan，Kazakhstan，Kenya ，Morocco，New Zealand，Russian Federation，Saudi，Arabia，South Africa United，Kingdom，United States。（2）指标的选取1,2地球是一个大生态系统，我们根据文献3选取了11个地球健康的指标，如下表所示：表1 地球健康指标地球健康社会因素城镇人口比重（%）人均耕地面积（公顷）高等院校入学率（%）经济因素人均GDP（美元）每千美元

4、的能源消耗量（千克石油当量）环境因素CO2排放量（人均公吨数）颗粒污染程度（每立方米微克）人均能源使用量（千克石油当量）人的因素人口密度（每平方公里）人口增长（%）五岁以下婴儿死亡率（%）选取21个国家11个指标的时间区间:1991-2010年；数据的来源:世界银行数据库45.2 主成分分析模型的建立与求解5.2.1主成分分析法的步骤1）对原始数据进行标准化处理假设进行主成分分析的指标变量有m个：，共有n个评价对象，第i个评价对象的第 j 个指标的取值为。将各指标值转换成标准化指标,其中， 即，为第 j 个指标的样本均值和样本标准差。对应地，称为标准化指标变量。2）计算相关系数矩阵 R式中=1

5、 ， =，是第i个指标与第 j 个指标的相关系数。3）计算特征值和特征向量计算相关系数矩阵R的特征值，及对应的特征向量，其中，由特征向量组成m个新的指标变量式中，是第1主成分，是第2主成分，是第 m 主成分。4）选择 p（）个主成分，计算综合评价值计算特征值的信息贡献率和累积贡献率。称为主成分的信息贡献率；为主成分的累积贡献率，当接近于1时，则选择前p个指标变量,作为p个主成分，代替原来m个指标变量，从而可对p个主成分进行综合分析。计算综合得分定义：Z为健康指数，表示健康状况的综合评价值5.2.2模型的求解以2000年的数据为例，我们选取5个主成分即令p=5，利用MATLAB软件对11个评价指

6、标进行主成分分析,相关系数矩阵的前几个特征根及其贡献率如表2。表2主成分分析结果序号特征根贡献率累积贡献率15.221847.470522.199519.995867.466331.267511.522778.989040.73696.699185.688050.69866.350592.038560.33963.087395.1258可以看出,前四个特征根的累计贡献率就达到85%以上,主成分分析效果很好。下面我们选取五个主成分（累计贡献率就达到92%）进行综合评价。前五个特征根对应的特征向量见表3。表3 标准化变量的前4个主成分对应的特征向量第1特征向量第2特征向量第3特征向量第4特征向量第

7、5特征向量0.37750.35660.3763-0.1293-0.2234-0.2237-0.04960.20210.4985-0.05980.04200.21460.1964-0.21260.55920.3196-0.44430.32210.4191-0.03340.1640-0.0233-0.0179-0.3379-0.6483-0.26540.5122-0.04110.4827-0.29197-0.07680.1891-0.33170.23330.073880.0144-0.3248-0.15670.07670.209290.44090.4558-0.15740.10300.260710

8、0.4813-0.00750.23450.27860.030211-0.41730.12230.6774-0.13470.0955竖列的1,2，11表示标准化后的11个数据，即，由此得到五个主成分：从主成分的系数可以看出,第一主成分主要反映了人均GDP、人口密度、高等院校入学率、人均耕地面积的信息；第二主成分主要反映了人均GDP、每千美元的能源消耗量、CO2排放量、人口增长、人口密度；第三主成分主要反映了人均GDP、每千美元的能源消耗量、五岁以下婴儿死亡率、人均耕地面积；第四主成分主要反映了城镇人口比重、每千美元的能源消耗量、颗粒污染程度、CO2排放量；第五主成分主要反映了人均能源使用量和颗粒

9、污染程度。综上可以看出人均GDP、每千美元的能源消耗量、五岁以下婴儿死亡率的影响相对较大一些。把各节点的原始11个指标的标准化数据代入五个主成分的表达式,就可以得到各地区的五个主成分值。分别以五个主成分的贡献率为权重,构建主成分综合评价模型:（1）把各节点的五个主成分值代入式（1）,可以得到各地区高教发展水平的综合评价值以及排序结果如表4表4 2000年的21个国家和世界的健康指数节点USAAUSCANNZLRUSJPN排名2.0542.0391.9680.7240.5390.529GBRFRADEUISRSAUITA120.4950.4290.4180.2910.1950.075ARGKAZ

10、ZAFBRAIRNWLD131415161718-0.075-0.094-0.558-0.788-0.85685-1.11022MAREGYCHNKEN19202122-1.34431-1.48728-1.49224-1.94967上述的USA等为各节点的缩写。由表4以及附录中表5的数据可以看出美国、澳大利亚、加拿大的健康指数在20年内都很高，说明这几个国家的健康状况较好；Morocco、Egypt Arab Rep. China、Kenya等国家的健康指数在20年后都很低，说明这几个国家的健康状况较差。同时可以世界平均水平的健康指数也较低，说明地球面临严峻的健康问题。运用同样的方法，可以求出

11、1991-2010的各节点的综合评价，所得的综合值Z进行如下标准化：标准化后的1991-2010年的21个国家和世界的健康指数见附录。5.3有向复杂网络模型的建立与求解针对上文已经得到的21个国家的生态健康指数值，考虑在同一个地球生物圈内各国的生态健康是相互作用的，不能独立于周边的国家的影响，而且随着全球化进程的不断加快，全球间的各个国家势必会联系越来越频繁。我们在考虑地球生态健康的时候，将各个国家看成一个个节点，不同节点由网络边连接起来，构成整体的生态系统网络。在整个生态系统大网络中各个节点间是否作用，或作用的强弱程度我们是通过对选取节点的健康值进行逐步回归得到的。本文运用逐步回归来分析来确

12、定有向复杂网络的边和对应权重。由于涉及国家数目过多，根据地理特性分别在各大洲共选取8个国家：澳大利亚、巴西、中国、日本、俄罗斯联邦、南非、英国、美国。依次将其中一个国家的生态健康值当做y值，其余各国指标为x值来进行逐步回归观察各国家节点在网络中是否关联及相应权重。5.3.1 逐步回归逐步回归的基本思想是：对全部因子按其对y影响程度大小（偏回归平方的大小），从大到小地依次逐个地引入回归方程，并随时对回归方程当时所含的全部变量进行检验，看其是否仍然显著，如不显著就将其剔除，知道回归方程中所含的所有变量对y的作用都显著是，才考虑引入新的变量。再在剩下的未选因子中，选出对y作用最大者，检验其显著性，显

13、著着，引入方程，不显著，则不引入。直到最后再没有显著因子可以引入，也没有不显著的变量需要剔除为止。逐步回归分析时在考虑的全部自变量中按其对y的贡献程度大小，由大到小地逐个引入回归方程，而对那些对y作用不显著的变量可能是中不被引入回归方程。另外，已被引入回归方程的变量在引入新变量进行F检验后失去重要性时，需要从回归方程中剔除出去。逐步回归算法步骤：Step 1 计算变量均值：差平方和：记各自的标准化变量为：Step 2 计算的相关系数矩阵Step 3 设已经选上了个变量：且互不相同，经过变换后为：对逐一计算标准化变量的偏回归平方和：记：作F检验：对给定的显著性水平，拒绝域为：Step 4 转到S

14、tep 3 循环，直至最终选上了个变量，且则对应的回归方程为：通过代数运算可得：逐步回归是用Matlab统计工具箱中的命令stepwise，它提供了一个交互式画面，通过这个工具你可以自由地选择变量，进行统计分析，其方法通常为：stepwise（x,y,inmodel,alpha）其中x是自变量数据，y是因变量数据，分别为 nm和 n1矩阵，inmodel是矩阵x的列数的指标，给出初始模型中包括的子集（缺省时设定为空），alpha为显著性水平。Stepwise Regression 窗口，显示回归系数及其置信区间，和其它一些统计量的信息。绿色表明在模型中的变量，红色表明从模型中移去的变量。在这个

15、窗口中有Export按钮，点击Export产生一个菜单，表明了要传送给Matlab工作区的参数，它们给出了统计计算的一些结果。下面以美国的指标为y值，其他国家的指标为x值进行逐步回归的结果：，由此认为美国和巴西，中国，日本，俄罗斯之间的节点存在边，各边权重分别为-0.4606,-0.4921,0.5677,0.191.具体运算界面见图1图1 逐步回归结果图（USA）同样方法可以得到其他各国节点在网络中的边和权重。各国数据如表6表 6网络中国家节点边和权重USA 常数项-0.4606-0.49210.56770.1910.5118-0.7831-0.56252.3883-0.4973-0.340

16、70.3206-0.4133-0.66250.7180.3529-0.6411-0.8937-0.49510.78070.8548-0.50610.43720.69330.2677-0.49530.2168-0.66110.36522.133-0.97610.45750.2597-1.3095-1.5355-0.11-0.65190.364-0.4864-0.6103经过matlab运算后的逐步回归结果整理后得到有向复杂网络图如图2，其中各节点的作用的方向和作用权重如图：图2 八国家的复杂网络图及权重5.3.2有向复杂网络模型计算分析随着复杂网络研究的不断深入，在原有的基于无权网络的ER随机图

17、模型、WS小世界模型和无标度网络模型BA模型基础上，文献9给出了一个由边和节点权值驱动网络动态演化的BBV模型。所谓的BBV网络模型6,7是2004年Barrat Barthelemy和vesPignani提出的一个无向加权的网络演化模型。首先定义一个有向加权网络:设图DWG=（V，E）表示一个复杂网络。其中V表示节点集，E表示边集。总节点数一记为N，边数记为M，E中每条边都有V中的一对点与之相对应。任意两点i和j对应的边（i,j）和（j,i）不同。入强度：出强度：其中所有Vi的所有邻接点Aij表示i和j之间有连接的边，且边的指向性是节点i指向节点j。有向复杂网络的平均最短路径长度：平均距离定

18、义：实际中网络的平均最短距离：在有向加权网络中，i的度有出入度之分。因此，需要分别考虑邻接点与i之间的关系。设节点i的邻接点为Ki，其中Kin，表示Ki作为源节点与i相连的节点,Kout表示Ki，作为目的节点与i相连的节点。在有向加权网络的分析中，为了更好的统一属性，主要考虑以i为源节点的相邻节点。设Ki个节点之间的实际边数为Ei，则一记节点i的聚类系数为Ci:网络的聚类系数记为：各表现网络节点的基本参数值，入强度，出强度，聚类系数，最短路径和权值的数据如表7所示：表7 网络节点的基本参数值入强度出强度最短路径权值美国1.71141.94050.1250.100027澳大利亚1.34560.6

19、0530.714290.078647巴西1.15860.46060.067717中国3.26823.51230.196428570.191018日本3.27213.08890.191246俄罗斯1.73842.09280.160714290.101605南非3.00281.92990.142857140.175506英国1.61233.4590.094235总计17.109417.08932.3928657189其中本复杂网络的实际网络平均最短距离：显然，的值是较小的，平均各节点仅需要1.5893的最短距离就可以建立联系，这点体现了复杂网络里的“小世界”特性。网络聚类系数：值在0,1间，更接近

20、与0，说明网络的的聚集程度较差，凝聚力较差。所有国家节点中，中国和日本的权值最大，对于整个系统健康状况有较大的影响。原因可能是这两个国家分别在生活用品和高科技产品有很重要的地位和交流的中心。对有向复杂网络的基本参数的讨论后，决定通过网络边上的权值来得到网络中各节点的权值。接下来用网络反馈加权计算八个国家组成的小生态系统的健康来和没有通过网络反馈而直接得到的综合值进行比较。因此，通过网络可以得到八个国家组成的小世界生态健康指标的总值：表8 网络加权得分和直接综合得分的比较时间网络反馈加权得分直接综合得分19910.3644471840.22931450220010.3978294660.2594

21、0183319920.3524591250.2157808820020.404084860.25595221319930.372457210.25221024520030.391744830.24429583819940.3739176070.25217372820040.3754495430.22110073919950.36951620.25194543620050.3821708670.23491205119960.3681700970.23782021220060.3408987020.20789179519970.3850665690.26287513520070.377473801

22、0.23769735219980.3869636680.28717443820080.3484959810.22216428319990.4143267660.29372693420090.3607179950.23795984720000.4130907510.28499716120100.3792223610.237210636图3 网络预测值与直接加权值的对比 图4 修正后模型与直接预测模型的对比由两者的对比可以观察出网络加权健康得分教直接预测得分高，根据两者特性和趋势对网络加权计算健康的分的模型做如下修正：修正模型后与直接预测值比较：见上图4通过图像可以看出二者总体的趋势和得分差距并不大，可以看出两种综合评价体系是相对合理的。由于网络加权评价模型考虑了网络中不同节点影响和反馈，因此相信有向复杂网络模型得到的评级得分时更具说明性的，在条件容许的条件下，把全球200多个国家全作为网络节点，运用有向复杂网络来考虑他们间的影响和联系来判断地球总体生态健康水平会更合理可靠。本模型有很大的推广和扩大价值。5.3