【正文】 endscatter(cc(:,1),cc(:,2),39。kkk=kkk39。n=size(ii,1)。 end。amp。amp。amp。amp。)。z=B(:,1)。Zn39。Ni39。Cu39。Cd39。)。B=xlsread(39。F:\A\39。)。
[ii,jj]=find(z==1000)。
z(i,j)=1000。amp。amp。amp。%插值
z=Z39。
[X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。}。,39。,39。,39。ss={39。F:\A\39。,1,39。
subplot(2,2,k),scatter(ii,jj,39。
disp(ii39。
end。(z(i,j)z(i+1,j1))amp。(z(i,j)z(i1,j1))amp。(z(i,j)z(i,j+1))amp。
for i=2:99
for j=2:99
if (z(i,j)z(i1,j))amp。,linspace(0,20000),39。
x=A(:,2)。Pb39。Hg39。Cr 39。As 39。,2,39。A4:E32239。
for i=1:5
for j=1:n。)
figure,surf(X,Y,Z)%三維曲面end
附錄二：單項污染指數(shù)求解的MATLAB程序：a=[



]function f=fun(a)
m=size(a,1)。shading interp%偽彩色圖
title(39。,linspace(0,20000),39。x=A(:,2)。F:\A\39。,1,39。,39。,39。)
legend(39。
scatter(x,y,20,39。y=A(:,3)。
scatter(x,y,20,39。y=A(:,3)。
scatter(x,y,20,39。y=A(:,3)。
scatter(x,y,20,39。y=A(:,3)。
scatter(x,y,20,39。
z=A(:,5)。A4:E32239。功能區(qū)39。%插值
pcolor(X,Y,Z)。
scatter(x,y,5,z)%散點圖
figure
[X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。)。八、附錄附錄一： 城區(qū)地形分布圖的MATLAB程序：A=xlsread(39。[2] HJ/T1662004，《土壤環(huán)境監(jiān)測技術規(guī)范》[S]. 北京：中國標準出版社，2004。要更好的分析出重金屬污染的主要原因，我們還需對該城市周邊農(nóng)田中農(nóng)藥的使用等因素造成的污染進行調(diào)查分析。為建立城市的地質演進模型，還應搜集巖石、土壤、大氣、水和生物等因素的相關信息，進而建立城市地質環(huán)境的綜合評價預測模型（參考文獻[5]）。解決問題二時，我們忽略了該市風向、天氣等因素對重金屬污染的影響。該模型能夠有效且快速的找出空間極大值，即可能的污染源。在規(guī)整的濃度分布矩陣基礎上分別使用pcolor，contourf，contour，surf等函數(shù)繪出了各種重金屬元素在城區(qū)的空間分布。表6：八種金屬元素空間分布極大值個數(shù)元素As Cd Cr Cu Hg Ni Pb Zn 個數(shù)6160576263605358 運用scatter函數(shù)畫出各重金屬元素空間分布極大值點的散點圖使數(shù)據(jù)可視化。通過搜索模型可以求出重金屬空間分布中的極值點即可能的污染源，再結合國家土壤環(huán)境質量標準[3]（下表5）通過MATLAB軟件對極值點進行篩選出，求出重金屬的主要污染源。并且污染源所在格子元素濃度大于周圍格子的元素濃度。綜上所述，可以認為工業(yè)區(qū)、主干道路區(qū)和生活區(qū)的活動是造成該城區(qū)表層土壤重金屬污染的主要原因。對于Hg, 其高含量點主要分布在交通繁忙的主干道路區(qū)周邊和工業(yè)區(qū)周邊, Hg污染的一個主要原因是由于燃煤造成的，無論是工業(yè)用煤還是居民用煤，而且燃燒方式落后。 根據(jù)空間分布圖、區(qū)域散點圖和主要重金屬含量土壤單項污染的指數(shù)進行分析：對于Cr和Ni，在來源上關聯(lián)較密切，該市表層土壤Cr和Ni基本未污染，只有個別點富集程度較高，污染達到中度污染，該富集中心的位置主要分布在生活區(qū)周邊和主干道區(qū)周邊，這可能是由于生活廢水的排放和交通源汽車尾氣的排放等原因造成的。二、問題二的求解： 模型建立：用SPSS11．0統(tǒng)計軟件對各種重金屬元素濃度和海拔做相關性分析，得出各種元素與元素之間和元素與海拔之間的相關系數(shù)矩陣及其相關性，結合第一問得出的空間分布圖和區(qū)域散點圖，參照主要重金屬含量土壤單項污染的指數(shù)，分析得出各重金屬污染的主要原因。根據(jù)HJ/T 1662004,內(nèi)梅羅綜合污染指數(shù)的分級標準（見表1），得出各個區(qū)域的污染等級。本文根據(jù)我國《土壤監(jiān)測技術規(guī)范》(HJ/T 1662004)[2]中推薦的內(nèi)梅羅綜合污染指數(shù)法進行評價。 模型建立：土壤環(huán)境質量單項污染指數(shù)主要用來評價某一污染物的污染程度，指數(shù)小污染輕，指數(shù)大污染則重。四、符號說明符號設定符號說明Pij區(qū)域i中第j個重金屬的污染分指數(shù)Cj第j個重金屬的實測濃度Sj第j元素的評價標準PN綜合污染指數(shù)Pj,ave平均單項污染指數(shù)Pj,max最大單項污染指數(shù)z濃度分布矩陣注：在此沒有設定的符號在下文中會具體說明。搜索極值的現(xiàn)代算法有模擬退火，遺傳算法，魚群算法等多種。(二) 問題2的分析：問題2要求通過數(shù)據(jù)分析來說明重金屬污染的主要原因。 假設問題二中只考慮題目中所給的8中重金屬，不考慮其它重金屬。(2) 通過數(shù)據(jù)分析，說明重金屬污染的主要原因。應用專門儀器測試分析，獲得了每個樣本所含的多種化學元素的濃度數(shù)據(jù)。對城市土壤地質環(huán)境異常的查證，以及如何應用查證獲得的海量數(shù)據(jù)資料開展城市環(huán)境質量評價，研究人類活動影響下城市地質環(huán)境的演變模式，日益成為人們關注的焦點。在模型一數(shù)據(jù)處理基礎上建立遍歷搜索模型，結合MATLAB軟件求出重金屬空間分布中的極值點即可能的污染源，得出極值點后再結合《國家土壤環(huán)境質量標準》通過MATLAB軟件對極值點進行篩選，得出8種重金屬元素的主要污染源。而重金屬污染對城市環(huán)境和人類健康造成了嚴重的威脅，因此對城市表層土壤重金屬污染的研究具有重大意義。我們鄭重承諾，嚴格遵守競賽規(guī)則，以保證競賽的公正、公平性。2011高教社杯全國大學生數(shù)學建模競賽承 諾 書我們仔細閱讀了中國大學生數(shù)學建模競賽的競賽規(guī)則。如有違反競賽規(guī)則的行為，我們將受到嚴肅處理。對于問題1，先用MATLAB軟件對所給數(shù)據(jù)進行處理，插值擬合得出8種主要重金屬元素在該城區(qū)的空間分布圖；再用內(nèi)梅羅綜合污染指數(shù)評價法建立模型進行求解。對于問題4，對所建立的模型進行分析，找出了各個模型的優(yōu)缺點。按照功能劃分，城區(qū)一般可分為生活區(qū)、工業(yè)區(qū)、山區(qū)、主干道路區(qū)及公園綠地區(qū)等，分別記為1類區(qū)、