[논문 리뷰] An application of Bell polynomials in numerical solving of nonlinear differential equations

Josef Rebenda|arXiv (Cornell University)|2019. 01. 27.

Iterative Methods for Nonlinear Equations인용 수 4

한 줄 요약

이 논문은 비선형 미분방정식의 수치적 해법을 단순화하기 위해 부분 일반 벨 다항식의 새로운 응용을 제안한다. 파아 디 브루노의 공식을 벨 다항식을 사용해 재구성함으로써, 기호 도함수 없이도 비선형 항을 효율적으로 계산할 수 있게 되었으며, 이는 재귀적 산술 연산을 통해 정확한 급수 해를 도출한다. 이는 정확한 해가 패턴 인식을 통해 확인된 두 개의 초기값 문제에서 입증되었다.

ABSTRACT

Partial ordinary Bell polynomials are used to formulate and prove a version of the Fa\\`{a} di Bruno's formula which is convenient for handling nonlinear terms in the differential transformation. Applicability of the result is shown in two examples of solving the initial value problem for differential equations which are nonlinear with respect to the dependent variable.

연구 동기 및 목표

• 비선형 미분방정식의 해법에서 미지 해의 복합 함수를 포함하는 비선형 항을 처리하는 데 발생하는 계산 복잡도를 해결하기 위해.
• 벨 다항식의 조합적 성질을 활용하여 비선형 항을 변환하는 체계적이고 계산 효율적인 접근법을 개발하기 위해.
• 기호 도함수 또는 복잡한 적분 계산이 필요한 기존의 반해석적 방법(예: 아도미안 분해법)에 대한 실용적인 대안을 제공하기 위해.
• 제안된 방법이 재귀적 산술 연산을 통해 정확한 급수 해를 도출할 수 있음을 보여주며, 닫힌 형태의 해 식별 가능성을 제시하기 위해.

제안 방법

• 논문은 복합 함수의 고차 도함수를 표현하기 위해 부분 일반 벨 다항식을 사용한 파아 디 브루노의 공식의 수정 버전을 도입한다.
• 함수 $ f(u(t)) $ 의 미분 변환은 기호 도함수 $ f $ 가 필요 없이 변환된 계수 $ U(k) $ 를 포함하는 다항식의 합으로 표현된다.
• 이 방법은 이전에 계산된 $ U(i) $ ( $ i < k $ ) 를 사용하여 $ U(k) $ 를 재귀적으로 계산함으로써 계산 효율성을 확보한다.
• 알고리즘은 역 미분 변환을 사용하여 계수들이 벨 다항식 전개로부터 유도된 테일러 유사 급수 형태로 해를 재구성한다.
• 이 방법은 수치적 적분, 미분 또는 반복 해법을 모두 피하고, 오직 기본 산술 연산에 의존한다.
• 이 방법은 두 개의 초기값 문제에 대해 검증되었으며, 식별 가능한 계수 패턴을 통해 정확한 해로 수렴하는 것으로 나타났다.

실험 결과

연구 질문

• RQ1 부분 벨 다항식을 사용하여 기호 미분 없이도 비선형 상미분방정식에서 비선형 항의 미분 변환을 간소화할 수 있는가?
• RQ2 벨 다항식을 사용해 파아 디 브루노의 공식을 재구성하면 반해석적 방법에서 계산 효율성이 어떻게 향상되는가?
• RQ3 이 방법을 비선형 초기값 문제에 적용했을 때, 얼마나 정확한 닫힌 형태의 해를 도출할 수 있는가?
• RQ4 기존 방법(예: 아도미안 분해법 또는 호모토피 편중법)에 비해 이 방법이 요구하는 연산 수에서 어떤 계산적 이점이 있는가?

주요 결과

• 비선형 항 $ f(u(t)) $ 의 미분 변환은 부분 일반 벨 다항식을 사용해 성공적으로 표현되었으며, 이는 기호 도함수 없이도 변환 계수 $ U(k) $ 에서 직접 계산 가능하게 하였다.
• 첫 번째 예제에서는 계수 패턴이 $ U(k)[0] = 0 $ ( $ k $ 가 짝수일 때), $ U(k)[0] = (-1)^{(k-1)/2} \frac{2^k}{k!} $ ( $ k $ 가 홀수일 때) 로 나타나, 정확한 해 $ u(t) = \sin 2t $ 를 도출하였다.
• 두 번째 예제에서는 계수들이 $ \sin 2t $ 의 테일러 급수와 정확히 일치하여, 패턴 인식을 통해 정확한 해가 확인되었다.
• 알고리즘은 산술 연산과 재귀적 계산에만 의존하여 계산 효율적이며, 적분, 도함수 또는 초기 추정치가 필요 없었다.
• 계수 패턴이 알려진 함수와 일치할 경우, 이 방법은 표준 수치 근사보다 더 정확한 닫힌 형태의 해 식별이 가능하게 하였다.
• 이 방법은 경계값 문제를 포함한 다른 문제들로 일반화 가능하며, 전개점 근처에서 해석성을 유지하는 함수 해석적 근사 제공한다.