Pythonでルジャンドル級数を微分する方法

Pythonでルジャンドル(Legendre)級数を微分するには、numpy.polynomial.legendreモジュールのlegder()メソッドを使用します。このメソッドは、指定された軸に沿ってm回微分した後のルジャンドル級数の係数cを返します。各反復処理のたびに、結果にはスカラー値sclが乗算されます。

legder()メソッドのパラメータ

第1パラメータ:c — ルジャンドル級数の係数を格納した配列です。cが多次元配列の場合、異なる軸はそれぞれ異なる変数に対応し、各軸における次数は対応するインデックスによって決まります。

第2パラメータ:m — 微分を行う回数です。非負の整数である必要があります(デフォルト:1)。

第3パラメータ:scl — スカラー値です。各微分のたびにこの値が乗算され、最終的な結果はscl**m倍になります。これは主に線形の変数変換を行う際に使用されます(デフォルト:1)。

第4パラメータ:axis — 微分を実行する対象となる軸です(デフォルト:0)。

実行手順

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

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

次に、係数の配列を作成します。

c = np.array([1,2,3,4])

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

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

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

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

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

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

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

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

最後に、legder()メソッドを呼び出してルジャンドル級数を微分します。

print("\nResult...\n",L.legder(c))

完全なコード例

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

# 係数の配列を作成
c = np.array([1,2,3,4])

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

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

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

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

# legder()メソッドでルジャンドル級数を微分
print("\nResult...\n",L.legder(c))

実行結果

Our Array...
   [1 2 3 4]

Dimensions of our Array...
1

Datatype of our Array object...
int64

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

Result...
   [ 6. 9. 20.]

このように、legder()メソッドを使うことで、ルジャンドル級数の微分係数を簡単に計算できます。出力結果から、元の4つの係数[1, 2, 3, 4]に対して1回微分を行うと、3つの係数[6., 9., 20.]が得られることがわかります。