Go Gonum 插值的横坐标为什么必须严格递增

Go Gonum 插值的横坐标 xs 必须严格递增，意思是对每个相邻元素都满足 xs[i] > xs[i-1]。重复值、倒序值都不符合要求；这不是格式偏好，而是插值区间、斜率计算和预测定位共同依赖的输入契约。

官方地址：https://pkg.go.dev/gonum.org/v1/gonum/interp

如果直接把不满足条件的数据传给 PiecewiseLinear、NaturalCubic、AkimaSpline 或其他可拟合预测器，Gonum 的 Fit 可能以 panic 暴露问题，而不是把它当作普通的可恢复业务错误。因此，外部数据进入插值器前，应该先完成长度、有限数值和顺序检查。

严格递增到底检查什么

严格递增只看相邻横坐标，不是“整体看起来大致从小到大”。例如 {0, 0.5, 1.2, 2} 合法，{0, 0.5, 0.5, 2} 因为出现重复值不合法，{0, 2, 1, 3} 因为倒序也不合法。

同时，xs 和 ys 是成对样本：xs[i] 对应 ys[i]。所以检查顺序时不能只把 xs 单独排序，否则会把测量值绑定到错误的横坐标。长度也必须一致，并且至少有两个点；某些三次插值器还需要更多节点。

原因一：每个相邻点必须形成正宽度区间

插值器把相邻采样点视为区间端点。最简单的线性插值会使用：

// slope 表示相邻采样点之间的平均变化率。
dx := xs[i+1] - xs[i]
slope := (ys[i+1] - ys[i]) / dx

当 xs[i+1] == xs[i] 时，dx 为零，两个不同的样本可能被挤在同一个横坐标上，斜率没有唯一含义；当 xs[i+1] 时，区间方向反转，后续基于相邻区间的定位和系数计算也失去统一语义。

Gonum 的实现会在计算区间差值时拒绝 dx 。这就是为什么“非递减”还不够：允许相等会让区间退化。

原因二：预测要在有序横坐标中定位区间

拟合完成后，预测值需要找到包含查询点 x 的区间，逻辑可以抽象成 xs[i] 。Gonum 的内部定位使用二分查找，这要求 xs 已经有序；否则查找返回的索引不能稳定代表正确的相邻节点。

这也解释了为什么只在预测阶段排序不是好办法：预测器内部已经按照拟合时的横坐标建立了区间数据，事后更换顺序会让已保存的 ys、斜率或多项式系数与区间不再对应。

原因三：三次插值还要使用相邻间距建立系数

线性插值只直接依赖相邻斜率，三次插值还会使用多个区间的宽度和边界条件。以自然三次样条为例，区间宽度会进入带状线性方程组；如果某一段宽度为零，方程中的分母或系数就失去正常的数值意义。

因此，严格递增不是只针对某一个算法的限制。不同预测器的最低点数、边界条件和失败方式可能不同，但它们都需要非空、方向一致的横坐标区间。

在调用 Fit 前做一次可恢复检查

推荐把来自文件、数据库或网络的数据先转换为样本对，再检查有限数值和严格递增。下面的函数只负责输入契约，不调用 Gonum，也不会改变调用方的数据：

package main

import (
    "fmt"
    "math"
)

// validateSamples 检查样本数量、有限数值和横坐标的严格递增关系。
// 它返回可读错误，让边界层可以拒绝坏数据，而不是让 Fit 触发 panic。
func validateSamples(xs, ys []float64) error {
    if len(xs) != len(ys) {
        return fmt.Errorf("xs 与 ys 长度不一致: %d != %d", len(xs), len(ys))
    }
    if len(xs)  0 && xs[i] 

这里显式排除了 NaN 和无穷值。仅用 xs[i] 不能把 NaN 当作正常横坐标处理，因为浮点比较遇到 NaN 时不会提供可靠的顺序关系。

排序时必须把 xs 和 ys 绑在一起

如果数据只是乱序，但每个样本对本身可信，可以按横坐标对样本结构排序。重复横坐标则不能靠排序解决：排序只能把相同值放到一起，仍然需要业务规则决定保留哪个值、求平均，还是直接拒绝。

type sample struct {
    x float64
    y float64
}

// sortSamples 只重排完整样本对，避免 xs 与 ys 失去对应关系。
func sortSamples(samples []sample) {
    sort.SliceStable(samples, func(i, j int) bool {
        return samples[i].x 

上面的示例需要在完整程序中额外导入 sort。真实项目中更稳妥的做法是：先校验数值，再稳定排序，再重新校验严格递增，最后才把结果交给 Fit。

一个安全的 Gonum 调用边界

import "gonum.org/v1/gonum/interp"

// fitLinear 只在输入契约满足时创建预测器。
func fitLinear(xs, ys []float64) (interp.Predictor, error) {
    if err := validateSamples(xs, ys); err != nil {
        return nil, err
    }

    var model interp.PiecewiseLinear
    if err := model.Fit(xs, ys); err != nil {
        return nil, fmt.Errorf("拟合插值器失败: %w", err)
    }
    return model, nil
}

PiecewiseLinear.Fit 本身在满足数据契约时通常返回 nil；前置检查的价值是把外部输入错误变成调用方可以处理的 error。若改用 NaturalCubic 或其他三次模型，仍然保留同一个检查层，再根据模型的点数要求处理拟合错误。

常见误区与结论

1. 把“有序”写成“非递减”：重复横坐标会令 dx=0，必须使用严格大于。
2. 只排序 xs：这会破坏 (x, y) 样本对，必须连同 ys 一起重排。
3. 把 Fit 的 panic 当作普通 error：边界层应在 Fit 前主动校验，并在需要时增加恢复策略。
4. 只检查横坐标：NaN、无穷值、长度不一致和点数不足也应在数据入口处理。

一句话总结：Gonum 要求横坐标严格递增，是因为插值器必须把样本组织成方向一致的非空区间，并据此计算斜率、建立三次系数和用二分查找定位预测区间。修复数据时，先保持样本对，再决定排序或合并重复点，最后才调用 Fit。

相关问题

重复横坐标可以直接取平均吗？

只有当业务确认同一横坐标的多个观测值可以聚合时才可以。平均、取最后一次或按权重聚合都会改变数据含义，不能由插值器自动猜测。

为什么 Fit 失败却没有返回 error？

部分 Gonum 插值器把长度、点数和横坐标顺序视为调用前必须满足的契约，违反时会 panic；线性方程组等算法性失败才可能通过 error 返回。

横坐标间距必须相等吗？

不必须。严格递增允许间距不相等，关键是每个相邻差值都大于零。三次插值会把不同间距纳入系数计算。

可以在 Fit 后修改原始 xs 吗？

不要依赖这种做法。预测器已经按拟合时的数据建立内部结构；若数据需要变更，应重新校验并重新 Fit。

资料来源：https://pkg.go.dev/gonum.org/v1/gonum/interp；https://github.com/gonum/gonum/blob/v0.17.0/interp/interp.go；https://github.com/gonum/gonum/blob/v0.17.0/interp/cubic.go