2011年A题的作图程序.doc

clear,clc count=319;%count±íê?è??ùμ?μ???êy %fidheng=fopen(hengzuobiao.txt,r); %heng=zeros(1,count); %[heng,count]=fread(fidheng,inf); x=load(xzuobiao.txt); x=x; y=load(yzuobiao.txt); y=y; z=load(zzuobiao.txt); z=z; as=load(as.txt); as=as; cd=load(cd.txt); cd=cd; cr=load(cr.txt); cr=cr; cu=load(cu.txt); cu=cu; hg=load(hg.txt); hg=hg; ni=load(ni.txt); ni=ni; pb=load(pb.txt); pb=pb; zn=load(zn.txt); zn=zn; gongneng=load(gongnengqu.txt); gongneng=gongneng; %plot3(x,y,as); hold on As=ones(size(as))*as; [X,Y]=meshgrid(x,y); x1=linspace(min(x),max(x)); y1=linspace(min(y),max(y)); [X1,Y1]=meshgrid(x1,y1); %As1=griddata(X,Y,As,X1,Y1) %As1=interp2(x,y,As,X,Y); bianhao=1:1:319; %?-ê?êy?Y′|àí maxas=max(as); minas=min(as); as=(as-minas)/(maxas-minas); maxcd=max(cd); mincd=min(cd); cd=(cd-mincd)/(maxcd-mincd); maxcr=max(cr); mincr=min(cr); cr=(cr-mincr)/(maxcr-mincr); maxcu=max(cu); mincu=min(cu); cu=(cu-mincu)/(maxcu-mincu); maxhg=max(hg); minhg=min(hg); hg=(hg-minhg)/(maxhg-minhg); maxni=max(ni); minni=min(ni); ni=(ni-minni)/(maxni-minni); maxpb=max(pb); minpb=min(pb); pb=(pb-minpb)/(maxpb-minpb); maxzn=max(zn) minzn=min(zn) zn=(zn-minzn)/(maxzn-minzn); %plot(bianhao,as,bianhao,cd,bianhao,cr,bianhao,cu,bianhao,hg,bianhao,ni,bianhao,pb,bianhao,zn) %mesh(X1,Y1,As1); subplot(4,2,1); plot(bianhao,as); subplot(4,2,2); plot(bianhao,cd); subplot(4,2,3); plot(bianhao,cr); subplot(4,2,4); plot(bianhao,cu); subplot(4,2,5); plot(bianhao,hg); subplot(4,2,6); plot(bianhao,ni); subplot(4,2,7); plot(bianhao,pb); subplot(4,2,8); plot(bianhao,zn); data=[x;y;as]; data1=[x;y;gongneng]; x=data(:,1); y=data(:,2); as=data(:,3); gongz=data1(:,3); %scatter(x,y,z);%é￠μ?í? figure [xx,yy,zz]=griddata(x,y,as,linspace(min(x),max(x)),linspace(min(y),max(y)),v4);%2??μ [xx,yy,zz]=griddata(x,y,gongz,linspace(min(x),max(x)),linspace(min(y),max(y)),v4);%2??μ daxiao=size(zz); zz=zz-3.6; for i=1:100 for j=1:100 if zz(i,j)0 zz(i,j)=0; end end end xi=linspac