[論文レビュー] Rational Hermite-Padé Approximants vs Padé Approximants. Numerical Results

Dobrolyubov, Egor O., Nikolay Ikonomov+2 more|arXiv (Cornell University)|Jun 12, 2023

Mathematical functions and polynomials被引用数 6

ひとこと要約

この論文は、数値例を通じて、多価な解析関数を近似する際、有理関数のハーモニック・ペーア近似（HPA）が対角ペーア近似（PA）を著しく上回ることを示している。特にリーマン面の複数のシートにわたる値の再構成において顕著である。主な結果は、同じ数の級数係数を用いる場合、HPAはPAよりも収束が速く、自由なヴァン・デル・ポール方程式の周波数関数のような解析的構造をより正確に捉えることができることである。これは、物理的でないシートでも同様に成り立つ。

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.

研究の動機と目的

• 多価な解析関数を近似する有理関数のハーモニック・ペーア近似（HPA）が対角ペーア近似（PA）を上回ることを、数値的に示すこと。
• 分岐点構造が複雑な関数、たとえば自由なヴァン・デル・ポール方程式の周波数関数に対して、HPAの収束性および解析的接続の性質を分析すること。
• HPAがリーマン面の複数のシートにわたる関数値を回復する理論的優位性を、数値的証拠として提示すること。特に物理的でない（ゼロでない）シートにおいて顕著である。
• Nuttallの領域および均衡測度を用いた代数関数の解析的構造をよりよく捉える方法として、HPAの有効性を示すことで、理論的第二部の基盤を築くこと。

提案手法

• 研究では、関数 $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$ 以下の多項式である。
• ハーモニック・ペーア近似では、$f(z)$ に対するタイプIIの有理関数近似 $P_{2m,1}/P_{2m,0}$ を用いる。これは最初の $3m+2$ 個の係数を含む方程式系から構成され、$N = 2n+1 = 3m+2$ である。
• 収束性は、零点数え測度 $\frac{1}{n}\chi(p_n)$ および $\frac{1}{2m}\chi(p_{2m,0})$ の弱-*収束を、区間 $[-1,1]$ および $[a,b]$ 上の均衡測度 $\tau_E^\vphantom{p}$ および $\lambda_E^\vphantom{p}$ に分析する。
• 漸近的収束速度は、$|f_*(\zeta) - R(1/\zeta)|^{1/N} \to e^{-g_D(\zeta,\infty)}$（PAの場合）および $e^{-\frac{1}{3}G_E^{\lambda_F}(\zeta) - g_D(\zeta,\infty)}$（HPAの場合）の極限を用いて評価する。ここで $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)$ に対して、その級数展開を用い、PAとHPAの両方がリーマン面の各シートにわたる $\nu(z)$ の再構成性能を比較する数値実験を実施している。

実験結果

リサーチクエスチョン

• RQ1 有理関数のハーモニック・ペーオ近似（HPA）は、対角ペーオ近似（PA）と比較して、多価な解析関数を近似する際にどのように数値的に優れているか？
• RQ2 ハーモニック・ペーオ近似（HPA）は、リーマン面の物理的でないシートにおける関数値を、ペーオ近似（PA）よりも正確に再構成できるか？
• RQ3 複素領域におけるペーオ近似およびハーモニック・ペーオ近似の漸近的誤差減衰率は何か？また、グリーン関数および均衡測度とどのように関係しているか？
• RQ4 同じ数の級数係数が与えられた場合、高次のハーモニック・ペーオ近似（タイプII）は、対角ペーオ近似よりも収束が速く、より良い解析的接続を実現できるか？
• RQ5 ペーオ近似およびハーモニック・ペーオ近似の分母多項式の零点分布は、リーマン面上でどのように均衡測度に弱-*収束するか？

主な発見

• 同じ数の級数係数 $N = 2n+1 = 3m+2$ を用いる場合、漸近的上限 $|f_* - R|^{1/N} \to e^{-\frac{1}{3}G_E^{\lambda_F}(\zeta) - g_D(\zeta,\infty)} < e^{-g_D(\zeta,\infty)}$ により、有理関数のハーモニック・ペーオ近似 $P_{2m,1}/P_{2m,0}$ は対角ペーオ近似 $[n/n]_f$ よりも誤差の減衰が速いことが示された。
• HPAの分母多項式の零点数え測度は、$[-1,1]$ 上の均衡測度 $\lambda_E^\vphantom{p}$ に弱-*収束するが、PAのそれは $\tau_E^\vphantom{p}$ に収束する。これは、解析的構造の近似がより洗練されていることを示している。
• 領域 $D = \widehat{\mathbb{C}} \setminus \operatorname{supp}(\varkappa)$ において、HPAの誤差はPAの誤差よりも厳密に速く減衰する。これは指数部に追加の $\frac{1}{3}G_E^{\lambda_F}(\zeta)$ 項があるためであり、リーマン面上での収束性の向上を反映している。
• 自由なヴァン・デル・ポール方程式の周波数関数 $\nu(z)$ に対する数値結果は、HPAが複数のシート、特に物理的でないシートにおける関数値をPAよりも正確に再構成できることを確認している。
• HPAの収束は、Nuttallのシートのうち最高のものののぞき、Komlovの理論的結果と整合的である。一方、PAは初期（ゼロ）シートでのみ関数値を回復している。
• 代数関数に対しては、HPA法がPA法よりもリーマン面全体の構造を捉えるのに効果的であることが確認された。特に、複数の分岐点と非自明なモノドロミーを持つ関数において顕著である。