Saya memiliki masalah nilai eigen kubik besar:
Aku bisa mengatasi ini dengan mengkonversi ke masalah nilai eigen linear tetapi akan menghasilkan sistem sebagai besar:
di mana dan z = λ y . Apa teknik lain yang tersedia untuk memecahkan masalah nilai eigen kubik? Saya pernah mendengar bahwa ada versi Jacobi-Davidson yang akan menyelesaikannya tetapi belum menemukan implementasi.
Juga, saya harus dapat menargetkan nilai eigen spesifik yang serupa dengan metode shift-and-invert ARPACK dan menemukan vektor eigen terkait.
Jawaban:
Dengan protokol komunikasi terbalik ARPACK, Anda tidak perlu menyimpan matriks secara eksplisit: Anda hanya perlu menyediakan dua fungsi yang menghitung:3 n × 3 n
dan [ x y z ] → [ A 1 x + A 2 y + A 3 z y z ]⎡⎣⎢xyz⎤⎦⎥→ ⎡⎣⎢- A0xyz⎤⎦⎥ ⎡⎣⎢xyz⎤⎦⎥→ ⎡⎣⎢SEBUAH1x + A2y+ A3zyz⎤⎦⎥
Mengenai transformasi invert, Anda dapat melakukan hal yang sama, yaitu mengimplementasikannya sendiri dengan menggunakan panggilan balik yang menghitungx ↦ M- 1x x ↦ Mx λ′s λ- 1 M.- 1x M. SEBUAH0
sumber