Python NumPyでラゲール級数を加算する方法：lagadd()メソッドの使い方

PythonのNumPyライブラリには、多項式演算を簡単に行える便利なモジュールが数多く用意されています。その中でもnumpy.polynomial.laguerreモジュールは、ラゲール多項式（Laguerre polynomial）に関連する各種操作を提供しています。本記事では、2つのラゲール級数を加算するためのpolynomial.laguerre.lagadd()メソッドの使い方を、具体的なコード例とともに解説します。

lagadd()メソッドとは

lagadd()メソッドは、引数として渡された2つのラゲール級数c1 + c2の和を表す配列を返します。係数は低次の項から高次の項へと順に並べて指定する必要があります。たとえば、[1, 2, 3]という配列は、P_0 + 2*P_1 + 3*P_2 という級数を意味します。

• c1：ラゲール級数の係数を格納した1次元配列（低次→高次の順）
• c2：ラゲール級数の係数を格納した1次元配列（低次→高次の順）

実装手順

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

まず、NumPyとラゲール多項式モジュールをインポートします。

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

ステップ2：係数の1次元配列を作成する

続いて、ラゲール級数の係数を格納した1次元配列を2つ作成します。

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

ステップ3：配列の内容・データ型・次元・形状を確認する

作成した配列が正しく定義されているか、内容や属性を表示して確認しましょう。

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)

ステップ4：lagadd()メソッドで級数を加算する

最後に、polynomial.laguerre.lagadd()メソッドを使って2つのラゲール級数を加算します。このメソッドは、両者の和を表すラゲール級数の配列を返します。

print("\nResult (sum)....\n",L.lagadd(c1, c2))

完全なコード例

ここまでの手順をまとめたサンプルコード全体は以下の通りです。

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

# ラゲール級数の係数を持つ1次元配列を作成
c1 = np.array([1,2,3])
c2 = np.array([3,2,1])

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

# lagadd()メソッドで2つのラゲール級数を加算
# 戻り値は和を表すラゲール級数の配列
print("\nResult (sum)....\n",L.lagadd(c1, c2))

実行結果

上記のコードを実行すると、以下のような出力が得られます。

Array1...
[1 2 3]

Array2...
[3 2 1]

Array1 datatype...
int64

Array2 datatype...
int64

Dimensions of Array1...
1

Dimensions of Array2...
1

Shape of Array1...
(3,)

Shape of Array2...
(3,)

Result (sum)....
[4. 4. 4.]

出力を見ると、対応する係数同士が足し合わされ、[4. 4. 4.]という結果が返されていることがわかります。このように、lagadd()メソッドを使えば、ラゲール級数の加算をわずか1行で実現できます。減算を行いたい場合は、同様にlagsub()メソッドを使用することも可能です。