Go Gonum interp 怎样在插值位置计算导数

在 Go 里用 Gonum interp 计算插值位置的导数，关键不是对 Predict(x) 的结果再做一次差分，而是选择实现了 interp.DerivativePredictor 的插值器。以 interp.NaturalCubic 为例，先用严格递增的 xs 和等长的 ys 完成 Fit，然后在同一个位置分别调用 Predict(x) 与 PredictDerivative(x)。

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

函数值回答“曲线在这里是多少”，导数回答“曲线在这里朝哪个方向、以多快的速度变化”。两者都来自已经拟合好的插值器，但语义和单位不同；查询位置最好限制在采样区间内。

先选能提供导数的插值器

Gonum interp 的基础接口 interp.Predictor 只承诺有 Predict(x)。如果文章标题中的“导数”是业务结果的一部分，就要选择同时实现 PredictDerivative(x) 的类型，也就是 interp.DerivativePredictor。官方文档列出的 AkimaSpline、ClampedCubic、FritschButland、NaturalCubic、NotAKnotCubic 和 PiecewiseCubic 都提供了导数查询能力，但它们的边界条件和使用前提并不相同。

一个实用的选择方式是：没有额外导数数据、需要平滑曲线时先看 NaturalCubic 或其他三次插值器；需要保持单调性时再考虑 FritschButland；已经有每个节点的导数时，才考虑 PiecewiseCubic.FitWithDerivatives。不要因为类型名称里有“Cubic”就默认所有边界行为一致。

准备满足 Fit 契约的采样数据

Fit 是查询前的接口边界。xs 与 ys 必须长度一致，横坐标至少有两个点，并且 xs 必须严格递增。文档明确说明，输入不满足这些条件时部分实现会 panic；因此来自文件、网络或用户输入的数据，应在进入插值器前由调用方完成检查。

package main

import (
    "fmt"

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

func main() {
    // 横坐标必须严格递增，纵坐标与它保持一一对应。
    xs := []float64{0, 1, 2, 3, 4}
    ys := []float64{0, 0.8, 2.7, 4.1, 5.0}

    var spline interp.NaturalCubic
    // Fit 失败时不要继续查询未完成的插值器。
    if err := spline.Fit(xs, ys); err != nil {
        panic(fmt.Errorf("拟合插值曲线失败: %w", err))
    }

    queryX := 2.5
    // 同一个位置分别得到曲线高度和局部变化率。
    value := spline.Predict(queryX)
    slope := spline.PredictDerivative(queryX)
    fmt.Printf("x=%.2f, y=%.6f, dy/dx=%.6f\n", queryX, value, slope)
}

这个示例没有把某个固定数值写成“标准输出”，因为插值算法、采样点和边界条件都会影响结果。真正要记录的是：拟合是否成功、查询位置属于哪个区间，以及函数值和导数分别代表什么。

在同一个位置读取函数值和导数

Predict(x) 返回插值函数在 x 处的估计值，PredictDerivative(x) 返回该位置的一阶导数估计。若原始数据的横坐标单位是秒、纵坐标单位是米，那么导数的单位就是米/秒；这也是为什么不能只看数值大小而忽略单位。

如果调用方只关心导数，可以让业务函数接收一个 interp.DerivativePredictor，这样接口会明确表达依赖：

package main

import (
    "fmt"

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

func reportSlope(p interp.DerivativePredictor, x float64) {
    // 先读取函数值，便于把趋势判断和当前状态一起记录。
    value := p.Predict(x)
    slope := p.PredictDerivative(x)
    fmt.Printf("位置 %.3f 的函数值=%.6f，导数=%.6f\n", x, value, slope)
}

func main() {
    // 使用同一套采样点建立可提供导数的插值器。
    xs := []float64{0, 1, 2, 3}
    ys := []float64{1, 1.4, 2.2, 3.6}
    var p interp.NaturalCubic
    if err := p.Fit(xs, ys); err != nil {
        panic(fmt.Errorf("拟合失败: %w", err))
    }

    reportSlope(&p, 1.5)
}

使用接口的好处是把“必须能够给出导数”写进函数签名，后续替换为其他实现时，调用方不需要依赖具体插值器名称。

处理区间边界与算法取舍

Gonum interp 的包文档说明，插值区间之外的返回值未定义，库会尽量返回合理结果，但调用方不能把这种结果当成可靠外推。生产代码应在查询前判断 xs[0] ，或者明确把区间外策略交给上层。

func inInterpolationRange(xs []float64, x float64) bool {
    // 空切片没有有效区间，避免访问首尾元素越界。
    if len(xs) == 0 {
        return false
    }
    // 这里只接受闭区间内的查询，区间外不假装是可靠外推。
    return xs[0] 

还要把算法选择和业务目标放在一起考虑：

常见问题

Predict 的结果能直接当导数吗？

不能。Predict 是函数值，PredictDerivative 才是插值器提供的一阶导数估计；两者可能数值接近，但物理含义和单位不同。

为什么 Fit 之前不能调用 PredictDerivative？

插值器需要先从采样点建立内部曲线。应先完成 Fit 并处理返回的错误，再把拟合后的对象交给查询逻辑。

区间外查询返回了一个数字，可以继续使用吗？

不建议直接使用。官方包文档把插值区间之外的结果定义为未定义；如果业务确实需要外推，应选择明确支持外推的方案并单独记录边界策略。

小结

Go Gonum interp 计算导数的最小路径是：选择实现 interp.DerivativePredictor 的类型，准备严格递增且长度匹配的 xs、ys，成功 Fit 后在目标位置同时调用 Predict 和 PredictDerivative。把区间检查、错误处理和算法取舍写在接口边界上，才能让导数结果真正适合后续的趋势判断或控制逻辑。