Existe uma generalização da Lei de Inércia de Sylvester para o problema simétrico de autovalor generalizado?

9

Eu sei que, para resolver o problema simétrico de autovalor , podemos usar a Lei de Inércia de Sylvester, que é o número de autovalores de A menor que a igual ao número de entradas negativas de D de onde a matriz diagonal D vem da fatoração LDL de um - um I = G D G T . Então, pelo método de bissecção, podemos encontrar todos ou alguns autovalores, conforme desejado. Desejo saber se existe uma generalização da Lei de Inércia de Sylvester para problemas simétricos de autovalores generalizados, que estão resolvendo A x =Ax=λxAaDDAaI=LDLT , onde A e B são matrizes simétricas. Obrigado.Ax=λBxAB

Willowbrook
fonte

Respostas:

5

Sim, se o lápis é definitivo, ou seja, se e B são hermitianos e B é positivo positivo. Em seguida, a assinatura de um - σ B tem a mesma interpretação para o problema de valores próprios ( A - λ B ) x = 0 , como no caso B = I . Um resultado mais geral desse tipo vale para qualquer problema de valor próprio não linear definido A ( λ ) x = 0 . Veja a Seção 5.3 do meu livroABBAσB(AλB)x=0B=IA(λ)x=0

Arnold Neumaier, Introdução à análise numérica, Cambridge Univ. Press, Cambridge 2001.

Para , a prova da minha afirmação pode ser deduzida do argumento dado por Jack Poulson ao observar que C - σ I e A - σ B são congruentes, portanto, têm a mesma inércia.(AλB)x=0CσIAσB

Em particular, pode-se calcular directamente a inércia de , e não necessita de um fatoração de Cholesky de B para formar C . De fato, se B está mal condicionado, a formação numérica de C degrada a qualidade do teste de inércia.AσBBCBC

Arnold Neumaier
fonte
Bom argumento sobre o mau condicionamento de B; Eu acho que sua abordagem é melhor se alguém realmente estiver interessado apenas em calcular a inércia. A abordagem que sugeri é típica para resolver o problema de autovalor (no caso em que está bem condicionado). B
Jack Poulson
@JackPoulson: O teste de inércia é geralmente aplicado para obter os autovalores em um intervalo específico quando e B são escassos e seu padrão de esparsidade articular gera pouco preenchimento. Mas seu C já estará denso quando B for tridiagonal, portanto, use nunca é adequado para encontrar os autovalores de um grande problema de autovalor generalizado esparso. (Considerando que, se o problema não é grande, há pouco ponto em usar a inércia, como encontrar todos os valores próprios é geralmente rápido o suficiente.)ABCB
Arnold Neumaier
Certamente; parece que deixei por engano a palavra "densa" fora do meu comentário.
Jack Poulson
3

BBB=LLH

Ax=LLHxλ,

e essa equação pode ser manipulada para mostrar que

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

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

SSHCAL1LHCACσIσA

Jack Poulson
fonte
Um voto negativo sem nenhuma crítica construtiva?
Jack Poulson
2
Não saí do computador do meu escritório e, por acaso, meu colega de escritório se deparou com essa guia no meu navegador e votou negativamente na resposta, peço desculpas pelo mal-entendido e perguntarei por que ele votou negativamente.
Shuhao Cao
B(A,B)AB
@ Jon: Suspiro. Não é para isso que serve o voto negativo.
22612 Jack Poulson
Eu sei! Eu já disse a ele "leia a regra" depois que descobri que ele usou minha conta para votar uma resposta relevante!
Shuhao Cao