罗德里格斯公式Rodrigues: 代码: import numpy as np def rodrigues_rotation(r, theta): # n旋转轴[3x1] # theta为旋转角度 # 旋转是过原点的,n是旋转轴 r = np.array(r).reshape(3, 1) rx, ry, rz = r[:, 0