自制无math.h的C++正弦计算程序算法异常求助

看起来你的正弦值计算程序遇到了典型的逻辑偏差，导致结果完全偏离预期——输入45度得到-0.868，而正确值应该是√2/2≈0.7071。结合你的描述和代码片段，我整理了几个最可能的问题点和排查步骤：

1. 首先排查角度转弧度的实现（最可能的原因） #

45度转换成弧度的正确值是π/4 ≈ 0.785398，如果你的deg_to_rad函数逻辑出错，后续的所有计算都会跑偏。常见错误包括：

• 搞反了转换公式：正确公式是 度数 × π ÷ 180，如果写成度数 × 180 ÷ π，45度会变成≈2578弧度，这个值远超出正弦函数的周期范围，结果自然完全错误；
• 使用了精度极低的π值（比如只写3.14），或者错误的符号（比如用了-3.14）；
• 用整数除法导致精度丢失：比如代码里写了deg * PI / 180但PI是int类型，或者180是整数导致截断小数。

快速验证方法：单独调用deg_to_rad(45)并打印结果，看是否接近0.7854。如果这里错了，先修复这个函数。

2. 检查正弦级数的符号与项计算逻辑 #

你提到用了特定算法计算正弦值，大概率是泰勒级数展开。正弦的泰勒展开式（在x趋近于0时收敛）是：

sin(x) = x - x³/3! + x⁵/5! - x⁷/7! + ... + (-1)^n × x^(2n+1)/(2n+1)!

如果你的代码犯了以下错误，会直接导致结果符号错误或数值偏差：

• 符号没有交替：比如每一项都用减号，或者都用加号，结果会完全偏离；
• 阶乘函数溢出：你的factorial函数返回int类型，而13!（6227020800）已经超过了int的最大值（通常是2^31-1=2147483647），溢出后阶乘会变成负数，导致后续项的符号和数值彻底混乱；
• 幂次或阶乘计算错误：比如误用了偶数次幂，或者阶乘计算少乘了数。

排查建议：

• 把factorial的返回类型改成long long或者double，避免整数溢出；
• 临时在代码中打印级数每一项的数值和符号，看是否符合“加、减、加、减”的交替规律；
• 手动计算前3项的和（x - x³/6 + x⁵/120），代入x=0.785，结果应该≈0.7068，接近正确值。

3. 优化参数类型避免精度丢失 #

你的sin函数参数是int input（角度）和int decimal_count，但在级数计算中，应该用double类型的弧度值参与运算，避免整数运算导致的精度截断。比如如果计算x的幂次时用了int类型，会直接丢掉小数部分，结果肯定不对。

示例修正代码片段 #

给你一个简化的正确实现参考，你可以对比自己的代码找差异：

#include <iostream>
using namespace std;

const double PI = 3.141592653589793;

// 角度转弧度
double deg_to_rad(int deg) {
    return deg * PI / 180.0; // 用double除法，保证精度
}

// 计算阶乘（用long long避免溢出）
long long factorial(int n) {
    long long result = 1;
    for (int i = 1; i <= n; ++i) {
        result *= i;
    }
    return result;
}

// 泰勒级数计算正弦值
double sin(double x, int terms) {
    double result = 0.0;
    int sign = 1; // 交替符号，初始为正
    double x_pow = x; // 递推计算x的幂次，提升效率
    for (int n = 0; n < terms; ++n) {
        int exponent = 2 * n + 1;
        // 累加项
        result += sign * x_pow / factorial(exponent);
        // 更新幂次和符号
        x_pow *= x * x;
        sign *= -1;
    }
    return result;
}

int main() {
    int input_deg = 45;
    double rad = deg_to_rad(input_deg);
    // 用5项足够得到接近0.707的结果
    double sin_val = sin(rad, 5);
    cout << "sin(" << input_deg << "°) = " << sin_val << endl;
    return 0;
}

运行这段代码，输出应该接近0.7071，符合预期。

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