如何计算3D向量夹角及圆锥朝向球体的三轴旋转角度

首先明确前提：圆锥局部坐标系中，锥尖沿Y轴向上（初始朝向向量为(0,1,0)），中心在原点；我们需要让锥尖始终指向球体，即圆锥的局部Y轴对齐从圆锥中心指向球体的向量。

你提到的单一夹角公式只能得到两个向量的夹角，但3D旋转需要明确旋转轴和角度，或者分解为三个欧拉角（绕X、Y、Z轴的旋转角度）。下面分两种方法说明：

方法一：直接分解为欧拉角（偏航角、俯仰角、滚转角） #

这种方法适合需要明确三个轴旋转角度的场景，注意欧拉角存在万向锁问题，但对于圆锥这类对称物体，滚转角可以灵活处理。

步骤如下：

1. 计算目标方向向量
   先算出从圆锥中心指向球体的向量，然后归一化：

   glm::vec3 cone_pos = ...; // 圆锥世界坐标
   glm::vec3 sphere_pos = ...; // 球体世界坐标
   glm::vec3 target_dir = glm::normalize(sphere_pos - cone_pos);
   

2. 计算偏航角（Yaw，绕Y轴旋转）
   偏航角控制圆锥在水平面上的朝向，先把目标方向投影到XZ平面，再计算与局部X轴的夹角：

   glm::vec3 xz_proj = glm::vec3(target_dir.x, 0.f, target_dir.z);
   if (glm::length(xz_proj) < 1e-6f) {
       // 目标在正上方或正下方，偏航角可设为0
       float yaw = 0.f;
   } else {
       glm::vec3 xz_normalized = glm::normalize(xz_proj);
       float yaw = std::atan2(xz_normalized.z, xz_normalized.x);
   }
   

3. 计算俯仰角（Pitch，绕X轴旋转）
   俯仰角控制圆锥上下方向的朝向，计算目标方向与XZ平面的夹角：

   float xz_length = glm::length(xz_proj);
   float pitch = std::atan2(target_dir.y, xz_length);
   

4. 计算滚转角（Roll，绕Z轴旋转）
   由于圆锥绕自身Y轴旋转不影响锥尖朝向球体，滚转角可直接设为0，或根据需求（比如让圆锥的局部X轴保持水平）调整，无需额外计算：

   float roll = 0.f;
   

方法二：轴角转欧拉角（更稳定，避免万向锁隐患） #

先通过轴角得到旋转，再转换为欧拉角，适合需要精确旋转的场景：

1. 计算旋转轴和角度
   初始朝向向量为up = glm::vec3(0.f,1.f,0.f)，目标方向为target_dir：

   glm::vec3 rotation_axis = glm::cross(up, target_dir);
   float rotation_angle = std::acos(glm::dot(up, target_dir));
   // 处理向量平行的特殊情况（夹角为0或180度）
   if (glm::length(rotation_axis) < 1e-6f) {
       rotation_axis = glm::vec3(1.f,0.f,0.f); // 随便选一个垂直轴
   } else {
       rotation_axis = glm::normalize(rotation_axis);
   }
   

2. 轴角转四元数
   四元数比欧拉角更稳定，避免万向锁：

   glm::quat rotation_quat = glm::angleAxis(rotation_angle, rotation_axis);
   

3. 四元数转欧拉角
   根据你的旋转顺序（比如X-Y-Z）转换为三个轴的旋转角度：

   glm::vec3 euler_angles = glm::eulerAngles(rotation_quat);
   float pitch = euler_angles.x;
   float yaw = euler_angles.y;
   float roll = euler_angles.z;
   

注意事项 #

• 注意坐标系的左右手系：不同引擎的坐标系规则不同，atan2的参数顺序、欧拉角的旋转方向可能需要调整。
• 角度单位：代码中计算的是弧度，若引擎使用角度，需转换为glm::degrees(angle)或std::rad2deg(angle)。
• 万向锁：当俯仰角接近±90度时，欧拉角会出现万向锁，此时偏航和滚转的区分会失效，这种情况下建议直接使用四元数控制旋转，而非欧拉角。

内容的提问来源于stack exchange，提问作者Summit