[论文解读] Nonlinear Sequence Transformations: Computational Tools for the Acceleration of Convergence and the Summation of Divergent Series

Ernst Joachim Weniger|ArXiv.org|Jul 10, 2001

Iterative Methods for Nonlinear Equations参考文献 135被引用 12

一句话总结

本文系统综述了非线性序列变换作为加速收敛和求和发散级数的计算工具。重点介绍了Wynn的epsilon和rho算法等非正则非线性方法，尽管缺乏理论正则性，但通常优于传统的Padé逼近，且在科学与工程的数值应用中展现出显著的实用性优势。

ABSTRACT

Convergence problems occur abundantly in all branches of mathematics or in the mathematical treatment of the sciences. Sequence transformations are principal tools to overcome convergence problems of the kind. They accomplish this by converting a slowly converging or diverging input sequence $\{s_n \}_{n=0}^{\infty}$ into another sequence $\{s^{\prime}_n \}_{n=0}^{\infty}$ with hopefully better numerical properties. Padé approximants, which convert the partial sums of a power series to a doubly indexed sequence of rational functions, are the best known sequence transformations, but the emphasis of the review will be on alternative sequence transformations which for some problems provide better results than Padé approximants.

研究动机与目标

• 提供非线性序列变换作为克服数值分析中收敛问题的工具的系统性概述。
• 证明尽管非线性变换不具备正则性，但通常在收敛加速方面优于正则线性方法。
• 突出非线性变换在现代计算数学中的实际主导地位，及其在物理学、工程学和数值积分中的应用。
• 探讨组合多个序列变换的局限性，原因在于后续步骤中输入数据的数值污染。
• 将序列变换的适用范围从标量序列扩展至向量和矩阵问题，并讨论其在预测高阶微扰级数系数中的应用。

提出的方法

• 回顾从Richardson外推法和Aitken的Δ²过程等早期方法到现代非线性算法的序列变换发展历程。
• 分析Wynn的epsilon和rho算法、Brezinski的theta算法等非线性变换的理论与实际特性。
• 将非线性变换与Padé逼近进行比较，表明尽管Padé是一类著名的序列变换，但非线性替代方法在特定问题上常表现更优。
• 讨论用于预测微扰级数高阶系数的递归方案，尤其针对非厄米特的PT对称哈密顿量。
• 研究在数值求积中的应用，特别是利用序列变换改善半无限区间上振荡积分的收敛性。
• 考虑向量和矩阵序列的序列变换公式化，引用近期关于其向多维问题扩展的文献。

实验结果

研究问题

• RQ1 为何非线性序列变换在缺乏理论正则性的情况下，仍常在实际收敛加速中优于正则线性变换？
• RQ2 Wynn的epsilon算法或Brezinski的theta算法等非线性序列变换在特定级数类中，与Padé逼近相比能在多大程度上提供更优的收敛行为？
• RQ3 组合多个序列变换存在哪些局限性？为何输入数据在后续步骤中的数值污染会阻碍进一步的收敛改进？
• RQ4 如何利用序列变换预测微扰级数中的未知高阶系数，特别是在理论物理中？
• RQ5 序列变换在何种方式下可推广至向量和矩阵序列？这对应用数学中的数值方法有何影响？

主要发现

• Wynn的epsilon和rho算法、Brezinski的theta算法等非线性序列变换，通常在实际应用中相较于Padé逼近实现更优的收敛加速效果，尤其在复杂问题中表现突出。
• 理论正则性虽具吸引力，但可能限制线性变换的效率；非正则的非线性方法尽管缺乏对正确极限的保证，却常表现更优。
• 通过Aitken的Δ²过程、Wynn的epsilon算法和theta算法的递归方案，成功预测了高阶微扰级数系数，为理解非厄米特的PT对称哈密顿量提供了洞见。
• 序列变换在数值求积中极为有效，尤其在半无限区间上的振荡积分中，能显著改善收敛行为。
• 尽管在实践中表现成功，非线性变换的理论性质仍理解不足，且组合多个变换常因输入数据的数值污染而失败。
• 将序列变换扩展至向量和矩阵序列是当前活跃的研究领域，其在迭代求解器和多维数值问题中具有重要应用。