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

《matlab湖泊水质富营养化的综合评价》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、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表二湖泊水质评价标准评价参数极贫营养贫营养中营养富营养极富营养总磷<1423110>660耗氧量<0.090

2、.361.807.10>27.1透明度>37122.40.55<0.17总氨<0.020.060.311.20>4.6要求（1）试用以上数据，分析总磷，耗氧量，透明度，总氨这4个指标对湖泊水质评价富营养化的作用（2）对这5个湖泊的水质综合评价，确定水质等级。2问题的分析和模型的建立在进行综合评价之前，首先要对评价的指标进行分析。通常评价指标分成效益型，成本型和固定型指标。效益型指标是指那些数值越大影响力越大的统计指标（也称正向型指标）；成本型指标是指数值越小越好的指标（也称逆向型指标）；而固定型指标是指数值越接近于某个常数越好的指标（也称适度型指标）

3、。如果每个评价指标的属性不一样，则在综合评价时就容易发生偏差，必须先对各评价指标统一属性。建模步骤（1）建立无量纲化实测数据矩阵和评价标准矩阵4实测数据矩阵等级标准矩阵，然后建立无量纲化实测数据矩阵和无量纲化等级标准矩阵，其中利用Matlab，我们得到X=[13010.300.352.76;10510.700.402.0;201.44.50.22;306.260.251.67;2010.130.500.23];Y=[1423110660;0.090.361.807.1027.1;37122.40.550.17;0.020.060.311.204.6]

4、;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.00000.96260.71431.00000.80771.00000.62500.72460.15380.13080.05560.07970.23080.58501.00000.60510.15380.94670.50000.08

5、33B=0.00150.00610.03480.16671.00000.00330.01330.06640.26201.00000.00460.01420.07080.30911.00000.00430.01300.06740.26091.00004（2）计算各评价指标的权重首先计算矩阵B的各行向量的均值和标准差然后计算变异系数，最后对变异系数归一化得到各指标的权重为w=0.27670.24440.23470.2442程序命令：b=B';t=std(b)./mean(b);w=t/sum(t);（3）建立各湖泊水质的综合评价模型通常我们可以利用向量之

6、间的距离来衡量两个向量之间的接近程度，在Matlab中，有以下的函数命令来计算向量之间的距离；Disp(w,p):计算w中的每个行向量和p中每个列向量之间的欧式距离；Mandist(w,p):绝对值距离计算A中各行向量到B中各列向量之间的欧氏距离若，则第个湖泊属于第级计算A中各行向量到B中各列向量之间的绝对值距离若，则第个湖泊属于第级结果及程序：jd=dist(A,B)mjd=mandist(A,B)jd=1.84721.83121.73741.37690.28811.59591.57981.48591.12710.503440.21850.2045

7、0.13670.33831.79171.32011.30381.20820.83920.95911.07931.06500.98670.73281.3450mjd=3.66313.63033.43742.67830.32313.14363.11082.91782.15870.84270.40620.37340.21100.57873.58002.40712.37432.18141.42231.57911.67011.63741.44441.06602.3161所以，从上面的计算可以看出，尽管欧氏距离和绝对值距离的意义完全不一样，但对湖泊水质的评价等级

8、是一样的，这表明了我们方法的稳定性。问题：利用其他的方法来计算个评价指标的权重，比较结果是否有所改变。4