Draw sphere with slats

Question

viet le 2016-9-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/301697-draw-sphere-with-slats

评论： Thorsten 2016-9-2

I would like to draw a hemisphere with slats as attached figures. 12 identical slits created in hemisphere, and width slit and width slat are equal.

2 个评论
显示无隐藏无

michio 2016-9-2

Duplicate question? http://jp.mathworks.com/matlabcentral/answers/301578-create-a-hemisphere-with-slats

Stephen23 2016-9-2

Duplicate:

https://www.mathworks.com/matlabcentral/answers/301578-create-a-hemisphere-with-slats

@viet le: please do not post duplicate questions. This actually makes it harder for us to help you, because it makes it hard for us to keep track of what information and explanations you have given, and also what answer and advice you have already been given.

If you want to add information to a question than please add comments to your original question.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Thorsten 2016-9-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/301697-draw-sphere-with-slats#answer_233603

编辑：Thorsten 2016-9-2

在 MATLAB Online 中打开

function hemispherewithstripes
clf
Nstripes = 12;
stripecolor = 'r';
Nhemicircles = 2*Nstripes;
r = 1; % radius of half-sphere
Npoints = 100; % number of points along each hemicircle
angle = linspace(0, pi, Npoints);
% xyz coordinates of a single hemicircle
x = r*cos(angle) + 1;
y = zeros(size(angle));
z = r*sin(angle);
% rotate hemicirclec
%rotation angle:
theta = linspace(-pi/2, pi/2, Nhemicircles);
% xyz coordinates of rotated hemicircles; preallocate for speed
XYZ(Npoints, 3, Nstripes) = 0; 
for i = 1:Nhemicircles
  % rotation matrix along x-direction:
  Rotx = [1  0             0; 
          0  cos(theta(i)) sin(theta(i)); 
          0 -sin(theta(i)) cos(theta(i))];
  XYZ(:,:,i) = [x(:) y(:) z(:)]*Rotx;
  % uncomment to plot hemicircles
  % plot3(XYZ(:,1,i), XYZ(:,2,i), XYZ(:,3,i), 'k'), hold on
end
hold on
% plot patch between every two hemicircles 
for i=1:2:Nhemicircles
  patch([XYZ(:,1,i)' flipud(XYZ(:,1,i+1))'],...
        [XYZ(:,2,i)' flipud(XYZ(:,2,i+1))'],...
        [XYZ(:,3,i)' flipud(XYZ(:,3,i+1))'],...
        stripecolor, 'EdgeColor', stripecolor)
end
view(3) % set default 3D viewing axes