我试图在Scipy中插值后获得一个中间点。但是似乎没有任何直接方法可以做到这一点。谁能帮忙?
示例代码:
from scipy.interpolate import CubicSpline
import matplotlib.pyplot as plt
x = np.arange(10)
y = np.sin(x)
cs = CubicSpline(x, y)
xs = np.arange(-0.5, 9.6, 0.1)
plt.figure(figsize=(6.5, 4))
plt.plot(x, y, 'o', label='data')
plt.plot(xs, np.sin(xs), label='true')
plt.plot(xs, cs(xs), label="S")
plt.plot(xs, cs(xs, 1), label="S'")
plt.plot(xs, cs(xs, 2), label="S''")
plt.plot(xs, cs(xs, 3), label="S'''")
plt.xlim(-0.5, 9.5)
plt.legend(loc='lower left', ncol=2)
plt.show()
正如CubicsPline文档所说,它返回的对象是可召唤的,带有参数
- x:array_like
评估interpolant的点。- nu:int,可选
衍生物的评估顺序。- 外推:{bool,'周期性',none},可选
so, cs(5)
给出了在5处样条的值,而 cs(5, 2)
是其第二个衍生物,等等,等等。