PythonとNumPyを使ってルジャンドル級数同士の減算を行う方法

はじめに

Python NumPyでは、polynomial.legendre.legsub() メソッドを使うことで、あるルジャンドル（Legendre）級数から別のルジャンドル級数を減算することができます。このメソッドは、2つのルジャンドル級数の差を表す配列を返します。

具体的には、legsub() は2つのルジャンドル級数 c1 − c2 の差を計算します。引数として渡すのは係数のシーケンスで、低次の項から高次の項へと順番に並べます。たとえば [1, 2, 3] は P₀ + 2·P₁ + 3·P₂ という級数を表します。パラメータ c1 および c2 は、係数が低次から高次の順に並んだルジャンドル級数の1次元配列です。

手順

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

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

次に、ルジャンドル級数の係数を持つ1次元配列を作成します。

c1 = np.array([2,3,4])
c2 = np.array([4,3,2])

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

print("Array1...\n",c1)
print("\nArray2...\n",c2)

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

print("\nArray1 datatype...\n",c1.dtype)
print("\nArray2 datatype...\n",c2.dtype)

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

print("\nDimensions of Array1...\n",c1.ndim)
print("\nDimensions of Array2...\n",c2.ndim)

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

print("\nShape of Array1...\n",c1.shape)
print("\nShape of Array2...\n",c2.shape)

最後に、あるルジャンドル級数から別の級数を減算するために、polynomial.legendre.legsub() メソッドを実行します。

print("\nResult (difference)....\n",L.legsub(c1, c2))

完全なコード例

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

# ルジャンドル級数の係数を持つ1次元配列を作成
c1 = np.array([2,3,4])
c2 = np.array([4,3,2])

# 係数の配列を表示
print("Array1...\n",c1)
print("\nArray2...\n",c2)

# データ型を表示
print("\nArray1 datatype...\n",c1.dtype)
print("\nArray2 datatype...\n",c2.dtype)

# 次元数を確認
print("\nDimensions of Array1...\n",c1.ndim)
print("\nDimensions of Array2...\n",c2.ndim)

# 形状を確認
print("\nShape of Array1...\n",c1.shape)
print("\nShape of Array2...\n",c2.shape)

# legsub()メソッドで2つのルジャンドル級数の差を計算
print("\nResult (difference)....\n",L.legsub(c1, c2))

出力結果

Array1...
   [2 3 4]

Array2...
   [4 3 2]

Array1 datatype...
int64

Array2 datatype...
int64

Dimensions of Array1...
1

Dimensions of Array2...
1

Shape of Array1...
(3,)

Shape of Array2...
(3,)

Result (difference)....
   [-2. 0. 2.]

結果の解説

出力結果 [-2. 0. 2.] は、2つの級数の対応する係数同士の差を表しています。つまり、c1 の各係数から c2 の対応する係数を引いた値が、新しいルジャンドル級数の係数になっています。このように legsub() を使えば、多項式そのものを展開して手動で計算する必要なく、ルジャンドル級数の減算を簡単かつ効率的に行うことができます。