2011数学建模全国赛a题代码

clc,clearA=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）functionR=byhs(w,v);R(1)=0;fori=1:318N=fix(log(abs(v(i)-v(i+1))))forj=1:N-1ZZ=w(j)-w(j+1);pf=ZZ*ZZ;s=sum(pf);endr=s/(2*N);R(i+1)=r;endclc,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,cleara=load('adata1.txt');%把附件1后4列数据保存到adata1.txtb=load('adata2.txt');%把采样点8种元素数据保存到adata2.txtx0=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)]);fori=1:8Fv{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