Go Gonum FritschButland 与普通三次样条有什么区别

我第一次在 Go 里给一组单调测量值做曲线插值时，直觉上会把“更平滑”当成更好的答案。真正换到 Gonum 的 interp 包后，选择点其实很明确：如果数据的上升、下降方向不能被插值曲线破坏，优先看 FritschButland；如果你更关心函数值、一阶导数和二阶导数都连续，并且能接受自然边界条件，再考虑 NaturalCubic。

官方资料：https://pkg.go.dev/gonum.org/v1/gonum/interp

FritschButland 不是“低配三次样条”，而是把单调性和局部性放在前面、只保证 C¹ 的分段三次插值；NaturalCubic 则用两端二阶导数为零换取 C² 连续。两者没有绝对优劣，关键是业务数据是否允许曲线在节点之间出现额外趋势。

先把“普通三次样条”对应到 Gonum 类型

Gonum 并没有一个名为“普通三次样条”的类型。本文把这个说法具体化为 interp.NaturalCubic：它是分段三次插值器，函数值、一阶导数和二阶导数连续，并使用左端、右端二阶导数都为零的自然边界条件。

interp 里还有 ClampedCubic、NotAKnotCubic 等类型。如果你的“普通样条”实际指其中一种，边界条件就不一样，不能直接套用下面的结论。

两者真正不同的是约束目标

FritschButland 的官方描述有三个关键词：monotone、local、C¹。它不要求调用方提供导数，会根据节点值构造局部斜率；当相邻数据保持同一方向时，曲线会尽量保持该方向。代价是曲线张力较高，节点之间可能显得偏平，不会为了追求“圆润”而随意鼓出额外的峰谷。

NaturalCubic 关注的是更高阶的光滑性：相邻分段不仅函数值和一阶导数接上，二阶导数也连续。它通过整体的线性方程组求解分段系数，端点默认采用自然边界。这个模型适合把数据看成一条需要连续曲率的曲线，但它本身不是单调性约束器。

最小代码：同一组数据切换两种插值器

package main

import (
	"fmt"

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

func main() {
	// X 必须严格递增，Y 与 X 一一对应；这里使用持续上升的数据。
	xs := []float64{0, 1, 2, 3, 4}
	ys := []float64{10, 12, 15, 19, 24}
	x := 2.5

	var fb interp.FritschButland
	if err := fb.Fit(xs, ys); err != nil {
		// FritschButland 当前 Fit 通常返回 nil，但保留错误分支便于替换实现。
		panic(err)
	}

	var nc interp.NaturalCubic
	if err := nc.Fit(xs, ys); err != nil {
		// NaturalCubic 需要解方程组，失败时不能继续使用未完成的拟合结果。
		panic(err)
	}

	fbValue := fb.Predict(x)
	ncValue := nc.Predict(x)
	fbSlope := fb.PredictDerivative(x)
	ncSlope := nc.PredictDerivative(x)
	fmt.Printf("x=%.1f FB=%.4f slope=%.4f NC=%.4f slope=%.4f\n", x, fbValue, fbSlope, ncValue, ncSlope)
}

这段代码的重点不是比较某一个小数结果，而是展示 API 形状相同、内部目标不同：两者都实现 Fit、Predict 和 PredictDerivative，因此可以放在相同的业务接口后面；真正的差别要回到数据约束和边界条件。

如果只是要做单调报表、累计指标或校准曲线，我会先试 FritschButland。如果曲率连续本身会影响后续计算，例如需要稳定观察导数变化，则再评估 NaturalCubic，并用业务数据检查它是否满足单调性要求。

边界条件比类型名更容易踩坑

• X 不能重复或倒序。两种类型的 Fit 都要求 xs 严格递增，长度至少为 2，并且与 ys 等长。
• FritschButland 的 nil 不等于输入可以随便传。文档说明合法输入下它总是返回 nil，但长度、顺序和切片对应关系不满足要求时仍会触发 panic。
• NaturalCubic 的错误要保留。它还可能因为求解所需线性方程组失败而返回 error，生产代码不要忽略。
• 区间外不要把结果当成可靠外推。interp 文档明确说，插值区间之外的返回值未定义，只是尽量给出合理结果。调用前应限制查询范围，或自行定义区间外策略。

另外，单调性是数据属性，不是每一组数据都能被“曲线算法”自动修复。如果原始节点本身先升后降，FritschButland 会围绕节点趋势工作，而不是把一组非单调数据强行变成全局单调曲线。

我会怎样做最后选择

可以按下面的顺序判断：

1. 数据是否必须保持节点之间的方向？是，就优先 FritschButland。
2. 是否需要二阶导数连续，或者后续算法依赖连续曲率？是，再评估 NaturalCubic。
3. 是否能接受两端采用二阶导数为零的自然边界？不能，就不要把 NaturalCubic 当成默认答案，应查看 ClampedCubic 或 NotAKnotCubic。
4. 查询点是否可能落在节点区间外？是，先设计外推策略，再决定使用哪种插值器。

我的经验是：监测值、累计量、按区间增长的业务指标，先用 FritschButland 把“不要凭空造峰谷”作为底线；几何曲线或需要连续曲率的数值处理，再把 NaturalCubic 放进候选，并针对端点和单调性做单独测试。

相关问题

FritschButland 是不是线性插值？

不是。它仍是分段三次插值，只是通过局部斜率控制保持单调，并连续到一阶导数。

NaturalCubic 一定比 FritschButland 更平滑吗？

如果“平滑”指二阶导数连续，NaturalCubic 的定义更强；如果指不出现业务上不允许的额外趋势，则不能只看连续阶数，必须检查单调性和边界行为。

两种类型能否共用一个接口？

可以。它们都能按 Fitter 和 Predictor 的形状使用；只是调用方必须保留 Fit 的错误处理，并在接口文档中写清区间外行为。

参考资料

Gonum interp 官方文档：https://pkg.go.dev/gonum.org/v1/gonum/interp

Gonum 官方仓库：https://github.com/gonum/gonum