聚类分析 熵方法加权 方差越大,权重越大(论文资料).docVIP

李 郭 0809501024 应数081 实验目的：学习利用距离函数建模的方法，掌握客观性权重的变异系数以及综合评价的基本方法，熟练掌握MATLAB处理矩阵的各种方法。 实验内容及要求： 问题:近年来我国淡水水质富营养化的污染日趋严重，如何对湖泊水质的富营养化进行综合评价是摆在我们面前的一项重要任务。表1和表2分别是我国5个湖泊实测数据和水质评价标准。 表3.4 全国5个主要湖泊品和技术参数的实测数据 　总磷 　（mg/L） 　耗氧量 　（mg/L） 透明度 （m） 总氮 （mg/L） 杭州西湖 130 10.30 0.35 2.76 武汉东湖 105 10.70 0.40 2.0 青海湖 20 1.4 4.5 0.22 巢湖 30 6.26 0.25 1.67 滇池 20 10.13 0.50 0.23 表3.5 湖泊水质评价标准 评价指标 极贫营养 贫营养 中营养 富营养 极富营养 总磷 1 4 23 110 660 耗氧量 0.09 0.36 1.80 7.10 27.1 透明度 37 12 2.4 0.55 0.17 总氮 0.02 0.06 0.31 1.20 4.6 是利用以上数据，分析总磷，耗氧量，透明度和总氮这4种指标对湖泊水质富营养化所起作用。 对上述5个湖泊的水质进行综合评估，确定水质等级。 问题的分析 在进行综合评价之前，首先要对评价的指标进行分析。通常评价指标分为效益型，成本型和固定型指标（另外还有区间型指标，偏离型指标等等）。在进行综合评定之前首先对各评定指标同一属性。由生物学知识可知，本题除了透明度指标为成本型指标以外，其余指标都为效益型指标。 在指标清楚之后，我们还要进行赋权和标准化。在初次建立模型的时候我们选择了：变异系数赋权法和线性比例变换法的标准化方法 模型的建立 根据表3.4和表3.5，我们得到实测数据矩阵X=（）和等级标准矩阵Y=（） 建立无量纲化实测数据矩阵：A=（）和无量纲化等级标准矩阵B=（）。 其中我们用线性比例变换的标准化方法得到： 利用MATLAB计算得 B= 计算评价指标的权重 首先计算矩阵B的各行向量的均值与标准差， 然后计算变异系数 （i=1，2，3，4） 利用MATLAB计算得：W=[0.2767 0.2444 0.2347 0.2442] W表示个指标的权重。根据权重的大小，即可说明总磷，耗氧量，透明度和总氮4种指标对湖泊水质富营养化所起作用。由上可知个指标的作用很接近，比较而言总磷所起作用最大，耗氧量，总氮次之，透明度的额作用最小。 建立个湖泊水质的综合评价模型，利用欧氏距离和绝对值距离进行建模 计算A中各行向量到B中各列向量的欧氏距离 若，则第i个湖泊属于第k级（i=1,2,3,4,5）。 计算A中各行向量到B中各列向量的绝对值距离 若，则第i个湖泊属于第k级（i=1，2，3，4，5）。 计算结果如下表所示：表3.6 欧氏距离判别表 级别 杭州西湖 1.8473 1.8312 1.7374 1.3768 0.2881 5 武汉东湖 1.5960 1.5798 1.4859 1.1271 0.5034 5 青海湖 0.2185 0.2044 0.1367 0.3383 1.7918 3 果湖 1.3201 1.3037 1.2083 0.8392 0.9591 4 滇池 1.0793 1.0649 0.9867 0.7328 1.3451 4 表3.7 绝对值距离判别表 级别 杭州西湖 3.6632 3.6303 3.4375 2.6782 0.3231 5 武汉东湖 3.1436 3.1107 2.9179 2.1586 0.8427 5 青海湖 0.4062 0.3733 0.2109 0.5788 3.5801 3 果湖 2.4072 2.3743 2.1815 1.4222 1.5791 4 滇池 1.6701 1.6372 1.4444 1.0661 2.3162 4 从上面的计算可知，尽管欧氏距离与绝对值距离意义不同，但是对个湖泊水质的富营养化评价等级是一样的，表明我们给出的方法具有稳定性。 附录：用MATLAB编程 X=[130 10.3 0.35 2.76;105 10.7 0.4 2.0;20 1.4 4.5 0.22;30 6.26 0.25 1.67;20 10.13 0.5 0.23]; %输入实测数据矩阵 A1=X(:,1)./max(X(:,1));A2=X(:,2)./ max(X(:,2));A3= min(X(:,3))./X(:,3); A4=X(:,4)./ m