PythonでChebyshev級数を積分する：chebint()メソッドによる積分次数の設定方法

PythonでChebyshev（チェビシェフ）級数を積分するには、numpy.polynomial.chebyshevモジュールに含まれるchebint()メソッドを使用します。このメソッドは、Chebyshev級数の係数配列cを指定した回数だけ積分した結果を返します。各反復処理のたびに、結果の級数はscl倍され、積分定数kが加算されます。

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

chebint()メソッドには以下の6つのパラメータがあります。

• c（第1引数）：Chebyshev級数の係数を格納した配列です。多次元配列の場合、異なる軸がそれぞれ異なる変数に対応し、各軸の次数は対応するインデックスによって決まります。
• m（第2引数）：積分の次数（何回積分するか）を指定します。正の整数である必要があります。デフォルト値は1です。
• k（第3引数）：積分定数を指定します。1回目の積分結果のゼロにおける値がリストの最初の値、2回目の積分結果のゼロにおける値が2番目の値、というように対応します。k == []（デフォルト）の場合、すべての定数はゼロに設定されます。m == 1の場合は、リストの代わりに単一のスカラー値を指定することもできます。
• lbnd（第4引数）：積分の下限です。デフォルト値は0です。
• scl（第5引数）：各積分の後、積分定数を加算する前に、結果に乗算されるスケール係数です。デフォルト値は1です。
• axis（第6引数）：積分を行う軸を指定します。デフォルト値は0です。

実装手順

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

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

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

ステップ2：Chebyshev級数の係数配列を作成する

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

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

print("Our coefficient 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：chebint()メソッドで積分を実行する

Chebyshev級数を積分するには、chebyshev.chebint()メソッドを使用します。ここでは、積分次数m = 3として3回積分しています。

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

完全なコード例

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

# Chebyshev級数の係数配列を作成
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()メソッドでChebyshev級数を3回積分
print("\nResult...\n",C.chebint(c, m = 3))

実行結果

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.03125 -0.1875 0.04166667 -0.05208333 0.01041667 0.00625 ]

まとめ

この記事では、NumPyのchebyshev.chebint()メソッドを使用してChebyshev級数を積分する方法を解説しました。積分次数mを指定することで、任意の回数だけ繰り返し積分が可能です。また、積分定数kや積分下限lbnd、スケール係数sclなどを組み合わせることで、柔軟な数値積分処理を実現できます。多項式補間や近似計算においてChebyshev級数を扱う際に、ぜひ活用してみてください。