如何用Scipy CubicSpline控制一二阶导数并设置边界条件？

问题描述 #

能否使用CubicSpline连接两点并为其设置特定的一阶、二阶导数？我尝试了以下代码：

boundary_conditions = (((1, dxdy1_spline), (2, 0.0)), ((1, dxdy2_spline), (2, 0.0)))  
spline = CubicSpline(x_spline_values, y_spline_values, bc_type = boundary_conditions) 

运行时出现错误：ValueError: 指定的导数阶数必须是1或2。

问题原因 #

代码里的边界条件格式错误，多余的嵌套括号导致CubicSpline无法正确解析导数阶数和对应值的配对关系，从而抛出非法阶数的错误。

解决方案 #

调整边界条件的格式，去掉多余的嵌套，用列表或字典明确每个端点的导数条件：

写法一（列表形式）

boundary_conditions = [
    [(1, dxdy1_spline), (2, 0.0)],
    [(1, dxdy2_spline), (2, 0.0)]
]
spline = CubicSpline(x_spline_values, y_spline_values, bc_type=boundary_conditions)

写法二（字典形式）

boundary_conditions = (
    {'derivatives': [(1, dxdy1_spline), (2, 0.0)]},
    {'derivatives': [(1, dxdy2_spline), (2, 0.0)]}
)
spline = CubicSpline(x_spline_values, y_spline_values, bc_type=boundary_conditions)

注意事项 #

当仅用两个点构造三次样条时，样条会退化为三次多项式，此时每个端点最多可指定2个导数条件（一阶+二阶），只要格式正确就能正常生成样条，但要确保条件本身符合三次多项式的数学约束，避免出现无解情况。

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