[논문 리뷰] Rational Hermite-Padé Approximants vs Padé Approximants. Numerical Results

Dobrolyubov, Egor O., Nikolay Ikonomov외 2명|arXiv (Cornell University)|2023. 06. 12.

Mathematical functions and polynomials인용 수 6

한 줄 요약

이 논문은 수치 예제를 통해 다가치 해석 함수를 근사할 때 유리 함수 허미트-파데 근사법(Hermite-Padé approximants, HPAs)이 대각 파데 근사법(Padé approximants, PAs)보다 유의미하게 뛰어나다는 것을 보여준다. 특히 리만 곡면의 여러 겹에서의 값 복원에 있어서 그렇다. 주요 결과는 동일한 수의 급수 계수를 사용할 경우, HPAs는 PAs보다 수렴 속도가 빠르며, 자유 바나드 폴 방정식의 주파수 함수와 같은 분석적 구조를 더 정확히 포착한다는 것이다. 이는 비물리적 겹에서 특히 그렇다.

ABSTRACT

The main purpose of the paper is to present some powerful data on the advantage of the rational approximation procedure based on Hermite-Padé polynomials over the Padé approximation procedure. The first part of the paper is devoted to some numerical examples in this direction. The second part will be devoted to some theoretical results. In particular, we demonstrate our ideas about the advantage of rational Hermite-Padé approximants over Padé approximants analyzing the analytical structure of the frequency function $ν$ of the free Van der Pol equation.

연구 동기 및 목표

• 다중가치 해석 함수를 근사하는 데 있어, 유리 함수 허미트-파데 근사법(Hermite-Padé approximants, HPAs)이 대각 파데 근사법(diagonal Padé approximants, PAs)보다 수치적으로 뛰어나다는 것을 입증하는 것.
• 복잡한 분지점 구조를 가진 함수, 예를 들어 자유 바나드 폴 방정식의 주파수 함수와 같은 함수에 대해, HPAs의 수렴성 및 해석적 연장 성질을 분석하는 것.
• HPAs가 리만 곡면의 여러 겹에서 함수 값을 복원하는 데 이론적 우수성을 보여주는 수치적 증거를 제공하는 것, 특히 물리적(영) 겹을 초월한 영역에서의 복원에 초점 맞춤.
• Nuttall의 영역과 평형 측도를 통한 대수적 함수의 분석적 구조 포착 능력 향상을 보여주는 수치적 예시를 제시함으로써, 이론적 제2부의 기초를 마련하는 것.

제안 방법

• 연구는 급수 $f(z) = \sum_{k=0}^\infty c_k z^k$의 첫 $2n+1$ 계수를 일치시켜 정의된 대각 파데 근사법 $[n/n]_f$를 사용한다. 이는 $(Q_n f - P_n)(z) = O(z^{2n+1})$ 를 만족시키며, $Q_n$과 $P_n$은 각각 차수 $n$ 이하의 다항식이다.
• 허미트-파데 근사법의 경우, 첫 $3m+2$ 계수를 포함하는 방정식 시스템을 바탕으로 유형 II 유리 근사 $P_{2m,1}/P_{2m,0}$ 를 $f(z)$에 대해 구성한다. 여기서 $N = 2n+1 = 3m+2$이다.
• 수렴 성질은 영역 $[-1,1]$ 과 $[a,b]$ 상에서 각각의 균형 측도 $\tau_E^\vphantom{p}$ 와 $\lambda_E^\vphantom{p}$ 로 수렴하는 영역 수세기 측도 $\frac{1}{n}\chi(p_n)$ 와 $\frac{1}{2m}\chi(p_{2m,0})$ 의 약한-* 수렴을 통해 분석된다.
• 점 渐近 수렴 속도는 파데 근사법의 경우 $|f_*(\zeta) - R(1/\zeta)|^{1/N} \to e^{-g_D(\zeta,\infty)}$ 이고, 허미트-파데 근사법의 경우 $e^{-\frac{1}{3}G_E^{\lambda_F}(\zeta) - g_D(\zeta,\infty)}$ 로 평가된다. 여기서 $g_D$ 는 그린 함수이고, $G_E^{\lambda_F}$ 는 그린의 잠재력이다.
• 복소수 영역 $D = \widehat{\mathbb{C}} \setminus \operatorname{supp}(\varkappa)$ 에서 오차의 감쇠 속도를 비교함으로써, HPA 오차는 추가적인 $\frac{1}{3}G_E^{\lambda_F}(\zeta)$ 항 덕분에 PA 오차보다 더 빠르게 감쇠됨을 보여준다.
• 자유 바나드 폴 방정식의 주파수 함수 $\nu(z)$ 에 대해 수치 실험을 수행하였으며, 그 급수 전개를 사용하여 PAs와 HPAs가 리만 겹 전반에 걸쳐 $\nu(z)$ 를 얼마나 잘 복원하는지 비교하였다.

실험 결과

연구 질문

• RQ1 유리 함수 허미트-파데 근사법(Hermite-Padé approximants)은 다가치 해석 함수를 근사할 때 대각 파데 근사법(Padé approximants)보다 수치적으로 어떻게 비교되는가?
• RQ2 허미트-파데 근사법은 리만 곡면의 비물리적 겹에서 파데 근사법보다 함수 값을 더 정확히 복원할 수 있는가?
• RQ3 복소수 영역에서 파데 근사법과 허미트-파데 근사법의 점 渐近 오차 감쇠 속도는 무엇이며, 그는 그린 함수와 평형 측도와 어떻게 관련되는가?
• RQ4 같은 수의 급수 계수를 사용할 때, 고차수의 허미트-파데 근사법(유형 II)은 대각 파데 근사법보다 더 빠른 수렴과 더 나은 해석적 연장 성질을 가지는가?
• RQ5 파데 근사법과 허미트-파데 근사법의 분모 다항식의 영점 분포는 리만 곡면 상에서 평형 측도로 수렴하는가?

주요 결과

• 같은 수의 급수 계수 $N = 2n+1 = 3m+2$ 를 사용할 경우, 유리 함수 허미트-파데 근사법 $P_{2m,1}/P_{2m,0}$ 는 대각 파데 근사법 $[n/n]_f$ 보다 더 빠른 오차 감쇠 속도를 보이며, 점 渐近 경계 $|f_* - R|^{1/N} \to e^{-\frac{1}{3}G_E^{\lambda_F}(\zeta) - g_D(\zeta,\infty)} < e^{-g_D(\zeta,\infty)}$ 를 통해 이를 입증하였다.
• HPA의 분모 다항식의 영점 수세기 측도는 $[-1,1]$ 상에서 평형 측도 $\lambda_E^\vphantom{p}$ 로 약한-* 수렴하는 반면, PA의 경우 $\tau_E^\vphantom{p}$ 로 수렴하여, 분석적 구조를 더 정교하게 포착한다는 것을 시사한다.
• 영역 $D = \widehat{\mathbb{C}} \setminus \operatorname{supp}(\varkappa)$ 에서 HPA 오차는 추가적인 $\frac{1}{3}G_E^{\lambda_F}(\zeta)$ 항 덕분에 PA 오차보다 엄격히 더 빠르게 감쇠되며, 이는 리만 곡면 상에서의 수렴 향상 정도를 반영한다.
• 자유 바나드 폴 방정식의 주파수 함수 $\nu(z)$ 에 대한 수치 결과는 HPAs가 여러 겹—특히 비물리적 겹—에서 PAs보다 더 정확하게 값을 복원한다는 것을 확인한다.
• HPA의 수렴은 Nuttall의 모든 겹에서 함수 값을 회복함을 보여주며, 이는 Komlov의 이론적 결과와 일치한다. 반면 PA는 오직 초기(영) 겹에서만 값을 회복한다.
• 대수적 함수의 경우, HPA 방법이 PA 방법보다 리만 곡면의 전체적인 구조를 더 효과적으로 포착함을 확인하였으며, 이는 함수에 다수의 분지점과 비자명한 단일화 성질이 존재할 경우 더욱 두드러진다.