高教社杯全国大学生数学建模竞赛全国一等奖A题城市表层土壤重金属污染分析

各重金属分布浓度等高线及采样点坐标综合分布图绘制程序：

d=data; % 只需从excel输入三列数据，格式为：[xi yi zi]，i表示行数，xi为采样点坐标x值，yi为采样点坐标y值，zi为某重金属浓度值（此程序需将第三列的值更换8次运行8次得到论文中8幅各重金属浓度等高线及采样点坐标综合分布图）。 x=d(:,1); % 采样点坐标x值 y=d(:,2); % 采样点坐标y值 z=d(:,3); % 重金属浓度值 nx=linspace(min(x),max(x),40); ny=linspace(min(y),max(y),40); [xx,yy]=meshgrid(nx,ny); zz=griddata(x,y,z,xx,yy,'v4'); contourf(xx,yy,zz,10); hold on for i=1:319 for i=1:44

plot(d(i,1),d(i,2),'bs'); end hold on for i=45:80

plot(d(i,1),d(i,2),'b+'); end hold on for i=81:146

plot(d(i,1),d(i,2),'bp'); end hold on for i=147:284

plot(d(i,1),d(i,2),'ko'); end

hold on for i=285:319

plot(d(i,1),d(i,2),'r<'); end end

8种重金属的浓度等高线在该地区三维地形图曲面的投影图程序：

d=data; % 只需从excel输入三列数据，格式为：[xi yi zi]，i表示行数，xi为采样点坐标x值，yi为采样点坐标y值，zi为某重金属浓度值（此程序需将第三列的值更换8次运行8次得到论文中8幅各重金属浓度等高线在该地区三维地形图曲面的投影图）。

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,'v4'); 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,'v4'); [cc,hh]=contour(xx,yy,zz,10); p=clabel(cc,hh) for n=p'

set(n,'string',sprintf('%.0f',get(n,'userdata'))) end