clc,clear
A=xlsread('cumcm2011A附件_数据.xls','附件1','A4:E322');
C=xlsread('cumcm2011A附件_数据.xls','附件2','A4:I322');
x=A(:,2);y=A(:,3);z=A(:,4);qy=A(:,5);
As=C(:,2);Cd=C(:,3);Cr=C(:,4);Cu=C(:,5);
Hg=C(:,6);Ni=C(:,7);Pb=C(:,8);Zn=C(:,9);
[X,Y,Z]=griddata(x,y,z,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,Z)
title('采样地形图');
[X,Y,AS]=griddata(x,y,As,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,AS)
title('As浓度随平面地形分布图');
[X,Y,CD]=griddata(x,y,Cd,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,CD)
title('Cd浓度随平面地形分布图');
[X,Y,CR]=griddata(x,y,Cr,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,CR)
title('Cr浓度随平面地形分布图');
[X,Y,CU]=griddata(x,y,Cu,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,CU)
title('Cu浓度随平面地形分布图');
[X,Y,NI]=griddata(x,y,Ni,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,NI)
title('Ni浓度随平面地形分布图');
[X,Y,HG]=griddata(x,y,Hg,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,HG)
title('Hg浓度随平面地形分布图');
[X,Y,PB]=griddata(x,y,Pb,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,PB)
title('Pb浓度随平面地形分布图');
[X,Y,ZN]=griddata(x,y,Zn,linspace(0,28654)',linspace(0,18449),'v4'); figure,contourf(X,Y,ZN)
title('Zn浓度随平面地形分布图');
[X,Y,QY]=griddata(x,y,qy,linspace(0,28654)',linspace(0,18449),'v4');
figure,contourf(X,Y,QY)
title('区域分布图');
(2)
function R=byhs(w,v) ;
R(1)=0;
for i=1:318
N=fix(log(abs(v(i)-v(i+1))))
for j=1:N-1
ZZ=w(j)-w(j+1);
pf=ZZ*ZZ;
s=sum(pf);
end
r=s/(2*N);
R(i+1)=r;
end
clc,clear%将附件2中重金属元素浓度粘贴到紧邻附件1的右侧
B=xlsread('cumcm2011A附件_数据','附件1','A4:M322'); A=sortrows(B,2);
x=A(:,2);y=A(:,3);
As=A(:,6);Cd=A(:,7);Cr=A(:,8);Cu=A(:,9);
Hg=A(:,10);Ni=A(:,11);Pb=A(:,12);Zn=A(:,13);
Asx=byhs(As,x);
Cdx=byhs(Cd,x);
Crx=byhs(Cr,x);
Cux=byhs(Cu,x);
Hgx=byhs(Hg,x);
Nix=byhs(Ni,x);
Pbx=byhs(Pb,x);
Znx=byhs(Zn,x);
clc,clear%将附件2中重金属元素浓度粘贴到紧邻附件1的右侧
B=xlsread('cumcm2011A附件_数据','附件1','A4:M322'); A=sortrows(B,3);
x=A(:,2);y=A(:,3);
As=A(:,6);Cd=A(:,7);Cr=A(:,8);Cu=A(:,9);
Hg=A(:,10);Ni=A(:,11);Pb=A(:,12);Zn=A(:,13);
Asy=byhs(As,y);
Cdy=byhs(Cd,y);
Cry=byhs(Cr,y);
Cuy=byhs(Cu,y);
Hgy=byhs(Hg,y);
Niy=byhs(Ni,y);
Pby=byhs(Pb,y);
Zny=byhs(Zn,y);
(3)
clc, clear
a=load('adata1.txt'); %把附件1后4列数据保存到adata1.txt
b=load('adata2.txt'); %把采样点8种元素数据保存到adata2.txt
x0=a(:,1); y0=a(:,2); z0=a(:,3); %分别提取x,y,z的坐标
xmm=minmax(x0') %提取x的最大值和最小值
ymm=minmax(y0')
zmm=minmax(z0')
[xi,yi]=meshgrid([xmm(1):100:xmm(2)],[ymm(1):100:ymm(2)]); for i=1:8
Fv{i}=TriScatteredInterp(x0,y0,b(:,i));
vi{i}=Fv{i}(xi,yi);
figure(i), subplot(121), c{i}=contour(xi,yi,vi{i})
subplot(122),contourf(xi,yi,vi{i}),clabel(c{i})
end