JavaScript实现拉格朗日插值多项式计算结果异常问题排查

结论 #

你的拉格朗日插值代码实现没有任何逻辑错误，结果不符合预期是对二次多项式插值的曲线走向存在错误预估导致的。

问题原因 #

拉格朗日插值生成的是穿过所有给定点的N次多项式（N=给定点数-1），二次多项式本质是抛物线，并不一定是单调曲线：
我们可以手动计算你给出的异常用例的正确结果验证：
给定点为 [[0, 100], [2.5, 66], [10, 33]]，求x=9时的y值，代入二次拉格朗日插值公式：

L(x) = y0*(x-x1)(x-x2)/[(x0-x1)(x0-x2)] + y1*(x-x0)(x-x2)/[(x1-x0)(x1-x2)] + y2*(x-x0)(x-x1)/[(x2-x0)(x2-x1)]

代入数值计算后得到结果约为31.42，和你的代码输出完全一致。
进一步推导该二次多项式的表达式为 y = 0.92x² -15.9x + 100，这是开口向上的抛物线，极值点（最低点）在x≈8.64的位置：

• x < 8.64 时，y随x增大递减
• x > 8.64 时，y随x增大递增
  所以x=9时y值低于x=10时的33是完全符合数学规律的结果。

替代方案 #

如果你需要插值结果符合「x越大、分数y越低」的单调递减预期，不建议使用高次拉格朗日插值，高阶插值容易出现振荡的龙格现象，推荐改用分段线性插值，实现简单且保证区间内单调：

function getScore (thresholds, macro) {
  // 按x坐标升序排序，避免输入顺序异常
  const sortedThresholds = [...thresholds].sort((a, b) => a[0] - b[0])
  const len = sortedThresholds.length
  // 超出左边界返回最左端点y值
  if (macro <= sortedThresholds[0][0]) {
    return sortedThresholds[0][1]
  }
  // 超出右边界返回最右端点y值
  if (macro >= sortedThresholds[len - 1][0]) {
    return sortedThresholds[len - 1][1]
  }
  // 查找x所在区间做线性插值
  for (let i = 0; i < len - 1; i++) {
    const [x0, y0] = sortedThresholds[i]
    const [x1, y1] = sortedThresholds[i + 1]
    if (macro >= x0 && macro <= x1) {
      return y0 + (y1 - y0) * (macro - x0) / (x1 - x0)
    }
  }
}

用该函数测试你的异常用例，x=9时输出约为37.4，符合你高于33的预期。

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