使用 scipy sqlrep 进行定期插值

Question

我无法让 scipy.interpolate.splrep 处理周期性数据点... 我想我没有正确处理“per”参数。到目前为止，我使用此代码：

from scipy.interpolate import splrep
import numpy as np
x = np.concatenate((np.linspace(-180, 180, 720),np.linspace(180, -180, 720)))
y = np.concatenate((np.linspace(-180, 180, 720),np.linspace(-180, 180, 720)))
spl = splrep(x, y, per=True)

但它不起作用，引发“输入数据错误”

Answer 1

正如the documentation 所述，数据点x 和y 需要定义一条曲线y = f(x)。这意味着：x 中的每个点都需要在y 中恰好有一个 点。

举个例子sin(x)的这个数据：

x 中的每个点在y 中都有一个点，并且一维样条拟合效果很好。

现在比较（类似于）您提出的数据：

x 中的每个点出现两次，更糟糕的是，y 中出现两个不同的值。因此，该数据无效。