通过四元数求机器人本体坐标旋转量

是的，通过两次姿态数据（以四元数表示）的差值，可以确定机器人在两个时刻之间的旋转角度变化。具体步骤如下：

1. 获取四元数：假设两个时刻的四元数分别为 ( q_1 ) 和 ( q_2 )。
2. 计算四元数的差值 ：
   • 将四元数 ( q_1 ) 的逆（反转）表示为 ( q_1^{-1} )。
   • 进行四元数乘法：( q_{\Delta} = q_2 \cdot q_1^{-1} )，得到差值四元数 ( q_{\Delta} )。
3. 转换为旋转角度 ：
   • 从差值四元数 ( q_{\Delta} ) 中提取旋转轴 ( \mathbf{v} ) 和旋转角度 ( \theta )。
   • 四元数 ( q_{\Delta} = \left(\cos\left(\frac{\theta}{2}\right), \sin\left(\frac{\theta}{2}\right)\mathbf{v}\right) )。

对于具体的计算，你可以使用以下公式：

• 四元数的逆：( q_1^{-1} = (q_{w1}, -q_{x1}, -q_{y1}, -q_{z1}) )（假设 ( q_1 = (q_{w1}, q_{x1}, q_{y1}, q_{z1}) )）。
• 四元数乘法：( q_{\Delta} = q_2 \cdot q_1^{-1} )。

最后，从 ( q_{\Delta} ) 提取旋转轴和旋转角度。

import numpy as np

def quaternion_conjugate(q):
    q = np.array(q)
    return np.array([q[0], -q[1], -q[2], -q[3]])

def quaternion_multiply(q1, q2):
    w1, x1, y1, z1 = q1
    w2, x2, y2, z2 = q2
    w = w1*w2 - x1*x2 - y1*y2 - z1*z2
    x = w1*x2 + x1*w2 + y1*z2 - z1*y2
    y = w1*y2 - x1*z2 + y1*w2 + z1*x2
    z = w1*z2 + x1*y2 - y1*x2 + z1*w2
    return np.array([w, x, y, z])

def quaternion_to_axis_angle(q):
    if q[0] > 1:
        q = q / np.linalg.norm(q)  # normalize the quaternion if needed
    angle = 2 * np.arccos(q[0])
    s = np.sqrt(1 - q[0]**2)  # assuming q[0] is the scalar part
    if s < 0.001:
        x = q[1]
        y = q[2]
        z = q[3]
    else: 
        x = q[1] / s
        y = q[2] / s
        z = q[3] / s
    return angle, np.array([x, y, z])

# 示例四元数
q1 = [0.9659, 0, 0.2588, 0]  # 初始姿态
q2 = [1, 0, 0, 0]            # 最终姿态

# 计算 q1 的逆
q1_inv = quaternion_conjugate(q1)

# 计算差值四元数 q_Δ
q_delta = quaternion_multiply(q2, q1_inv)

# 提取旋转轴和角度
angle, axis = quaternion_to_axis_angle(q_delta)

print(f"旋转角度: {np.degrees(angle)} 度")
print(f"旋转轴: {axis}")

通过以上计算，你可以得到机器人在两个时刻之间的旋转角度和旋转轴。🤖✨