[論文レビュー] Geometric involutive bases for positive dimensional polynomial ideals and SDP methods

Greg Reid, Fei Wang+1 more|arXiv (Cornell University)|May 28, 2014

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

ひとこと要約

本稿は、正次元多項式イデアルにおける記号的・数値的幾何的インボリューション基底を、半定値計画法（SDP）手法と統合することで進展させ、実根の効率的計算を可能にする。次数低減技術と包含テストを導入し、中間系を単純化することで、モーメント行列のランク安定化と射影に基づく刈り込みによって正確な実根基底を達成する。

ABSTRACT

Geometric involutive bases for polynomial systems of equations have their origin in the prolongation and projection methods of the geometers Cartan and Kuranishi for systems of PDE. They are useful for numerical ideal membership testing and the solution of polynomial systems. In this paper we further develop our symbolic-numeric methods for such bases. We give methods to explicitly extract and decrease the degree of intermediate systems and the output basis. Algorithms for the numerical computation of involutivity criteria for positive dimensional ideals are also discussed. We were also motivated by some remarkable recent work by Lasserre and collaborators who employed our prolongation projection involutive criteria as a part of their semi-definite based programming (SDP) method for identifying the real radical of zero dimensional polynomial ideals. Consequently in this paper we begin an exploration of the interaction between geometric involutive bases and these methods particularly in the positive dimensional case. Motivated by the extension of these methods to the positive dimensional case we explore the interplay between geometric involutive bases and the new SDP methods.

研究の動機と目的

• 正次元多項式イデアルへの幾何的インボリューション基底法の拡張を図り、従来の記号的・数値的アプローチの限界を克服すること。
• 中間系を単純化し、出力基底を改善することで、数値的安定性と効率性を高める次数低減技術の開発。
• 多項式系の実根を計算するための最新のSDPに基づく手法と幾何的インボリューション基底の相乗効果の解明。
• 数値的設定においてカルタンのインボリューション基準を正次元イデアルに適応し、実装可能性を確保すること。
• 正次元ケースにおける解の計算に、延長-射影アルゴリズムとモーメント行列補完を統合するフレームワークの提供。

提案手法

• カルタンの基準を正次元イデアルに一般化し、幾何的延長-射影アルゴリズムを用いてインボリューション基底を計算する。
• 反復処理中に高次で冗長な系を除外するため、アルゴリズム4.1における包含テストを適用し、行列サイズと計算コストを低減する。
• M をモーメント行列、ker M をその核とするサブルーチン gen(ker M) を用いてモーメント行列を構築し、生成子を計算する。
• ランク安定化（rank(M) = dim ker gif(Q)）を終了基準として用い、ラッセールの零次元法に類似した条件を満たす。
• 次元表（例：図5）を用いて記号空間からの射影系を抽出し、次数低減と単純化を実現する。
• 正規化と代入（例：z = 2）を適用し、高次生成子を単純化することで、最小で解釈可能な基底を得る。

実験結果

リサーチクエスチョン

• RQ1 伝統的なグレブナー基底法がスケーラビリティの問題を抱える正次元多項式イデアルにおいて、幾何的インボリューション基底を効果的に拡張する方法は何か？
• RQ2 射影と次数低減技術は、中間系の単純化と数値的条件の改善にどのような役割を果たすか？
• RQ3 幾何的インボリューション基底とSDPに基づくモーメント行列手法の相互作用を、実根計算のために形式化・活用するにはどうすればよいか？
• RQ4 モーメント行列のランク安定化条件が、正次元ケースにおけるインボリューション基底計算の終了基準としてどの程度有効に機能するか？
• RQ5 アルゴリズム4.1における包含テストは、基底の完全性を損なわずに、冗長で高次の系を信頼性高く除外できるか？

主な発見

• アルゴリズムは、正次元系の実根の幾何的インボリューション基底を正常に計算でき、出力 {z−2, x(z−2), y(z−2), z(z−2), x²−y²+3} は手計算で検証済み。
• モーメント行列のランクは7であり、gen(ker M) の核の次元は13であった。これは良好に条件付けられた中間系であることを示している。
• 包含テストを適用した後、図5の第1列における ℓ=1 のみが保持され、システムサイズが縮小され、冗長成分が排除された。
• ℓ=2、次数1の射影系から、正規化と小項の除去を経て、簡略化された生成子 z−2 が得られた。
• z=2 を2次生成子に代入することで、x²−y²+3 に低減され、一貫性と最小形式が確認された。
• d=5=r で終了が達成され、出力に5つの生成子が含まれ、ランク安定化条件は零次元ケースにおけるラッセールの基準と一致した。