numpy.polynomial.legendre.legvander3d#

polynomial.legendre.legvander3d(x, y, z, deg)[source]#

給定次數的偽范德蒙矩陣。

返回給定次數 deg 和樣本點 (x, y, z) 的偽范德蒙矩陣。如果 lmn 是在 xyz 中給定的次數,則偽范德蒙矩陣定義為

\[V[..., (m+1)(n+1)i + (n+1)j + k] = L_i(x)*L_j(y)*L_k(z),\]

其中 0 <= i <= l0 <= j <= m0 <= j <= nV 的前導索引指向點 (x, y, z),而最後一個索引編碼勒讓德多項式的次數。

如果 V = legvander3d(x, y, z, [xdeg, ydeg, zdeg]),則 V 的列對應於形狀為 (xdeg + 1, ydeg + 1, zdeg + 1) 的 3-D 係數陣列 c 中的元素,順序為

\[c_{000}, c_{001}, c_{002},... , c_{010}, c_{011}, c_{012},...\]

np.dot(V, c.flat)legval3d(x, y, z, c) 在捨入誤差範圍內將會相同。這種等效性對於最小平方擬合以及評估大量相同次數和樣本點的 3-D 勒讓德級數都很有用。

參數:
x, y, z類陣列

點座標的陣列,所有陣列的形狀都相同。dtype 將根據是否有任何元素是複數而轉換為 float64 或 complex128。純量會轉換為 1-D 陣列。

deg整數列表

最大次數的列表,形式為 [x_deg, y_deg, z_deg]。

返回:
vander3dndarray

返回矩陣的形狀為 x.shape + (order,),其中 \(order = (deg[0]+1)*(deg[1]+1)*(deg[2]+1)\)。dtype 將與轉換後的 xyz 相同。