Go Gonum quad 怎样使用 Hermite 积分节点

如果你想用 Go Gonum quad 处理从负无穷到正无穷的积分，且被积表达式外面带有 e-x² 权重，应该选择 quad.Hermite{}。它不是把有限区间自动切成若干段，而是提供一组适用于 Hermite 求积的节点和权重。

官方文档：https://pkg.go.dev/gonum.org/v1/gonum/integrate/quad

最小可用写法是把边界设为 math.Inf(-1) 和 math.Inf(1)，传入正数采样点数，并让回调只负责计算 f(x)。Hermite 实际近似的是 ∫e-x²f(x)dx，边界和权重模型不能照搬有限区间的 Legendre 用法。

下面按“现象—验证—修复”的顺序，把 quad.Fixed、quad.Hermite 与手动读取节点的关系拆开。

先确认 Hermite 到底在计算哪类积分

Gonum 文档给出的 Hermite 规则对应全实数范围上的带权积分：∫-∞∞ e-x²f(x)dx。因此，f 回调里应该放权重之外的函数部分，而不是再把同一个平方指数权重重复乘一次。

这也是最容易误判的地方：如果业务目标是有限区间上的普通积分，优先考虑 quad.Legendre{}；如果目标是正态分布期望，也要先把概率密度的常数因子和指数形式整理清楚，再决定回调返回什么。

用 quad.Fixed 传入 Hermite 节点规则

quad.Fixed 的参数顺序是函数、积分下界、积分上界、点数、规则和并发度。Hermite 规则要求边界是两个无穷值，点数 n 必须为正数；并发度小于等于零时按串行方式调用回调。

package main

import (
    "fmt"
    "math"

    "gonum.org/v1/gonum/integrate/quad"
)

func main() {
    // 这里的 f 只表示权重之外的函数部分，权重由 Hermite 规则处理。
    f := func(x float64) float64 {
        return x * x
    }

    // Hermite 规则的边界必须覆盖整个实数轴，n 是求积节点数量。
    value := quad.Fixed(
        f,
        math.Inf(-1),
        math.Inf(1),
        16,
        quad.Hermite{},
        0, // 并发度为 0，表示串行调用 f。
    )

    // x^2 的解析结果是 sqrt(pi)/2，可用于解释结果量级。
    fmt.Println(value)
}

这个例子表达的是 ∫e-x²x²dx，解析结果为 √π/2。示例里的 16 只是一个点数选择，不代表所有函数都能用同一个点数达到相同精度；实际项目应根据函数变化速度和可接受误差安排点数。

发现结果异常时，先查边界和权重

如果把 -5、5 传给 Hermite，问题不在“点数太少”，而在规则契约不匹配。Gonum 的 Hermite 实现会检查边界是否分别为负无穷和正无穷；有限区间应换用适合有限边界的规则。

另一个常见错误是把 e-x² 同时写进回调，又把 Hermite 当作普通采样器。这样会改变被积函数，导致结果多乘一次权重。排查时可以把目标积分先写成“规则权重 × 回调函数”的形式，再逐项对照代码。

回调还应尽量保持纯计算特征。若把计数器、共享切片或非并发安全的缓存放进 f，随后将并发度改为正数，回调可能同时被调用，结果就会依赖共享状态。

需要节点时直接读取并保持配对

如果除了积分值，还要查看采样节点，可以创建等长的 x 和 weight 切片，再调用 Hermite{}.FixedLocations。它会原地写入节点与权重；同一索引的两个元素必须始终配对。

package main

import (
    "fmt"
    "math"

    "gonum.org/v1/gonum/integrate/quad"
)

func main() {
    const n = 8
    x := make([]float64, n)
    weight := make([]float64, n)

    // FixedLocations 会原地填充 Hermite 节点和对应权重。
    quad.Hermite{}.FixedLocations(x, weight, math.Inf(-1), math.Inf(1))

    var sum float64
    for k := range x {
        // x[k] 与 weight[k] 属于同一次采样，不能交叉或单独排序。
        sum += weight[k] * (x[k] * x[k])
    }

    // sum 与 quad.Fixed 对同一规则的加权求和含义一致。
    fmt.Println(sum)
}

这里的手动求和只是为了观察规则内部的“节点—权重—函数值”关系。正式业务若只需要积分值，直接调用 quad.Fixed 更省事；只有在需要复用节点、记录采样或组合多个函数时，才适合读取切片。

点数、并发和边界的速查

最后可以用三个问题复查：积分是否真的带平方指数权重？边界是否是整个实数轴？回调是否只表达权重之外的函数，并且能安全地被重复调用？这三项都成立时，Hermite 节点才是在正确的问题上发挥作用。

相关问题

Hermite 能不能直接计算有限区间积分？

不适合。Hermite 的实现要求负无穷到正无穷的边界；有限区间应选择与有限边界假设匹配的规则。

为什么节点数增加后结果还会变化？

固定点求积是近似。函数变化越快、尾部贡献越复杂，通常越需要比较不同点数下的结果，而不能只凭一次计算判断精度。