[論文レビュー] Eigendecomposition of Block Tridiagonal Matrices

Aliaksei Sandryhaila, José M. F. Moura|arXiv (Cornell University)|Jun 2, 2013

Matrix Theory and Algorithms参考文献 41被引用数 11

ひとこと要約

本稿では、固有値を行列多項式の零点に結びつけることで、ブロック三重対角行列に対する閉形式の固有ベクトル表現を提示する。この手法により、固有ベクトルの計算がより小さい $K\times K$ 行列の零空間計算に還元され、高速な固有ベクトル展開アルゴリズムの実現が可能となる。具体的には、'スパイダー'グラフのスペクトル解析において、$O(N\log^2 L)$ の計算量を達成し、直接計算の $O(N^2)$ より顕著に改善される。

ABSTRACT

Block tridiagonal matrices arise in applied mathematics, physics, and signal processing. Many applications require knowledge of eigenvalues and eigenvectors of block tridiagonal matrices, which can be prohibitively expensive for large matrix sizes. In this paper, we address the problem of the eigendecomposition of block tridiagonal matrices by studying a connection between their eigenvalues and zeros of appropriate matrix polynomials. We use this connection with matrix polynomials to derive a closed-form expression for the eigenvectors of block tridiagonal matrices, which eliminates the need for their direct calculation and can lead to a faster calculation of eigenvalues. We also demonstrate with an example that our work can lead to fast algorithms for the eigenvector expansion for block tridiagonal matrices.

研究の動機と目的

• 信号処理、物理学、データサイエンスの分野において、大規模ブロック三重対角行列の固有値分解の高コストを軽減すること。
• $N\times N$ 行列演算を直接行わずに、行列多項式の構造を活用して固有値と固有ベクトルを計算するフレームワークの構築。
• $K \ll N$ である $K\times K$ 行列の零空間に基づく、固有ベクトルの閉形式表現の確立。
• 再帰的構造と分割統治技法を活用して、固有ベクトル展開の高速化を実現すること。
• 「スパイダー」グラフと呼ばれるブロック三重対角隣接行列のクラスにおける、本手法の効率性をスペクトル解析を通じて実証すること。

提案手法

• ブロック三重対角行列の構造に由来する三項再帰を満たす $K\times K$ 行列多項式 $\mathbf{P}_n(x)$ の族を定義する。
• ブロック三重対角行列の固有値が $\det \mathbf{P}_L(x)$ の零点に一致することを証明し、スペクトル特性を多項式の根に結びつける。
• 固有ベクトルを、$\mathbf{P}_L(\lambda)\mathbf{v} = \mathbf{0}$ の解であるブロックベクトルとして構成し、$K\times K$ 行列の零空間計算に問題を還元する。
• 再帰的因数分解と構造的行列分解（例：DST-1、近隣変換）を用いて、固有ベクトル基底における高速な行列-ベクトル積を実現する。
• FFT や DCT アルゴリズムにインspiredされた分割統治戦略を適用し、固有ベクトル展開の計算量を $O(N^2)$ から $O(N\log^2 L)$ に低減する。
• 隣接行列がブロック三重対角形式をとる「スパイダー」グラフに対してフレームワークを検証し、導出された固有ベクトル構造を用いた高速変換を導出する。

実験結果

リサーチクエスチョン

• RQ1 一般のブロック三重対角行列の固有値は、どのように行列多項式を用いて特徴付けられるか？
• RQ2 大規模な $N\times N$ システムを解かずに、ブロック三重対角行列の固有ベクトルを閉形式で表現できるか？
• RQ3 ブロック三重対角行列の固有ベクトル展開の計算量は何か？ また、これを高速化できるか？
• RQ4 本手法は、スケーラビリティと数値的効率性の観点から、標準的固有値分解と比較してどのように差を示すか？
• RQ5 本フレームワークは、「スパイダー」グラフのような構造的グラフに応用可能か？ これにより高速スペクトル変換が可能になるか？

主な発見

• ブロック三重対角行列の固有値は、行列のブロックを含む再帰的構成によって得られる行列多項式 $\mathbf{P}_L(x)$ の行列式の零点に一致する。
• 固有ベクトルは、各固有値 $\lambda$ において $\mathbf{P}_L(\lambda)$ の零空間に属するブロックをもつベクトルとして閉形式で表現され、計算が $K\times K$ 零空間問題に還元される。
• 「スパイダー」グラフの隣接行列に対する固有ベクトル展開は、$O(N\log^2 L)$ の演算で計算可能であり、標準的な $O(N^2)$ のコストと比べ顕著に改善される。
• DST-1 や近隣変換などの構造的行列分解を活用することで、本手法は $O(N\log^2 L)$ の計算量スケーリングを実現し、高速アルゴリズムを可能にする。
• 本フレームワークは、対称行列に限らない非対称ブロック三重対角行列へも直交多項式法を一般化し、適用範囲を拡大する。
• $K=2$ の場合、本手法は標準的な DCT-1 変換に還元され、$O(N\log N)$ の計算量を示し、既知の高速変換理論と整合する。