scipy.interpolate.UnivariateSpline.
roots#
- UnivariateSpline.roots()[source]#
返回样条的零点。
另请参阅
注释
限制:FITPACK 只支持三次样条。对于非三次样条,请使用 PPoly.root(有关示例,请参见下文)。
示例
对于某些数据,此方法可能会错过根。当其中一个样条节点(FITPACK 自动放置)恰好与真实根重合时,就会发生这种情况。解决方法是转换为
PPoly
,它使用不同的根查找算法。例如,
>>> x = [1.96, 1.97, 1.98, 1.99, 2.00, 2.01, 2.02, 2.03, 2.04, 2.05] >>> y = [-6.365470e-03, -4.790580e-03, -3.204320e-03, -1.607270e-03, ... 4.440892e-16, 1.616930e-03, 3.243000e-03, 4.877670e-03, ... 6.520430e-03, 8.170770e-03] >>> from scipy.interpolate import UnivariateSpline >>> spl = UnivariateSpline(x, y, s=0) >>> spl.roots() array([], dtype=float64)
转换为 PPoly 对象确实找到了 x=2 处的根
>>> from scipy.interpolate import splrep, PPoly >>> tck = splrep(x, y, s=0) >>> ppoly = PPoly.from_spline(tck) >>> ppoly.roots(extrapolate=False) array([2.])