Python:SciPy.interpolate 分段多项式

Question

import numpy as np
from scipy.interpolate import PiecewisePolynomial

xi = np.array([1,10])
yi = np.array([10,1])

p = PiecewisePolynomial(xi,yi)

不会产生两点的线性插值，但是

ZeroDivisionError: integer division or modulo by zero

哪里出了问题？

Answer 1

将 yi 替换为

yi = np.array([[10], [1]])

分段多项式要求y数组是类数组或数组列表结构。 y 的每个元素都可以是 x 及其后续导数的函数值。上面对 y 的更改创建了正确的线性插值

p = PiecewisePolynomial(xi,yi)
p.__call__([5.])
>> array([6.])
p.__call__([2.])
>> array([9.])