三次样条插值中途退出

Question

我正在尝试进行三次样条插值，由于某种原因，插值在其中间下降。这是非常神秘的，我无法在网上任何地方找到类似的事件。

这是我的论文所以我已经排除了一些标签等以使其故意模糊，但所有相关代码如下。对于背景，这是一个与天文学相关的情节。

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

W = np.array([0.435,0.606,0.814,1.05,1.25,1.40,1.60])
sum_all = np.array([sum435,sum606,sum814,sum105,sum125,sum140,sum160])
sum_can = np.array([sumc435,sumc606,sumc814,sumc105,sumc125,sumc140,sumc160])

fall = CubicSpline(W,sum_all)
newallx=np.arange(0.435,1.6,0.001)
newally=fall(newallx)

fcan = CubicSpline(W,sum_can)
newcanx=np.arange(0.435,1.6,0.001)
newcany=fcan(newcanx)

#----plot

plt.plot(newallx,newally)
plt.plot(newcanx,newcany)
plt.plot(W,sum_all,marker='o',color='r',linestyle='')
plt.plot(W,sum_can,marker='o',color='b',linestyle='')
plt.yscale("log")
plt.ylabel("Flux S$_v$ [erg s$^-$$^1$ cm$^-$$^2$ Hz$^-$$^1$]")
plt.xlabel("Wavelength [n$\lambda$]")
plt.show()

我从中获得的情节就是这样，插值中有明显的差距：

如果您想知道，这些是sum_all和sum_can数组中的值（我认为这无关紧要，但以防您希望数字自己绘制）：

sum_all:
[  3.87282732e+32   8.79993191e+32   1.74866333e+33   1.59946687e+33
   9.08556547e+33   6.70458731e+33   9.84832359e+33]
can_all:
[  2.98381061e+28   1.26194810e+28   3.30328780e+28   2.90254609e+29
   3.65117723e+29   3.46256846e+29   3.64483736e+29]

差距发生在[0.606,1.26194810e + 28]和[0.814,3.30328780e + 28]之间。如果我将间隔从0.001更改为更高的值，很明显，图表实际上并未中断，而是仅在y轴上低于0（但图表是连续的）。那为什么这样做呢？当然这不是一个正确的插值？只是用我们的眼睛看，这显然不是这两点之间的良好内插联系。

任何提示或评论都将非常感激。非常感谢你！

Answer 1

在线性范围内可以更好地观察到击穿的原因。

我们看到样条实际上传递到0以下，这在对数刻度上是未定义的。

所以我建议首先取数据的对数，对对数缩放的数据进行样条插值，然后按第10次幂缩放。

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

W = np.array([0.435,0.606,0.814,1.05,1.25,1.40,1.60])
sum_all = np.array([  3.87282732e+32,   8.79993191e+32,   1.74866333e+33,   1.59946687e+33,
   9.08556547e+33,   6.70458731e+33,   9.84832359e+33])
sum_can = np.array([  2.98381061e+28,   1.26194810e+28,   3.30328780e+28,   2.90254609e+29,
   3.65117723e+29,   3.46256846e+29,   3.64483736e+29])

fall = CubicSpline(W,np.log10(sum_all))
newallx=np.arange(0.435,1.6,0.001)
newally=fall(newallx)

fcan = CubicSpline(W,np.log10(sum_can))
newcanx=np.arange(0.435,1.6,0.01)
newcany=fcan(newcanx)


plt.plot(newallx,10**newally)
plt.plot(newcanx,10**newcany)
plt.plot(W,sum_all,marker='o',color='r',linestyle='')
plt.plot(W,sum_can,marker='o',color='b',linestyle='')
plt.yscale("log")

plt.ylabel("Flux S$_v$ [erg s$^-$$^1$ cm$^-$$^2$ Hz$^-$$^1$]")
plt.xlabel("Wavelength [n$\lambda$]")
plt.show()