欧拉角转四元数及四元数乘向量计算出现精度异常问题

你遇到的-8.74228e-08这类极小值是单精度浮点数的固有计算误差，并非代码逻辑错误：

1. 欧拉角转弧度时，π（对应180度）是无理数，无法用浮点数精确表示，导致cos(π/2)的计算结果不是严格的0，而是接近0的极小值；
2. 这个极小值会传递到四元数的w分量，进而在四元数与向量的乘法运算中，产生x分量的微小误差；
3. glm出现同样问题是因为其底层也是基于浮点运算，这类误差是通用现象，不是库的bug。

1. 阈值过滤（最简便） #

直接将绝对值小于极小阈值（如1e-6）的分量置为0，消除浮点误差的影响：

Vector3 sanitizeVector(const Vector3& vec) {
    const float epsilon = 1e-6f;
    return Vector3(
        fabs(vec.x) < epsilon ? 0.0f : vec.x,
        fabs(vec.y) < epsilon ? 0.0f : vec.y,
        fabs(vec.z) < epsilon ? 0.0f : vec.z
    );
}

// 使用示例
Vector3 forward = quaternion * Vector3(0.0f, 0.0f, 1.0f);
forward = sanitizeVector(forward);

2. 特殊角度硬编码优化 #

针对180度这类特殊欧拉角，提前判断并生成精确的四元数，避免浮点计算误差：

Quaternion Quaternion::Euler(const Vector3& euler) {
    float xRad = Radians(euler.x);
    float yRad = Radians(euler.y);
    float zRad = Radians(euler.z);

    const float pi = 3.14159265358979323846f;
    const float eps = 1e-6f;

    // 处理Y轴180度旋转的特殊情况
    if (fabs(yRad - pi) < eps) {
        return Quaternion(0.0f, 1.0f, 0.0f, 0.0f);
    }
    if (fabs(yRad + pi) < eps) {
        return Quaternion(0.0f, -1.0f, 0.0f, 0.0f);
    }
    // 可扩展添加X/Z轴90、180、270度的特殊处理

    xRad /= 2;
    yRad /= 2;
    zRad /= 2;

    return Quaternion(
        cos(zRad) * cos(yRad) * sin(xRad) - sin(zRad) * sin(yRad) * cos(xRad),
        cos(zRad) * sin(yRad) * cos(xRad) + sin(zRad) * cos(yRad) * sin(xRad),
        sin(zRad) * cos(yRad) * cos(xRad) - cos(zRad) * sin(yRad) * sin(xRad),
        cos(zRad) * cos(yRad) * cos(xRad) + sin(zRad) * sin(yRad) * sin(xRad)
    );
}

3. 四元数归一化 #

浮点运算可能导致四元数模长偏离1，归一化后可减少后续乘法的误差累积：

Quaternion Quaternion::Normalize() const {
    float magnitude = sqrt(x*x + y*y + z*z + w*w);
    if (magnitude < 1e-6f) {
        return Quaternion(0.0f, 0.0f, 0.0f, 1.0f);
    }
    float invMag = 1.0f / magnitude;
    return Quaternion(x*invMag, y*invMag, z*invMag, w*invMag);
}

// 使用示例
Quaternion quat = Quaternion::Euler(rotation).Normalize();
Vector3 forward = quat * Vector3(0.0f, 0.0f, 1.0f);

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