使用SciPy插值数据

Question

我有两个分别对应于x和y值的数据数组，我想用三次样条插值。

我试图这样做，但是我的插值函数没有通过我的数据点。

    import numpy as np
    import matplotlib.pyplot as plt
    from scipy.interpolate import interp1d

    re = np.array([0.2,2,20,200,2000,20000],dtype = float)
    cd = np.array([103,13.0,2.72,0.800,0.401,0.433],dtype = float)

    plt.yscale('log')
    plt.xscale('log')
    plt.xlabel( "Reynold's number" )
    plt.ylabel( "Drag coefficient" )
    plt.plot(re,cd,'x', label='Data')

    x = np.linspace(0.2,20000,200000)
    f = interp1d(re,cd,kind='cubic')
    plt.plot(x,f(x))

    plt.legend()

    plt.show()

我最终看到的是这样的；

这显然是我功能的糟糕表现。 我在这里想念什么？

谢谢。

Answer 1

通过执行以下操作，您可以获得预期的结果（对数轴上的平滑样条）：

f = interp1d(np.log(re),np.log(cd), kind='cubic')
plt.plot(x,np.exp(f(np.log(x))))

这将在日志空间中构建插值并正确绘制它。 以线性比例绘制数据，以查看三次方必须如何翻转才能使尾部位于左侧。

Answer 2

您缺少的主要内容是轴上的对log缩放。 对于您输入的数据，显示的样条曲线并不是不合理的结果。 尝试使用plt.xscale('linear')而不是plt.xscale('log')绘制图。 至少在原始数据上，三次样条也许不是最佳插值技术。 更好的选择可能是在数据导入日志中进行插值。