Python NumPyのchebmul()メソッドでチェビシェフ級数同士を乗算する方法

Pythonで1つのチェビシェフ級数を別のチェビシェフ級数と乗算するには、numpy.polynomial.chebyshev.chebmul() メソッドを使用します。このメソッドは、2つのチェビシェフ級数の積を表す係数の配列を返します。

引数として渡すのは、低次の項から高次の項へと順に並べた係数のシーケンスです。例えば、[1, 2, 3] は T_0 + 2*T_1 + 3*T_2 という級数を表します。パラメータ c1 と c2 は、低次から高次の順に並べられたチェビシェフ級数の係数を持つ1次元配列です。

実行手順

1. 必要なライブラリをインポートする

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

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

2. チェビシェフ級数の係数配列を作成する

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

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

3. 配列の内容を表示する

print("Array1...\n",c1)
print("\nArray2...\n",c2)

4. データ型を確認する

print("\nArray1 datatype...\n",c1.dtype)
print("\nArray2 datatype...\n",c2.dtype)

5. 配列の次元を確認する

print("\nDimensions of Array1...\n",c1.ndim)
print("\nDimensions of Array2...\n",c2.ndim)

6. 配列の形状を確認する

print("\nShape of Array1...\n",c1.shape)
print("\nShape of Array2...\n",c2.shape)

7. chebmul()メソッドで乗算を実行する

最後に、numpy.polynomial.chebyshev.chebmul() メソッドを使って、2つのチェビシェフ級数を乗算します。

print("\nResult (multiplication)....\n",C.chebmul(c1,c2))

完全なコード例

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

# チェビシェフ級数の係数を持つ1次元配列を作成
c1 = np.array([1,2,3])
c2 = np.array([3,2,1])

# 係数配列を表示
print("Array1...\n",c1)
print("\nArray2...\n",c2)

# データ型を表示
print("\nArray1 datatype...\n",c1.dtype)
print("\nArray2 datatype...\n",c2.dtype)

# 両方の配列の次元を確認
print("\nDimensions of Array1...\n",c1.ndim)
print("\nDimensions of Array2...\n",c2.ndim)

# 両方の配列の形状を確認
print("\nShape of Array1...\n",c1.shape)
print("\nShape of Array2...\n",c2.shape)

# chebmul()メソッドでチェビシェフ級数を乗算
print("\nResult (multiplication)....\n",C.chebmul(c1,c2))

出力結果

Array1...
[1 2 3]

Array2...
[3 2 1]

Array1 datatype...
int64

Array2 datatype...
int64

Dimensions of Array1...
1

Dimensions of Array2...
1

Shape of Array1...
(3,)

Shape of Array2...
(3,)

Result (multiplication)....
[ 6.5 12. 12. 4. 1.5]

このように、chebmul() メソッドを使用することで、2つのチェビシェフ級数の積を簡単に計算できます。結果は新しいチェビシェフ級数の係数配列として返され、元の級数よりも次数が高くなる点に注意してください。