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

1、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,linspa

2、ce(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&#

3、39;);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(&#

4、39;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=gridda

5、ta(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)',lin

6、space(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&#

7、39;,'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=xlsr

8、ead('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:8F