[論文レビュー] A Probabilistic Proof of the Perron-Frobenius Theorem

Peter W. Glynn, Paritosh Y. Desai|arXiv (Cornell University)|Aug 15, 2018

Markov Chains and Monte Carlo Methods参考文献 9被引用数 4

ひとこと要約

本稿は、マルコフ連鎖のダイナミクスを用いてペロン＝フロベニウス固有値と固有ベクトルを期待値として表現することで、ペロン＝フロベニウスの定理に対する確率的証明を提示する。有限次元および無限次元の設定において、支配的固有ペアの存在と一意性を確立し、古典的定理を連続状態空間へと拡張し、離散化を伴わないモンテカルロ計算を可能にする。

ABSTRACT

The Perron-Frobenius theorem plays an important role in many areas of management science and operations research. This paper provides a probabilistic perspective on the theorem, by discussing a proof that exploits a probabilistic representation of the Perron-Frobenius eigenvalue and eigenvectors in terms of the dynamics of a Markov chain. The proof provides conditions in both the finite-dimensional and infinite-dimensional settings under which the Perron-Frobenius eigenvalue and eigenvectors exist. Furthermore, the probabilistic representations that arise can be used to produce a Monte Carlo algorithm for computing the Perron-Frobenius eigenvalue and eigenvectors that will be explored elsewhere.

研究の動機と目的

• マルコフ連鎖の長期的挙動を用いて、ペロン＝フロベニウス固有値および固有ベクトルの確率的解釈を提供すること。
• ペロン＝フロベニウスの定理を非負の無限次元行列および連続状態空間作用素へと拡張すること。
• これらの一般化された設定において、支配的固有値および正の固有ベクトルが存在するための条件を確立すること。
• 状態空間の離散化を伴わないペロン＝フロベニウス固有ペアを計算するモンテカルロアルゴリズムの基盤を構築すること。
• 特に準定常分布および再生構造の文脈において、定理の確率的側面を統一的かつ明確にすること。

提案手法

• マルコフ連鎖の再生サイクル上の期待値として、ペロン＝フロベニウス固有値および固有ベクトルを表現する。
• 出生死滅過程の遷移核を用いて、準定常挙動をモデル化する非確率的行列 $ B $ を定義する。
• 再生時間 $ \tau $ を導入し、$ B^m $ から導かれる遷移核 $ P_* $ を持つマルコフ連鎖 $ X $ を定義する。これにより一様エルゴディック性が保証される。
• 不変測度 $ \pi_* $ を定義し、それによって一意な正の左固有ベクトル $ u_* $ および右固有ベクトル $ \eta_* $ を特徴付ける。
• $ B^m u_* = \lambda_* u_* $ および $ \eta_* B^m = \lambda_* \eta_* $ を確立し、これにより $ Bu_* = \lambda_*^{1/m} u_* $ および $ \eta_* B = \lambda_*^{1/m} \eta_* $ が得られる。
• 一様エルゴディック性および劣マルティンゲールの議論を活用し、正の固有ベクトルがスケーリングを除き一意であることを証明する。

実験結果

リサーチクエスチョン

• RQ1 どのようにしてマルコフ連鎖のダイナミクスを用いてペロン＝フロベニウス固有値および固有ベクトルを確率的に表現できるか？
• RQ2 無限次元の非負行列において、正の固有値および対応する正の固有ベクトルが存在するための条件は何か？
• RQ3 ペロン＝フロベニウス理論は、連続状態空間の非負作用素へとどのように拡張できるか？
• RQ4 確率的表現は、状態空間の離散化を伴わずに固有ペアを効率的にモンテカルロ計算可能にするか？
• RQ5 なぜ無限次元設定において、固有値 $ \lambda_*^{1/m} $ がペロン＝フロベニウス固有値の正しい一般化となるのか？

主な発見

• ペロン＝フロベニウス固有値 $ \lambda_*^{1/m} $ は $ B^n(x,C) $ の指数的成長率を支配し、$ n \to \infty $ のとき $ B^n(x,C) \sim \lambda_*^{n/m} \frac{u_*(x)\eta_*(C)}{\int_S u_*(w)\eta_*(dw)} $ となる。
• 右固有ベクトル $ u_* $ は、一様エルゴディック性および劣マルティンゲール収束のおかげで、正のスケーリングを除き $ B^m u = \lambda_* u $ の唯一の正の解である。
• 左固有ベクトル $ \eta_* $ は、同じ条件下で、正のスケーリングを除き $ \eta B^m = \lambda_* \eta $ の唯一の正の解である。
• $ P_* $ および $ \tilde{P}_* $ を用いた確率的表現により一様エルゴディック性が保証され、これにより不変測度 $ \pi_* $ への一意性と収束が保証される。
• 条件 A4 が満たされれば、再生連鎖の最小化およびドリフト条件が保証され、正の固有値および正の固有ベクトルの存在が確立される。
• この枠組みは自然に連続状態空間の非負作用素へと拡張され、有限次元ペロン＝フロベニウスの定理の連続的アナロジーを提供する。