如何在Julia中执行B样条插值？求解示例数据的yq值

下面提供两种经过验证的B样条插值实现方案，解决你获取yq的问题：

方案一：使用Interpolations包 #

这是Julia生态中常用的插值工具包，正确处理边界外插即可解决问题：

1. 安装并导入包

using Pkg
Pkg.add("Interpolations")
using Interpolations

2. 定义原始数据与查询点

x = collect(300:50:900)
y = rand(13)
xq = collect(290:4:910)

3. 创建带边界处理的B样条插值对象
   这里采用三次B样条，边界用线性外插处理超出原始x范围的查询点：

# 构建插值对象，Cubic表示三次B样条，Line(OnGrid())定义边界处理逻辑
itp = interpolate((x,), y, BSpline(Cubic(Line(OnGrid()))))
# 包装为支持外插的对象
etp = extrapolate(itp, Line())

4. 生成查询点对应的y值

yq = etp.(xq)

方案二：使用Diercks包 #

Diercks专注于B样条实现，调用方式更简洁：

1. 安装并导入包

using Pkg
Pkg.add("Diercks")
using Diercks

2. 定义数据并创建B样条曲线

x = collect(300:50:900)
y = rand(13)
xq = collect(290:4:910)

# 默认创建三次B样条曲线
bs = BSpline(x, y)

3. 计算查询点的y值

yq = evaluate(bs, xq)

常见失败原因说明 #

你之前尝试失败可能是以下原因：

• 未处理xq中超出原始x范围的外插点，导致报错
• 插值阶数与节点配置不匹配（比如误用OnCell代替OnGrid）
• 未正确调用包中的核心插值/评估函数

内容的提问来源于stack exchange，提问作者Manu Han