Python NumPy로 점 x 목록에서 라게르(Laguerre) 급수 평가하기

Python NumPy에서 특정 점 x에서 라게르(Laguerre) 급수를 평가하려면 polynomial.laguerre.lagval() 메서드를 사용합니다. 이 메서드는 주어진 계수 배열에 해당하는 라게르 급수를 지정한 점들에서 계산하여 결과를 반환합니다.

lagval() 메서드의 매개변수

1. x

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

2. c (계수 배열)

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

3. tensor

tensor가 True이면 계수 배열의 형태가 오른쪽에 1을 추가하여 확장됩니다. 확장 횟수는 x의 차원 수만큼이며, 이때 스칼라는 차원 0으로 간주됩니다. 그 결과 c의 모든 계수 열이 x의 모든 요소에 대해 평가됩니다. False이면 x가 c의 열에 걸쳐 브로드캐스트(broadcast)되어 평가됩니다. 이 키워드는 c가 다차원일 때 유용하며, 기본값은 True입니다.

단계별 구현 방법

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

import numpy as np
from numpy.polynomial import laguerre 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]

lagval() 메서드를 사용해 점 x에서 라게르 급수를 평가합니다.

print("\nResult...\n",L.lagval(x,c))

전체 예제 코드

import numpy as np
from numpy.polynomial import laguerre 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)

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

# lagval() 메서드로 라게르 급수 평가
print("\nResult...\n",L.lagval(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...
    [ 3.5 76. 223.5]

위 결과에서 볼 수 있듯이, 계수 [1, 2, 3]으로 정의된 라게르 급수가 점 [5, 10, 15]에서 각각 3.5, 76.0, 223.5로 평가되었습니다. lagval() 메서드는 입력이 리스트여도 자동으로 ndarray로 변환해 처리하므로 편리하게 사용할 수 있습니다.