3D cylinder code debugging

Question

Miles Hao Peng Su 2022-12-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1870122-3d-cylinder-code-debugging

评论： Miles Hao Peng Su 2022-12-4

Hi, I am trying to implement a 3D cylinder code, and here is the code:

nCS = 2; %number of cross section of cylinder
nNodes = 50; %number of nodes per cross section
rSC = size * ones(1, nNodes);       %set length of each vertices
theta = linspace(0, 2* pi, nNodes); %set angle of each vertices
polar(theta, rSC) 
zSC = linspace(0, hSC, nCS)';   %create z coordinates, with height hSC
[xSC,ySC] = pol2cart(theta, rSC);   %transform from polar to cartisen
XSC = repmat(xSC, nCS, 1);      %repeat for array of number of cross sections
YSC = repmat(ySC, nCS, 1);
ZSC = repmat(zSC, 1, nNodes);   %repeat for number of z-layers
xSC_lid = zeros(2, nNodes);     %creating the top and bottom lid
ySC_lid = zeros(2, nNodes);
zSC_lid = repmat([0, hSC], 1, nNodes);
X = [xSC_lid(1,:); XSC; xSC_lid(2,:)];
Y = [ySC_lid(1,:); YSC; ySC_lid(2,:)];
Z = [zSC_lid(1,:); ZSC; zSC_lid(2,:)]; 
surf(X, Y, Z)

The error message is

What is wrong and how to solve this?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Matt J 2022-12-3

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1870122-3d-cylinder-code-debugging#answer_1119122

编辑：Matt J 2022-12-3

在 MATLAB Online 中打开

Much of your code seems to be an unnecessary reinvention of the cylinder() command.

nCS = 2; %number of cross section of cylinder

nNodes = 50; %number of nodes per cross section

[X,Y,Z]=cylinder(ones(1,nCS),nNodes);

surf(X,Y,Z,'FaceAlpha',0.7)

patch('XData',X(1,:),'YData',Y(1,:), 'ZData',Z(1,:),'FaceColor','r');

patch('XData',X(nCS,:),'YData',Y(nCS,:), 'ZData',Z(nCS,:),'FaceColor','r');

6 个评论
显示 4更早的评论隐藏 4更早的评论

Miles Hao Peng Su 2022-12-4

编辑：Miles Hao Peng Su 2022-12-4

在 MATLAB Online 中打开

Hi! May I ask what is the function of the ones(1,nSC) in the code?

[X,Y,Z]=cylinder(ones(1,nCS),nNodes);

I am curious, because when I pass just 1 through returns me the same result. Like so:

[XSC,YSC,ZSC]=cylinder(1,nNodes);

Why do we need to pass radius 2 times? Oh and anyway I only need the top and bottom 2 slices, so I don't really know if I need the nCS there.

Miles Hao Peng Su 2022-12-4

Hi, I got any other qusetion >< with this method is it possible to generate a hollow cylinder? like the walls have a certain thickness and the middle is void. Or would you suggest to use multicylinder for this?

请先登录，再进行评论。