关于利用matlab制作三维图形的经典例子.docVIP

附录 附录一： 程序 1． 城区各功能区的地形二维等高线污染浓度散点图分布图程序 clc clear close all D=load(坐标数据．dat); M=load(浓度含量shuju．dat); x=D(:,2); y=D(:,3); z=D(:,4); c=D(:,5); xi=linspace(min(x),max(x),300); yi=linspace(min(y),max(y),300); [xi,yi]=meshgrid(xi,yi); zi=griddata(x,y,z,xi,yi); ci=griddata(x,y,c,xi,yi); marker={d,^,s,o,V}; color={b,y,c,r,m}; mat={As,Cd,Cr,Cu,Hg,Ni,Pb,Zn}; str={等高线,生活区,工业区,山林区,交通区,绿地区}; for j=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]); legend( 生活区,工业区,山林区,交通区,绿地区); %legend( {k,r,y,c,b}); title([重金属 ,mat{j}, 空间二维等高线分布图]) xlabel(X---轴) ylabel(Y---轴) colormap summer colorbar grid on hold on for i=1:5 loc=c==i; plot(x(loc),y(loc),marker{i},markerfacecolor,color{i},MarkerEdgeColor,color{i}); end % legend(str,location,best) for k=1:length(x) text(x(k)-200,y(k)+200,num2str(M(k,j+1)),fontsize,8); end legend( 功能一区——生活区,功能二区——工业区,功能三区——山林区,功能四区——交通区,功能五区——绿地区); end 2． 重金属 ‘**’含量程序 clc; clear; load shuju1．txt A=shuju1; x=A(:,1);y=A(:,2);z=A(:,3); scatter(x,y,5,z)%散点图 figure(1) [X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200),linspace(min(y),max(y),200),v4);%插值 pcolor(X,Y,Z); shading interp%伪彩色图 figure(2),contourf(X,Y,Z) %等高线图 legend(污染程度分布); title(重金属As含量空间分布平面图) xlabel(x—轴),ylabel(y—轴);%zlabel(z—轴) figure(3),surf(X,Y,Z)%三维曲面 legend(污染程度分布); title(重金属As含量空间分布平面图) xlabel(x—轴),ylabel(y—轴);zlabel(z—轴) 3 重金属综合污染指数程序 %Cd综合污染指数 Grandvalue=190; maxvalue=0; sum=0; C=[ 352．60 556．20 113．70 240．50 120．30 514．10 246．30 158．00 211．30 236．50 ]; % 可以输入任意长度的污染浓度数据 C=C n=length(C); for i=1:n; if C(i)maxvalue maxvalue= C(i); else sum=sum+ C(i)/Grandvalue; end end p=sqrt(((sum/n)^2+(maxvalue/Grandvalue)^2)/2) 附录二： 污染状况空间分布图