与えられたベクトルの固有ベクトル成分を計算する


8

エルミートスパース演算子Mの固有空間に分解できるベクトルVあります。M

V=ivim^i

最大の(大きさ)に対応する(固有ベクトル自体)を見つける方法はありますか?m^ivi

基本的に、事前にわからないの固有ベクトルを含む、合計の最大の数項が必要です。M

具体的には、最大の対応するの固有ベクトルを同時に見つけたい 、最大のを見つける。最初にスペクトル全体を見つけることはできません。M|vi|viM

私が考えていたいくつかの可能性:

「Wieldant's Deflation」の反対を使用して、マトリックスを「膨張」させることができます。

M1=M+σ[Σivim^i]VH=M+σVVH

異なるの固有値は、シフトされます。固有ベクトルは変化しないため、と抽出できると思います。問題は、外積が密であることです。m^iλi+σ|vi|2σviV

別の可能性:

べき乗法(収束までにベクトル掛け続ける)は、最大の固有値を持つの成分を見つけます。この方法の欠点は、の大きさを制御できないため、すべてのコンポーネントを見つけ、最大のものを見つけることです。MVVvi

最大のコンポーネントのみに収束するようにこれを制御する方法はありますか?


私は現実世界を考えているので、確率的シミュレーションが必要だと考えています。ベクトルがあり、部分空間でその場所を見つけたいとします。部分空間は固有ベクトルによって広がります。秘訣は、固有ベクトルがすべてではなく、一部だけであるということです。ドット積を使用して、固有ベクトルの1つまたは複数に寄与する固有ベクトルを確認できます。運が良ければ、固有ベクトルのサブセットでベクトルを広げることができます。そうでない場合、あなたはあなたが持っている固有ベクトルに関してあなたのベクトルを分解し、あなたが持っていないベクトルを見つけることができます。
EngrStudent 2015年

解決策については、scicomp.stackexchange.com / questions / 28111 /…をご覧ください。解決策の詳細を質問に追加しました。
as2457

回答:


1

行列はエルミート行列であるため、これをハミルトニアンとして使用して、想像上の時間で伝播させることができます。つまり、次の微分方程式系を解きます。

idVdt=MV

これに対する一般的な解決策は次のとおりです。

V(t)=V0eiMt

次に、をフーリエ変換します。ピークの高さと配置により、さまざまな固有ベクトルに沿った成分とそれに関連する固有値がわかります。これは、超高速原子物理学では「スペクトル法」と呼ばれることもあります。V(t)V(0)

固有値を取得したら、任意の固有固有値ソルバーを使用して固有ベクトルを見つけます。


教科書、講義ノート、またはこの「スペクトル法」の紹介が含まれているソースをお勧めできますか グーグル検索は本当に私に満足のいく結果を提供しませんでした。
Marco Breitig、2015

@Marco Breitig:どちらも見つけることができませんでした。それは、原子物理学の口承伝統のほんの一部だと思います。
Dan

0

承知しました。ましょう、各列はエルミート演算子の別個の固有ベクトルであるような行列である。の隠蔽性から、は可逆であることがわかります。を取得するには、線形システム解きます。次に、ようなインデックス選択します; の番目の列は、最大のに対応する固有ベクトルになります。M^MMM^viM^v=Vimaxj|vj|=|vi|iM^vi


1
すみません、よくわかりませんでした。私はすべての固有ベクトルを計算してその方法で答えを得ることができることを知っていますが、最大成分に対応する固有ベクトルが必要なだけです。可能であれば、すべての固有ベクトルの計算を避けて、どの成分が最大かを調べたいと思います。V
Andrew Spott 14

2
@AndrewSpott:その要件を質問に追加してください。固有値を与えられた固有ベクトルを見つける(または最大の大きさの固有値を見つける)のは簡単です。与えられた基底の観点から最高の低ランクの近似を見つけることも簡単です。ただし、ここでは、固有基底を解き、同時に最高の低ランク近似を見つける必要があります。それは可能かもしれませんが、すぐに頭に浮かぶアルゴリズムはありません。
Geoff Oxberry

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