如果安装了scipy version>;=0.18.0,则可以使用scipy.interpolate中的cubic spline函数进行三次样条曲线插值。
您可以通过在python中运行以下命令来检查scipy版本:#!/usr/bin/env python3
import scipy
scipy.version.version
如果scipy版本为>;=0.18.0,则可以运行以下三次样条曲线插值的示例代码:#!/usr/bin/env python3
import numpy as np
from scipy.interpolate import CubicSpline
# calculate 5 natural cubic spline polynomials for 6 points
# (x,y) = (0,12) (1,14) (2,22) (3,39) (4,58) (5,77)
x = np.array([0, 1, 2, 3, 4, 5])
y = np.array([12,14,22,39,58,77])
# calculate natural cubic spline polynomials
cs = CubicSpline(x,y,bc_type='natural')
# show values of interpolation function at x=1.25
print('S(1.25) = ', cs(1.25))
## Aditional - find polynomial coefficients for different x regions
# if you wan