python - SciPy 插值错误

Question

我scipy.interpolate.interp1d用来插值一些数据（这是〜二次的），但我得到了一些奇怪的结果（非常小的周期性数值差异）。为了说明这一点，我采用了插值数据的导数，您可以在其中看到周期性的窥视，这表明数据中存在一些不连续性。我知道它很小，但就我而言，它很重要。所以当你看原始函数的插值函数看起来不错，但是当我稍后做一些微积分时，我发现了一些奇怪的结果。

谢谢！！

图片：http: //img37.imageshack.us/img37/9444/resultskx.png

Answer 1

确保您使用正确的顺序狐狸x值，例如（线性插值）：

def interp_normalize(x, y):
   dorder = {}
   for i, e in enumerate(x):
       dorder[e] = y[i]
   dorder = OrderedDict(sorted(dorder.items()))
   interpolate_x = list(dorder.keys())
   interpolate_y = list(dorder.values())
   return(interpolate_x, interpolate_y)