論文「Quaternion Matrices の左固有値の計算」の技術的サマリー

Michael Sebek（プラハ工科大学）によるこの論文は、四元数行列の**左固有値（left eigenvalues）**を計算するための実用的で信頼性の高いニュートン法に基づく手法を提案しています。非可換な四元数体の特性により、複素数行列とは異なる数学的・数値的課題が存在しますが、本論文は標準的な実数/複素数線形代数カーネルのみを用いて、任意のサイズの行列に対して再現性のある計算を実現しています。

以下に、問題定義、手法、主要な貢献、結果、および意義について詳細にまとめます。

---

1. 問題定義と背景

四元数行列の左固有値問題
四元数行列 $A \in \mathbb{H}^{n \times n}$ に対して、非ゼロの四元数ベクトル $v \in \mathbb{H}^n$ と四元数 $\lambda \in \mathbb{H}$ が存在し、以下の関係式を満たすとき、$\lambda$ を $A$ の左固有値と呼びます。
$$Av = \lambda v$$

既存の課題

• 非可換性: 四元数の乗法は非可換であるため、右固有値（$Ay = y\mu$）とは異なり、左固有値は四元数相似変換に対して不変ではありません。
• スペクトルの複雑さ: 左固有スペクトル $\sigma_\ell(A)$ は、孤立点だけでなく、連続的な集合（特に「球面成分」と呼ばれる 2 次元球面）を含む可能性があります。
• 数値計算の難しさ:
  • 固有ベクトルの右スケーリング（$v \to vq$）による自由度（ゲージ）を明示的に固定する必要があります。
  • 一般に $n$ 個の孤立固有値が存在するとは限らず、$n$ 個を超える孤立固有値や、固有スペクトルが不足する（deficient）現象が発生します。
  • 連続的な解集合（球面）が存在する場合、ヤコビアンが特異となり、標準的なニュートン法の局所収束理論が適用できません。
  • 従来の式変形（行列式の消去など）は、$n$ が増えると非線形問題の次数が急増し、数値的に不安定になります。

2. 提案手法：ゲージ制約付きニュートン法

本論文の核心は、固有値と固有ベクトルを同時に求めるための直接ニュートン法の定式化にあります。

2.1. ゲージ制約の導入

固有ベクトル $v$ は右から任意の四元数 $q$ を掛けた $vq$ も解となるため、ニュートン法を適用するにはこの自由度を固定する必要があります。

• 制約条件: 特定のインデックス $j$ をピボットとして選び、以下の 4 つの実数制約を課します。
  1. 正規化: $\|v\|_2 = 1$
  2. ピボット成分の実数化: $v_j \in \mathbb{R}$ （かつ $v_j > 0$ で符号を固定）
• これにより、解空間を局所的に切断（cross-section）し、実数座標系における正方系の方程式を構築します。

2.2. 実数埋め込み（Real Embedding）

四元数演算を直接扱うのではなく、四元数行列を実数行列に埋め込むことで、標準的な線形代数ソルバーを使用可能にします。

• 四元数 $q = a + bi + cj + dk$を実数ベクトル$(a, b, c, d)^T$として表現し、四元数左乗法を$4 \times 4$実数行列$L(q)$ で表現します。
• 行列 $A \in \mathbb{H}^{n \times n}$ を $4n \times 4n$ の実数行列 $\rho(A)$ に変換します。
• 固有値問題 $Av - \lambda v = 0$ は、この埋め込みを用いて実数線形方程式系として定式化され、ニュートンステップ（修正項の計算）が標準的な実数線形ソルバーで解かれます。

2.3. アルゴリズムの流れ

1. マルチスタート戦略: 複数の初期値からニュートン反復を開始し、解空間を探索します。
2. ニュートン修正: 欠陥方程式 $Av - \lambda v = 0$ とゲージ制約を同時に満たすように、$(4n+4) \times (4n+4)$ の実数線形システムを解いて修正項 $(\Delta \lambda, \Delta v)$ を求めます。
3. リゲージング（Regauging）: 各ステップでベクトルを再正規化し、ゲージ条件を厳密に満たします。
4. 収束判定と重複排除: 残差（residual）が閾値以下であれば収束とみなし、発見された解が既存の解と重複していないか（距離と残差に基づき）を判定して重複を排除します。
5. 球面成分の検出: 収束した点が球面上に分布している場合、クラスタリングと球面フィッティングを行い、球面成分として検出・記述します。

3. 主要な貢献

1. 数学的定式化の確立:

   • 左固有値対 $(\lambda, v)$ に対するニュートン反復を導出し、ゲージ制約を明示的に固定することで、実解析的な正方系を構築しました。
   • 孤立した単純な固有値対に対して、局所的な二次収束が保証されることを証明しました（ヤコビアンが正則である場合）。

2. 頑健なソルバーの構築:

   • 試行ベースのマルチスタート戦略、受け入れ/リスタート/重複排除ロジックを組み合わせたグローバルソルバーを開発しました。
   • 特異な状況（多重固有値、球面成分、固有スペクトル不足）に対しても、再現性のある候補と診断情報を提供します。

3. 実装と検証:

   • 標準的な実数/複素数線形代数ライブラリのみを使用する、コンパクトなスタンドアロン MATLAB 実装を提供しました。
   • 文献の既知例（$2 \times 2$ や $3 \times 3$ のケース）およびベンチマークセット（ランダム行列、エルミート行列、疎行列など）で広範に検証しました。
   • 最大 $64 \times 64$ のサイズまで、再現性のある計算と証明（残差に基づく証明書）を達成しました。

4. 結果と発見

• 精度: 孤立した固有値に対して、機械精度に近い残差（$10^{-16}$ レベル）を達成し、二次収束を確認しました。
• 非一般的な現象の検出:
  • $n$ を超える孤立固有値: $3 \times 3$ 行列で 5 つの異なる孤立左固有値を持つ例を特定しました。
  • 固有スペクトル不足（Left-spectrum deficiency）: $4 \times 4$ 行列で、期待される 4 つではなく 2 つの孤立固有値しか存在しないケースを再現しました。
  • 球面成分と孤立点の共存: 球面成分と孤立点が混在する構成を特定し、ニュートン法が球面上の異なる点に収束することを利用し、球面モデル（中心・半径）をフィッティングして検出しました。
• 性能: 文献の例やベンチマークセットにおいて、高い成功率（100% の固有値発見）を記録しました。

5. 意義と結論

この研究は、四元数行列の左固有値計算において、理論的な複雑さ（非可換性、連続スペクトル）を克服するための実用的な計算フレームワークを提供しました。

• 実用性: 四元数演算を直接実装する必要がなく、既存の高性能な実数線形代数ソルバーを流用できるため、大規模な行列計算にも拡張可能です。
• 信頼性: 単なる近似値の提示にとどまらず、残差に基づく「証明書（certificate）」を提供することで、数値的に不安定な状況や特異なスペクトル構造においても結果の信頼性を保証します。
• 学術的貢献: 左固有スペクトルの多様な構造（球面、不足、多重化など）を数値的に探求するための基盤となり、制御工学、信号処理、幾何学などにおける四元数行列の応用をさらに促進する可能性があります。

著者は、この手法の MATLAB 実装と補足資料（実験プロトコル、追加の例、診断データ）を公開し、研究の再現性と再利用性を担保しています。