基于CORDIC算法实现log10的C语言代码精度问题求助

针对你遇到的问题，核心原因大概率出在CORDIC迭代的方向逻辑错误或预计算常数表的准确性上，以下是具体分析和排查步骤：

1. 优先检查预计算的CORDIC常数表 #

log10的CORDIC迭代依赖预计算的d_k = log10(1 + 2^-k)（或对应负项），你遇到的两个特殊值刚好是1+2^-2=1.25和1+2^-3=1.125，对应的精确log10值为：

• log10(1.25) ≈ 0.09691001301
• log10(1.125) ≈ 0.05115252248

如果你的常数表中这两个值计算错误（比如误用自然对数ln()而非log10()，或精度不足），会直接导致输入这两个值时出现大误差，且迭代次数增加无法修正。

排查动作：打印你的d_k数组，对比上述精确值，确认是否一致。

2. 检查迭代方向与x缩放的对应关系 #

CORDIC计算log的核心逻辑是通过迭代缩放x趋近于1，同时累加对应log增量。正确的逻辑应为：

// 伪代码示例
float x = input;
float y = 0.0f;
for (int k = 0; k < ITER_COUNT; k++) {
    float factor = 1.0f + ldexpf(1.0f, -k); // 1 + 2^-k
    if (x > 1.0f) {
        // x大于1，除以factor让它趋近于1，同时累加log10(factor)
        y += d_k;
        x /= factor;
    } else {
        // x小于1，乘以factor让它趋近于1，同时减去log10(factor)
        y -= d_k;
        x *= factor;
    }
}
// 最终y即为log10(input)

如果你的代码中缩放方向搞反（比如x>1时乘以factor而非除以），会导致x越来越偏离1，y的累加值完全错误，且增加迭代次数只会让误差越来越大——这完全符合你描述的现象（1.125/1.25误差大，迭代次数无效）。

排查动作：单步调试x=1.25的第一次迭代，观察x的变化：

• 正确情况：x=1.25 → 除以1.5 → 约0.8333
• 错误情况：x=1.25 → 乘以1.5 → 1.875（方向完全错误）

3. 验证ldexpf的使用是否正确 #

你提到用ldexpf计算x·2ⁿ，需确认在迭代中是否用对了指数符号：

• 计算2^-k应使用ldexpf(1.0f, -k)，而非ldexpf(1.0f, k)；
• 如果是在输入预处理阶段使用ldexpf，由于你限定输入1≤x≤2（在收敛区间内），这部分应该没问题，但仍需确认是否误修改了x的初始值。

4. 排除浮点数精度以外的逻辑错误 #

由于1.125和1.25是float可精确表示的值，若迭代逻辑正确，即使迭代次数较少，误差也会逐步收敛。如果增加迭代次数后误差未减小，说明迭代过程根本没有向正确方向收敛，大概率是上述方向或常数表的问题。

总结 #

优先排查常数表的准确性和迭代时x的缩放方向，这两个点是导致你遇到的特殊值大误差、迭代无效的最常见原因。

内容的提问来源于stack exchange，提问作者emacs drives me nuts