清空筛选条件
清空筛选条件

How to do a 3D circle in matlab

71 次查看（过去 30 天）
john pag
john pag 2013-6-17
评论： Rudranarayan Kandi 2018-9-12
Hello,
I need a code for a 3D circle. I have the code but i can not to do 3D.
Thank you

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

采纳的回答

Sean de Wolski
Sean de Wolski 2013-6-17

更多回答（1 个）

Kye Taylor
Kye Taylor 2013-6-17
编辑：Kye Taylor 2013-6-17
I assume you do not have a parametric form for the circle in 3D, but you do have the equation for the plane where the circle lives in 3D. In that case, i would solve this problem by creating the data in the xy-plane
% Original points, original plane
t = linspace(0,2*pi);
x = cos(t);
y = sin(t);
z = 0*t;
pnts = [x;y;z];
% unit normal for original plane
n0 = [0;0;1];
n0 = n0/norm(n0);
Then, I would rotate the circle data into a new plane
% unit normal for plane to rotate into
% plane is orthogonal to n1... given by equation
% n1(1)*x + n1(2)*y + n1(3)*z = 0
n1 = [1;1;1];
n1 = n1/norm(n1);
% theta is the angle between normals
c = dot(n0,n1) / ( norm(n0)*norm(n1) ); % cos(theta)
s = sqrt(1-c*c); % sin(theta)
u = cross(n0,n1) / ( norm(n0)*norm(n1) ); % rotation axis...
u = u/norm(u); % ... as unit vector
C = 1-c;
% the rotation matrix
R = [u(1)^2*C+c, u(1)*u(2)*C-u(3)*s, u(1)*u(3)*C+u(2)*s
u(2)*u(1)*C+u(3)*s, u(2)^2*C+c, u(2)*u(3)*C-u(1)*s
u(3)*u(1)*C-u(2)*s, u(3)*u(2)*C+u(1)*s, u(3)^2*C+c];
% Rotated points
newPnts = R*pnts;
Visualize:
plot3(newPnts(1,:),newPnts(2,:),newPnts(3,:),'ko')
  3 个评论
显示 1更早的评论
Vince
Vince 2018-1-13
Fantastic answer for creating a circle located in arbitrary 3D space!
Rudranarayan Kandi
Rudranarayan Kandi 2018-9-12
very nice approach for plotting in a different plane with normal specified.!!!! Cheers!!!

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Data Distribution Plots 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by