[论文解读] Polynomial Homotopies for Dense, Sparse and Determinantal Systems

Jan Verschelde|ArXiv.org|Jul 12, 1999

Polynomial and algebraic computation参考文献 109被引用 24

一句话总结

本文提出了针对稠密、稀疏和行列式多项式系统的最优多项式同伦连续方法，利用Bézout、混合体积和Schubert演算的根数计数方法构建同伦，确保所有解路径收敛。关键贡献在于为这三类系统中的通用实例提供了具有线性计算复杂度的数值求解框架，通过PHC软件实现并验证了在代数几何与工程领域的应用。

ABSTRACT

Numerical homotopy continuation methods for three classes of polynomial systems are presented. For a generic instance of the class, every path leads to a solution and the homotopy is optimal. The counting of the roots mirrors the resolution of a generic system that is used to start up the deformations. Software and applications are discussed.

研究动机与目标

• 开发数值稳定且最优的同伦连续方法，确保对所有孤立解的收敛性。
• 通过同伦变形，将形式化的代数几何根数计数（Bézout、Bernshteïn、Schubert）与有效的数值计算相连接。
• 通过利用几何与组合结构，实现对机制设计、控制理论和枚举几何等应用中复杂系统的高效求解。
• 提供一个计算框架——在PHC软件中实现——能够以高精度和高效率求解多项式系统，包括完全实解的情况。
• 通过引入多面体同伦与SAGBI同伦，解决数值延续中路径发散的挑战，确保通用实例下最优的路径追踪。

提出的方法

• 基于系数-参数延续原理构建同伦，将起始系统的解变形为目标系统的解。
• 对稠密系统使用Bézout界，对稀疏系统使用牛顿多面体的混合体积，对行列式系统使用Schubert演算，以计算精确的根数。
• 采用基于牛顿多面体混合剖分的多面体同伦，实现无发散路径的最优路径追踪，适用于通用实例。
• 通过紧化与齐次坐标的应用，确保路径追踪过程中的数值稳定性和几何保真度。
• 引入嵌入方法，以计算正维解分支上的通用点，从而在数值代数几何中实现全局收敛。
• 利用PHC软件包实现并测试同伦方法，支持在实际应用中使用多面体与SAGBI同伦。

实验结果

研究问题

• RQ1 如何使同伦连续方法对稠密、稀疏和行列式多项式系统的通用实例达到最优，即无发散路径？
• RQ2 代数几何中的形式化根数计数（Bézout、Bernshteïn、Schubert）在多大程度上可转化为求解多项式系统的有效数值算法？
• RQ3 同伦方法能否系统性地应用于寻找多项式系统中的完全实解，特别是在力学与控制理论问题中？
• RQ4 当通过多面体与行列式形式化利用多项式系统的几何与组合结构时，求解其计算复杂度如何？
• RQ5 如何设计数值软件以支持灵活且结构感知的同伦方法，使其在效率与精度上优于黑箱求解器？

主要发现

• 对于n=7的循环系统（cyclic-n），多面体同伦是最优的，所有924条路径均收敛至有限解；而对于n≥8，混合体积会高估，且出现分量。
• Stewart-Gough平台问题（stewgou40）恰好有40个实解，通过最优同伦连续方法确认，解决了机器人学中长期存在的开放问题。
• 极点配置问题（pole28sys）实现了完全实解实例，所有1,430个解均为实解且孤立，经SAGBI同伦验证。
• cassou系统混合体积为24，但有8条路径发散至无穷远，需借助多面体终点策略分离病态的有限根。
• 当n=8时，循环系统有34,940个根，包含1,747个分量，表明仅依赖根数计数不足以应对，需采用分量感知的同伦方法。
• PHC软件成功求解了超过100个应用驱动的系统，涵盖计算机图形学、控制理论与机械设计，验证了该框架的通用性与鲁棒性。