聚类分析在数学建模中的应用-课件.pptxVIP

聚类分析在数学建模中的应用 hjyahead@ 2011大学生数学建模夏令营B题 水资源短缺风险，泛指在特定的时空环境条件下，由于来水和用水两方面存在不确定性，使区域水资源系统发生供水短缺的可能性以及由此产生的损失。 北京是世界上水资源严重缺乏的大都市之一，其人均水资源占有量不足300m3，为全国人均的1/8，世界人均的1/30，属重度缺水地区，附表中所列的数据给出了1979年至2000年北京市水资源短缺的状况。北京市水资源短缺已经成为影响和制约首都社会和经济发展的主要因素。如何对水资源风险的主要因子进行识别，对风险造成的危害等级进行划分，对不同风险因子采取相应的有效措施规避风险或减少其造成的危害，这对社会经济的稳定、可持续发展战略的实施具有重要的意义。 2011大学生数学建模夏令营B题 考虑以下问题： 1 评价判定北京市水资源短缺风险的主要风险因子； 2 建立一个数学模型对北京市水资源短缺风险进行综合评价，作出风险等级划分并陈述理由。 3 对北京市未来两年水资源的短缺风险进行预测，并提出应对措施。 4 以北京市水行政主管部门为报告对象，写一份建议报告。 年份 缺水量（x亿立方米） 1979 4.69 1980 24.54 1981 24.11 1982 10.62 1983 12.86 1984 0.74 1985 -6.29 1986 9.52 1987 -7.71 1988 3.25 1989 23.09 1990 5.26 1991 -0.26 1992 23.99 1993 25.55 1994 0.45 1995 14.54 1996 -5.86 1997 18.07 1998 2.73 1999 27.49 2000 23.54 2001 19.7 2002 18.5 2003 17.4 2004 13.2 2005 11.3 2006 9.8 2007 11 2008 0.9 定义模糊集合 u 表示风险程度 x 表示缺水量 构造隶属函数 将1979-2008年的缺水量换算成对应的风险程度用0到1区间的数字表示； 表示缺水量的最小值， 表示缺水量的最大值 年份 缺水量（x亿立方米） 风险程度（u（x）） 1979 4.69 0.1241 1980 24.54 0.8394 1981 24.11 0.8172 1982 10.62 0.2712 1983 12.86 0.3415 …… …… ……. 2003 17.4 0.5089 2004 13.2 0.3529 2005 11.3 0.2917 2006 9.8 0.2474 2007 11 0.2825 2008 0.9 0.0598 通过以上两图的比较可以说明隶属度真正反映了30年缺水的状况，很好的反应了风险的模糊性。 通过对北京30年间风险程度的不同进行分类采用系统聚类的方法 开始时每个对象自成一类， 然后每次将最相似的两类合并，合并后重新计算新类与其他类的距离或相似度， 这一过程一直继续，一直到所有对象归为一类为止，并类的过程用一张普类系谱图表示。 层次（系统）聚类法（Hierarchical Clustering） ——计算步骤 ①计算n个样本两两间的距离{dij}，记D ②构造n个类，每个类只包含一个样本； ③合并距离最近的两类为一新类； ④计算新类与当前各类的距离；若类的个数等于1，转到5）；否则回3）； ⑤画聚类图； ⑥决定类的个数和类； 由spss软件求对北京30年风险的分类树状图 对于这30年的数据共可分为30组、5组、4组、2组、1组。 根据系统聚类的结果，然后采用模糊聚类的方法分别分出5组、4组、2组 鉴于30组和1组情况太过极端对风险的分类太过分散和统一 舍去这两种情况，分别采用模糊聚类的方法对风险分成5组、4组和2组 模糊聚类的步骤 第一步：建立数据观测矩阵sample； sample=[4.69 24.54 24.11 10.62 12.86 0.74 -6.29 9.52 -7.71 3.25 23.09 5.26 -0.26 23.99 25.55 0.45 14.54 -5.86 18.07 2.73 27.49 23.54 19.7 18.5 17.4 13.2 11.3 9.8 11 0.9]; 第二步：标定模糊相似矩阵data=(rij),矩阵元素rij 表示向量xi和xj的相似程度，由隶属函数W可求得data； data=[ 0.1241 0.8394 0.8172 0.2712 0.3415 0.0576 0.0016 0.2396 0