[論文レビュー] The Lyapunov matrix equation. Matrix analysis from a computational perspective

Valeria Simoncini|arXiv (Cornell University)|Jan 29, 2015

Matrix Theory and Algorithms参考文献 22被引用数 8

ひとこと要約

本稿は、Lyapunov行列方程式 $AX + XA^T = D$ の解における崩壊（decay）およびスパarsity（疎性）のパターンを調査し、$D$ がスパースである場合、解 $X$ が $D$ の非ゼロ要素の周囲に構造的な崩壊パターンを有することを示している。著者らは、これらのパターンを活用して、$D$ が低ランクでない場合にも効率的かつスパースな近似解を得る手法を提案し、行列構造と崩壊行動に関する事前知識に基づく新たな計算戦略を提示する。

ABSTRACT

Decay properties of the solution $X$ to the Lyapunov matrix equation $AX + X A^T = D$ are investigated. Their exploitation in the understanding of equation matrix properties, and in the development of new numerical solution strategies when $D$ is not low rank but possibly sparse is also briefly discussed.

研究の動機と目的

• スパースであるが必ずしも低ランクでない $D$ を持つLyapunov方程式 $AX + XA^T = D$ の解行列 $X$ における崩壊およびスパarsityパターンを理解・特徴づけること。
• これらの崩壊特性が、大規模なLyapunov方程式を解くための新しい計算戦略の設計にどのように寄与するかを調査すること。
• 従来、逆行列や帯行列構造に対して知られていた古典的な崩壊結果を、関数に基づく行列解析を用いてLyapunov方程式の解へと拡張すること。
• $D$ がフルランクまたは高ランクであっても、$D$ がスパースである場合に限って、$X$ が崩壊構造を有することを示し、その構造を活用してスパース行列として効率的に近似できることを示すこと。
• 収束速度と $X$ のスパarsityおよび崩壊パターンとの関連を明確にすることで、投影法およびKrylov型反復法の使用を促進すること。

提案手法

• 本稿は、Kronecker積表現を用いてLyapunov方程式を大規模な線形系 $\mathcal{A}\mathbf{x} = \mathbf{b}$ として表現し、スペクトル理論および行列関数の性質を用いた解析を可能としている。
• 特に、解析的関数 $f$ に対する $f(A)$ の成分に適用可能な境界を用いて、解行列 $X$ の崩壊推定値を導出している。
• 解 $X$ を右辺 $D$ の関数として分析し、$D$ の各非ゼロ要素が $X$ に局所的かつ崩壊する成分を寄与させることを示しており、その崩壊パターンは $A$ の固有値特性に依存する。
• 対称正定値帯行列 $A$ の場合、古典的な逆行列崩壊境界（例：Demkoら）をLyapunov解へと一般化し、$X$ における指数的非対角成分の崩壊を示している。
• 非正規だが対角化可能な行列に対しても、[8] の結果を応用することで、対称でない場合への適用範囲を拡大している。
• 数値実験を通じて、特に $D$ が対角行列の場合に予測された崩壊パターンが妥当であることを検証しており、$D$ の非ゼロ対角成分に対応する位置に局所的ピークが現れることが確認されている。

実験結果

リサーチクエスチョン

• RQ1 Lyapunov方程式 $AX + XA^T = D$ の解 $X$ における崩壊およびスパarsityパターンは、$D$ のスパarsity構造にどのように依存するか？
• RQ2 $A$ の固有値特性と $D$ の非ゼロ要素の位置に基づいて、事前に $X$ の崩壊行動を予測できるか？
• RQ3 $D$ が低ランクでない場合に、$X$ のスパarsityおよび崩壊構造をどれだけ活用して、効率的かつスパースな近似解法を設計できるか？
• RQ4 Krylov部分空間法の収束速度は、$X$ のスパarsityおよび崩壊パターンとどのように関連するか？
• RQ5 Sylvester方程式 $A_1X + XA_2 = D$ の解に対しても、同様に崩壊に基づく戦略を用いて解析および近似可能か？

主な発見

• 対称正定値帯行列 $A$ の場合、Lyapunov方程式の解 $X$ は指数的非対角成分の崩壊を示し、その崩壊率は $A$ の条件数およびバンド幅に依存する。これは逆行列の崩壊と類似した性質である。
• 対角行列 $D$ の場合、解 $X$ は $D$ の非ゼロ対角成分に対応する位置に局所的ピークを示し、それらの周囲に向かって崩壊する。これは解の加法的構造を裏付ける。
• $D$ の各非ゼロ要素に対応する成分 $X_j$ は、$(j,j)$ 位置を中心として崩壊パターンを有しており、スパース近似のための切り捨てが可能となる。
• Krylov部分空間法によるLyapunov方程式の解法の収束速度は、$X$ のスパarsityおよび崩壊構造と関連しており、両者とも同じ多項式近似特性に支配されている。
• $X$ の崩壊パターンにより、$D$ がフルランクであっても、$D$ がスパースである限り、ある崩壊閾値を超える成分を切り捨てることで、高精度なスパース近似が可能となる。
• 関数計算を用いることで、非正規だが対角化可能な行列に対しても理論的崩壊境界を拡張でき、このフレームワークの適用範囲を対称な場合に限らないように拡大できる。