第三问模型建立与求解.doc.doc

基于最小加权距离和求污染源模型 因为每种元素的污染源应该在以该污染物浓度最大的点为圆心的圆内,且到其周围污染程度较大的点的总加权距离最小,所以本文建立模型(1),其中Wi为加权因子. 具体实现步骤: (1)观察每个重金属元素的浓度分布图 由观察可得出,As的污染源个数为3个,Cd的污染源个数为4个Cr为1个Cu2个,Hg3个,Ni1个,Pb2个,Zn2个.并将其周围污染程度大的点与其归为一组. 因为权重因子主要与距离有关,离污染源越大,浓度梯度下降越大,所以本文以R=1000m为半径,画出每组的同心圆,确定权重因子Wi,(2). AS Cd Cr Cu 利用Lingo软件进行最优化求解,找出污染源坐标. 模型求解: 以As为例,其污染源个数为3个,具体分组如下表: 组一 编号 As x y 浓度变化 Wi 84 30.13 18134 10046 0 82 10.99 18738 10921 19.14 0.018002257 77 8.50 17087 11933 21.63 0 206 7.56 20282 8590 22.57 0.008697495 256 7.34 15198 10100 22.79 0 196 8.67 17005 7212 21.46 0.007034682 101 7.63 19968 12961 22.5 0.006533291 195 7.56 17203 6218 22.57 0.005729008 250 9.13 14065 10987 21 0.005028254 307 7.34 14482 12692 22.79 0.005053439 72 8.06 14298 7418 22.07 0.004746339 284 7.12 16629 14481 23.01 0.004913097 215 7.56 21018 5764 22.57 0.004371743 308 8.23 14318 13569 21.9 0.004216728 251 7.34 12734 10344 22.79 0.004213971 292 7.12 17414 15476 23.01 0.004200822 138 9.17 15007 5535 20.96 0.003818685 247 9.58 12641 9560 20.55 0.003726539 92 7.41 23664 9790 22.72 0 267 8.67 12153 12336 21.46 0.003350821 287 8.00 20983 15862 22.13 0.003417075 组二 编号 As x y 浓度变化 Wi 29 21.87 4742 7293 30 18.38 4948 7293 3.49 0.016941748 31 10.53 5567 6782 11.34 0.011685421 19 7.41 5481 6004 14.46 0.009732131 221 8.67 5006 8846 13.2 0.008379356 23 9.39 3573 6213 12.48 0.007841659 22 9.84 3299 6018 12.03 0.006247403 41 16.58 6869 7286 5.29 0.002487071 16 11.45 4777 4897 10.42 0 18 7.84 6534 5641 14.03 0 42 7.41 7056 8348 14.46 0 20 8.50 4592 4603 13.37 0.004962512 237 8.90 7106 9467 12.97 0 152 8.23 5101 4080 13.64 0.004218992 153 9.35 5438 3994 12.52 0.003713371 223 7.12 6395 10443 14.75 0.004146287 7 8.94 2883 3617 12.93 0.003138883 8 9.62 2383 3692 12.25 0.002845595 154 8.90 5382 3012 12.97 0 150 8.00 4020 2990 13.87 0 233 10.27 7776 10613 11.6 0.00