行列の


8

与えられた大行列Aと固有値σ1≥σ2≥…、私は、これらの値のサブセットのみを決定言いたいσ5,σ8およびσ19。これを実行できるアルゴリズムはありますか、または上位19の固有値を実行できる最高のものを見つけていますか?


マトリックスには特別な特性がありますか?あなたはおおよその値を知っています、σ 8およびσ 19を?σ5σ8σ19
— nicoguaro

いいえ。ただし、特別な特性を想定した結果が出た場合は、遠慮なくお知らせください。
— dexter04

関連:stackoverflow.com/questions/12167654/…おそらく、その裏でメソッドを決定できるでしょう。
— carlosvalderrama

クリロフ部分空間、アーノルディ、ランチョスの反復が必要な方向に何かあると思います。
— carlosvalderrama

2
行列がエルミート行列である場合、固有値のすべてが実数であるため、順序付けられた固有値のセットがあることが保証されます。また、おおよその値がわかっている場合は、shift-invertメソッドを使用してみてください。結局のところ、最初の固有値を見つける方が簡単かもしれません。n
— nicoguaro

回答:


6

いいえ、私が知る限り、これらの固有値の場所をおおよそ知らない限り、何もありません。行列のスペクトルのサブセットを計算できるメソッドについては、以下を生成できるメソッドのみを知っています。

  • 「スペクトルの極値」からの固有値。たとえば、絶対値が最大のもの、または最も負の実数/虚数部分を持つもの。たとえば、Arnoldi /部分空間の反復(適切にシフト)。
  • 複素平面の特定の点に最も近い固有値。たとえば、係数の最小値(ゼロに最も近い)。たとえば、(有理)シフトおよび反転Arnoldiまたは逆べき乗反復。
  • 例えば、飛行機の特定の領域における固有値、[−1,1]×[−1,1]。FEASTタイプのメソッドはこれを行うと思いますが、それらについてはほとんど知りません。

私が理解していることから、FEASTタイプのメソッドを使用すると、(等高線積分によって)複雑な平面の指定された領域の固有値の数を取得することもできるため、スペクトルの中央から指定された固有値が必要な場合あなたの行列のスペクトルがであること(シフトおよびスケーリング)を前提としています、例えば300分のために、あなたは二分検索のソートを実行することができ[−1,1](それは簡単にするため、実際のです)。1に最も近い固有値μ計算する12。[1の固有値の数を計算します[12,1]。300より大きい場合は、[μ,1]中点に最も近い固有値を探し、それ以外の場合は[−1,μ]の中点に最も近い固有値を探します。繰り返す。

LAあなたの例のように、300を19に置き換えた場合、これらのどれもがバニラアーノルディと競合することはないと思います。そして、たぶん300でさえも。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.