[論文レビュー] An algorithm to find the spectral radius of nonnegative tensors and its convergence analysis

Yuning Yang, Qingzhi Yang+1 more|arXiv (Cornell University)|Feb 14, 2011

Tensor decomposition and applications参考文献 8被引用数 7

ひとこと要約

本稿では、非負行列に用いられるスムージング法を拡張することで、非負テンソルのスペクトル半径を計算する反復的アルゴリズムを提案する。可約でない非負テンソルに対しては、アルゴリズムがスペクトル半径およびその対応する正の固有ベクトルにグローバルに収束することが証明されており、従来の手法に対する信頼性が高く収束性に優れた代替手法を提供する。

ABSTRACT

In this paper we propose an iterative algorithm to find out the spectral radius of nonnegative tensors. This algorithm is an extension of the smoothing method for finding the largest eigenvalue of a nonnegative matrix \cite{s14}. For nonnegative irreducible tensors, we establish the converges of the algorithm. Finally we report some numerical results and conclude this paper with some remarks.

研究の動機と目的

• 非負テンソルのスペクトル半径を計算する信頼性が高く収束性に優れたアルゴリズムの開発。
• 非負行列に用いられるスムージング法を高次テンソルへ拡張すること。
• 可約でない非負テンソルについて理論的収束性の確立。
• 合成的およびランダムなテンソル例におけるアルゴリズムの性能に関する数値的検証の提供。

提案手法

• アルゴリズムは、対角変換の系列と反復的精緻化を用いてスペクトル半径を近似する。
• 高次テンソルの固有値問題における非線形性に対処するためにスムージング技術を適用する。
• 上界 $ R^{(k)} $ と下界 $ r^{(k)} $ の間を交互に計算し、これらがスペクトル半径に収束する。
• 固有ベクトルおよび固有値の推定値を更新するためのベクトル正規化とパワー型反復の利用が鍵となる。
• 収束性は背理法を用いて分析され、境界が収束しないと仮定し、可約でないテンソルのペロン＝フロベニウス性質を用いて矛盾を導出する。
• 収束判定条件として、差分 $ R^{(k)} - r^{(k)} \leq 10^{-7} $ に基づく終了条件を実装する。

実験結果

リサーチクエスチョン

• RQ1 非負行列に用いられるスムージング法を、非負テンソルのスペクトル半径計算に一般化できるか？
• RQ2 提案されたアルゴリズムは、可約でない非負テンソルに対して真のスペクトル半径に収束するか？
• RQ3 次数や次元が異なるテンソルに対して、アルゴリズムの数値的性能はいかがなっているか？
• RQ4 対角変換および正規化の選択が収束挙動に与える影響は何か？
• RQ5 アルゴリズムは可約テンソルに対しても効果的に適用可能か？その場合、どのような条件下で成立するか？

主な発見

• 定理3.2で示されるように、任意の非負で可約でないテンソルについて、アルゴリズムはスペクトル半径 $ \rho(\mathcal{A}) $ に収束する。
• 例1の3次3次元テンソルでは、$ \rho(\mathcal{B}) = 5.79262 $ が計算され、50回目の反復で収束に到達した。
• 表2に示すように、テストされたすべてのランダムテンソル（可約を含む）についても、アルゴリズムはスペクトル半径に収束した。
• 収束速度は $ \mathcal{A}(\alpha) = \mathcal{B} + \alpha\mathcal{I} $ の対角シフトパラメータ $ \alpha $ の選択に影響を受けるが、最適値の特定は今後の課題である。
• 数値結果から、$ R^{(k)} - r^{(k)} \leq 10^{-7} $ および固有方程式の残差誤差が $ 10^{-8} $ 未満であるなど、高い精度が得られ、収束性が強く保証されていることが示された。