为何SciPy的Cubic Spline生成不符合预期的曲线形状？

问题背景 #

尝试用SciPy的CubicSpline绘制类拱形曲线，但插值结果（黑色线条）在控制点间出现不符合逻辑的形状，预期是类似点4到点8的平滑拱形走势（绿色线条）。

复现代码 #

from scipy.interpolate import CubicSpline
import numpy as np
import matplotlib.pyplot as plt

# 点2、4、8、13、15的坐标
x = [-0.0243890844, -0.0188174509, -0.00021640210000000056, 0.0202699043, 0.0239562802]
y = [-0.0117638968, 0.00469300617, 0.0177650191, 0.00215831073, -0.0154924048]

cs = CubicSpline(x, y)
dsX = np.linspace(x[0], x[-1], num=1000)

plt.plot(dsX, cs(dsX), 'k')
plt.plot(x, y, 'mo')
plt.show()

效果对比 #

问题成因 #

CubicSpline默认使用自然边界条件（bc_type='natural'），该条件强制曲线在首尾端点的二阶导数为0。你的控制点首尾Y值均为最低点，自然边界会让曲线在两端平缓延伸，为满足二阶导数连续的要求，中间区域会出现不符合预期的扭曲，与你想要的单峰拱形走势冲突。

此外，控制点中峰值点（第三个点）后Y值快速下降，自然三次样条为了适配边界条件，会产生过度弯曲的形状。

修复方法 #

修改CubicSpline的边界条件，或改用保形插值方法：

1. 修改边界条件为夹紧型 #

夹紧边界（bc_type='clamped'）强制首尾端点的一阶导数为0，让曲线在端点处水平收尾，更贴合拱形的形态：

cs = CubicSpline(x, y, bc_type='clamped')

2. 手动指定首尾一阶导数 #

若你有明确的端点斜率需求，可直接传入导数参数。例如指定首尾斜率为0：

# 格式为((导数类型, 导数数值), (导数类型, 导数数值))，1代表一阶导数
cs = CubicSpline(x, y, bc_type=((1, 0), (1, 0)))

3. 使用保形插值（Pchip） #

PchipInterpolator会保留控制点的单调性和趋势，避免过度扭曲，适合需要贴合直观形状的场景：

from scipy.interpolate import PchipInterpolator
cs = PchipInterpolator(x, y)

是否遗漏配置参数？ #

是的，你未指定bc_type参数。默认的自然边界条件不匹配你的拱形需求，修改该参数是核心修复点。

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