重金属污染的附录.docx

经过对数据和题目的分析，直接使用matlab（附录一）使用二次插值法可以画出8种主要重金属元素在该城区的空间分布图，由图形来反映各种的分布情况，附录一：plot(xm(1),ym(1),ow) hold for i=1:319 switch gongneng(i) case 1 plot(xm(i),ym(i),or) case 2 plot(xm(i),ym(i),*g) case 3 plot(xm(i),ym(i),*b)case 4 plot(xm(i),ym(i),*r) case 5 plot(xm(i),ym(i),ob) end end load data.txt load data2.txt x=data(:,1); y=data(:,2); z=data2(:,8); xi=0:100:max(x); yi=min(y):100:max(y); zi=griddata(x,y,z,xi,yi,cubic); contour(xi,yi,zi) title(Zn的空间分布图)由图形可以得出：某种重金属元素对某个区域具有非常显著的影响，因此我们再次利用软件matlab（附录二）求出各种重金属在各个区域最大的浓度和所在的区域，具体的信息可以在下表体现：附录二：for i=1:8; [a,b]=max(data(:i)); End 由matlab（附录三、四）软件可以得到各种重金属元素在各个区域的污染指数。附录三：%划分区域for i=1:size(qu,1) if qu(i)==1 shenghuo=[shenghuo;i]; end if qu(i)==2 gongye=[gongye;i]; end if qu(i)==3 shanqu=[shanqu;i]; end if qu(i)==4 daolu=[daolu;i]; end if qu(i)==5 gongyuan=[gongyuan;i]; end end shenghuo,gongye,shanqu,daolu,gongyuan,结果：1）生活区为以下采样点：12 13 16 18 20 21 23 25 33 36 39 42 63 67 68 72 94 106 107 152 154 156 157 158 176 180 183 184 186 187 195 198 199 243 254 258 259 267 268 269 272 273 276 306 2）工业区为以下采样点：4 6 8 29 30 31 103 153 177 185 197 221 223 224 225 226 228 230 231 232 233 237 238 239 246 247 248 249 255 262 270 271 274 275 277 278 3）山区为以下采样点：62 76 80 81 82 92 100 101 104 108 109 110 111 112 114 115 116 117 119 120 121 122 123 124 125 126 127 128 129 131 132 133 134 135 136 137 138 189 194 207 209 210 211 217 218 219 220 251 280 281 282 286 287 288 290 291 295 296 297 298 299 300 301 302 304 305 4）主干道路区为以下采样点：1 2 3 5 7 9 10 14 15 17 19 22 26 27 28 32 34 35 37 38 40 41 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 64 65 66 71 73 75 77 78 79 83 84 85 86 87 88 89 90 91 93 95 96 97 98 99 102 105 113 118 130 150 151 155 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 178 179 181 182 188 191 192 196 201 203 204 205 206 208 212 214 215 216 222 227 229 234 235 236 240 241 242 244 245 250 252 256 257 260 261 263 264 265 266 279 283 284 285 289 292 293 303 318 5）公园绿地区为以下采样点：11 24 69 70 74 139 140 141 142