[논문 리뷰] Computing the Inverse Mellin Transform of Holonomic Sequences using Kovacic's Algorithm

Jakob Ablinger|arXiv (Cornell University)|2018. 01. 02.

Numerical methods for differential equations인용 수 12

한 줄 요약

이 논문은 해석학적 해를 초월하여 리우빌리안 해를 포함하는 선형 미분방정식을 해결하는 코바치치 알고리즘의 통합을 통해, 해석학적 수열의 역 멀티노 변환을 계산할 수 있도록 HarmonicSums 컴퓨터 대수 패키지를 확장한다. 이 방법은 대수적 및 로그 함수를 포함하는 반복적 적분 형태로 역 멀티노 변환을 계산할 수 있게 하여 기호 적분 및 특수 함수 이론에서 해결 가능한 문제의 범위를 크게 넓힌다.

ABSTRACT

We describe how the extension of a solver for linear differential equations by Kovacic's algorithm helps to improve a method to compute the inverse Mellin transform of holonomic sequences. The method is implemented in the computer algebra package HarmonicSums.

연구 동기 및 목표

• 기존 방법이 d’Alembertian 해를 갖는 해석학적 수열에 대해서만 역 멀티노 변환을 계산할 수 있는 한계를 극복하기 위해.
• HarmonicSums 내의 선형 미분방정식 해법을 확장하여 d’Alembertian 해가 아닌 해를 처리할 수 있도록 코바치치 알고리즘을 적용하기 위해.
• 리우빌리안 해를 갖는 미분방정식을 유도하는 해석학적 수열의 역 멀티노 변환 계산을 가능하게 하기 위해.
• 대수적 및 로그 함수를 포함하는 반복적 적분 형태로 역 멀티노 변환을 체계적으로 표현하는 방법을 제공하기 위해.
• 일반화된 조화합수, 이항계수, 중첩 합을 포함한 명시적 예제를 통해 방법의 효과성을 입증하기 위해.

제안 방법

• 수열과 그 생성함수의 해석학적 성질을 활용하여, 역 멀티노 변환 문제를 선형 미분방정식의 해법 문제로 환원하기 위해.
• 해석학적 닫힘 성질을 이용하여 멀티노 변환 $ F(n) $의 해석학적 재귀식으로부터 미지의 함수 $ f(x) $에 대한 선형 미분방정식을 유도하기 위해.
• 코바치치 알고리즘을 적용하여 유도된 선형 미분방정식을 해결하고, d’Alembertian 해를 초월하여 리우빌리안 해를 포함한 해를 확장하기 위해.
• 반복적 적분과 대수적 함수로 구성된 미분방정식의 일반해를 구성하기 위해.
• 수치적 또는 기호적 적분을 통해 계산된 멀티노 변환 $ F(n) $의 초기 조건을 일치시켜 적분 상수를 결정하기 위해.
• 일반해와 초기 조건을 조합하여 $ f(x) $의 닫힌 형태 표현, 즉 역 멀티노 변환을 얻기 위해.

실험 결과

연구 질문

• RQ1 코바치치 알고리즘이 컴퓨터 대수 시스템에 효과적으로 통합되어 d’Alembertian 해를 초월하는 역 멀티노 변환 문제의 해를 확장할 수 있는가?
• RQ2 어떤 종류의 해석학적 수열이 d’Alembertian 해는 아니지만 리우빌리안 해를 갖는 미분방정식을 유도하는가?
• RQ3 d’Alembertian 해가 부족할 경우, 역 멀티노 변환은 어떻게 반복적 적분 형태로 표현할 수 있는가?
• RQ4 대수적 및 로그 함수를 포함하는 해를 포함하도록 해법기의 확장을 계산적으로 어떻게 영향을 미치는가?
• RQ5 이 방법은 일반화된 조화합수, 이항계수, 비d’Alembertian 구조를 갖는 중첩 합에 대해 체계적으로 적용될 수 있는가?

주요 결과

• 코바치치 알고리즘의 HarmonicSums 통합은 관련 미분방정식의 해가 리우빌리안일 경우, d’Alembertian 클래스를 초월하여 해석학적 수열의 역 멀티노 변환 계산을 가능하게 하였다.
• 수열 $ F(n) = \frac{4}{9} \frac{(3n)!}{(2n)!n!} $에 대해, 역 멀티노 변환은 $ f(x) = \frac{\text{const}}{\text{대수적 및 루트 표현식}} $ 형태로 표현되며, 모멘트를 일치시켜 상수를 결정하였다.
• 해당 방법은 $ \binom{3i}{i} \frac{1}{i} $의 역 멀티노 변환을 성공적으로 계산하였으며, 삼중 적분과 대수적 함수를 포함하는 해를 도출하였고, 최종 표현식은 $ \frac{\text{const}}{27x - 4} $와 중첩된 적분을 포함하였다.
• $ \binom{4n}{2n} $의 역 멀티노 변환은 $ \frac{1}{\text{루트 표현식}} $을 포함하는 적분 형태로 표현되었으며, 이는 대칭적인 이항계수를 다룰 수 있는 능력을 보여주었다.
• $ \frac{1}{n \binom{3n}{n}} $의 역 멀티노 변환은 $ \frac{1}{\text{제곱근의 합의 세제곱근}} $을 포함하는 적분 형태로 주어졌으며, 이는 특이적이고 다항식이 아닌 생성함수를 다룰 수 있는 능력을 보여주었다.
• 이 방법은 $ \binom{3i}{i} \frac{1}{i^2} $를 포함하는 수열에 대해 기호적으로 역 멀티노 변환을 계산할 수 있게 하였으며, 삼중 적분과 로그 항을 포함하는 해를 도출하여 고차 특이성을 갖는 경우에도 방법의 강건성을 확인하였다.