20xx年數(shù)學(xué)建模獲獎(jiǎng)?wù)撐腶題城市表層土壤重金屬污染分析-資料下載頁(yè)
2025-08-08 21:51本頁(yè)面
【正文】 )
figure,contourf(X,Y,Z) %等高線圖
figure,contour(X,Y,Z)
title(39。功能區(qū)39。)
figure,surf(X,Y,Z)%三維曲面end
附錄二：?jiǎn)雾?xiàng)污染指數(shù)求解的MATLAB程序：a=[
]function f=fun(a)
m=size(a,1)。
n=size(a,2)。
c=[]
b=[
130
31
35
31
69
]。
b=b39。
for i=1:5
for j=1:n。
c(i,j)=a(i,j)/b(j)end
end
附錄三：由各區(qū)的平均單項(xiàng)污染指數(shù)Pj,ave和最大單項(xiàng)污染指數(shù)Pj,max求各區(qū)的綜合污染指數(shù)的MATLAB程序：Pj,avePj,maxfunction f=fun4(x)
a=((x(1)^2+x(2)^2)/2)^(1/2)附錄四： 重金屬元素砷、鎘、鉻、銅的污染源分布圖的MATLAB程序：A=xlsread(39。F:\A\39。,1,39。A4:E32239。)。B=xlsread(39。F:\A\39。,2,39。B4:I32239。)。ss={39。As 39。,39。Cd39。,39。Cr 39。,39。Cu39。,39。Hg39。,39。Ni39。,39。Pb39。,39。Zn39。}。
x=A(:,2)。y=A(:,3)。
for k=1:4
z=B(:,k)。
[X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。,linspace(0,20000),39。v439。)。%插值
z=Z39。
for i=2:99
for j=2:99
if (z(i,j)z(i1,j))amp。amp。(z(i,j)z(i+1,j))amp。amp。(z(i,j)z(i,j+1))amp。amp。(z(i,j)z(i,j1))amp。amp。(z(i,j)z(i1,j1))amp。amp。(z(i,j)z(i1,j+1))amp。amp。(z(i,j)z(i+1,j1))amp。amp。(z(i,j)z(i+1,j+1))。
z(i,j)=1000。
end。
end。
end。
[ii,jj]=find(z==1000)。
disp(ii39。)。disp(jj39。)。
subplot(2,2,k),scatter(ii,jj,39。*39。),title(ss{k})
end
重金屬元素汞、鎳、鉛、鋅的污染源分布圖的MATLAB程序：A=xlsread(39。F:\A\39。,1,39。A4:E32239。)。B=xlsread(39。F:\A\39。,2,39。B4:I32239。)。ss={39。As 39。,39。Cd39。,39。Cr 39。,39。Cu39。,39。Hg39。,39。Ni39。,39。Pb39。,39。Zn39。}。
x=A(:,2)。y=A(:,3)。
for k=5:8
z=B(:,k)。
[X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。,linspace(0,20000),39。v439。)。%插值
z=Z39。
for i=2:99
for j=2:99
if (z(i,j)z(i1,j))amp。amp。(z(i,j)z(i+1,j))amp。amp。(z(i,j)z(i,j+1))amp。amp。(z(i,j)z(i,j1))amp。amp。(z(i,j)z(i1,j1))amp。amp。(z(i,j)z(i1,j+1))amp。amp。(z(i,j)z(i+1,j1))amp。amp。(z(i,j)z(i+1,j+1))。
z(i,j)=1000。
end。
end。
end。
[ii,jj]=find(z==1000)。
disp(ii39。)。disp(jj39。)。
subplot(2,2,k4),scatter(ii,jj,39。*39。),title(ss{k})
end
附錄五：篩選污染源的MATLAB程序（以砷為例，其它元素的程序只需調(diào)整參數(shù)即可）：A=xlsread(39。F:\A\39。,1,39。A4:E32239。)。B=xlsread(39。F:\A\39。,2,39。B4:I32239。)。 ss={39。As 39。,39。Cd39。,39。Cr 39。,39。Cu39。,39。Hg39。,39。Ni39。,39。Pb39。,39。Zn39。}。x=A(:,2)。y=A(:,3)。z=B(:,1)。[X,Y,Z]=griddata(x,y,z,linspace(0,30000)39。,linspace(0,20000),39。v439。)。%插值z(mì)=Z39。k=z。for i=2:99 for j=2:99 if (z(i,j)z(i1,j))amp。amp。(z(i,j)z(i+1,j))amp。amp。(z(i,j)z(i,j+1))amp。amp。(z(i,j)z(i,j1))amp。amp。(z(i,j)z(i1,j1))amp。amp。(z(i,j)z(i1,j+1))amp。amp。(z(i,j)z(i+1,j1))amp。amp。(z(i,j)z(i+1,j+1))。 z(i,j)=1000。 end。 end。end。 [ii,jj]=find(z==1000)。 c=[ii,jj]。n=size(ii,1)。for i=1:n kk(i)=k(ii(i),jj(i))。endpp=kk39。kkk=find(kk500)。kkk=kkk39。n=size(kkk,1)for i=1:n rr(i)=pp(kkk(i))。endrrfor i=1:n cc(i,1)=c(kkk(i),1)。 cc(i,2)=c(kkk(i),2)。endscatter(cc(:,1),cc(:,2),39。*39。)axis([0,100,0,100])title(ss{1})
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1