Go Gonum 插值器在数据区间之外怎样处理预测值

不能把Gonum所有插值器的区间外行为概括为自动线性外推。interp包说明将样本区间之外的返回值视为未定义的外推语义；具体实现可能仍返回一个数。当前PiecewiseLinear实现对左侧越界返回首个y，对右侧越界返回最后一个y，应用应明确选择接受这种端点保持、拒绝越界，还是采用另行定义的外推模型。

官方地址：https://www.gonum.org/

假设校准表只有x=0、2、5三个点，对应y=10、14、20。业务查询x=6时，即使Predict返回20，也不说明这个位置的真实值就是20。问题在于结果属于哪个数据范围，而不是函数是否“有返回值”。下面限定PiecewiseLinear，给出拒绝与夹取两种显式策略。

先区分包级合同与具体实现的端点行为

样本区间为[xs[0], xs[n-1]]，左右端点本身仍在区间内。区间内插值依据相邻样本，区间外已经离开样本覆盖范围。包级说明不提供统一的外推可靠性承诺，因此切换插值器或升级版本时，不能只沿用之前观察到的一个端点结果。

当前PiecewiseLinear在区间外使用端点值。它不会因为名字里有Linear，就自动沿着首尾线段延伸直线。当前PiecewiseCubic的Predict也有端点值处理，但不能由此推断所有Predictor或所有导数方法都遵循同一业务策略。

第一张图按“样本合同”和“预测对象”分组：xs、ys、有效区间属于样本合同；PiecewiseLinear、Predict与端点值属于预测对象。xs与ys组成拟合数据，xs关联有效区间，插值器关联两组样本，Predict关联插值器和端点值。普通连线仅说明静态依赖。

拟合之前，保存范围并检查输入合同

PiecewiseLinear.Fit要求至少两个点、xs与ys等长、横坐标严格递增。不满足相应前提会panic，不能只写一个if err!=nil就认为覆盖了所有输入问题。重复横坐标、NaN和无穷值也不应含糊地交给业务预测路径。

本文包装函数在构造时检查输入和斜率，保存最小与最大横坐标。以后查询只需要先判断是否在这个范围内，无需从预测值倒猜越界。极端浮点值仍可能使计算溢出，因此结果也检查有限性；这个例子不是任意数值范围下的稳定性证明。

显式选择拒绝或夹取，并保留结果来源

Reject表示越界即返回错误，不提供业务预测值；Clamp表示先把查询x限制到边界，再调用插值器，并返回clamped标记。区间内返回inside。两种策略由业务选择，不能把夹取值静默标成普通插值结果。

下面是完整main示例，需要模块已加入 gonum.org/v1/gonum 依赖。代码只使用PiecewiseLinear，不为了配图运行模型。

package main

import (
    "fmt"
    "log"
    "math"

    "gonum.org/v1/gonum/interp"
)

type Mode int
const (
    Reject Mode = iota
    Clamp
)
type BoundedLinear struct {
    pl interp.PiecewiseLinear
    min, max float64
    mode Mode
}

func finite(x float64) bool { return !math.IsNaN(x) && !math.IsInf(x, 0) }

// 构造函数提前检查数据，调用方不得在构造期间并发修改输入。
func NewBoundedLinear(xs, ys []float64, mode Mode) (*BoundedLinear, error) {
    if mode != Reject && mode != Clamp {
        return nil, fmt.Errorf("未知边界策略")
    }
    if len(xs)  0 {
            dx, dy := xs[i]-xs[i-1], ys[i]-ys[i-1]
            // 同时防止极端数值相减或斜率计算溢出。
            if dx  m.max {
        if m.mode == Reject {
            return 0, "outside", fmt.Errorf("查询超出[%g,%g]", m.min, m.max)
        }
        source = "clamped"
        if x 

BoundedLinear应通过NewBoundedLinear构造，不应手工创建零值实例后调用Predict。模型与输入切片不能在拟合期间被其他协程修改；上线使用还要明确是否允许重新拟合以及模型替换方式。示例中构造后只读，不提供动态更新入口。

第二张图按“查询合同”和“结果合同”分组：查询x、有效区间、边界策略属于前者，错误、预测值与来源标记属于后者。查询x关联有效区间与边界策略，边界策略关联三种结果字段。它不是执行流程图，重点是调用方拿到哪些信息。

对照理论结果，确认你真正想要的行为

对本文样本，区间内相邻段的斜率均为2。按公式手工推导，各种策略的结果如下；代码未在本文制作过程中实际执行，表格不是实测日志。

若沿斜率2真正做线性外推，x=-1与6会分别得到8和22，与端点保持的10和20不同。本文没有实现这种外推，因为它引入了区间外趋势继续成立的假设。不能只为了消除越界错误，就把这个假设默认为业务事实。

Clamp包装看似与当前原生结果相同，价值在于策略和标记由应用自己掌握。未来更换预测器时，只要包装层合同不变，就仍能明确拒绝或夹取，避免依赖某个实现恰好返回端值。

集成时处理缺失、过期与策略变更

越界不等于数据缺失，也不等于时间过期。区间判断只覆盖x数值范围；如果横坐标表示时间，还需要单独判断采样时间与模型是否过期。如果数据存在大空洞，查询可能在整体区间内，却离最近样本很远，应另定最大间距合同。

拒绝策略适合要求查询必须有校准依据的接口；夹取可用于明确允许边缘保持的展示或近似场景，但调用方应看见clamped。需要区间外预测时，另写允许距离、趋势模型和误差假设，不要把它藏在一个无状态的默认分支里。

日志建议记录模型版本、样本范围、查询位置、策略和来源标记。不要只看“没有NaN”就判断正确，也不要返回零值后让调用方误以为预测就是0。出现error时，应先处理错误，再读取y。

边界清单与相关疑问

• 拟合输入等长、有限且横坐标严格递增。
• 保留样本有效区间，区分端点与真正越界。
• 选择Reject或Clamp，返回来源标记。
• 覆盖左越界、端点、内部、右越界与非有限查询值。
• 更换插值器时重新核对区间外值与导数合同。

区间外一定会返回NaN吗？不是，当前PiecewiseLinear会返回端点值。包级“未定义”不应被解读成保证NaN。

ClampedCubic是不是自动夹取查询区间？不能由名字这样理解。该名称涉及样条端点条件，应用的查询边界策略应独立表达。

三次插值器更平滑，外推就更可靠吗？平滑描述样本区间内的曲线性质，不自动提供区间外数据依据。外推仍需要额外模型假设。

技术依据：官方interp包合同、PiecewiseLinear实现、三次插值器实现。校准表、边界包装函数与理论对照为本文原创设计。