[论文解读] A new algorithm for computing Pade approximants

В. М. Адуков, Ibryaeva O. L|arXiv (Cornell University)|Dec 24, 2011

Statistical and numerical algorithms参考文献 7被引用 6

一句话总结

该论文提出了一种新算法，通过基于托普利茨矩阵核结构导出的最小次数分母，定义了约化帕德逼近。该方法通过确保分母在次数上唯一最小，避免了由数值舍入引起的弗鲁萨特双点（Froissart doublets）——即虚假的极点-零点对，从而在非正规情况下也能保证稳定且精确的有理逼近。

ABSTRACT

In paper a new definition of reduced Pade approximant and algorithm for its computing is proposed. Our approach is based on the investigation of the kernel structure of the Toeplitz matrix. It is shown that the reduced Pade approximant always has nice properties which classical Pade approximant possesses only in the normal case. The new algorithm allows us to avoid Froissart doublets appearance induced by computer roundoff in the non-normal Pade table.

研究动机与目标

• 解决在信号处理等应用中，帕德逼近因数值舍入而产生虚假极点-零点对（弗鲁萨特双点）的问题。
• 开发一种对经典帕德逼近的鲁棒替代方法，以克服在非正规情况下可能出现的非唯一性和不稳定性。
• 确保所得有理逼近保持精确的次数和亏值指标，这些性质通常仅在正规情况下才能保证。
• 提供一种数值稳定的算法，可直接从与幂级数系数相关的托普利茨矩阵的核中计算最小次数分母。

提出的方法

• 将约化解析帕德逼近定义为所有由托普利茨矩阵核导出的可能分母中，次数最小的有理函数。
• 利用托普利茨矩阵的核结构参数化所有可能的分母，并通过秩分析识别最小次数的分母。
• 应用定理 3.2，基于从截断托普利茨子矩阵的秩确定的基本指数 μ₁，计算第一个基本多项式（即最小次数分母）。
• 利用定理 4.2 和定理 4.3 检测并消除分子和分母中的虚假零系数，确保最终有理函数处于最简形式。
• 在 Maple 中实现该算法，使用过程 `ReducedPade(f, m, n, a, d)`，其中 d 控制系数删除，并通过根求解和与精确有理函数的比较验证结果。
• 通过数值实验表明，约化解析逼近在机器精度范围内与原函数完全匹配，且无虚假根。

实验结果

研究问题

• RQ1 是否可以定义一种帕德逼近，使其分母在非正规情况下也具有唯一最小的次数？
• RQ2 如何利用托普利茨矩阵的核结构计算最小次数分母，同时避免引入数值伪影？
• RQ3 在实际计算中，该算法在多大程度上能防止因浮点数舍入而产生的弗鲁萨特双点？
• RQ4 即使经典帕德逼近失效，约化解析帕德逼近仍保留哪些性质，如次数唯一性和无公共根？
• RQ5 该算法能否在 Maple 等符号计算系统中高效实现，同时保持精度并避免虚假极点-零点对？

主要发现

• 约化解析帕德逼近在常数因子意义下唯一确定，其分子和分母的次数精确，且亏值指标已知。
• 该算法通过分析托普利茨矩阵的核，成功计算出最小次数分母，避免了非唯一或虚假解。
• 在例 5.1 中，函数 $ f(z) = \frac{(z+1)(z-2)}{(z+2.1)(z-1)} $ 的 (4,4) 逼近准确恢复了极点 $ z = -2.1, 1 $ 和零点 $ z = -1, 2 $，无虚假根。
• 在例 5.2 中，函数 $ f(z) = \frac{z+1.01}{(z+2)(z-2.01)} $ 的 (2,3) 逼近准确恢复了零点 $ z = -1.01 $ 和极点 $ z = -2, 2.01 $，无弗鲁萨特双点。
• 在例 5.3 中，当 $ d = 1 $ 时保留零系数，导致出现虚假根（如 $ \sim 10^5 $）；但当 $ d = 0 $ 时删除零系数，可得到正确的极点和零点，证明了系数删除的必要性。
• 当 $ d = 0 $ 时，该算法始终能精确匹配原始有理函数的根结构，证明其对舍入引起的双点具有强鲁棒性。