Pythonでラゲール級数を微分する方法：lagder()による導関数の計算とスカラー倍

Pythonでラゲール（Laguerre）級数を微分するには、numpy.polynomial.laguerreモジュールのlaguerre.lagder()メソッドを使用します。このメソッドは、指定した軸に沿って係数配列cをm回微分した結果のラゲール級数係数を返します。さらに、各微分のたびに結果へスカラーsclが乗算され、最終的にはscl**m倍された値が得られます。

引数cは、各軸に沿って次数の低いものから高いものへと並べられた係数の配列です。例えば、[1,2,3]は 1*L_0 + 2*L_1 + 3*L_2 という級数を表します。また、[[1,2],[1,2]]は axis=0 をx、axis=1 をyとした場合、1*L_0(x)*L_0(y) + 1*L_1(x)*L_0(y) + 2*L_0(x)*L_1(y) + 2*L_1(x)*L_1(y) を意味します。

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

• c（第1引数）： ラゲール級数の係数配列。多次元配列の場合、異なる軸が異なる変数に対応し、各軸の次数は対応するインデックスによって決まります。
• m（第2引数）： 微分の回数。非負の整数である必要があります（デフォルト：1）。
• scl（第3引数）： 各微分の結果に乗算されるスカラー値。最終的な結果は scl**m 倍になります。これは変数の線形変換を行う際に便利です（デフォルト：1）。
• axis（第4引数）： 微分を適用する軸（デフォルト：0）。

実行手順

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

import numpy as np
from numpy.polynomial import laguerre 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)

ラゲール級数を微分するために、laguerre.lagder()メソッドを使用します。ここでは、2回微分し、各微分にスカラー「-1」を乗算しています。

print("\nResult...\n",L.lagder(c, 2, scl = -1))

コード例

import numpy as np
from numpy.polynomial import laguerre 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)

# lagder()メソッドでラゲール級数を微分（2回微分、スカラー-1を乗算）
print("\nResult...\n",L.lagder(c, 2, scl = -1))

出力結果

Our Array...
    [1 2 3 4]

Dimensions of our Array...
1

Datatype of our Array object...
int64

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

Result...
    [11. 4.]

このように、lagder()メソッドを使うことで、ラゲール級数の微分を簡単に実行できます。m回の微分とスカラー倍を組み合わせることで、変数の線形変換を伴う微分計算にも柔軟に対応できます。