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

1问题近年来我国淡水湖水质富营养化的污染日益严重，如何对湖泊水质的富营养化进行综合评价与治理是摆在我们面前的任务，下面两个表格分别为我国5个湖泊的实测数据和湖泊水质评价标准。表二湖泊水质评价标准评价参数极贫营养贫营养中营养富营养极富营养总磷14耗氧量66027.1透明度37122.40.55总氨0.020.060.311.204.6要求（1）试用以上数据，分析总磷，耗氧量，透明度，总氨这营养4个指标对湖泊水质评价富化的作用（2）对这5个湖泊的水质综合评价，确定水质等级。表一全国五个主要湖泊评价参数的实测数据总磷（mg/L）耗氧量（mg/L）透明度（m）总氨(mg/L)杭州西湖13010.300.352.76武汉东湖10510.700.402.0青海湖201.44.50.22巢湖306.260.251.67滇池20pn TFII h/-10.130.500.232问题的分析己卜模型的建在进行综合评价之前，首先要对评价的指标进行分析。通常评价指标分成效益型，成本型和固定型 指标。效益型指标是指那些数值越大影响力越大的统计指标（也称正向型指标）；成本型指标是指 数值越小越好的指标（也称逆向型指标）；而固定型指标是指数值越接近于某个常数越好的指标（也称适度型指标）。如果每个评价指标的属性不一样，则在综合评价时就容易发生偏差，必须先对各评价指标统一属性。 建模步骤（i）建立无量纲化实测数据矩阵和评价标准矩阵-13010.30.352.76 10510.70.42.0实测数据矩阵X201.44.50.22306.260.251.672010.130.50.23-14231106600.090.361.87.1027.1等级标准矩阵Y=37122.40.550.170.020.0610.321.4.6然后建立无量纲化实测专数据矩阵A和无量纲化等级标准矩阵B,其中Xij / m ax Xijj 壬 31yj mktk -3a=bkt=kk 二ij =利用| m in Xij/ Xqj =3Matlab，我们得到1m 1 nyktVkt、 k3X=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;B 仁丫(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.0833B =B =0.00150.00610.03480.16671.00000.00330.01330.06640.26201.00000.00460.01420.07080.30911.00000.00430.01300.06740.26091.0000(2 )计算各评价指标的权重首先计算矩阵B的各行向量的均值和标准差1 (加一侦J 5 厢,s =丫 =1,2, 3, 4然后计算变异系数Wj = /顶，最后对变异系数归一化得到各指标的权重为W =0.27670.24440.23470.2442程序命令：b=B;t=std(b)./mea n(b); w=t/sum(t);(3)建立各湖泊水质的综合评价模型通常我们可以利用向量之间的距离来衡量两个向量之间的接近程度，在Matlab中，有以下的函数命令来计算向量之间的距离；Disp(w,p):计算w中的每个行向量和p中每个列向量之间的欧式距离；Mandist(w,p):绝对值距离计算A中各行向量到B中各列向量之间的欧氏距离ikdij -(日若dik嗖肿，则第i个湖泊属于第k级 计算A中各行向量到B中各列向量之间的绝对值距离Dij|ak_ bJi 若Dik =min DJ，则第i个湖泊属于第k级结果及程序：jd=dist(A,B) mjd=ma ndist(A,B)jd =1.37690.28811.12710.50341.84721.83121.73741.59591.57981.48590.21850.20450.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所以，从上面的计算可以看出，尽管欧氏距离和绝对值距离的意义完全不一样，质的评但对湖泊水价等级是一样的，这表明了我们方法的稳定性。问题：利用其他的方法来计算个评价指标的权重，比较结果是否有所改变。