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

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); v4 v4 ); title( 采样地形图 ); v4 v4 ); title( As 浓度随平面地形分布图 ); v4 v4 ); title( Cd 浓度随平面地形分布图 ); v4 v4 ); title( Cr 浓度随平面地形分布图 ); v4 v4 ); title( Cu 浓度随平面地形分布图 ); v4 v4 ); title( Ni 浓度随平面地形分布图 ); v4 v4 ); title( Hg 浓度随平面地形分布图 ); v4 v4 ); title( Pb 浓度随平面地形分布图 ); v4 v4 ); title( Zn 浓度随平面地形分布图 ); v4 );[X,Y,QY]=griddata(x,y,qy,linspace(0,28654),linspace(0,18449), figure,contourf(X,Y,QY) title( v4 ); (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