matlab湖泊水质富营养化的综合评价.doc

1 问题近年来我国淡水湖水质富营养化的污染日益严重，如何对湖泊水质的富营养化进行综合评价与治理是摆在我们面前的任务，下面两个表格分别为我国5个湖泊的实测数据和湖泊水质评价标准。表一 全国五个主要湖泊评价参数的实测数据总磷（mg/L）耗氧量(mg/L)透明度(m)总氨(mg/L)杭州西湖13010.300.352.76武汉东湖10510.700.402.0青海湖201.44.50.22巢湖306.260.251.67滇池2010.130.500.23 表二 湖泊水质评价标准评价参数极贫营养贫营养中营养富营养极富营养总磷660耗氧量27.1透明度37122.40.550.17总氨4.6要求（1）试用以上数据，分析总磷，耗氧量，透明度，总氨这4个指标对湖泊水质评价富营养化的作用（2）对这5个湖泊的水质综合评价，确定水质等级。2 问题的分析和模型的建立在进行综合评价之前，首先要对评价的指标进行分析。通常评价指标分成效益型，成本型和固定型指标。效益型指标是指那些数值越大影响力越大的统计指标（也称正向型指标）；成本型指标是指数值越小越好的指标（也称逆向型指标）；而固定型指标是指数值越接近于某个常数越好的指标（也称适度型指标）。如果每个评价指标的属性不一样，则在综合评价时就容易发生偏差，必须先对各评价指标统一属性。建模步骤（1） 建立无量纲化实测数据矩阵和评价标准矩阵实测数据矩阵等级标准矩阵，然后建立无量纲化实测数据矩阵和无量纲化等级标准矩阵，其中 利用Matlab，我们得到X=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;Y=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;B1=Y(1,:)./660;B2=Y(2,:)./27.1;B3=0.17./Y(3,:);B4=Y(4,:)./4.6;B=B1;B2;B3;B4;A1=X(:,1)./130;A2=X(:,2)./10.7;A3=0.25./X(:,3);A4=X(:,4)./2.76;A=A1,A2,A3,A4;A = 1.0000 0.9626 0.7143 1.0000 0.8077 1.0000 0.6250 0.7246 0.1538 0.1308 0.0556 0.0797 0.2308 0.5850 1.0000 0.6051 0.1538 0.9467 0.5000 0.0833B = 0.0015 0.0061 0.0348 0.1667 1.0000 0.0033 0.0133 0.0664 0.2620 1.0000 0.0046 0.0142 0.0708 0.3091 1.00000.0043 0.0130 0.0674 0.2609 1.0000（2）计算各评价指标的权重首先计算矩阵B的各行向量的均值和标准差然后计算变异系数，最后对变异系数归一化得到各指标的权重为w =0.2767 0.2444 0.2347 0.2442程序命令：b=B;t=std(b)./mean(b);w=t/sum(t);（3）建立各湖泊水质的综合评价模型通常我们可以利用向量之间的距离来衡量两个向量之间的接近程度，在Matlab中，有以下的函数命令来计算向量之间的距离；Disp(w,p): 计算w中的每个行向量和p中每个列向量之间的欧式距离；Mandist(w,p):绝对值距离计算A中各行向量到B中各列向量之间的欧氏距离若，则第个湖泊属于第级计算A中各行向量到B中各列向量之间的绝对值距离若，则第个湖泊属于第级结果及程序：jd=dist(A,B)mjd=mandist(A,B)jd = 1.8472 1.8312 1.7374 1.3769 0.2881 1.5959 1.5798 1.4859 1.1271 0.5034 0.2185 0.2045 0.1367 0.3383 1.7917 1.3201 1.3038 1.2082 0.8392 0.95911.0793 1.0650 0.9867 0.7328 1.3450mjd = 3.6631 3.6303 3.4374 2.6783 0.3231 3.1436 3.1108 2.9178 2.1587 0.8427 0.4062 0.3734 0.2110 0.5787 3.5800 2.4071 2.3743 2.1814 1.4223 1.5791 1.6701 1.6374 1.4444 1.0660 2