Pythonで多次元係数を持つルジャンドル級数を微分する方法

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

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

• 第1パラメータ（c）：ルジャンドル級数の係数を格納した配列です。cが多次元の場合、各軸はそれぞれ異なる変数に対応し、各軸の次数は対応するインデックスによって決まります。
• 第2パラメータ（m）：微分を行う回数です。非負の整数である必要があります。（デフォルト：1）
• 第3パラメータ（scl）：スカラー値です。各微分の結果はscl倍され、最終的な結果はscl**m倍になります。これは変数の線形変換を行う際に役立ちます。（デフォルト：1）
• 第4パラメータ（axis）：微分を実行する対象となる軸です。（デフォルト：0）

実行手順

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

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

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

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

c = np.arange(4).reshape(2,2)

ステップ3：配列を表示する

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

ステップ4：次元数を確認する

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

ステップ5：データ型を確認する

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

ステップ6：配列の形状を確認する

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

ステップ7：legder()メソッドでルジャンドル級数を微分する

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

完全なコード例

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

# 多次元の係数配列を作成
c = np.arange(4).reshape(2,2)

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

Dimensions of our Array...
2

Datatype of our Array object...
int64

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

Result...
    [[2. 3.]]

この例では、形状(2, 2)の2次元配列に対してlegder()メソッドを適用しています。デフォルトではaxis=0（最初の軸）に沿って1回微分が行われ、その結果として形状(1, 2)の配列[[2. 3.]]が出力されています。