Pythonでチェビシェフ級数を積分する方法【chebint()徹底解説】

Pythonでチェビシェフ級数を積分するには、NumPyの numpy.polynomial.chebyshev モジュールが提供する chebint() メソッドを使用します。このメソッドは、チェビシェフ級数の係数配列 c を、下限 lbnd から指定した軸に沿って m 回積分した結果の係数を返します。各反復処理のたびに、結果の級数は scl 倍され、積分定数 k が加算されます。

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

第1パラメータ「c」： チェビシェフ級数の係数を格納した配列です。c が多次元の場合、異なる軸はそれぞれ異なる変数に対応し、各軸の次数は対応するインデックスによって決まります。

第2パラメータ「m」： 積分の階数（積分を繰り返す回数）です。正の整数である必要があります。（デフォルト：1）

第3パラメータ「k」： 積分定数です。リストの1番目の値が1回目の積分におけるゼロでの値、2番目の値が2回目の積分におけるゼロでの値、というように対応します。k == []（デフォルト）の場合、すべての積分定数はゼロに設定されます。m == 1 の場合は、リストの代わりに単一のスカラー値を指定することも可能です。

第4パラメータ「lbnd」： 積分の下限です。（デフォルト：0）

第5パラメータ「scl」： 各積分の後、積分定数を加算する前に結果に乗算されるスカラー値です。（デフォルト：1）

第6パラメータ「axis」： 積分を実行する対象の軸です。（デフォルト：0）

実行手順

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

import numpy as np
from numpy.polynomial import chebyshev as C

次に、チェビシェフ級数の係数配列を作成します。

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

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

print("Our coefficient 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)

chebint() メソッドを使ってチェビシェフ級数を積分します。

print("\nResult...\n",C.chebint(c))

サンプルコード

import numpy as np
from numpy.polynomial import chebyshev as C

# チェビシェフ級数の係数配列を作成
c = np.array([1,2,3])

# 係数配列を表示
print("Our coefficient 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)

# chebint()メソッドでチェビシェフ級数を積分
print("\nResult...\n",C.chebint(c))

実行結果

Our coefficient Array...
[1 2 3]

Dimensions of our Array...
1

Datatype of our Array object...
int64

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

Result...
[ 0.5 -0.5 0.5 0.5]

このように、chebint() メソッドを使えば、引数を適切に指定するだけでチェビシェフ級数の積分を簡単に実行できます。積分回数や積分定数、下限などを変更することで、さまざまな条件に応じた計算にも柔軟に対応可能です。