[论文解读] Numerically validating the completeness of the real solution set of a system of polynomial equations

Daniel A. Brake, Jonathan D. Hauenstein+1 more|arXiv (Cornell University)|Feb 1, 2016

Polynomial and algebraic computation参考文献 38被引用 7

一句话总结

本文提出一种结合数值代数几何与平方和规划的数值方法，用于验证多项式方程组计算得到的实解集是否完整。通过检验解集的Zariski闭包是否与理想的真实根式匹配，该方法利用半定规划实现完整性认证，成功验证了具有有限和无限实解的系统（包括含不等式的系统）的完整性。

ABSTRACT

Computing the real solutions to a system of polynomial equations is a challenging problem, particularly verifying that all solutions have been computed. We describe an approach that combines numerical algebraic geometry and sums of squares programming to test whether a given set is "complete" with respect to the real solution set. Specifically, we test whether the Zariski closure of that given set is indeed equal to the solution set of the real radical of the ideal generated by the given polynomials. Examples with finitely and infinitely many real solutions are provided, along with an example having polynomial inequalities.

研究动机与目标

• 解决验证多项式方程组所有实解是否已被计算的问题，特别是在存在无穷多解或不等式的情况下。
• 为数值求解器开发一种全局停止准则，以确认解集的完备性。
• 提供一种结合数值计算、代数几何与优化技术的认证方法。
• 通过确保捕获完整的实解集，实现对实代数簇的可靠分解。

提出的方法

• 该方法首先计算一个数值近似集 $ S $，它是实解集Zariski闭包的子集。
• 然后在 $ S $ 上进行多项式插值，生成真实根式理想的候选生成元。
• 对每个候选生成元，应用平方和（SOS）规划，通过半定规划验证其是否属于真实根式。
• 验证过程独立处理每个多项式，支持并行化处理并提升数值鲁棒性。
• 该方法利用如下事实：若所有在 $ S $ 上消失的多项式均属于真实根式，则 $ S $ 在实解集中是Zariski稠密的。
• 使用数值同伦延续法与梯度下降同伦法对实解进行采样，尤其适用于正维数分量。

实验结果

研究问题

• RQ1 能否对数值计算得到的多项式系统实解集进行认证，以确认无实解遗漏？
• RQ2 如何利用平方和规划验证一组多项式是否生成理想的真实根式？
• RQ3 Zariski闭包在认证实解集完备性中起什么作用？
• RQ4 该方法能否处理具有无穷多实解的系统或涉及多项式不等式的系统？
• RQ5 如何将数值方法与代数几何结合，以实现解集完备性的全局有效认证？

主要发现

• 该方法成功验证了六变量系统（含12个实解）的实解集是完备的，确认计算集的Zariski闭包与真实根式一致。
• 对于Seiler系统（一条含实点的一维曲线），该方法通过验证复曲线上存在光滑实点，确认理想等于其真实根式。
• 在 $3\times3$ $\phi^4$ 能量景观中，该方法认证了仅存在三个实临界点，与精确解集一致。
• 通过将验证框架扩展至半代数集，该方法成功认证了含多项式不等式的系统的完备性。
• 所有真实根式候选生成元均通过SOS规划在 $\alpha \leq 2$ 条件下得到验证，证明了该方法的数值可靠性。
• 该方法支持多项式的并行验证，且无需符号域扩张，提升了效率与可扩展性。