Python SciPy interpolate.NdPPoly用法及代碼示例

本文簡要介紹 python 語言中 scipy.interpolate.NdPPoly 的用法。

用法:

class  scipy.interpolate.NdPPoly(c, x, extrapolate=None)#

分段張量積多項式

點處的值xp = (x', y', z', ...)通過首先計算區間索引來評估i這樣：

x[0][i[0]] <= x' < x[0][i[0]+1]
x[1][i[1]] <= y' < x[1][i[1]+1]
...

然後計算：

S = sum(c[k0-m0-1,...,kn-mn-1,i[0],...,i[n]]
        * (xp[0] - x[0][i[0]])**m0
        * ...
        * (xp[n] - x[n][i[n]])**mn
        for m0 in range(k[0]+1)
        ...
        for mn in range(k[n]+1))

其中k[j] 是多項式在維度 j 中的次數。這種表示是分段多元冪基礎。

參數 ：：

c： ndarray，形狀(k0，...，kn，m0，...，mn，...)

多項式係數，每個維度 j 具有多項式階 kj 和 mj+1 間隔。

x： ndim-tuple of ndarrays, 形狀 (mj+1,)

每個維度的多項式斷點。這些必須按升序排列。

extrapolate： 布爾型，可選

是否根據第一個和最後一個間隔推斷出界點，或者返回 NaN。默認值：真。

注意：

冪基中的高階多項式可能在數值上不穩定。

屬性 ：：

x： ndarray 的元組

斷點。

c： ndarray

多項式的係數。