python - Python/Scipy 中的分段常数或 0 度样条插值

Question

尽管 scipy 的文档表明scipy.interpolate.UnivariateSpline将以 k<=5 的顺序运行，但实际上，对于 >=1 的顺序还有一个额外的约束。有谁知道用 Numpy/Scipy 实现 0 度样条或分段常数插值的方法？

Answer 1

您可以interp1d使用kind='zero'或kind='nearest'。

Answer 2

我尝试了一个简单的技巧1）对逐步常数函数进行数值积分，你会得到虚线2）使用1.度样条拟合3）计算样条的导数（方法.derivative（））

它已经奏效了，但它可能会更好。