Gibt es eine Verallgemeinerung des Sylvester-Trägheitsgesetzes für das symmetrische verallgemeinerte Eigenwertproblem?

9

Ich weiß, dass wir zur Lösung des symmetrischen Eigenwertproblems das Sylvester-Trägheitsgesetz verwenden können, dh die Anzahl der Eigenwerte von A kleiner als a entspricht der Anzahl der negativen Einträge von D, wobei die Diagonalmatrix D von der stammt LDL - Faktorisierung von A - a I = L D L T . Dann können wir durch die Halbierungsmethode alle oder einige Eigenwerte nach Wunsch finden. Ich möchte wissen, ob es eine Verallgemeinerung des Sylvester-Trägheitsgesetzes für symmetrische verallgemeinerte Eigenwertprobleme gibt, nämlich die Lösung von A x =Ax=λxAaDDAaI=LDLT , wobei A und B symmetrische Matrizen sind. Vielen Dank.Ax=λBxAB

Willowbrook
quelle

Antworten:

5

Ja, wenn der Bleistift eindeutig ist, dh wenn und B hermitisch sind und B eindeutig positiv ist. Dann wird die Signatur von A - σ B hat die gleiche Auslegung für das Eigenwertproblem ( A - λ B ) x = 0 , wie im Fall B = I . Ein allgemeineres Ergebnis dieser Art gilt für jedes bestimmte nichtlineare Eigenwertproblem A ( λ ) x = 0 . Siehe Abschnitt 5.3 meines BuchesABBAσB(AλB)x=0B=IA(λ)x=0

Arnold Neumaier, Einführung in die numerische Analyse, Cambridge Univ. Press, Cambridge 2001.

Für kann der Beweis meiner Behauptung aus dem Argument von Jack Poulson abgeleitet werden, dass C - σ I und A - σ B kongruent sind und daher die gleiche Trägheit haben.(AλB)x=0CσIAσB

Insbesondere kann man die Trägheit von direkt berechnen und benötigt keine Cholesky-Faktorisierung von B , um C zu bilden . Wenn B schlecht konditioniert ist, verschlechtert die numerische Bildung von C die Qualität des Trägheitstests.AσBBCBC

Arnold Neumaier
quelle
Guter Punkt über die schlechte Konditionierung von B; Ich denke, dass Ihr Ansatz besser ist, wenn man wirklich nur daran interessiert ist, die Trägheit zu berechnen. Der von mir vorgeschlagene Ansatz ist typisch für die tatsächliche Lösung des Eigenwertproblems (in dem Fall, in dem gut konditioniert ist). B
Jack Poulson
@JackPoulson: Der Trägheitstest wird normalerweise angewendet, um die Eigenwerte in einem bestimmten Intervall zu erhalten, wenn und B spärlich sind und ihr gemeinsames Sparsity-Muster nicht zu viel Füllung erzeugt. Ihr C ist jedoch bereits dicht, wenn B tridiagonal ist, und verwendet es daher ist niemals geeignet, um die Eigenwerte eines großen, spärlich verallgemeinerten Eigenwertproblems zu finden. (Wenn das Problem nicht groß ist, macht es wenig Sinn, Trägheit zu verwenden, da das Finden aller Eigenwerte normalerweise schnell genug ist.)ABCB
Arnold Neumaier
Bestimmt; es scheint, dass ich fälschlicherweise das Wort "dicht" aus meinem Kommentar herausgelassen habe.
Jack Poulson
3

BBB=LLH

Ax=LLHxλ,

und diese Gleichung kann manipuliert werden, um dies zu zeigen

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

CL1ALHA(A,B)CC(A,B)

SSHCAL1LHCACσIσA

Jack Poulson
quelle
Eine Ablehnung ohne konstruktive Kritik?
Jack Poulson
2
Ich habe mich nicht am Computer meines Büros abgemeldet, und mein Amtskollege ist zufällig auf diese Registerkarte in meinem Browser gestoßen und hat die Antwort abgelehnt. Ich entschuldige mich für das Missverständnis und werde ihn fragen, warum er dies abgelehnt hat.
Shuhao Cao
B(A,B)AB
@ Jon: Seufz. Dafür ist die Abstimmung nicht gedacht.
Jack Poulson
Ich kenne! Ich sagte ihm bereits "Bitte lies die Regel", nachdem ich festgestellt hatte, dass er mein Konto verwendet hatte, um eine relevante Antwort abzustimmen!
Shuhao Cao