[论文解读] 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

一句话总结

本文通过数值例子表明，有理赫尔mite-帕德逼近（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.

研究动机与目标

• 通过数值方法证明有理赫尔mite-帕德逼近（HPA）在逼近多值解析函数方面优于对角帕德逼近（PA）。
• 分析具有复杂分支点结构的函数（如自由范德波尔方程的频率函数）的HPA的收敛性与解析延拓性质。
• 提供数值证据，支持HPA在恢复黎曼面多个面片上函数值（尤其是物理面片之外）方面的理论优势。
• 通过展示HPA如何更优地捕捉代数函数的解析结构（借助Nuttall区域与平衡测度），为理论第二部分奠定基础。

提出的方法

• 本研究采用对角帕德逼近$[n/n]_f$，通过匹配幂级数$f(z) = \sum_{k=0}^\infty c_k z^k$的前$2n+1$个系数来定义，即$(Q_n f - P_n)(z) = O(z^{2n+1})$，其中$Q_n$与$P_n$为次数不超过$n$的多项式。
• 对于赫尔mite-帕德逼近，采用II型有理逼近$P_{2m,1}/P_{2m,0}$来逼近$f(z)$，其构造基于前$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)$项。
• 在自由范德波尔方程的频率函数$\nu(z)$上进行数值实验，利用其级数展开，并比较PA与HPA在重构$\nu(z)$于黎曼面各面片上的性能。

实验结果

研究问题

• RQ1 有理赫尔mite-帕德逼近与对角帕德逼近在数值逼近多值解析函数方面有何差异？
• RQ2 赫尔mite-帕德逼近是否能比帕德逼近更准确地重构黎曼面非物理面片上的函数值？
• RQ3 帕德与赫尔mite-帕德逼近在复域中的渐近误差衰减速率如何？其与格林函数及平衡测度有何关联？
• RQ4 在给定相同数量级数系数的情况下，使用高阶赫尔mite-帕德逼近（II型）是否比对角帕德逼近实现更快收敛与更优的解析延拓？
• RQ5 帕德与赫尔mite-帕德逼近的分母多项式零点分布如何收敛至黎曼面上的平衡测度？

主要发现

• 在相同级数系数数量$N = 2n+1 = 3m+2}$下，有理赫尔mite-帕德逼近$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误差严格快于PA误差，因其指数中多出$\frac{1}{3}G_E^{\lambda_F}(\zeta)$项，反映了在黎曼面上更优的收敛性。
• 对自由范德波尔方程频率函数$\nu(z)$的数值结果确认，HPA能比PA更准确地重构多个面片（尤其是非物理面片）上的函数值。
• HPA的收敛性被证明可恢复除最高Nuttall面片外的所有Nuttall面片上的函数值，与Komlov的理论结果一致；而PA仅能恢复初始（零）面片上的值。
• 研究证实，对于代数函数，HPA方法在捕捉完整的黎曼面结构方面比PA方法更有效，尤其当函数具有多个分支点与非平凡单值性时。