python - 在 Scipy 中使用 splev 评估样条的导数

Question

我使用 splprep 从一组点创建了一个 bspline，如下所示:

tck,uout = splprep([x,y],s=0.,k=2,per=False)

现在，我正在尝试使用以下方法评估样条曲线的导数:

dx,dy = splev(uout,tck,der=1)

我发现 splev 返回两个导数列表。

鉴于 Spline 已参数化(比如在 u 中)，它是否返回 dx/du 和 dy/du ？

如果不是，如何正确评估导数 (dy/dx)？

Answer 1

是的，如果 der = 1，则列表是每个点的 dx/du 和 dy/du 的值。则梯度为 dy/dx = dy/du/dx/du。

我有点担心 splprep 调用:s 是可选的，但如果已定义，它的值应该与点数大致相同(越大意味着越平滑)。 per 是一个整数值，而不是 bool 值。三次样条比二次样条表现得更好。 http://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.splprep.html