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）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;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);for i=1:8Fvi=TriScatteredInterp