Python 最小二乘自然样条

Question

我正在尝试找到一个适合自然样条的数值包，该样条最小化加权最小二乘法。

scipy 中有一个程序包可以实现我想要的非自然样条曲线。

import numpy as np
import matplotlib.pyplot as plt
from scipy import interpolate, randn

x = np.arange(0,5,1.0/6)
xs = np.arange(0,5,1.0/500)

y = np.sin(x+1) + .2*np.random.rand(len(x)) -.1

knots = np.array([1,2,3,4])
tck = interpolate.splrep(x,y,s=0,k=3,t=knots,task=-1)
ys = interpolate.splev(xs,tck,der=0)

plt.figure()
plt.plot(xs,ys,x,y,'x')

Answer 1

这个 tar file 里面的 spline.py 文件来自 this page默认情况下进行自然样条拟合。 this page 上也有一些代码声称主要是你想要的。 pyD3D package还有一个natural spline function在其 pyDataUtils模块。最后一个对我来说看起来是最有希望的。但是，它似乎没有设置您自己的结的选项。也许如果您查看源代码，您可以找到纠正该问题的方法。

另外，我找到了this message在 Scipy 邮件列表上说使用 s=0.0(如在您给定的代码中)使样条曲线使用您的上述程序自然地根据消息的作者进行拟合。我确实找到了 this splmake可以选择进行自然样条拟合的函数，但在查看源代码后我发现它尚未实现。