[论文解读] Solving Parameterized Polynomial Systems with Decomposable Projections

Carlos Améndola, Lindberg, Julia+1 more|arXiv (Cornell University)|Dec 28, 2016

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

一句话总结

本文提出一种基于单值性的新颖算法，用于求解参数多项式系统，通过利用可分解投影来降低计算复杂度，尤其在伽罗瓦群为非本原时效果显著。该方法采用迹测试作为停止准则，并通过证明解的数量为1442而非此前推测的1443，解决了高斯混合模型中的一个长期悬而未决的猜想。

ABSTRACT

The Galois group of a parameterized polynomial system of equations encodes the structure of the solutions. This monodromy group acts on the set of solutions for a general set of parameters, that is, on the fiber of a projection from the incidence variety of parameters and solutions onto the space of parameters. When this projection is decomposable, the Galois group is imprimitive, and we show that the structure can be exploited for computational improvements. Furthermore, we develop a new algorithm for solving these systems based on a suitable trace test. We illustrate our method on examples in statistics, kinematics, and benchmark problems in computational algebra. In particular, we resolve a conjecture on the number of solutions of the moment system associated to a mixture of Gaussian distributions.

研究动机与目标

• 通过利用伽罗瓦群中的结构对称性，开发一种高效求解参数多项式系统的算法。
• 通过避免符号消去，采用数值同伦延续方法，解决高次系统中消去带来的计算瓶颈。
• 通过新型迹测试算法，为基于单值性的求解器提供可靠的停止准则。
• 解决代数统计中的开放猜想，特别是关于高斯混合模型中解的数量问题。
• 在基准问题（如循环n次根系统和运动链机制）上展示该方法的有效性。

提出的方法

• 该方法利用可分解投影，其中投影 π: V(F) → ℂᵏ 可分解为 α∘β，从而将解空间划分为更小且具有对称性的子组件。
• 提出适配于可分解投影的见证集，通过数值同伦延续实现对解分量的高效追踪。
• 开发一种新型迹测试算法，以验证解集的完备性，替代依赖于次数估计或总路径计数的方法。
• 利用单值群作用探索解纤维，通过群不变坐标实现对称性约化，以最小化路径追踪量。
• 对于具有二面体或其他群对称性的系统，使用雷诺olds算子构造中间不变量，以定义分解结构。
• 该算法使用数值代数几何工具实现，路径追踪在β层投影上执行，从而降低计算负载。

实验结果

研究问题

• RQ1 在参数多项式系统中，可分解投影是否可用于减少数值同伦延续中的路径追踪数量？
• RQ2 伽罗瓦群的非本原性是否意味着可被算法利用的结构性简化？
• RQ3 能否构建一种迹测试，以可靠地确认在可分解设置下解集的完备性？
• RQ4 高斯混合模型的矩系统中真实解的数量是多少？是否与先前的猜想一致？
• RQ5 与总次数法和多面体同伦方法相比，该新型单值性算法在基准系统上的效率和准确性如何？

主要发现

• 本文通过证明高斯混合模型矩系统中的解数为1442而非先前推测的1443，解决了代数统计中的一个猜想。
• 在四杆连杆机构运动学问题中，该方法将路径追踪数从 3! × 1442 = 8652 条减少至仅5028条，显著提升了效率。
• 对于循环5次根系统，基于迹测试的单值法平均追踪3664.2条路径，而多面体方法仅需70条，但后者实现了高可靠性的完整解恢复。
• 在运行时间方面，单值法（次数变体）优于总次数法和多面体同伦方法，n=9循环根问题中实现10倍加速（5.567秒 vs. 677.038秒）。
• 迹测试算法在所有测试案例中（包括高次系统）成功确认了解集的完备性，未出现误报或遗漏解。
• 该方法成功利用二面体对称性对循环n次根系统进行了分解，通过雷诺olds算子不变单项式构造了非平凡的因子分解 π = α∘β。