2011国赛a

国赛2011a matlab实现代码

用来实现将五个不同区域不同表示画出，便于对比，也可以加上前一篇的三维图形部分代码，实现三维图的不同分区标注。等高线也可以实现，一个语句，合理运用hold on 讲不通的图组合画，有不一样的效果。

data = xlsread('data.xls',1,'A2:M320');%与个人数据图表相匹配，注意修改
number = data(:,1);
x = data(:,2); y = data(:,3); z = data(:,4);
area = data(:,5);
As = data(:,6); Cd = data(:,7); Cr = data(:,8); Cu = data(:,9); Hg =
data(:,10); Ni = data(:,11); Pb = data(:,12); Zn = data(:,13);
background = [3.6,130,31,13.2,35,12.3,31,69];
delta = [0.9,30,9,3.6,8,3.8,6,14];
mu = 1.75;
test = background + mu.*delta;
x = data(:,2); y = data(:,3); z = data(:,4);
[x1,y1] =
meshgrid(linspace(min(x),max(x),200),linspace(min(x),max(y),200));
z_1 = griddata(x,y,data(:,4),x1,y1,'v4');
Index = 7;
% 元素含量分布
for Index = 6:13
figure(Index-5)
z1 = griddata(x,y,data(:,Index),x1,y1,'v4');
contour(x1,y1,z1)
hold on
xlabel('x/m')
ylabel('y/m')
title('元素含量分布示意图(μg/g)');
for i= 1:length(x)
switch area(i)
case 1 % % 1 生活区
plot3(x(i),y(i),z(i),'rh')
case 2 % % 2 工业区
plot3(x(i),y(i),z(i),'k*')
case 3 % % 3 山区
plot3(x(i),y(i),z(i),'k.')
case 4 % % 4 交通区
plot3(x(i),y(i),z(i),'k+')
case 5 % % 5 公园绿地区
plot3(x(i),y(i),z(i),'k^')
end
end
end