Misalkan adalah matriks pasti simetris, positif. cukup besar sehingga mahal untuk menyelesaikan secara langsung. A A x = b
Adakah algoritma iteratif untuk menemukan nilai eigen terkecil dari yang tidak melibatkan pembalikan di setiap iterasi?A
Artinya, saya harus menggunakan algoritma iteratif seperti gradien konjugasi untuk menyelesaikan , jadi berulang kali menerapkan sepertinya seperti "lingkaran dalam" yang mahal. Saya hanya perlu vektor eigen tunggal.A - 1
Terima kasih!
linear-algebra
eigensystem
eigenvalues
iterative-method
Justin Solomon
sumber
sumber
eigs
-routine. Ini adalah metode berulang. Ada beberapa opsi untuk menentukan nilai eigen yang Anda inginkan, mis . Real terkecil .Jawaban:
Hitung nilai eigen dengan magnitudo terbesar dari A (dengan, katakanlah, ).λmax A
eigs('lm')
Kemudian menghitung besarnya terbesar (negatif) nilai eigen λ m a x dari M = A - λ m a x I (sekali lagi, melalui panggilan standar untuk ).λ^max M=A−λmaxI
eigs('lm')
Cari eigenvector Anda oleh pemecahan .( A - λ m i n I ) v = 0v (A−λminI)v=0
sumber