用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=f（x）。这意味着：对于x中的每一个点，y中必须有一个点。在

以sin(x)的数据为例：

x中的每个点在y中正好有一个点，并且1d样条函数很好地工作。在

现在把你提出的数据进行比较：

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