20xx年數(shù)學(xué)建模a題優(yōu)秀論文-城市表層土壤重金屬污染分析(文件)
【正文】 umcm2020A,附件 _數(shù)據(jù) .xls39。 x=A(:,2)。 x=x(x1(:))。) hold on。 x=x(x2(:))。) hold on。 x=x(x3(:))。) hold on。 x=x(x4(:))。) hold on。 x=x(x5(:))。) title(39。,39。,39。) matlab程序 d=xlsread(39。A2:C32039。 z=d(:,3)。 zz=griddata(x,y,z,xx,yy,39。 shading interp colormap(gray)。 nx=linspace(min(x),max(x),100)。v439。 set(n,39。,get(n,39。 p3=[ ]。 w1=p1./b w2=p2./b w3=p3./b w4=p4./b w5=p5./b p1=[ ]。 p5=[ ]。 w3=p3./b。p1m=max(w1)。p3m=max(w3)。p5m=max(w5)。p4m=max(w4)。p2m=max(w2)。 w5=p5./b。 w1=p1./b。 p3=[ ]。 p5=[ ]。))) end 單項(xiàng)污染指數(shù) p1=[ ]。,sprintf(39。 [cc,hh]=contour(xx,yy,zz,10)。 [xx,yy]=meshgrid(nx,ny)。 y=d(:,2)。)。 ny=linspace(min(y),max(y),100)。 （只要改變 ’ A2:C320’，就可以得到八種重金屬的分布圖） x=d(:,1)。附件 _數(shù)據(jù) .xls39。,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。)。 [4] 韓中庚，數(shù)學(xué)建模競賽 —— 獲獎(jiǎng)?wù)撐木x與點(diǎn)評，北京：科學(xué)出版社， 2020 八．附錄 城區(qū)地形分布圖的 MATLAB程序： A=xlsread(39。 本模型的建立對于研究減少城市污染問題和保護(hù)環(huán)境具有重要意義，尤其在當(dāng)今以高能耗高污染的生產(chǎn)模式為主的工業(yè)時(shí)代，該模型的建立對于研究城市規(guī)劃，和工廠位置選擇以及交通建設(shè)時(shí)具有重要的才能考價(jià)值，同時(shí)利用該模型也可以研究物質(zhì)擴(kuò)散現(xiàn)象的規(guī)律 。 所以要想得出各區(qū)在不同位置的污染等級，需進(jìn)一步求出各種重金屬的空間分布函數(shù)。 ,amp。首先我們需要求出空間中任意一點(diǎn)對應(yīng)的污染物濃度與時(shí)間的關(guān)系表達(dá)式；然后再利用問題一中求解不同區(qū)域重金屬的污染程度的模型對城市進(jìn)行污染程度分析；最后對由統(tǒng)計(jì)得到的各個(gè)樣點(diǎn)在多個(gè)時(shí)間內(nèi)的污染程度進(jìn)行分析，即可研究出城市地質(zhì)環(huán)境的演變模式。 （ 3）在問題三上，我們應(yīng)用 MATLAB 軟件對數(shù)據(jù)進(jìn)行處理，得到了污染源的可能分布地區(qū)，共有三處；然后再用 MATLAB 進(jìn)行三次擬合，比較簡便準(zhǔn)確地找到了污染源的位置。由于所給數(shù)據(jù)采樣點(diǎn)的不規(guī)則性 ， 首先使用 griddata函數(shù)對所給數(shù)據(jù)進(jìn)行插值規(guī)整得出一個(gè) X,Y分別等步長的某種元素的濃度分布矩陣。 問題三的求解 設(shè)第 m 個(gè) 污染源的污染值為 omX ，利用 MATLAB 軟件對數(shù)據(jù)進(jìn)行處理即可得到 omX 以及污染源的位置（見圖 12），了解重金屬污染物的傳播特征，考慮第 i 個(gè)采樣點(diǎn)受到污染源的影響與該采樣點(diǎn)和污染源的海拔差 Hj、空間距離Rj、以及第 m 個(gè) 污染源的污染值 omX 三個(gè)因素有關(guān)，所以構(gòu)建目標(biāo)函數(shù) （ 3） 其中 （ 4） （ 5） 最后再利用 MATLAB 對構(gòu)建的目標(biāo)函數(shù)（ 4）進(jìn)行數(shù)據(jù)擬合，得出污染源的傳播范圍及其傳播特征。工業(yè)排放也是表層土壤 Hg 污染的另一個(gè)重要來源，主要在大面積污染的幾個(gè)工業(yè)濃集中心。 對于 Pb 和 Cd ，在來源上關(guān)聯(lián)較密切， Pb 和 Cd 的高含量點(diǎn)主要分布在 交通 區(qū)和工業(yè)區(qū)周邊 ,這可能是因?yàn)?Pb 和 Cd 來自該市中心交通源汽車尾氣的排放、汽車輪胎的磨損和冶煉廠的廢水、塵埃和廢渣 ,以及電鍍、電池、顏料、塑料穩(wěn)定劑、涂料工業(yè)的廢水等。 參照分布圖 可見各重金屬濃度均和海拔成負(fù)相關(guān)，即海拔越高，其含各種重金屬濃度越低； Cr 和 Ni 的相關(guān)性最好，相關(guān)系數(shù)最大，為 ，其次為 Pb 和Cd ，相關(guān)系數(shù)為 0． 660，以下是 Cr 和 Cu 的相關(guān)性較好，相關(guān)系數(shù)是 ，其它元素之間的相關(guān)性并不是很好。然后得出各個(gè)區(qū)里面主要重金屬含量的平均值， 用 EXCEL做 出了各區(qū)域各重金屬平均值的 二維折線 圖 （見附錄 ） ， 可看作各區(qū)中主要重金屬含量值。單項(xiàng)污染指數(shù)和綜合污染指數(shù)的計(jì) 算公式如下： ij j jP C S? (1) 22, , m a x( ) 2N j a ve jP P P?? (2) 圖（ 9）： 鋅在該城區(qū)的空間分布圖 圖（ 10）：銅 在該城區(qū)的空間分布圖 當(dāng)ijP≤ 1時(shí)，表示土壤未受該因子污染，當(dāng)ijP1時(shí)，表示土壤受該因子污染。綜合污染指數(shù)可以綜合判斷某土壤多種污染物的聯(lián)合污染效應(yīng)。 圖 （ 2） 的 X、 Y軸為地理坐標(biāo)值（單位： m）。 （四）對問題四的分析： 首先對我們建立的模型進(jìn)行優(yōu)缺點(diǎn)分析，然后通過查詢資料了解影響 城市地質(zhì)環(huán)境演變 的其他一些因素，從而選擇模型解決問題。 （二） 對問題 2 的分析： 問題 2要求 通過數(shù)據(jù)分析來說明重金
公司管理相關(guān)推薦
鄂ICP備17016276號-1