Python NumPy로 1D 계수 배열을 활용해 점 (x, y)에서 2D 에르미트 급수 평가하기

Python NumPy에서 hermite.hermval2d() 메서드를 사용하면 점 (x, y)에서 2D 에르미트(Hermite) 급수를 손쉽게 평가할 수 있습니다. 이 메서드는 x와 y의 대응하는 값들이 짝을 이루어 만든 점들에서 이차원 다항식의 값을 반환합니다.

hermval2d() 메서드의 매개변수 이해하기

첫 번째 매개변수: x, y — 이차원 급수는 점 (x, y)에서 평가되며, x와 y는 반드시 동일한 형태(shape)를 가져야 합니다. x 또는 y가 리스트나 튜플이면 먼저 ndarray로 변환되며, 이미 ndarray라면 그대로 사용됩니다. 만약 ndarray가 아니라면 스칼라 값으로 취급됩니다.

두 번째 매개변수: C(계수 배열) — 다중 차수(multi-degree) i, j에 해당하는 항의 계수가 c[i, j]에 저장되도록 정렬된 배열입니다. c의 차원이 2보다 큰 경우, 남은 인덱스들은 여러 개의 계수 집합을 나타내는 데 사용됩니다.

단계별 구현 방법

먼저 필요한 라이브러리를 임포트합니다.

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

1차원 계수 배열을 생성합니다.

c = np.array([3, 5])

생성한 배열을 화면에 출력합니다.

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

배열의 차원을 확인합니다.

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

배열의 데이터 타입을 확인합니다.

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

배열의 형태(shape)를 확인합니다.

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

이제 hermite.hermval2d() 메서드를 사용하여 점 (x, y)에서 2D 에르미트 급수를 평가합니다.

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

전체 예제 코드

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

# 1차원 계수 배열 생성
c = np.array([3, 5])

# 배열 출력
print("Our Array...\n",c)

# 차원 확인
print("\nDimensions of our Array...\n",c.ndim)

# 데이터 타입 확인
print("\nDatatype of our Array object...\n",c.dtype)

# 형태(shape) 확인
print("\nShape of our Array object...\n",c.shape)

# hermite.hermval2d() 메서드로 점 (x, y)에서 2D 에르미트 급수 평가
print("\nResult...\n",H.hermval2d([1,2],[1,2],c))

실행 결과

Our Array...
   [3 5]

Dimensions of our Array...
1

Datatype of our Array object...
int64

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

Result...
   [ 59. 105.]

실행 결과를 보면 계수 배열 [3, 5]를 사용했을 때, 점 (1, 1)에서는 59, 점 (2, 2)에서는 105라는 값이 계산됩니다. 이처럼 hermval2d() 메서드를 활용하면 복잡한 수학적 연산 없이도 다양한 좌표에서 에르미트 다항식의 값을 빠르게 구할 수 있습니다.