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

PythonのNumPyライブラリでは、numpy.polynomialモジュールのlegvander()メソッドを使うことで、ルジャンドル多項式に基づく擬似ファンデルモンド行列を簡単に生成できます。ファンデルモンド行列は、最小二乗法による多項式近似や補間など、数値計算のさまざまな場面で活用される重要な行列です。

legvander()メソッドの概要

このメソッドは、指定した点における各次数のルジャンドル多項式の値を列として並べた、擬似ファンデルモンド行列を返します。返される行列の形状は x.shape + (deg + 1,) となり、最後のインデックスが対応するルジャンドル多項式の次数を表します。データ型（dtype）は、変換後のxと同じものが使用されます。

主なパラメータ

• x: 評価点の配列を指定します。要素に複素数が含まれる場合は complex128、含まれない場合は float64 に変換されます。スカラー値を渡した場合は、自動的に1次元配列へ変換されます。
• deg: 結果として得られる行列の最大次数を指定します。

実行手順

ステップ1：必要なライブラリをインポートする

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

ステップ2：配列を作成する

x = np.array([0, 1, -1, 2])

ステップ3：配列の内容と属性を確認する

# 配列を表示
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)

ステップ4：legvander()メソッドで行列を生成する

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

完全なコード例

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

# 配列を作成
x = np.array([0, 1, -1, 2])

# 配列を表示
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 1 -1 2]

Dimensions of our Array...
1

Datatype of our Array object...
int64

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

Result...
    [[ 1. 0. -0.5]
    [ 1. 1. 1. ]
    [ 1. -1. 1. ]
    [ 1. 2. 5.5]]

解説

出力結果を見ると、入力配列 [0, 1, -1, 2] の各要素に対して、0次・1次・2次のルジャンドル多項式の値が計算されていることがわかります。各行は1つの評価点に対応し、左から順に P₀(x)=1、P₁(x)=x、P₂(x)=(3x²−1)/2 の値が格納されています。

たとえば x=2 の場合、P₁(2)=2、P₂(2)=(3×4−1)/2=5.5 となるため、最終行は [1. 2. 5.5] となります。このようにして得られた行列は、ルジャンドル多項式によるデータフィッティング（legfit()）や多項式評価（legval()）などの処理で利用できます。