我有以下程序：

nknots = 4
x_i = [0, 1, 2, 3]
y_i = [1, np.exp(1), np.exp(2), np.exp(3)]
coeff = interpolate.make_interp_spline(x_i, y_i, bc_type="natural")

y_i(x) = a + b*(x - x_i) + c*(x - x_i)^2 + d*(x - x_i)^3

当我这样做时

print(coeff(x_i))

我只得到a值的数组：

[ 1.          2.71828183  7.3890561  20.08553692]

但是，我缺少b、c和d系数的数组。我怎么才能提取出来呢？或者是我遗漏了什么步骤？我阅读了有关make_interp_Spline的Scipy文档，但我不了解如何获取b、c和d系数。

推荐答案

我建议您查看interpolate.CubicSpline。如果你想要的是多项式系数，那就方便多了。使用您的变量：

spl = interpolate.CubicSpline( x_i, y_i )
spl.c

array([[-1.57973952e-01,  2.93118310e-01, -1.35144359e-01],
       [ 1.11022302e-16, -4.73921855e-01,  4.05433076e-01],
       [-3.01723742e-01, -7.75645598e-01, -8.44134377e-01],
       [ 1.00000000e+00,  5.40302306e-01, -4.16146837e-01]])

有关分段多项式系数的存储方式，请参阅PPoly文档。