Python NumPy를 활용한 점 x에서의 르장드르 급수 평가 방법

Python NumPy에서 특정 점 x에서 르장드르(Legendre) 급수를 평가하려면 polynomial.legendre.legval() 메서드를 사용하면 됩니다. 이 메서드는 르장드르 다항식 계수 배열과 평가하고자 하는 점을 입력받아 해당 지점에서의 급수 값을 계산해 줍니다.

legval() 메서드의 주요 매개변수

첫 번째 매개변수 x: 평가할 지점입니다. x가 리스트(list)나 튜플(tuple)이면 ndarray로 변환되며, 그렇지 않으면 변경되지 않은 채 스칼라로 취급됩니다. 어떤 경우든 x 또는 그 요소는 자기 자신뿐 아니라 c의 요소와도 덧셈과 곱셈 연산을 지원해야 합니다.

두 번째 매개변수 C: 계수 배열로, 차수 n에 해당하는 계수들이 c[n]에 저장되어 있는 형태로 정렬됩니다. c가 다차원 배열인 경우 나머지 인덱스들은 여러 개의 다항식을 나타냅니다. 2차원 배열의 경우 계수들은 c의 열(column)에 저장되어 있다고 생각할 수 있습니다.

세 번째 매개변수 tensor: True로 설정하면 계수 배열의 형태(shape)가 오른쪽에 1이 추가되어 확장되는데, 이는 x의 각 차원마다 하나씩 추가되는 것입니다(스칼라의 경우 차원은 0으로 간주). 그 결과 c의 모든 계수 열이 x의 모든 요소에 대해 평가됩니다. False로 설정하면 x가 c의 열 위로 브로드캐스트(broadcast)되어 평가됩니다. 이 키워드는 c가 다차원일 때 유용하며, 기본값은 True입니다.

단계별 구현 방법

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

import numpy as np
from numpy.polynomial import legendre as L

계수 배열을 생성합니다.

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)

점 x에서 르장드르 급수를 평가하기 위해 polynomial.legendre.legval() 메서드를 사용합니다.

print("\nResult...\n",L.legval(1,c))

전체 예제 코드

import numpy as np
from numpy.polynomial import legendre as L

# 계수 배열 생성
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)

# legval() 메서드로 점 x=1에서 르장드르 급수 평가
print("\nResult...\n",L.legval(1,c))

실행 결과

Our Array...
    [1 2 3]

Dimensions of our Array...
1

Datatype of our Array object...
int64

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

Result...
6.0

위 예제에서 계수 [1, 2, 3]에 해당하는 르장드르 급수를 x = 1 지점에서 평가한 결과 값은 6.0입니다. 이처럼 legval() 메서드를 사용하면 복잡한 수학적 계산 없이도 르장드르 급수를 손쉽게 평가할 수 있습니다.