在 3D 空间中绘制一个旋转的平移圆

Question

Here is a little image 从我目前的情况来看。目前我正在尝试绘制一个以前面的红色矢量为中心的圆圈。圆圈应该接触蓝色矢量。 但我可以自己调整参数（当然还有其他圈子）。

所以我想要一个通用的方法来在 3d 空间中绘制一个具有任意中心点并围绕任意轴旋转的圆。

下面的代码是我目前的草稿。它围绕任意中心点绘制一个水平圆圈，但旋转部分仍然缺失。

set parametric
set urange [0:pi/2]
set vrange [0:2*pi]
 [..]
# radius
r = 0.1
# center (d,e,k) 
d = 0.1
e = 0.2
k = 0.3
fTx(u,v) = d+r*cos(v)
fTy(u,v) = e+r*sin(v)
fTz(u,v) = k         

splot [..]  ,\
      fTx(u,v), fTy(u,v), fTz(u,v)

Answer 1

您可以使用rotation matrices 旋转圆的点。例如，对于围绕 x 轴旋转角度 th，点 (fTx, fTy, fTz) 与 Rx(th) 相乘：

         ( 1     0         0    )
Rx(th) = ( 0  cos(th)  -sin(th) )
         ( 0  sin(th)   cos(th) )

对于旋转 th=24°，您的 Gnuplot 脚本可以扩展如下：

th = 24.0*pi/180
rotx(u,v) = fTx(u,v)
roty(u,v) = cos(th)*fTy(u,v) - sin(th)*fTz(u,v)
rotz(u,v) = sin(th)*fTy(u,v) + cos(th)*fTz(u,v)

splot rotx(u,v), roty(u,v), rotz(u,v)

绕另一个轴旋转是直线前进。

再想一想这个主题，在将旋转的圆圈移动到(d, e, k)之前先旋转可能会更容易：

set parametric
set urange [0:pi/2]
set vrange [0:2*pi]
# radius
r = 0.1
# center (d,e,k)
d = 0.1
e = 0.2
k = 0.3
# rotation angles (th, ..)
th = 24.0*pi/180.0

fTx(u,v) = r*cos(v)
fTy(u,v) = r*sin(v)
fTz(u,v) = 0      

rotx(u,v) = fTx(u,v)
roty(u,v) = cos(th)*fTy(u,v) - sin(th)*fTz(u,v)
rotz(u,v) = sin(th)*fTy(u,v) + cos(th)*fTz(u,v)

movx(u,v) = d + rotx(u, v)
movy(u,v) = e + roty(u, v)
movz(u,v) = k + rotz(u, v)

splot movx(u,v), movy(u,v), movz(u,v)