C++实现均匀1D网格线性插值的优化方案（适配VS2017）

实现方案 #

因为你用的是均匀一维网格，完全不需要二分查找，直接通过数学计算就能O(1)复杂度拿到插值区间的左右下标，效率远高于二分查找的O(logN)。

核心逻辑

均匀网格的位置和数组下标是线性对应的：

• 先做边界校验，查询位置小于最小位置直接返回第一个值，大于最大位置直接返回最后一个值
• 计算单位步长：步长 = (最大位置 - 最小位置) / (数值数组长度 - 1)，你当前场景步长为1
• 把查询位置标准化到下标区间：标准化位置 = (查询位置 - 最小位置) / 步长
• 左边界下标就是标准化位置向下取整的结果，插值系数t就是标准化位置减去左边界下标
• 带入线性插值公式计算结果即可

C++标准库相关支持

• 线性插值计算：C17开始<cmath>头文件提供了std::lerp(a, b, t)函数，直接实现了a*(1-t) + b*t的逻辑，精度和稳定性比自己手写的公式更好，VS2017 15.8及以上版本支持该特性，只要在项目属性中开启C17编译标准即可使用
• 如果你后续需要适配非均匀网格的场景，可以用<algorithm>头文件的std::lower_bound查找左边界下标，该函数内部就是优化过的二分实现，不需要自己手写二分逻辑，减少出错概率

示例代码

#include <iostream>
#include <vector>
#include <cmath>

int main() {
    // 预定义的数值数组
    const std::vector<double> values = {0.1, 0.22, 0.35, 0.5, 0.78, 0.92};
    const double min_pos = 0.0;
    const double max_pos = 5.0;
    const int interval_count = values.size() - 1;
    const double step = (max_pos - min_pos) / interval_count;

    double target_pos = 2.3;
    // 边界截断处理
    if (target_pos <= min_pos) {
        std::cout << values[0] << std::endl;
        return 0;
    }
    if (target_pos >= max_pos) {
        std::cout << values.back() << std::endl;
        return 0;
    }

    // 均匀网格直接计算下标，O(1)复杂度
    double normalized_pos = (target_pos - min_pos) / step;
    int left_idx = static_cast<int>(std::floor(normalized_pos));
    double t = normalized_pos - left_idx;

    // 插值计算二选一：
    // 1. 自写公式，兼容所有C++版本
    double result = values[left_idx] * (1 - t) + values[left_idx + 1] * t;
    // 2. C++17标准库函数，更稳定
    // double result = std::lerp(values[left_idx], values[left_idx + 1], t);

    std::cout.precision(4);
    std::cout << result << std::endl; // 输出0.395，和示例计算结果一致
    return 0;
}

内容的提问来源于stack exchange，提问作者Nick X Tsui