長江水质的评价和预测数据处理.docVIP

长江水质的评价和预测1）构建n个观测站的m个评价指标的判断矩阵 R=（i=1,2…n, j=1,2…m） 2）将判断矩阵R归一化处理，得到矩阵 (i=1,2,3…n, j=1,2,…m） 在上式中，有 计算评价指标的熵权矩阵公式W： 其中矩阵元素的计算公式为 ，且满足 相对差： 对相对差的求和： 定义衡量标准： 定义程度系数： 表示水质某时期内长江总体污染的程度； k表示要调查的月份数(本题k=28)； 表示第q个月的综合评测指数； 横坐标1代表2003.6,2代表2003.7,3代表2003.8依次类推； 纵坐标代表值；线是值 目标函数为： (i=1,2;n=1,2,3…7) 其中为n观测站污染源的排污量；为n观测站观测到的i污染物的质量；为来自n观测站上一个观测站的排污量。 对于长江干流上的七个主要观测站：攀枝花、重庆、?宜昌、岳阳、九江、安庆、南京，可分别称其为、、、、、、。 对于攀枝花至重庆、重庆至宜昌、宜昌至岳阳、岳阳至九江、九江至安庆、安庆至南京的长江干流流域段，可分别定义为：、、、、、。 染物的浓度为： 其中c为上游起点的污染物浓度；k为污染物的降解系数；x为两站点间的距离；v为该河段的平均流速。 由质量与浓度、体积的关系式：，可得下游x处的污染物质量： 质量变为： 其中x为两观测站的距离，v为该河段的流速。 其中为n站点i污染物的浓度，为观测站的水流量。 本地区的污染物源产量为： 其中 为两站点之间的平均距离； 式中i=1,2；n=1,2,3…,7 。 长江干流的江水做匀加速运动(a≤0或a≥0)。定义、、、、、这6段干流流域段的平均流速矩阵如下： 注：1）其中为流域段上游的观测站在i时段的水流速。 2）其中为流域段下游的观测站在i时段的水流速。 3）其中i为从2004年4月至2005年4月的13个时段。 由原始数据经分析处理，结合公式，运用MATLAB编程（程序见附录）得到速度矩阵如下表： 流域段 时段 2004．04 2.9 1.5 0.9 0.95 1.05 1.15 2004．05 2.8 1.35 0.85 1 1.1 1.15 2004．06 3 1.65 1.25 1.4 1.5 1.55 2004．07 3.2 1.85 1.45 1.5 1.55 1.65 2004．08 2.95 1.75 1.4 1.45 1.6 1.7 2004．09 4.95 3.25 1.8 2 2.75 3.4 2004．10 2.7 1.9 1.55 1.6 1.65 1.8 2004．11 2.3 1.3 0.75 0.85 0.9 0.95 2004．12 2.3 1.1 0.75 0.8 0.8 0.85 2005．01 1.8 1 0.55 0.65 0.7 0.75 2005．02 1.5 0.7 0.5 0.65 0.7 0.75 2005．03 1.4 0.65 0.5 0.7 0.8 0.85 2005．04 1.65 0.8 0.45 0.6 0.75 0.8 速度矩阵表 单位：m/s 定义、、、、、的流域段长度分别为：、、、、、， 由问题原始数据累积差可得： 通过L与矩阵的计算( 过程见附录MATLAB程序)可以得到时间矩阵： 其中为i时刻，江水完全流过水区域段所需的时间(单位：天)。 时间矩阵表 流域段 时段 2004．04 3.7915 6.0031 5.0797 6.0916 1.8078 4.6699 2004．05 3.9269 6.6701 5.3785 5.7870 1.7256 4.6699 2004．06 3.6651 5.4574 3.6574 4.