【正文】 素在該城區(qū)的空間分布，并分析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度.(2) 通過數(shù)據(jù)分析，說明重金屬污染的主要原因.(3) 分析重金屬污染物的傳播特征，由此建立模型，確定污染源的位置.(4) 分析你所建立模型的優(yōu)缺點(diǎn)，為更好地研究城市地質(zhì)環(huán)境的演變模式，還應(yīng)收集什么信息？有了這些信息，如何建立模型解決問題？2 基本假設(shè)(1)假設(shè)采樣數(shù)據(jù)確定可信，可以真實反映出當(dāng)?shù)刂亟饘俚奈廴厩闆r.(2)假設(shè)暫不考慮影響土壤金屬分布的物理性質(zhì).(3)假設(shè)暫不考慮自然環(huán)境對金屬濃度的影響.(4)假設(shè)重金屬元素的背景值可以真實地反映出當(dāng)?shù)夭皇芪廴緯r的真實值.(5)污染物濃度在y、z軸上的分布符合高斯分布（正態(tài)分布）. (6)源強(qiáng)是連續(xù)均勻的.(7)在擴(kuò)散過程中污染物質(zhì)量是守恒的.3 符號說明225。:輕度污染。figure,surf(X,Y,Z)附錄2利用excel 把采樣點(diǎn)位置數(shù)據(jù)按不同區(qū)域分開，分別存為data1，data2，data3，data4，,y,z與附錄1中表示的數(shù)據(jù)相同.x1= data1(:,1)。 z3=data3(:,3)。)。ks39。 =。附件239。V4V7FACTORS(8)fontfamily:Monospaced。V5V8/ROTATIONv439。for k=1:nc=x0(1)au(2)/au(1)。endfor i=1:n1b(i,1)=z(i+1)。fontweight:normal。FACTORS(3)V2D:\課件\畢業(yè)論文\A\39。VARIMAXSORTV1WINDOW=FRONT.FACTOR/FILE=39。v39。scatter3(x2,y2,z2,39。 z5=data5(:,3)。 y2=data2 (:,2)。 y=data (:,2)。 Gaussian diffusion model。 :中度污染。 y1=data1 (:,2)。scatter3(x3,y3,z3,39。contourf(X,Y,A)。) ,hold onx4= data4(:,1)。 =sqrt((^2+^2)/2)交通區(qū)：= mean(a(4,:))。V5V8ITERATE(25)fontsize:13pt。V6NOROTATE)。ycl(k+1)=c*exp(au(1) *k)+au(2)/au(1)。endfor k=2:nz(k)=*x1(k)+*x1(k1)。fontstyle:normal。/CRITERIAV1WINDOW=FRONT.FACTOR/FILE=39。/ROTATION/FORMAT/ANALYSISDataSet1scatter3(x5,y5,z5,39。 z2=data2 (:,3)。 y5=data5 (:,2)。),hold onx2= data2(:,1)。三次多項式插值法對8種重金屬進(jìn)行插值，比較形象的觀察出各種金屬在該城區(qū)空間分布的極值點(diǎn)，從而得到污染源的位置.(1)模型一中，8種重金屬元素空間分布圖的點(diǎn)都是離散的點(diǎn)，利用MATLAB的多項式插值法畫圖，得出的數(shù)據(jù)有些誤差.(2) 模型三的分布中，忽略了污染源對風(fēng)向、水流的考慮，而且直接從圖中觀察得出污染源的位置，得出的數(shù)據(jù)不夠準(zhǔn)確.為了更好地研究城市地質(zhì)環(huán)境的演變模式，還應(yīng)收集不同時間、不同季節(jié)中重金屬的污染濃度；更深層土壤中的污染濃度；風(fēng)向、水流、降雨量對重金屬污染的影響；生活區(qū)、工業(yè)區(qū)和交通區(qū)等重要污染區(qū)的污水、廢氣和廢棄物的排放量以及土壤的降解速度.[6]由于該題目中給出的信息是少量的、不完全，因此，在本問中，選用灰色預(yù)測模型，：對于給定的原始數(shù)據(jù)序列：，經(jīng)一次累加得到一個新序列：.設(shè)滿足一階常微分方程其中為待估計參數(shù)，為 發(fā)展灰數(shù)，.按最小二乘法求解得：其中 .把式代入式，并解微分方程，得到灰色預(yù)測模型為：通過灰色預(yù)測模型，可以得出以后不同時間各種重金屬在土壤中的濃度，通過分析上述方法得到的數(shù)據(jù)，就可以預(yù)測城市未來的地質(zhì)環(huán)境的演變模式. 編程實現(xiàn)(附錄8).7 模型的推廣本文中所建立的模型是金屬物在土壤中擴(kuò)散的的模型，該模型可以用在金屬物在大氣中擴(kuò)散的模型、石油在海洋中泄露的模型、放射性氣體在空間中擴(kuò)散的模型以及炸彈爆炸的煙霧濃度的擴(kuò)散模型等.8 結(jié)語隨著我國經(jīng)濟(jì)的發(fā)展和城市人口的增加，分析了城市的污染程度、確定了污染源，發(fā)現(xiàn)土壤重金屬污染已經(jīng)比較嚴(yán)重，而且重金屬的除去與清潔需要一個相當(dāng)漫長的過程，因此，人們必須要提高對城市土壤的保護(hù)意識，并且做好對土壤重金屬的清潔工作.參考文獻(xiàn)：[1] [M].北京：電子工業(yè)出版社，2012，10.[2] 王新智等. 基于單項指數(shù)法和模糊綜合評價法對松花江吉林市段水質(zhì)的評價[J]. 環(huán)境科學(xué)與管理，2012，9，9期，37卷： 184187.[3] 寇文杰等. 內(nèi)梅羅指數(shù)法在水質(zhì)評價中存在問題及修正[J].南水北調(diào)與水利科技，2012，8，4期，10卷：3847.[4] 呂振通，[M]. 北京：機(jī)械工業(yè)出版社，2009，6.[5] 石東偉，陳冬娜. 高斯擴(kuò)散模型在確定污染源位置中的應(yīng)用[J].河南科技學(xué)院學(xué)報， 2012，4，2期，40卷：5558.[6] 唐麗芳，賈冬青，孟慶鵬. 用MATLAB實現(xiàn)灰色預(yù)測GM(1，1)模型[J]. 滄州師范?？茖W(xué)校學(xué)報，2008，6，2期，24卷：3637.[7] [M].北京：中國水利水電出版社，2005，5.[8] 陳澤，[M].北京：電子工業(yè)出版社，2011，6.[9] Mark [M].北京：機(jī)械工業(yè)出版社，2005，6.[10] Frank R. Giordano，etc . First Course in Mathematical Modeling[M]. 北京：機(jī)械工業(yè)出版社，2009,10.[11] Schneider，Leo Modelling[J]. The Mathematics Teacher，1996， 89（7）： 608.附錄附錄1把附件1中的取樣點(diǎn)位置數(shù)據(jù)表命名為data.x=data(:,1)。 Principal ponent analysis。 :嚴(yán)重污染.下面用Excel 統(tǒng)計得出不同功能區(qū)中重金屬元素的含量平均值，并將得出的數(shù)據(jù)作為，8種重金屬的背景值分別為：，31，31，69 (mg/kg)，如表2.表1 各種金屬不同區(qū)域中的平均值(mg/kg)元素AsCdCrCuHgNiPbZn生活區(qū)工業(yè)區(qū)