插值在python中没有负值

Question

我一直试图从这些值创建一条平滑线,但我的结果中不能有负值.到目前为止,我尝试的所有方法都给出了负值.会爱一些帮助.

import matplotlib.pyplot as plt
from scipy.interpolate import UnivariateSpline
import numpy as np
y = np.asarray([0,5,80,10,1,10,40,30,80,5,0])
x = np.arange(len(y))

plt.plot(x, y, 'r', ms=5)
spl = UnivariateSpline(x, y)
xs = np.linspace(0,len(y)-1, 1000)
spl.set_smoothing_factor(2)

plt.plot(xs, spl(xs), 'g', lw=3)
plt.show()

Answer 1

已知样条拟合过冲.你好像在寻找一种所谓的单调插值器.例如,

In [10]: from scipy.interpolate import pchip

In [11]: pch = pchip(x, y)

产生

In [12]: xx = np.linspace(x[0], x[-1], 101)

In [13]: plt.plot(x, y, 'ro', label='points')
Out[13]: [<matplotlib.lines.Line2D at 0x7fce0a7fe390>]

In [14]: plt.plot(xx, pch(xx), 'g-', label='pchip')
Out[14]: [<matplotlib.lines.Line2D at 0x7fce0a834b10>]