【正文】 )z(i1,j1))amp。
for i=2:99
for j=2:99
if (z(i,j)z(i1,j))amp。
x=A(:,2)。Hg39。As 39。A4:E32239。)
figure,surf(X,Y,Z)%三維曲面end
附錄二：?jiǎn)雾?xiàng)污染指數(shù)求解的MATLAB程序：a=[



]function f=fun(a)
m=size(a,1)。,linspace(0,20000),39。F:\A\39。,39。)
legend(39。y=A(:,3)。y=A(:,3)。y=A(:,3)。y=A(:,3)。
z=A(:,5)。功能區(qū)39。
scatter(x,y,5,z)%散點(diǎn)圖
figure
[X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。八、附錄附錄一： 城區(qū)地形分布圖的MATLAB程序：A=xlsread(39。要更好的分析出重金屬污染的主要原因，我們還需對(duì)該城市周邊農(nóng)田中農(nóng)藥的使用等因素造成的污染進(jìn)行調(diào)查分析。解決問題二時(shí)，我們忽略了該市風(fēng)向、天氣等因素對(duì)重金屬污染的影響。在規(guī)整的濃度分布矩陣基礎(chǔ)上分別使用pcolor，contourf，contour，surf等函數(shù)繪出了各種重金屬元素在城區(qū)的空間分布。通過搜索模型可以求出重金屬空間分布中的極值點(diǎn)即可能的污染源，再結(jié)合國(guó)家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)[3]（下表5）通過MATLAB軟件對(duì)極值點(diǎn)進(jìn)行篩選出，求出重金屬的主要污染源。綜上所述，可以認(rèn)為工業(yè)區(qū)、主干道路區(qū)和生活區(qū)的活動(dòng)是造成該城區(qū)表層土壤重金屬污染的主要原因。 根據(jù)空間分布圖、區(qū)域散點(diǎn)圖和主要重金屬含量土壤單項(xiàng)污染的指數(shù)進(jìn)行分析：對(duì)于Cr和Ni，在來源上關(guān)聯(lián)較密切，該市表層土壤Cr和Ni基本未污染，只有個(gè)別點(diǎn)富集程度較高，污染達(dá)到中度污染，該富集中心的位置主要分布在生活區(qū)周邊和主干道區(qū)周邊，這可能是由于生活廢水的排放和交通源汽車尾氣的排放等原因造成的。根據(jù)HJ/T 1662004,內(nèi)梅羅綜合污染指數(shù)的分級(jí)標(biāo)準(zhǔn)（見表1），得出各個(gè)區(qū)域的污染等級(jí)。 模型建立：土壤環(huán)境質(zhì)量單項(xiàng)污染指數(shù)主要用來評(píng)價(jià)某一污染物的污染程度，指數(shù)小污染輕，指數(shù)大污染則重。搜索極值的現(xiàn)代算法有模擬退火，遺傳算法，魚群算法等多種。 假設(shè)問題二中只考慮題目中所給的8中重金屬，不考慮其它重金屬。應(yīng)用專門儀器測(cè)試分析，獲得了每個(gè)樣本所含的多種化學(xué)元素的濃度數(shù)據(jù)。在模型一數(shù)據(jù)處理基礎(chǔ)上建立遍歷搜索模型，結(jié)合MATLAB軟件求出重金屬空間分布中的極值點(diǎn)即可能的污染源，得出極值點(diǎn)后再結(jié)合《國(guó)家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》通過MATLAB軟件對(duì)極值點(diǎn)進(jìn)行篩選，得出8種重金屬元素的主要污染源。城市工業(yè)、經(jīng)濟(jì)的發(fā)展，污水排放和汽車尾氣排放等均能引起城市表層土壤重金屬污染。
關(guān)鍵詞：梅羅綜合污染指數(shù)評(píng)價(jià)法 污染等級(jí) 相關(guān)矩陣 遍歷搜索模型 污染源一 、問題重述隨著城市經(jīng)濟(jì)的快速發(fā)展和城市人口的不斷增加，人類活動(dòng)對(duì)城市環(huán)境質(zhì)量的影響日顯突出。現(xiàn)要求通過數(shù)學(xué)建模來完成以下任務(wù)：(1) 給出8種主要重金屬元素在該城區(qū)的空間分布，并分析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度。對(duì)于該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度的研究可以借助我國(guó)《土壤監(jiān)測(cè)技術(shù)規(guī)范》(HJ/T 1662004)[2]中推薦的內(nèi)梅羅綜合污染指數(shù)法進(jìn)行評(píng)價(jià)，求出不同區(qū)域重金屬的污染等級(jí)。(四) 問題4的分析：首先應(yīng)對(duì)問題一，二，三所建立的模型進(jìn)行優(yōu)缺點(diǎn)分析然后根據(jù)影響城市演化模型的因素，分析還應(yīng)搜集的數(shù)據(jù)以及模型如何建立的問題。目前土壤環(huán)境質(zhì)量評(píng)價(jià)方法有很多，各有優(yōu)點(diǎn)和缺點(diǎn)。如下表：表2：各區(qū)重金屬含量的平均值區(qū)域As (μg/g)Cd (ng/g)Cr (μg/g)Cu (μg/g)Hg (ng/g)Ni (μg/g)Pb (μg/g)Zn (μg/g)12345然后根據(jù)公式（1）、（2）結(jié)合MATLAB軟件算得各區(qū)重金屬單項(xiàng)污染指數(shù)和綜合污染指數(shù)，如下表：表3：各區(qū)重金屬單項(xiàng)污染指數(shù)和綜合污染指數(shù)區(qū)域單 項(xiàng) 污 染 指 數(shù)綜合污染指數(shù)AsCdCrCuHgNiPbZn12345再由內(nèi)梅羅綜合污染指數(shù)的分級(jí)標(biāo)準(zhǔn)得出各區(qū)的綜合污染等級(jí)，如下表：表4：各區(qū)綜合污染等級(jí)區(qū)域污染等級(jí)生活區(qū)5 重污染工業(yè)區(qū)5 重污染山區(qū)3 輕等污染主干道路區(qū)5 重污染公園綠地區(qū)4 中等污染從表中可以看出，該城區(qū)內(nèi)生活區(qū)、工業(yè)區(qū)、主干道路區(qū)屬于重污染區(qū)，公園綠地區(qū)屬于中等污染區(qū)，山區(qū)屬于輕度污染區(qū)。對(duì)于Cu,該市表層土壤Cu基本未污染，只有個(gè)別點(diǎn)富集程度較高，污染達(dá)到中度污染，該富集中心的位置主要分布在生活區(qū)周邊，這可能是由城市商業(yè)活動(dòng)、城市居民生活累加到土壤中的Cu。根據(jù)問題分析可知：污染源存在于二維網(wǎng)格中的某些格子中。用同樣方法得出其他7種重金屬在空間分布極大值個(gè)數(shù)（見表6）。模型三中依據(jù)模型一中建立的濃度分布矩陣建立了遍歷搜索模型。城市地質(zhì)環(huán)境是一個(gè)涉及到地球巖石圈表層的巖石、土壤、大氣、水和生物的復(fù)雜系統(tǒng)。
七、參考文獻(xiàn)[1] 張志涌，《精通MATLAB 》[M].北京：北京航天航空大學(xué)出版社， 234302，2003。A4:E32239。)。,1,39。
y=y(x1(:))。
y=y(x2(:))。
y=y(x3(:))。
y=y(x4(:))。
y=y(x5(:))。工業(yè)區(qū)39。F:\A\39。)。%插值
pcolor(X,Y,Z)。
b=b39。F:\A\39。,39。,39。
[X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。amp。amp。
[ii,jj]=find(z==1000)。F:\A\39。)。Cu39。Zn39。)。(z(i,j)z(i,j1))amp。(z(i,j)z(i+1,j+1))。disp(jj39。)。,39。,39。y=A(:,3)。for i=2:99 for j=2:99 if (z(i,j)z(i1,j))amp。(z(i,j)z(i1,j1))amp。 end。kkk=find(kk5