scipy - scipy.interpolate.splrep 为 b 样条提供不正确的插值

Question

我正在使用以下代码：

x = [80, 85, 90, 95, 100, 105, 110, 115]
y = [0.31976258043267208, 0.31718670335264199, 0.30116311412953689, 0.29687766899837931, 0.29802720516866577, 0.30707962224836099, 0.32328422706655086, 0.34467623793716823]

bSpline = scipy.interpolate.splrep(x, y, k=3)
yInterp = scipy.interpolate.splev(np.linspace(0,1,num=8), bSpline)

但是 bSpline 有奇怪的结（有些是 0），然后 yInterp 都搞砸了：

 bSpline: (array([  80,   80,   80,   80,   90,   95,  100,  105,  115, 115,  115,  115 ]), array([ 0.31976258,  0.33015199,  0.29797999,  0.29680536,  0.29635823,
    0.30911387,  0.32907255,  0.34467624,  0,  0, 0,  0 ]), 3)
 yInterp: [-29.46356522 -29.3160613  -29.16904367 -29.02251153 -28.87646408, -28.73090053 -28.58582008 -28.44122191]

yInterp 应该有一些值在正 0.30 左右，而不是 -28 和 -29

Answer 1

您传递的 x 值splrep范围从 80 到 115。您正在评估从 0 到 1 的值的样条曲线。如果您在输入范围内的值进行评估，输出将在 0.3 左右。