Python NumPyのhermval3d()メソッドで点（x, y, z）における3次元エルミート級数を評価する方法

PythonのNumPyライブラリでは、点（x, y, z）における3次元エルミート級数を評価するために、numpy.polynomial.hermiteモジュールのhermval3d()メソッドを使用します。このメソッドは、x、y、zの対応する値から構成される点の組に対して、多次元多項式の値を返します。

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

第1パラメータ（x, y, z）： 3次元級数を評価する点（x, y, z）を指定します。x、y、zは同じ形状である必要があります。いずれかがリストやタプルの場合は、まずndarrayに変換されます。ndarrayでない場合はスカラーとして扱われます。

第2パラメータ（C）： 係数の配列です。多項式次数i, j, kの項に対応する係数はc[i, j, k]に格納されている必要があります。cの维度が3より大きい場合、残りのインデックスは複数の係数セットを表します。

実装手順

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

import numpy as np
from numpy.polynomial import hermite as H

次に、係数用の3次元配列を作成します。

c = np.arange(24).reshape(2,2,6)

作成した配列を表示して確認しましょう。

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)

最後に、hermite.hermval3d()メソッドを使って、点（x, y, z）における3次元エルミート級数を評価します。このメソッドは、x、y、zの対応する値から構成される点の組に対する多次元多項式の値を返します。

print("\nResult...\n",H.hermval3d([1,2],[1,2],[1,2],c))

完全なコード例

import numpy as np
from numpy.polynomial import hermite as H

# 係数用の3次元配列を作成
c = np.arange(24).reshape(2,2,6)

# 配列を表示
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)

# hermval3d()メソッドで点（x, y, z）における3次元エルミート級数を評価
print("\nResult...\n",H.hermval3d([1,2],[1,2],[1,2],c))

実行結果

Our Array...
    [[[ 0 1 2 3 4 5]
    [ 6 7 8 9 10 11]]

    [[12 13 14 15 16 17]
    [18 19 20 21 22 23]]]

Dimensions of our Array...
3

Datatype of our Array object...
int64

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

Result...
    [-4050. 52240.]

このように、hermval3d()メソッドを使うことで、指定した点（x, y, z）における3次元エルミート級数の値を簡単に計算できます。エルミート多項式は量子力学や確率論などの分野で広く利用されており、NumPyの多項式モジュールを活用することで効率的に数値計算を行うことが可能です。