正负弧度下π计算及C++验证A与B相差π（精度0.01）的方法

1. 如何处理可为正负的弧度对应的π相关计算？ #

• 优先使用标准库π值：C++11及以上版本中，在定义_USE_MATH_DEFINES后，可通过<cmath>获取M_PI常量，避免硬编码3.14159等近似值带来的累积误差。
• 利用周期性简化运算：三角函数具有2π的周期性，正负弧度可通过加减2π的整数倍转换到[-π, π]或[0, 2π]区间，减少计算复杂度。例如用fmod(theta, 2*M_PI)获取等效周期内的弧度，再根据符号调整到目标区间。
• 保留符号的π偏移：计算与原弧度相差π的角度时，直接执行theta + M_PI或theta - M_PI即可，无需额外处理符号——两种结果在周期意义上等价，但保留符号能反映原始角度的方向信息。

2. C++中验证正负弧度角相差π（精度0.01） #

由于浮点数存在精度误差，直接比较差值等于±π不可靠，需结合周期特性与容差判断：

核心逻辑 #

两个角相差π弧度，等价于它们的差值在模2π后，结果的绝对值与π的差小于0.01。因为差值可能为π、-π、3π、-3π等（均为π加2π的整数倍），模2π后会落到[-π, π]区间，此时绝对值应接近π。

代码实现 #

#define _USE_MATH_DEFINES
#include <cmath>
#include <iostream>

const double TOLERANCE = 0.01;

bool isDifferencePi(double A, double B) {
    double diff = B - A;
    // 将差值映射到[-π, π]区间
    double wrapped_diff = fmod(diff, 2 * M_PI);
    if (wrapped_diff > M_PI) {
        wrapped_diff -= 2 * M_PI;
    } else if (wrapped_diff < -M_PI) {
        wrapped_diff += 2 * M_PI;
    }
    // 验证绝对值与π的误差在容差内
    return std::fabs(std::fabs(wrapped_diff) - M_PI) < TOLERANCE;
}

int main() {
    // 测试用例
    double A1 = M_PI / 2;
    double B1 = A1 + M_PI;
    std::cout << std::boolalpha << isDifferencePi(A1, B1) << "\n"; // true

    double A2 = -M_PI / 4;
    double B2 = A2 - M_PI;
    std::cout << std::boolalpha << isDifferencePi(A2, B2) << "\n"; // true

    double A3 = 3 * M_PI;
    double B3 = A3 + 3 * M_PI;
    std::cout << std::boolalpha << isDifferencePi(A3, B3) << "\n"; // true

    double A4 = 1.0;
    double B4 = 1.0 + 3.14;
    std::cout << std::boolalpha << isDifferencePi(A4, B4) << "\n"; // true

    double A5 = 0.0;
    double B5 = 0.0 + 3.13;
    std::cout << std::boolalpha << isDifferencePi(A5, B5) << "\n"; // false
    return 0;
}

注意事项 #

• 若编译器无法识别M_PI，需在包含<cmath>前定义_USE_MATH_DEFINES宏。
• 容差TOLERANCE可根据需求调整，此处设为0.01以满足精度要求。

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