[论文解读] alphaCertified: certifying solutions to polynomial systems

Jonathan D. Hauenstein, Frank Sottile|arXiv (Cornell University)|Nov 4, 2010

Polynomial and algebraic computation参考文献 29被引用 13

一句话总结

本文介绍了 alphaCertified，一种利用 Smale 的 α-理论来严格认证多项式组数值解的计算工具。该工具采用精确的有理数或任意精度浮点数算术，验证牛顿法的二次收敛性，区分不同解，并识别实解，适用于平方系统和超定系统，包括在施伯特问题中验证了 126 个实解。

ABSTRACT

Smale's alpha-theory uses estimates related to the convergence of Newton's method to give criteria implying that Newton iterations will converge quadratically to solutions to a square polynomial system. The program alphaCertified implements algorithms based on alpha-theory to certify solutions to polynomial systems using both exact rational arithmetic and arbitrary precision floating point arithmetic. It also implements an algorithm to certify whether a given point corresponds to a real solution to a real polynomial system, as well as algorithms to heuristically validate solutions to overdetermined systems. Examples are presented to demonstrate the algorithms.

研究动机与目标

• 解决数值同伦求解器在多项式组解的启发式验证中缺乏严格认证的问题。
• 实现 α-理论算法，利用精确有理数或任意精度浮点数算术，提供数值解的数学严格认证。
• 实现对近似解是否对应实解的认证，尤其针对实多项式组。
• 通过启发式验证随机平方子系统，将认证扩展至超定系统。
• 为数值代数几何领域的研究人员提供一种实用且用户友好的工具，以数学确定性验证解集。

提出的方法

• alphaCertified 计算 Smale 的 α-理论中的 α(f,x) 常数，以确定牛顿法在多项式组 f 的点 x 处是否能二次收敛至解。
• 对于平方多项式组的完全严格认证，使用精确有理数算术。
• 对于高精度认证，采用任意精度浮点数算术来评估 α(f,x)，并控制误差界。
• 通过分析系统的实结构以及牛顿迭代在复共轭下的行为，判断解是否为实解。
• 对于超定系统（N > n），生成多个随机平方子系统，对每个子系统应用 α-理论，并比较结果以启发式验证近似解在给定 δ 内彼此接近。
• 该软件与 Maple 集成，以协助构建复代数几何问题的输入文件。

实验结果

研究问题

• RQ1 α-理论能否有效实现，以提供多项式组数值解的严格认证？
• RQ2 alphaCertified 中的精确有理数算术能否保证平方系统解认证的数学正确性？
• RQ3 α-理论在多大程度上能区分多项式组中的不同解，并识别实解？
• RQ4 α-理论能否被调整以在牛顿法不动点可能并非实际解的超定系统中，启发式验证解？
• RQ5 在涉及瞬时曲线切平面的施伯特问题中，实解的个数存在何种下界？

主要发现

• alphaCertified 证实，在 1000 个所有实切平面的 Shapiro 猜想实例中，恰好有 126 个实解，从而在该设定下证实了该猜想。
• 在 42,000 个实例中（k=1，包含实与共轭复数切平面混合），程序识别出至少 100 个实解，支持存在非平凡下界的假设。
• 对于 k=2，24,000 个实例中至少有 100 个实解，其中 1,000 个实例得到 126 个实解，表明下界假设具有高度一致性。
• 在 23,500 个 k=3 的实例中，程序检测到至少 100 个实解；在 22,500 个 k=4 的实例中，最小值仍为 100，进一步支持存在 100 个解的下界。
• 在 1,000 个 Shapiro 猜想测试实例中，软件正确识别出全部 126 个解为实解，表明在 256 位精度下具有高度可靠性。
• 在 100,000 个实例（k=0 至 4）中，实解数量的频率表显示一致的下界，支持此类施伯特问题中实解最小数量的猜想。