【正文】 ???????????????? 5 污染程度的單項(xiàng)污染指數(shù)評(píng)價(jià) ?????????????? 7 綜合污染程度的評(píng)價(jià)模型 ???????????????? 9 問(wèn)題 2 的分析與求解 ????????????????????? 10 模型的建立 ?????????????????????? 10 模型的結(jié)果 ?????????????????????? 11 問(wèn)題 3 的分析與求解 ????????????????????? 12 土壤重金屬污染物的來(lái)源 ???????????????? 12 重金屬污染物的傳播特征 ???????????????? 13 基于 matlab 的數(shù)據(jù)處理與三維數(shù)據(jù)插值模型 ??????? ? 13 問(wèn)題 4 的分析 與求解 ????????????????????? 19 模型的優(yōu)點(diǎn) ?????????????????????? 19 數(shù)據(jù)的收集 ?????????????????????? 20 模型的建立與求解 ??????????????????? 20 7 模型的推廣 ?????????????????????????? 21 8 結(jié)語(yǔ) ????????????????????????????? 21 參考文獻(xiàn) ??????????????????????? ?????? 22 附錄 ??????????????????????????????? 23 1 城市表層土壤重金屬污染分析 摘 要： 本文主要對(duì)某城市表層土壤地質(zhì)環(huán)境重金屬污染的情況進(jìn)行分析 .首先，利用 MATLAB繪制出 8 種重金屬元素在該城區(qū)的空間分布，通過(guò)內(nèi)梅羅指數(shù)評(píng)價(jià)法來(lái)分析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度 .接著，通過(guò) SPSS，用相關(guān)性分析和主成分分析，來(lái)分析重金屬污染的主要原因 .然后，建立三次多項(xiàng)式插值模型和高斯模型，確定污染源的位置 .最后，運(yùn)用灰色預(yù)測(cè)模型對(duì)土壤中的重金屬的變化進(jìn)行預(yù)測(cè)，較好地研究了城市 地質(zhì)環(huán)境的演變模式 . 關(guān)鍵詞： 重金屬污染；內(nèi)梅羅指數(shù)評(píng)價(jià)法；相關(guān)性分析；主成分分析；高斯擴(kuò)散模型；灰色預(yù)測(cè)模型 中圖分類(lèi)號(hào)： X53 Analysis of Heavy Metal Pollution of Urban Topsoil Abstract: This paper analyses the heavy metal pollution situation of urban topsoil geological environment. First， using MATLAB to draw the spatial distribution of eight heavy metals in the city, and using nemerow pollution index method for the analysis of the heavy metal pollution of the city zone in different areas. Then， by using SPSS, pollution causes of heavy metals are analyzed with correlation analysis and principal ponent analysis. After that, establishing cubic polynomial interpolation model and Gaussian diffusion model to confirm the pollution sources. At last, grey prediction model is set up to forecast the changes of heavy metals in the soil, preferably forecasting the evolution of the geological environment model. Key words: Pollution of heavy metals。 y=data (:,2)。 y2=data2 (:,2)。 z5=data5(:,3)。 scatter3(x2,y2,z2,39。v39。fontfamily:Monospaced。 [X,Y,A]=griddata(x,y,a,linspace(0,27823)39。 au=b\y。 end b(:,2)=1。y= data(:,2)。附件 239。 z5=data5(:,3)。 y2=data2 (:,2)。) ,hold on x5= data5(:,1)。*39。v439。 Nemerow pollution index method 。 z=data (:,3)。 z2=data2 (:,3)。 scatter3(x5,y5,z5,39。ro39。) 25 附錄 4 矩陣 a 表示下表中各種重金屬在不同區(qū)域中的單項(xiàng)污染指數(shù) 元素 As Cd Cr Cu Hg Ni Pb Zn 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) 生活區(qū)： P = mean(a(1,:))。fontsize:13pt。,linspace(0,18449),39。 ycl(1)=x0(1)。 y(i)=x0(i+1)。textdecoration:none} 27 附錄 7 由于 8 種重金屬的空間分布 圖的程序都差不多，下面以 As 的空間分布圖為例 . x=data(:,1)。 /SHEET=name 39。 y5=data5 (:,2)。),hold on x2= data2(:,1)。mp39。 contourf(X,Y,Z) ,hold on scatter3(x1,y1,z1,39。三次多項(xiàng)式插值法分別對(duì) 8種重金屬污染物的濃度進(jìn)行數(shù)據(jù)插值，編程實(shí)現(xiàn) (附錄 7).畫(huà)出圖像后，觀察圖像的極值 點(diǎn)，根據(jù)附件 2的數(shù)據(jù)，利用 MATLAB的搜索函數(shù) find，先從重金屬的最高濃度開(kāi)始查找，然后找到該點(diǎn)對(duì)應(yīng)的 x， y坐標(biāo)，最后把找到的坐標(biāo)跟圖像對(duì)比，符合的就是圖像的極值點(diǎn)，找到第一點(diǎn)后，再用次高濃度查找，直到找到所有的極值點(diǎn)為止 . 2. 模型的結(jié)果 圖 5 As的濃度圖 從圖 5中，可以看出 As的濃度圖中分別有 2個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的位置為： 14 (18134， 10046)、 (12696， 3024). 圖 6 Cd的濃度圖 從圖 6中，可以看出 Cd的濃度圖中分別有 1個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的位置為：(21439， 11383). 圖 7 Cr的濃度圖 從圖 7中，可以看出 Cr的濃度圖中分別有 2個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的位置為：(3299， 6018)、 (4592， 4603). 15 圖 8 Cu的濃度圖 從圖 8中，可以看出 Cu的濃度圖中分別有 2個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的位置為：(2383， 3692)、 (3299， 6018). 圖 9 Hg的濃度圖 從圖 9中，可以看出 Hg的濃度圖中分別有 4個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的 位置為：(2383， 3692)、 (2708， 2295)、 (13694， 2357)、 (15248， 29106). 16 圖 10 Ni的濃度圖 從圖 10中，可以看出 Ni的濃度圖中分別有 1個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的位置為： (3299， 6018). 圖 11 Pb的濃度圖 從圖 11中，可以看出 Pb的濃度圖中分別有 2個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的位置為： (1647， 2728)、 (4777， 4897). 17 圖 12 Zn的濃度圖 從圖 11中，可以看出 Zn的濃度圖中分別有 5個(gè)極值點(diǎn)，可以大概確定極值點(diǎn)的位置為： (9328， 4311)、 (13797， 9621)、 (12696， 3024)、 (4948， 7293)、 (2383， 3692). 根據(jù)各種重金屬濃度圖中的極