clear all
clc
dat=[ 0 40 44 50 40 40 42 40 40 42 40 40 40
1 40 42 43 44 40 40 40 40 41 41 42 40
2 41 42 40 40 40 41 42 40 40 40 40 40
3 42 40 40 40 40 40 40 40 40 40 40 40
4 40 40 40 40 40 40 40 40 40 40 40 40
5 41 42 42 43 41 41 42 41 40 40 40 40
6 40 40 40 40 40 40 40 40 40 40 40 40
7 41 42 43 41 42 40 40 40 40 40 40 40
8 40 45 45 41 40 40 40 40 40 40 40 40
9 41 42 41 42 41 42 40 40 40 40 40 40
10 42 42 40 40 40 40 50 40 45 40 40 40];
hei=dat(:,1);
R_mat= dat(:,2:end);
n_lay=length(hei);
n_len= length(R_mat(1,:));
the = pi/6*[1:n_len];
col_mm=jet(n_lay);
hold on
for i=1:n_lay
z=hei(i);
for j = 1 :n_len;
r=R_mat(i,j);
the = pi/6 * (j-1) ;
x = r*cos(the);
y = r*sin(the);
X(i,j)=x;
Y(i,j)=y;
Z(i,j)=z;
plot3([0 x],[0 y],[z z],'-','linewidth',2+0.1*i,'color',col_mm(i,:));
end
end
X=[X';X(:,1)']';
Y=[Y';Y(:,1)']';
Z=[Z';Z(:,1)']';
h_su=surf(X,Y,Z)
alpha(h_su,0.2)
pbaspect([1 1 2])
box on
view(-5,12)