2013-01-22 31 views
5

実際の正方形の非負行列のPerron-Frobenius固有値を推定する最適化パッケージまたは方法はありますか?これは、Peron-Frobenius固有値が行列を反復することによって到達できるとすれば、正確な計算よりもはるかに高速である(特に大行列および/または希薄行列の場合)。私はこれを行う最適化されたパッケージが存在することを願っています。実数、正方行列、負でない行列のPerron-Frobenius固有値を推定するためのPythonパッケージ

答えて

4

scipy.sparse.linalgにはARPACK libraryを使用するeigseigshの機能があります。あなたはthis tutorialで続きを読むことができますが、aはおそらくスパース形式で、正方行列である場合、あなたはその最大の大きさの固有値を得ることができ、そのペロン・フロベニウス固有値すなわち、および対応する固有ベクトルとして:

val, vec = scipy.sparse.linalg.eigs(a, k=1, which='LM') 
+0

感謝!これはまさに私が望んでいたものです。 'which'フラグを忘れました。また、これは疎なパッケージに入れるのが面白い(そして素晴らしい)。ダブルポイント。 – gabe

関連する問題