numpy两四元数旋转累加

对于numpy的数组，可以使用Quaternions库来进行四元数的旋转操作。以下是使用numpy和Quaternions库进行两个四元数的旋转累加的示例代码：

import numpy as np
from Quaternions import Quaternions

# 创建两个四元数
quat1 = Quaternions(np.array([[1, 0, 0, 0]]))  # 第一个单位四元数
quat2 = Quaternions(np.array([[0, 1, 0, 0]]))  # 第二个单位四元数

# 对第一个四元数进行旋转
rot_quat1 = Quaternions.from_angle_axis(np.pi/4, np.array([0, 0, 1]))  # 绕z轴旋转π/4弧度
quat1_rotated = quat1 * rot_quat1  # 旋转后的四元数

# 对第二个四元数进行旋转
rot_quat2 = Quaternions.from_angle_axis(np.pi/2, np.array([1, 0, 0]))  # 绕x轴旋转π/2弧度
quat2_rotated = quat2 * rot_quat2  # 旋转后的四元数

# 两个旋转后的四元数累加
quat_sum = quat1_rotated + quat2_rotated

# 打印结果
print("quat_sum:", quat_sum)

运行上述代码，将得到两个旋转后的四元数分别为quat1_rotated和quat2_rotated，以及累加后的四元数quat_sum。