[論文レビュー] Polynomial Homotopies for Dense, Sparse and Determinantal Systems

Jan Verschelde|ArXiv.org|Jul 12, 1999

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

ひとこと要約

本稿では、密度的、スパース、行列式多項式系に対して最適な多項式ホモトピー続行法を提示する。Bézout、混合体積、Schubert計算による根数の計算を用いて、すべての解軌道が収束することを保証するホモトピーを構築する。主な貢献は、3つのクラスすべてにおいて一般系に対して線形計算複雑性を達成する数値的解法のフレームワークを提供することであり、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ホモトピーを導入し、一般系において最適な経路追跡を保証すること。

提案手法

• 係数パrameter継続の原則に基づくホモトピーを構築し、スタート系の解をターゲット系の解へ変形する。
• 密度的系にはBézoutの上限、スパース系にはニュートン多面体を用いた混合体積、行列式系にはSchubert計算を用いて正確な根数を計算する。
• ニュートン多面体の混合分割に基づく多面体的ホモトピーを採用し、一般系において発散路が存在しない最適な経路追跡を達成する。
• コンパactificationと斉次座標を用いて、経路追跡中に数値的安定性と幾何的忠実性を確保する。
• 正次元解成分上の一般点を計算する埋め込み法を導入し、数値代数幾何学におけるグローバル収束を可能にする。
• 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のcyclic系では34,940個の根と1,747個の成分を有しており、根数の計算だけでは不十分であり、成分を考慮したホモトピー法の必要性が示された。
• PHCソフトウェアは、コンピューターグラフィックス、制御理論、機械的設計からの100件以上の応用駆動系を正常に解き、本フレームワークの汎用性と頑健性を検証した。