[論文レビュー] GPGCD, an Iterative Method for Calculating Approximate GCD of Univariate Polynomials, with the Complex Coefficients

Akira Terui|arXiv (Cornell University)|Jul 12, 2010

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

ひとこと要約

本稿は、実数係数多項式の近似 GCD を計算するための反復的最適化手法 GPGCD を、複素数係数多項式に対しても適用可能にする拡張を行う。近似 GCD 問題を制約付き最適化問題に再定式化し、修正ニュートン法を用いて解くことで、STLN に基づく手法と同等の摂動を達成しつつ、中程度の次数の多項式では最大 30 倍の高速化を実現した。

ABSTRACT

We present an extension of our GPGCD method, an iterative method for calculating approximate greatest common divisor (GCD) of univariate polynomials, to polynomials with the complex coefficients. For a given pair of polynomials and a degree, our algorithm finds a pair of polynomials which has a GCD of the given degree and whose coefficients are perturbed from those in the original inputs, making the perturbations as small as possible, along with the GCD. In our GPGCD method, the problem of approximate GCD is transfered to a constrained minimization problem, then solved with a so-called modified Newton method, which is a generalization of the gradient-projection method, by searching the solution iteratively. While our original method is designed for polynomials with the real coefficients, we extend it to accept polynomials with the complex coefficients in this paper.

研究の動機と目的

• 実数係数多項式を対象として開発された GPGCD 法を、複素数係数多項式に対しても適用可能にする。
• 指定された次数の近似 GCD を計算する際、係数の摂動を最小限に抑えること。
• 同様の問題に対して、既存の手法（例：STLN に基づくアルゴリズム）と比較して計算効率を向上させること。
• 複素数多項式入力におけるノイズや悪条件性の影響に対しても、収束性と精度を保証すること。

提案手法

• 近似 GCD 問題を、係数摂動の 2-ノルムを最小化する制約付き最適化問題に変換する。
• 修正ニュートン法（勾配射影法の一般化）を用いて、反復的に制約付き最適化問題を解く。
• 解法は、結果的係数の (d−1)-番目の部分結果式がゼロであるという条件に依存しており、これにより次数 d の GCD が保証される。
• 複素数係数は、実部と虚部に分離し、最適化フレームワーク内で独立変数として扱う。
• アルゴリズムは、近似 GCD、摂動 ΔF と ΔG、および線形関係 A·F̃ + B·G̃ = 0 を満たす補因子 A と B を計算する。
• 収束制御のため、更新ベクトルのノルムに ε = 1.0×10⁻⁸ の閾値を設定している。

実験結果

リサーチクエスチョン

• RQ1 GPGCD 法は、精度を保持したまま、複素数係数多項式に対しても効果的に拡張可能か？
• RQ2 複素数係数多項式の近似 GCD 問題において、GPGCD 法の性能は STLN に基づく手法と比べて、摂動の大きさと計算時間の両面でどう異なるか？
• RQ3 修正ニュートン法は、複素数係数多項式の近似 GCD 問題に対し、収束性と安定性を保証するか？
• RQ4 多項式の次数とノイズレベルの変化が、複素数係数に対する GPGCD 法の精度と効率に与える影響は何か？

主な発見

• GPGCD 法は、STLN に基づく手法と同等の摂動を達成しており、全テストケースにおいて 3.72×10⁻³ から 4.84×10⁻³ の誤差を示した。
• GPGCD 法は、STLN に基づく手法と比較して 10～30 倍の高速化を達成しており、計算時間はそれぞれ 0.15～6.00 秒と 1.79～199.48 秒であった。
• 両手法とも 10 個のテストケースすべてで収束したため、次数 10～100 の多項式（近似 GCD の次数が 5～50）に対しても、頑健性と信頼性が確認された。
• GPGCD と STLN の反復回数はほぼ同一であり、実行時間の差にもかかわらず、収束行動は類似していることが示された。
• 多項式の次数が増加しても高い精度を維持しており、次数 30 以上では摂動レベルが約 4.5×10⁻³ に安定した。