python - 使用 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 中出现两个不同的值。因此，此数据无效。