Python NumPyでルジャンドル多項式の擬似ファンデルモンド行列を生成する方法

PythonのNumPyでは、numpy.polynomialモジュールに含まれる legvander() メソッドを使うことで、ルジャンドル多項式の擬似ファンデルモンド行列を簡単に生成できます。このメソッドは、指定した点の配列に対して各次数のルジャンドル多項式の値を計算し、それらを列として並べた行列を返します。

戻り値となる行列の形状は x.shape + (deg + 1,) であり、最後のインデックスが対応するルジャンドル多項式の次数を表します。データ型（dtype）は、入力配列 x を変換した後の型と同じになります。要素に複素数が含まれる場合は complex128、そうでない場合は float64 に変換されます。

パラメータについて

このメソッドには主に2つのパラメータがあります。

x： 点の座標を格納した配列です。スカラー値を渡した場合は自動的に1次元配列へ変換されます。また、dtypeは要素が複素数かどうかに応じて float64 または complex128 に変換されます。

deg： 生成される行列の最大次数を指定します。0から deg までの各次数のルジャンドル多項式が計算されます。

実行手順

まず、必要なライブラリをインポートします。

import numpy as np
from numpy.polynomial import legendre as L

点の配列を作成します。

x = np.array([0, 3.5, -1.4, 2.5])

作成した配列を表示します。

print("Our Array...\n",x)

配列の次元数を確認します。

print("\nDimensions of our Array...\n",x.ndim)

配列のデータ型を確認します。

print("\nDatatype of our Array object...\n",x.dtype)

配列の形状を確認します。

print("\nShape of our Array object...\n",x.shape)

legvander() メソッドを使って、ルジャンドル多項式の擬似ファンデルモンド行列を生成します。

print("\nResult...\n",L.legvander(x, 2))

サンプルコード全体

import numpy as np
from numpy.polynomial import legendre as L

# 配列を作成
x = np.array([0, 3.5, -1.4, 2.5])

# 配列を表示
print("Our Array...\n",x)

# 次元数を確認
print("\nDimensions of our Array...\n",x.ndim)

# データ型を確認
print("\nDatatype of our Array object...\n",x.dtype)

# 形状を確認
print("\nShape of our Array object...\n",x.shape)

# legvander() メソッドでルジャンドル多項式の擬似ファンデルモンド行列を生成
print("\nResult...\n",L.legvander(x, 2))

実行結果

Our Array...
    [ 0. 3.5 -1.4 2.5]

Dimensions of our Array...
1

Datatype of our Array object...
float64

Shape of our Array object...
(4,)

Result...
    [[ 1. 0. -0.5 ]
    [ 1. 3.5 17.875]
    [ 1. -1.4 2.44 ]
    [ 1. 2.5 8.875]]

このように、legvander(x, 2) を呼び出すことで、各点における0次・1次・2次のルジャンドル多項式の値が計算され、4×3の擬似ファンデルモンド行列が得られます。この手法は、ルジャンドル多項式による最小二乗近似や補間など、多項式フィッティングの基盤として広く活用されています。