数据集端点的数据点插值

不幸的是,我不熟悉scipy,但也许这会有帮助:

让我们假设,我们有一个数据集的点(如席,彝),I=1…N。从splerp函数中,我们可以在点'2'(df_second')和'n-1'(df_before_last)中找到一阶导数(实际上是有限差分)。我们需要两个线性方程组:

a1*x1^2 + b1*x1 + c1 = y1
a1*x2^2 + b1*x2 + c1 = y2
2*a1*x2 + b1 = df_second

和

a2*xn^2 + b2*xn + c2 = yn
a2*x(n-1)^2 + b2*x(n-1) + c2 = y(n-1)
2*a2*x(n-1) + b2 = df_before_last

解这个方程你会得到两个多边形:a1*x^2+b1*x+c1和a2*x^2+b2*x+c2,这两个多边形是插值的,并且可以很容易地访问样条函数。

实际上,只需在系统中加入附加方程,就可以提高多项式的阶数,并在第一点和最后一点设置插值图形的曲率和角度。