対称一般化固有値問題に対してシルベスター慣性法の一般化はありますか?


9

私は、対称固有値問題を解決するためにことを知っている、私たちはの固有値の数であるシルベスターの慣性法則、使用することができますAを未満の負のエントリ数と等しいD対角行列Dはから来ているのLDLの分解A - I = L D L T。次に、二分法により、必要に応じてすべてまたは一部の固有値を見つけることができます。私は解決され、対称一般化固有値問題のシルベスター慣性法の一般化、そこに存在するかどうかを知りたい、Ax=λxAaDDAaI=LDLTA BAx=λBxAとは対称行列です。ありがとう。B

回答:


5

はい、鉛筆が明確である場合、つまりとがエルミートで、が正定である場合。次に、のシグネチャは、固有値問題、場合と同じ解釈になります。この種のより一般的な結果は、明確な非線形固有値問題当てはまり。私の本のセクション5.3を参照B B A - σ B A - λ B X = 0 B = I A λ のx = 0ABBAσB(AλB)x=0B=IA(λ)x=0

アーノルドノイマイヤー、数値解析入門、ケンブリッジ大学 プレス、ケンブリッジ2001。

、私の主張の証拠があることを指摘すると、ジャックPoulsonによって与えられた引数から推論することができるおよびと同じ慣性を持っているため、一致しています。C - σ I A - σ B(AλB)x=0CσIAσB

特に、慣性を直接計算でき、を形成するためにコレスキー分解を必要としません。実際、が悪条件である場合、の数値形成は慣性テストの品質を低下させます。B C B CAσBBCBC


Bの悪条件についての良い点; 慣性の計算にのみ関心がある場合は、あなたのアプローチの方が良いと思います。私が提案したアプローチは、固有値問題を実際に解決するための典型的なものです(が条件が整っている場合)。B
Jack Poulson 2012

@JackPoulson:慣性テストは、通常時に特定の間隔で固有値を取得するために適用されとまばらで、その共同スパースパターンがであまり塗りつぶしをしませ生成しかし、あなたの。ときすでに密になり三重対角である、したがって、それを使用して大規模なスパース一般化固有値問題の固有値を見つけるのに適していることは決してありません。(問題が大きくない場合、すべての固有値を見つけるのは通常十分高速であるため、慣性を使用してもほとんど意味がありません。)B C BABCB
Arnold Neumaier 2012

もちろん; 私はコメントから誤って "dense"という単語を省いたようです。
Jack Poulson、2012

3

場合エルミート正定値、のコレスキー分解である、言う、与えることBBB=LLH

Ax=LLHxλ,

この方程式を操作して、

(L1ALH)(LHx)=(LHx)λ,

ここで、はの対称性を維持し、鉛筆と同じスペクトルを持つことは明らかです。したがって、を形成した後、コレスキー分解とそれに続く両側三角ソルブを使用すると、シルベスターの慣性則をに直接適用して、鉛筆固有値に関する情報を収集できます。CL1ALHA(A,B)CC(A,B)

シルベスターの慣性の法則は合同変換、たとえばに関して不変なので、行列は変換を通じて合同であることに注意してください、したがってはと同じ慣性を持っています。ただし、慣性が必要な場合、一部の非ゼロシフト場合、単純に考慮ことはできません。 C A L -SSHCA CAC-σIσAL1LHCACσIσA


建設的な批判のない反対票?
Jack Poulson 2012

2
私はオフィスのコンピューターからログアウトしていません。オフィスメイトがブラウザーのこのタブにたまたま出て、回答に反対票を投じました。誤解をお詫びし、なぜこれに反対票を投じたのか尋ねます。
Shuhao Cao 2012

B(A,B)AB

@ジョン:ため息。それは反対投票の目的ではありません。
Jack Poulson 2012

知っている!私は彼が私のアカウントを使用して関連する回答に反対票を投じたことに気付いた後、私はすでに「ルールを読んでください」と彼に言いました!
Shuhao Cao
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.