Pertanyaan yang diberi tag «eigenvalues»


1
Mengapa SciPy eigsh () menghasilkan nilai eigen yang salah dalam kasus osilator harmonik?
Saya sedang mengembangkan beberapa kode yang lebih besar untuk melakukan perhitungan nilai eigen dari matriks jarang besar, dalam konteks fisika komputasi. Saya menguji rutinitas saya terhadap osilator harmonik sederhana dalam satu dimensi, karena nilai eigen dikenal secara analitis. Melakukannya dan membandingkan rutinitas saya sendiri dengan pemecah bawaan SciPy, saya telah …


1
Nilai eigen terkecil tanpa terbalik
Misalkan adalah matriks pasti simetris, positif. cukup besar sehingga mahal untuk menyelesaikan secara langsung. A A x = bA ∈ Rn × nA∈Rn×nA\in\mathbb{R}^{n\times n}SEBUAHAAA x = bAx=bAx=b Adakah algoritma iteratif untuk menemukan nilai eigen terkecil dari yang tidak melibatkan pembalikan di setiap iterasi?ASEBUAHAASEBUAHAA Artinya, saya harus menggunakan algoritma iteratif seperti …

2
Menguji apakah sebuah matriks positif semi-pasti
Saya memiliki daftar dari matriks simetris yang saya perlu periksa untuk semi-definiteness positif (yaitu nilai eigennya non-negatif.)LL{\cal L} Komentar di atas menyiratkan bahwa seseorang dapat melakukannya dengan menghitung nilai eigen masing-masing dan memeriksa apakah nilainya tidak negatif (mungkin harus mengurus kesalahan pembulatan.) Menghitung nilai eigen cukup mahal dalam skenario saya, …


2
Vektor eigen dari penyesuaian norma kecil
Saya memiliki dataset yang perlahan berubah, dan saya harus melacak vektor eigen / nilai eigen dari matriks kovariansnya. Saya sudah menggunakan scipy.linalg.eigh, tapi terlalu mahal, dan tidak menggunakan fakta bahwa saya sudah memiliki dekomposisi yang hanya sedikit salah. Adakah yang bisa menyarankan pendekatan yang lebih baik untuk menangani masalah ini?

1
Implementasi metode Jacobi-Davidson untuk masalah nilai eigen kubik
Saya memiliki masalah nilai eigen kubik besar: ( A0+ λ A1+ λ2SEBUAH2+ λ3SEBUAH3) x =0.(SEBUAH0+λSEBUAH1+λ2SEBUAH2+λ3SEBUAH3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Aku bisa mengatasi ini dengan mengkonversi ke masalah nilai eigen linear tetapi akan menghasilkan sistem sebagai besar:32323^2 ⎡⎣⎢- A0000saya000saya⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥= λ ⎡⎣⎢SEBUAH1saya0SEBUAH20sayaSEBUAH300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[-SEBUAH0000saya000saya][xyz]=λ[SEBUAH1SEBUAH2SEBUAH3saya000saya0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ …
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.