20xx年數(shù)學(xué)建模a題優(yōu)秀論文-城市表層土壤重金屬污染分析-資料下載頁
2024-11-06 10:06本頁面
【導(dǎo)讀】的影響日顯突出。對城市土壤地質(zhì)環(huán)境異常的查證，以及如何應(yīng)用查證獲得的海。式，日益成為人們關(guān)注的焦點(diǎn)。對于問題二我們首先用SPSS軟件對八種重金屬的濃度與海拔做相關(guān)性分析，的主要原因主要來自工業(yè)區(qū)、交通區(qū)和生活區(qū)。對于問題三我們建立模型，建立目標(biāo)函數(shù)；三次擬合后，得到污染源的位置。濃度以及各污染源重金屬的產(chǎn)生率。為此，將所考察的城區(qū)劃分為間距。度）進(jìn)行取樣、編號(hào)，并用GPS記錄采樣點(diǎn)的位置。重金屬污染是迫在眉下的事。二次插值法可以畫出8種主要重金屬元素在該城區(qū)的空間分布圖。區(qū)域的污染指數(shù)來分析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度。問題2要求通過數(shù)據(jù)分析來說明重金屬污染的主要原因。問題3要求通過分析重金屬污染物的傳播特征，找出污染源的位置。時(shí)除用單項(xiàng)污染指數(shù)外，還常用綜合污染指數(shù)。壤多種污染物的聯(lián)合污染效應(yīng)。
【正文】 4] 韓中庚，數(shù)學(xué)建模競賽 —— 獲獎(jiǎng)?wù)撐木x與點(diǎn)評，北京：科學(xué)出版社， 2020 八．附錄 城區(qū)地形分布圖的 MATLAB程序： A=xlsread(39。F:\A\cumcm2020A,附件 _數(shù)據(jù) .xls39。,1,39。A4:E32239。)。 x=A(:,2)。y=A(:,3)。z=A(:,4)。 scatter(x,y,5,z) figure [X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。,linspace(0,20200),39。v439。)。 pcolor(X,Y,Z)。shading interp title(39。功能區(qū) 39。) figure,contourf(X,Y,Z) figure,contour(X,Y,Z) title(39。功能區(qū) 39。) figure,surf(X,Y,Z) 功能區(qū)分布散點(diǎn)圖的 MATLAB 程序： A=xlsread(39。F:\A\cumcm2020A,附件 _數(shù)據(jù) .xls39。,1,39。A4:E32239。)。 x=A(:,2)。y=A(:,3)。 z=A(:,5)。 x1=find(z==1)。 x=x(x1(:))。 y=y(x1(:))。 scatter(x,y,20,39。d39。) hold on。 x=A(:,2)。y=A(:,3)。 x2=find(z==2)。 x=x(x2(:))。 y=y(x2(:))。 scatter(x,y,20,39。h39。) hold on。 x=A(:,2)。y=A(:,3)。 x3=find(z==3)。 x=x(x3(:))。 y=y(x3(:))。 scatter(x,y,20,39。s39。) hold on。 x=A(:,2)。y=A(:,3)。 x4=find(z==4)。 x=x(x4(:))。 y=y(x4(:))。 scatter(x,y,20,39。p39。) hold on。 x=A(:,2)。y=A(:,3)。 x5=find(z==5)。 x=x(x5(:))。 y=y(x5(:))。 scatter(x,y,20,39。x39。) title(39。功能區(qū)分布 39。) legend(39。生活區(qū) 39。,39。工業(yè)區(qū) 39。,39。山區(qū) 39。,39。交通區(qū) 39。,39。公園綠地區(qū) 39。) matlab程序 d=xlsread(39。F:\A\cumcm2020A184。附件 _數(shù)據(jù) .xls39。,4,39。A2:C32039。)。 （只要改變 ’ A2:C320’，就可以得到八種重金屬的分布圖） x=d(:,1)。 y=d(:,2)。 z=d(:,3)。 nx=linspace(min(x),max(x),100)。 ny=linspace(min(y),max(y),100)。 [xx,yy]=meshgrid(nx,ny)。 zz=griddata(x,y,z,xx,yy,39。v439。)。 surfl(xx,yy,zz)。 shading interp colormap(gray)。 hold on d=data x=d(:,1)。 y=d(:,2)。 z=d(:,3)。 nx=linspace(min(x),max(x),100)。 ny=linspace(min(y),max(y),100)。 [xx,yy]=meshgrid(nx,ny)。 zz=griddata(x,y,z,xx,yy,39。v439。)。 [cc,hh]=contour(xx,yy,zz,10)。 p=clabel(cc,hh) for n=p39。 set(n,39。string39。,sprintf(39。%.0f39。,get(n,39。userdata39。))) end 單項(xiàng)污染指數(shù) p1=[ ]。 p2=[ ]。 p3=[ ]。 p4=[ ]。 p5=[ ]。 b=[ 130/1000 31 35 31 69 ]。 w1=p1./b w2=p2./b w3=p3./b w4=p4./b w5=p5./b p1=[ ]。 p2=[ ]。 p3=[ ]。 p4=[ ]。 p5=[ ]。 b=[ 130/1000 31 35 31 69 ]。 w1=p1./b。 w2=p2./b。 w3=p3./b。 w4=p4./b。 w5=p5./b。 p1p=sum(w1)/8。p1m=max(w1)。 p1z=sqrt((p1p^2+p1m^2)/2) p2p=sum(w2)/8。p2m=max(w2)。 p2z=sqrt((p2p^2+p2m^2)/2) p3p=sum(w3)/8。p3m=max(w3)。 p3z=sqrt((p3p^2+p3m^2)/2) p4p=sum(w4)/8。p4m=max(w4)。 p4z=sqrt((p4p^2+p4m^2)/2) p5p=sum(w5)/8。p5m=max(w5)。 p5z=sqrt((p5p^2+p5m^2)/2) As分布圖 Cd分布圖 0100200300400500600700As (μg/g) Cd(ng/g)Cr (μg/g) Cu (μg/g) Hg(ng/g)Ni (μg/g) Pb (μg/g) Zn (μg/g) 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) Cr分布圖 Cu分布圖 Hg 分布圖 Ni 分布圖 Pb 分布圖 Zn 分布圖
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1