MATLAB数学建模与仿真（第2版·微课视频版）程序代码 第12章-城市表层土壤重金属污染分析代码

实现调查区的地形图的MATLAB代码如下：A=xlsread('F:\A\cumcm2011A附件_数据.xls',1,'A4:E322');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)',linspace(0,20000),'v4');%插值pcolor(X,Y,Z);shadinginterp%伪彩色图title('功能区')figure,contourf(X,Y,Z)%等高线图figure,contour(X,Y,Z)title('功能区')figure,surf(X,Y,Z)%三维曲面图12-1调查地区的地形图实现功能区散点图的MATLAB代码如下：A=xlsread('F:\A\cumcm2011A附件_数据.xls',1,'A4:E322');x=A(:,2);y=A(:,3);z=A(:,5);x1=find(z==1);x=x(x1(:));y=y(x1(:));scatter(x,y,20,'d')holdon;x=A(:,2);y=A(:,3);x2=find(z==2);x=x(x2(:));y=y(x2(:));scatter(x,y,20,'h')holdon;x=A(:,2);y=A(:,3);x3=find(z==3);x=x(x3(:));y=y(x3(:));scatter(x,y,20,'s')holdon;x=A(:,2);y=A(:,3);x4=find(z==4);x=x(x4(:));y=y(x4(:));scatter(x,y,20,'p')holdon;x=A(:,2);y=A(:,3);x5=find(z==5);x=x(x5(:));y=y(x5(:));scatter(x,y,20,'x')title('功能区分布')legend('生活区','工业区','山区','主干道区','公园绿地区')A=xlsread('F:\A\cumcm2011A附件_数据.xls',1,'A4:E322');B=xlsread('F:\A\cumcm2011A附件_数据.xls',2,'B4:I322');x=A(:,2);y=A(:,3);fork=1:8z=B(:,k);scatter(x,y,5,z)%散点图figure[X,Y,Z]=griddata(x,y,z,linspace(0,30000)',linspace(0,20000),'v4');%插值pcolor(X,Y,Z);shadinginterp%伪彩色图title('功能区')figure,contourf(X,Y,Z)%等高线图figure,contour(X,Y,Z)title('功能区')figure,surf(X,Y,Z)%三维曲面end图12-2功能区散点图单项污染指数求解的MATLAB代码如下：a=[6.27289.96 69.02 49.4 93.04 18.34 69.11 237.017.25 393.11 53.41 127.54 642.36 19.81 93.04 277.934.04 152.32 38.96 17.32 40.96 15.45 36.56 73.295.71 360.01 58.05 62.21 446.82 17.62 63.53 242.856.26 280.5443.64 30.19 114.9915.29 60.71 154.24]m=size(a,1);n=size(a,2);c=[]b=[3.61303113.23512.33169];b=b';fori=1:5forj=1:n;c(i,j)=a(i,j)/b(j)endend表12-4各区重金属单项污染指数和综合污染指数区域单项污染指数综合污染指数AsCdCrCuHgNiPbZn11.74172.23052.22653.74242.65831.49112.22943.43493.170422.01393.02391.72299.662118.35311.61063.00134.02813.533131.12221.17171.25681.31211.17031.25611.17941.06221.253241.58612.76931.87264.712912.76631.43252.04943.51969.426451.73892.1581.40772.28713.28541.24311.95842.23542.7343实现重金属元素砷、镉、铬、铜的污染源分布图的MATLAB程序如下：A=xlsread('F:\A\cumcm2011A附件_数据.xls',1,'A4:E322');B=xlsread('F:\A\cumcm2011A附件_数据.xls',2,'B4:I322');ss={'As','Cd','Cr','Cu','Hg','Ni','Pb','Zn'};x=A(:,2);y=A(:,3);fork=1:4z=B(:,k);[X,Y,Z]=griddata(x,y,z,linspace(0,30000)',linspace(0,20000),'v4');%插值z=Z';fori=2:99forj=2:99if(z(i,j)>z(i-1,j))&&(z(i,j)>z(i+1,j))&&(z(i,j)>z(i,j+1))&&(z(i,j)>z(i,j-1))&&(z(i,j)>z(i-1,j-1))&&(z(i,j)>z(i-1,j+1))&&(z(i,j)>z(i+1,j-1))&&(z(i,j)>z(i+1,j+1));z(i,j)=1000;end;end;end;[ii,jj]=find(z==1000);disp(ii');disp(jj');subplot(2,2,k),scatter(ii,jj,'*'),title(ss{k})end图12-5砷、镉、铬、铜的空间分布极大值散点图实现重金属元素汞、镍、铅、锌的污染源分布图的MATLAB代码如下：A=xlsread('F:\A\cumcm2011A附件_数据.xls',1,'A4:E322');B=xlsread('F:\A\cumcm2011A附件_数据.xls',2,'B4:I322');ss={'As','Cd','Cr','Cu','Hg','Ni','Pb','Zn'};x=A(:,2);y=A(:,3);fork=5:8z=B(:,k);[X,Y,Z]=griddata(x,y,z,linspace(0,30000)',linspace(0,20000),'v4');%插值z=Z';fori=2:99forj=2:99if(z(i,j)>z(i-1,j))&&(z(i,j)>z(i+1,j))&&(z(i,j)>z(i,j+1))&&(z(i,j)>z(i,j-1))&&(z(i,j)>z(i-1,j-1))&&(z(i,j)>z(i-1,j+1))&&(z(i,j)>z(i+1,j-1))&&(z(i,j)>z(i+1,j+1));z