0

我比较了用于创建 3D 旋转矩阵的 scipy 和 MATLAB 函数的输出,这两个程序显示了不同的结果。

第一个 SciPys scipy.spatial.transform.Rotation

>>> import numpy as np
>>> from scipy.spatial.transform import Rotation as R

>>> R.from_euler('zyz', [np.pi/2, np.pi/2, 0]).as_dcm()
array([[ 0.00000000e+00, -2.22044605e-16,  1.00000000e+00],
       [ 1.00000000e+00,  2.22044605e-16,  0.00000000e+00],
       [-2.22044605e-16,  1.00000000e+00,  2.22044605e-16]])


>>> R.from_euler('zyz', [0, np.pi/2, np.pi/2]).as_dcm()
array([[ 0.00000000e+00, -1.00000000e+00,  2.22044605e-16],
       [ 2.22044605e-16,  2.22044605e-16,  1.00000000e+00],
       [-1.00000000e+00,  0.00000000e+00,  2.22044605e-16]])

现在matlabs eul2rotm函数:

>> eul2rotm([pi/2 pi/2 0], 'ZYZ')

ans =

    0.0000   -1.0000    0.0000
    0.0000    0.0000    1.0000
   -1.0000         0    0.0000


>> eul2rotm([0 pi/2 pi/2], 'ZYZ')

ans =

    0.0000   -0.0000    1.0000
    1.0000    0.0000         0
   -0.0000    1.0000    0.0000

所以最后,应用旋转的顺序似乎不同。但是哪一个是正确的呢?我很感激任何帮助。

使用 SciPy 1.3.1、NumPy 1.15.4 和 Matlab R2018b 进行测试。

4

0 回答 0