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

Python NumPy에서 특정 지점 x에서 르장드르(Legendre) 급수를 평가하려면 polynomial.legendre.legval() 메서드를 사용합니다. 이 메서드는 주어진 계수 배열에 해당하는 르장드르 급수를 입력된 지점들에서 계산해 줍니다.

legval() 메서드의 매개변수

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

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

세 번째 매개변수 tensor: True이면 계수 배열의 형태(shape) 오른쪽에 x의 각 차원마다 하나씩 1을 추가하여 확장합니다(스칼라의 경우 차원은 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는 리스트입니다.

x = [5, 10, 15]

이제 legval() 메서드를 사용하여 지점 x에서 르장드르 급수를 평가합니다.

print("\nResult...\n",L.legval(x,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)

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

# x는 리스트
x = [5, 10, 15]

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

실행 결과

Our Array...
    [1 2 3]

Dimensions of our Array...
1

Datatype of our Array object...
int64

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

Result...
    [ 122. 469.5 1042. ]

위 결과에서 볼 수 있듯이, 계수 [1, 2, 3]으로 정의된 르장드르 급수가 지점 5, 10, 15에서 각각 122, 469.5, 1042로 평가되었습니다. 이처럼 legval() 메서드를 활용하면 리스트 형태의 여러 지점에서도 손쉽게 르장드르 급수 값을 계산할 수 있습니다.