Apakah ada generalisasi Hukum Sylvester Inersia untuk masalah nilai eigen simetris yang digeneralisasi?


9

Saya tahu bahwa untuk memecahkan simetris eigen masalah , kita dapat menggunakan Sylvester Inersia Hukum, yaitu jumlah nilai eigen dari A kurang dari satu sama dengan jumlah entri negatif dari D mana matriks diagonal D berasal dari LDL faktorisasi A - sebuah I = L D L T . Kemudian, dengan metode pembagian dua, kita dapat menemukan semua atau beberapa nilai eigen yang diinginkan. Saya ingin tahu apakah ada generalisasi Hukum Sylvester Inersia untuk masalah nilai eigen umum simetris, yaitu menyelesaikan A x =Ax=λxAaDDAaI=LDLT , di mana A dan B adalah matriks simetris. Terima kasih.Ax=λBxAB

Jawaban:


5

Ya, jika pensilnya pasti, yaitu, jika dan B adalah Hermitian dan B adalah pasti positif. Maka tanda tangan dari A - σ B memiliki penafsiran yang sama untuk masalah nilai eigen ( A - λ B ) x = 0 seperti dalam kasus B = I . Hasil yang lebih umum dari jenis ini berlaku untuk masalah nilai eigen nonlinier yang pasti A ( λ ) x = 0 . Lihat Bagian 5.3 buku sayaABBAσB(AλB)x=0B=IA(λ)x=0

Arnold Neumaier, Pengantar analisis numerik, Cambridge Univ. Pers, Cambridge 2001.

Untuk , bukti pernyataan saya dapat disimpulkan dari argumen yang diberikan oleh Jack Poulson setelah mencatat bahwa C - σ I dan A - σ B adalah kongruen, maka memiliki inersia yang sama.(AλB)x=0CσIAσB

Secara khusus, satu langsung dapat menghitung inersia , dan tidak membutuhkan faktorisasi Cholesky dari B untuk membentuk C . Memang, jika B tidak dikondisikan maka pembentukan numerik C menurunkan kualitas tes inersia.AσBBCBC


Poin bagus tentang pengondisian B; Saya pikir pendekatan Anda lebih baik jika seseorang benar-benar hanya tertarik menghitung inersia. Pendekatan yang saya sarankan khas untuk benar-benar memecahkan masalah nilai eigen (dalam kasus di mana dikondisikan dengan baik). B
Jack Poulson

@JackPoulson: Tes inersia biasanya diterapkan untuk mendapatkan nilai eigen dalam interval tertentu ketika dan B jarang dan pola sparitas gabungannya menghasilkan pengisian yang tidak terlalu banyak. Tetapi C Anda akan menjadi padat ketika B adalah tridiagonal, karenanya menggunakannya tidak pernah cocok untuk menemukan nilai eigen dari masalah nilai eigen umum jarang yang besar. (Sedangkan jika masalahnya tidak besar, ada sedikit gunanya menggunakan inersia, karena menemukan semua nilai eigen biasanya cukup cepat.)ABCB
Arnold Neumaier

Pasti; sepertinya saya keliru meninggalkan kata "padat" dari komentar saya.
Jack Poulson

3

BBB=LLH

Ax=LLHxλ,

dan persamaan ini dapat dimanipulasi untuk menunjukkan hal itu

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

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

SSHCAL1LHCACσIσA


Sebuah downvote tanpa kritik yang membangun?
Jack Poulson

2
Saya belum logout di komputer kantor saya, dan teman kantor saya kebetulan menemukan tab ini di browser saya dan menurunkan jawabannya, saya minta maaf atas kesalahpahaman dan akan bertanya kepadanya mengapa dia menurunkan ini.
Shuhao Cao

B(A,B)AB

@ Jon: Huh. Bukan itu tujuan downvote.
Jack Poulson

Aku tahu! Saya sudah mengatakan kepadanya "tolong baca aturannya" setelah saya menemukan bahwa dia menggunakan akun saya untuk menurunkan jawaban yang relevan!
Shuhao Cao
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.