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

附录

附录一：程序

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) 附录二：污染状况空间分布图