Python에서 Hermite_e(확률론자 헤르미트) 급수를 적분하려면 numpy.polynomial.hermite_e 모듈의 hermeint() 메서드를 사용합니다. 이 메서드는 주어진 계수 배열에 대해 Hermite_e 급수의 부정적분을 수행하며, 다양한 옵션을 통해 적분 차수, 상수, 하한 등을 세밀하게 제어할 수 있습니다.
hermeint() 메서드의 매개변수
1. c (계수 배열)
Hermite_e 급수의 계수들을 담고 있는 배열입니다. 만약 c가 다차원 배열이라면, 각 축(axis)은 서로 다른 변수에 대응하며 각 축의 차수(degree)는 해당 인덱스 값으로 결정됩니다.
2. m (적분 차수)
적분을 수행할 횟수를 의미하며, 반드시 양수여야 합니다. 기본값은 1입니다.
3. k (적분 상수)
적분 상수를 지정합니다. 하한(lbnd)에서의 첫 번째 적분값이 리스트의 첫 번째 값이 되고, 두 번째 적분값이 두 번째 값이 되는 방식으로 순서대로 대응됩니다. k == [] (기본값)이면 모든 적분 상수는 0으로 설정됩니다. m == 1인 경우에는 리스트 대신 단일 스칼라 값을 전달할 수도 있습니다.
4. lbnd (적분 하한)
적분의 하한(lower bound)을 지정합니다. 기본값은 0입니다.
5. scl (스케일 팩터)
스칼라 값으로, 각 적분이 수행된 후 적분 상수가 더해지기 전에 그 결과에 곱해집니다. 기본값은 1입니다.
6. axis (축)
적분이 수행될 축을 지정합니다. 기본값은 0입니다.
구현 단계
먼저 필요한 라이브러리를 임포트합니다:
import numpy as np from numpy.polynomial import hermite_e as H
계수 배열을 생성합니다:
c = np.array([1,2,3])
배열을 출력합니다:
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)hermeint() 메서드를 호출하여 Hermite_e 급수를 적분합니다:
print("\nResult...\n",H.hermeint(c))전체 예제 코드
import numpy as np
from numpy.polynomial import hermite_e as H
# 계수 배열 생성
c = np.array([1,2,3])
# 배열 출력
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)
# hermeint() 메서드로 Hermite_e 급수 적분
print("\nResult...\n",H.hermeint(c))실행 결과
Our Array...
[1 2 3]
Dimensions of our Array...
1
Datatype of our Array object...
int64
Shape of our Array object...
(3,)
Result...
[1. 1. 1. 1.]위 실행 결과를 보면, 입력 배열 [1, 2, 3]에 대해 적분을 수행한 결과 [1. 1. 1. 1.]이 반환되었습니다. 적분 후 급수의 차수가 하나 증가하여 계수 배열의 길이도 하나 늘어난 것을 확인할 수 있습니다. 또한 정수형 입력에 대해서도 적분 결과는 실수(float)로 반환됩니다.