三维立体图matlab代码.doc

clc

clear

closeall

D=load(ZuoBiao.dat);

M=load(JinShu.dat);

x=D(:,2);

y=D(:,3);

z=D(:,4);

c=D(:,5);

xi=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);

ci=griddata(x,y,c,xi,yi);

marker={*,o,s,^,p};

color={k,r,y,c,b};

mat={As,Cd,Cr,Cu,Hg,Ni,Pb,Zn};

str={等高线,生活区,工业区,山林区,交通区,绿地区};

forj=1:8

%等高线图形

figure

contourf(xi,yi,zi,0:10:500);

%set(h,ShowText,on,TextStep,get(h,LevelStep)*2)

%clabel(C,h,[0:10:50,50:50:300,300:100:500]);

title([金属,mat{j},二维等高线分布图])

xlabel(X)

ylabel(Y)

colormapsummer

colorbar

gridon

holdon

fori=1:5

loc=c==i;

plot(x(loc),y(loc),marker{i},markerfacecolor,color{i},MarkerEdgeColor,color{i});

end

legend(str,location,best)

fork=1:length(x)

text(x(k)-200,y(k)+200,num2str(M(k,j+1)),fontsize,8);

end

end

%三维体图

figure

h=surf(xi,yi,zi);

set(h,cdata,ci);

colormaphsv

title(三维图立体（颜色条表示分类）)

xlabel(X)

ylabel(Y)

colorbar

hiddenoff

holdon

fori=1:5

loc=c==i;

plot3(x(loc),y(loc),z(loc),marker{i},markerfacecolor,color{i});

end

str{1}=三维图;

legend(str,location,best)