2011年A题的作图程序.doc

clear,clccount=319;%count%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 onAs=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;%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);%figurexx,yy,zz=griddata(x,y,as,linspace(min(x),max(x),linspace(min(y),max(y),v4);%xx,yy,zz=griddata(x,y,gongz,linspace(min(x),max(x),linspace(min(y),max(y),v4);%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 endendxi=linspace(min(x),max(x),100);yi=linspace(min(y),max(y),100);xi,yi=meshgrid(xi,yi);zi=griddata(x,y,z,xi,yi);%gongnengi=griddata(x,y,gongneng,xi,yi);marker=,*,s,p,o;color=r,b,k,c,y