[论文解读] Quadratic Interval Refinement for Real Roots

John Abbott|arXiv (Cornell University)|Mar 6, 2012

Numerical Methods and Algorithms被引用 4

一句话总结

本文提出了一种新型的根精化算法——二次区间精化（QIR），该算法结合了二分法的鲁棒性与牛顿法的二次收敛性。通过使用离散化线性插值和自适应精化因子，且无需计算导数，QIR在任意精度有理数算术中实现了具有竞争力的性能，同时保证收敛性并提供显式的误差界。

ABSTRACT

We present a new algorithm for refining a real interval containing a single real root: the new method combines characteristics of the classical Bisection algorithm and Newton's Iteration. Our method exhibits quadratic convergence when refining isolating intervals of simple roots of polynomials (and other well-behaved functions). We assume the use of arbitrary precision rational arithmetic. Unlike Newton's Iteration our method does not need to evaluate the derivative.

研究动机与目标

• 开发一种结合二分法鲁棒性与牛顿法快速收敛性的根精化算法。
• 消除牛顿法中因导数难以计算而带来的导数计算需求这一主要限制。
• 通过在整个精化过程中保持包含根的区间，确保收敛性并提供显式的误差界。
• 为任意精度算术提供一种实用且高效的牛顿迭代替代方案，尤其适用于多项式和性质良好的函数。
• 解决标准牛顿迭代在有理数算术中分子分母大小指数级增长的问题。

提出的方法

• QIR 使用一系列基于离散化线性插值的缩小步骤，对包含单个单重根的隔离区间进行精化。
• 每个缩小步骤将当前区间划分为 N 个相等的子区间，并利用线性插值预测包含根的子区间。
• 该方法在预测子区间的端点处计算函数值；若检测到符号变化，则更新区间；否则该步失败。
• 每次成功缩小后，精化因子 N 会增加，从而在区间缩小过程中实现二次收敛。
• 该算法使用任意精度有理数算术，确保中间值的精确表示并提供显式的误差界。
• 针对 N=4 的特殊情况进行了实现，以提高效率，并避免早期迭代中的不必要的精度损失。

实验结果

研究问题

• RQ1 能否设计一种根精化方法，使其兼具二分法的鲁棒性与牛顿法的二次收敛性？
• RQ2 是否可能在不计算函数导数的情况下实现二次收敛？
• RQ3 在任意精度算术中，无导数、基于区间的精化方法与牛顿迭代相比性能如何？
• RQ4 该方法能否避免有理数牛顿迭代中常见的分子分母大小指数级增长问题？
• RQ5 附近复根对新方法收敛行为有何影响？

主要发现

• 当函数在根附近具有有效的泰勒展开时，QIR 对单重根实现了二次收敛，收敛速度与牛顿法一致。
• 对于多项式 $x^5 - 2$，QIR 在六次迭代内将根的误差控制在 $2^{-32}$ 以内，且所有中间值的位数不超过 50 位，而牛顿法产生的近似值分母超过 17,000 位。
• 在实验中，QIR 在 6,600 秒内将 $f_2 = x^{50} + (10^{50}x - 1)^3$ 的根精化至 $10^{-10000}$ 精度，表现出即使在根接近共轭复根对时仍具鲁棒性。
• QIR 与 GMP 的专用根函数性能相当：对于 $\sqrt[3]{3}$，QIR 在 5.1 秒内达到 $10^6$ 位精度，而 GMP 耗时 17 秒。
• 对于 $\sqrt{5}$，QIR 比 GMP 的 mpz_sqrt 慢（35 秒 vs 8.2 秒，对应 $10^7$ 位精度），但在导数计算不切实际的一般函数场景下仍表现出极高的有效性。
• 该方法通过保持区间并使用受控精化，避免了有理数算术中分子分母无界增长的问题，使其适用于精确计算环境。