[论文解读] Two-dimensional generalization of the Muller root-finding algorithm and its applications

Plamen Fiziev, Denitsa Staicova|arXiv (Cornell University)|May 28, 2010

Model Reduction and Neural Networks参考文献 18被引用 7

一句话总结

本文提出了一种 Müller 根求解算法的二维推广，用于求解包含复变量的两个非线性超越方程组，尤其适用于涉及合流 Heun 函数的问题。在某些具有挑战性的情形下（如 Schwarzschild 与 Kerr 黑洞的准正规模频率计算），该方法的收敛速度优于牛顿法与 Broyden 法，即使导数计算失败或收敛缓慢。

ABSTRACT

We propose a new algorithm for solving a system of two nonlinear transcendental equations with two complex variables based on the Muller algorithm. The two-dimensional Muller algorithm is tested on systems of different type and is found to work comparably to Newton's method and Broyden's method in many cases. The new algorithm is particularly useful in systems featuring the Heun functions whose complexity may make the already known algorithms not efficient enough or not working at all. In those specific cases, the new algorithm gives distinctly better results than the other two methods. As an example for its application in physics, the new algorithm was used to find the quasi-normal modes (QNM) of Schwarzschild black hole described by the Regge-Wheeler equation. The numerical results obtained by our method are compared with the already published QNM frequencies and are found to coincide to a great extent with them. Also discussed are the QNM of the Kerr black hole, described by the Teukolsky Master equation.

研究动机与目标

• 解决牛顿法与 Broyden 法在求解涉及特殊函数（如 Heun 函数）的复数非线性系统时的局限性。
• 克服因合流 Heun 函数导致的导数计算困难与收敛缓慢等计算挑战。
• 开发一种无导数、快速收敛的算法，适用于理论物理中的谱问题，如黑洞的准正规模频率。
• 为标准迭代方法失效或收敛过慢的系统提供一种稳健的数值替代方案。
• 在实际物理问题（如描述黑洞准正规模频率的 Regge-Wheeler 方程与 Teukolsky 方程）上展示该算法的有效性。

提出的方法

• 将原始为一维设计、具有超线性收敛速度（约 1.84）的 Müller 算法推广至两个复变量。
• 基于三个初始点在二维空间中进行二次插值，以近似根，避免显式计算导数。
• 在二维框架内迭代调用一维 Müller 子程序，以逐步优化根的估计值。
• 通过调整一维子程序中的最大迭代次数 $ P $ 来调节收敛性，该参数显著影响整体性能与收敛速度。
• 将该方法应用于两个复数方程的方程组，特别关注函数排序与线性组合的使用，以提高收敛的可靠性。
• 在包含 Heun 函数的系统上进行数值测试，并将结果与牛顿法和 Broyden 法进行比较。

实验结果

研究问题

• RQ1 无导数的二维根求解算法是否能在包含复杂数、非多项式函数的系统中优于牛顿法与 Broyden 法？
• RQ2 在导数计算困难的合流 Heun 函数相关系统中，二维 Müller 算法的表现如何？
• RQ3 初始猜测值与函数排序对新算法收敛行为及根选择有何影响？
• RQ4 调整一维 Müller 子程序中的迭代次数 $ P $ 是否显著影响二维方法的收敛速度与精度？
• RQ5 该新算法能否可靠地计算由 Regge-Wheeler 方程与 Teukolsky 方程描述的 Schwarzschild 与 Kerr 黑洞的准正规模频率？

主要发现

• 在合流 Heun 函数相关的系统中，二维 Müller 算法的收敛速度优于牛顿法与 Broyden 法，即使导数计算失败或不准确。
• 对于 Schwarzschild 黑洞的准正规模频率，新方法计算得到的频率与已有文献结果高度一致，验证了其准确性。
• 在某些情况下，增加一维子程序中的最大迭代次数 $ P $ 反而减少了总计算时间，表明其依赖关系并非单调。
• 该算法对函数排序与初始条件敏感；通过重新排序或使用方程的线性组合，可有效防止发散或收敛至非期望根。
• 在 Heun 函数相关的系统中，尽管 Broyden 法通常具有鲁棒性，该方法的收敛速度仍优于 Broyden 法。
• 在具有实根或简单非线性的系统中，二维 Müller 方法的表现与牛顿法和 Broyden 法相当；但在涉及特殊函数的复杂数、奇异情形中表现尤为突出。