Convergence Analysis of Natural Power Method and Its Applications to Control

技術要約：ナチュラル・パワー法の収束解析とその制御への応用

問題提起
オンライン主成分アルゴリズムは、支配的な部分空間を抽出するために広く用いられている。連続時間におけるOjaフローおよびその離散時間オイラー近似は十分に研究されているが、これらは実部が最大の固有値を対象としている。対照的に、離散時間のナチュラル・パワー法（NPM）は、$U[k + 1] = AU[k](U[k]^\top A^\top AU[k])^{-1/2}$ で定義され、主に対称正定値行列に対して、絶対値が最大の固有値を抽出するものとして分析されてきた。一般の正方行列（非対称、不定値）に対するNPMの理論的理解、特にその収束特性と、支配的なモードが絶対値の最大の固有値によって定義される離散時間制御システムへの適用可能性については、理論的な空白が存在する。

手法
著者らは、一般の正方行列 $A \in \mathbb{R}^{n \times n}$ に対するNPMの厳密な収束解析を行う。

1. 収束解析： 本論文は、非対称行列に対して既存の結果（Hua et al., 1999）を一般化している。著者らはスティフェル多様体 $St(r, n)$を定義し、アルゴリズムの順方向不変集合を解析する。初期値が特定の集合$V^{(d)}{m,r}$内にあり、かつ絶対値にスペクトルギャップ（$|\lambda_m| > |\lambda{m+1}|$）が存在する場合、数列$U[k]$が、絶対値の大きい上位$r$個の固有値に対応する不変集合$U^{(d)}_{m,r}$ に収束することを証明する。
2. 定常性と修正： 標準的なNPMの反復は、一般に単一の定常点（不変集合の要素）に収束せず、集合内を振動することを示す。これに対処するため、著者らは追加の直交行列の乗算を含む修正アルゴリズム（式9）を提案し、不変集合の特定の要素への収束を保証する。
3. 次元削減： すでに計算された支配的部分空間から低次元の部分空間を効率的に抽出するために、低次化されたNPMを提案しており、これによりフルサイズの $n \times n$ 行列に対する繰り返しの演算を回避する。
4. 制御への応用： 理論的結果は、離散時間線形時不変（LTI）システムに対するモデル次数低減（MOR）、および線形時変（LTV）システムに対する低ランクコントローラ合成に適用される。著者らは、可観測性、可到達性、および安定性を保持する低次モデルを構築するためにNPMを利用する。

主な貢献と結果

• 収束特性の解明： 本論文は、NPMが絶対値の大きい上位 $r$ 個の固有値に対応する $r$-支配的部分空間を抽出することを確立する。これは、実部の最大値を対象とするOjaフローとは明確に異なる特性である。
• 一般行列への収束： 収束は、対称正定値行列だけでなく、一般の正方行列に対して証明されている。解析によれば、固有値の絶対値にスペクトルギャップが存在する場合、ほとんどすべての初期行列が不変集合に収束する。
• システム特性の保持： 離散時間LTIシステムにおいて、提案されたMORスキームは、元のシステムの最小性（可観測性と可到達性）および安定性を保持する。低次モデルは、支配的部分空間上に投影された元のシステムと同じ伝達関数を共有する。
• 特異摂動とLTVへの拡張： 本手法は、スローモードとファストモードを抽出することにより、特異摂動システムにも適用可能であることを示す。さらに、NPMの再帰的な性質により、時変の支配的部分空間を追跡するための、緩やかに変化するLTVシステムへの適用が可能となる。
• 数値検証： パラメータ化された行列 $A_\alpha$ を用いたシミュレーションにより、収束速度は固有値の絶対値の比（$|\lambda_{r+1}|/|\lambda_r|$）に依存することが示される。本論文では、この比が1に近い場合（例：$\alpha = 0.9$）、収束が遅くなり、それが派生したコントローラの追従性能に影響を与えることが指摘されている。

意義と主張
著者らは、本研究が一般の正方行列の文脈におけるNPMの基礎的な収束解析を提供し、Ojaフローと明確に区別するものであると主張している。主な意義は、支配的なダイナミクスが実部ではなく固有値の大きさによって支配される離散システムにおいて、この理論的空白を埋めることで、制御理論への実用的な適用を可能にすることにある。

著者らは、自らの貢献を以下のように控えめに位置づけている：

1. 一般行列に対するNPMの標準的な収束ベンチマークの提供。
2. モデル次数低減および低ランクコントローラ合成におけるNPMの有用性の実証。
3. 時変部分空間を追跡する潜在能力の提示。ただし、その性能は基礎となるNPMアルゴリズムの収束速度によって制限されることを認めている。

また、著者らは、誤差界の厳密な解析や、時変制御および推定問題へのより広範な適用にはさらなる研究が必要であると明記しており、現在のLTVシステム用コントローラ合成は、NPMの収束が十分に速い場合にのみ良好に機能する初期段階のステップであると述べている。