资源描述:
《2011数学建模全国赛a题代码.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
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,linspace(0,
2、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,c
3、ontourf(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浓度随平面地形分
4、布图');[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,li
5、nspace(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,1844
6、9),'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:M32
7、2');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中重
8、金属元素浓度粘贴到紧邻附件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)