Pertanyaan yang diberi tag «svd»

Apa itu algoritma sederhana untuk menghitung SVD matriks?2 × 22×22 \times 2 Idealnya, saya ingin algoritma yang kuat secara numerik, tetapi saya ingin melihat implementasi yang sederhana dan yang tidak terlalu sederhana. Kode C diterima. Adakah referensi untuk makalah atau kode?

26 linear-algebra  matrix  reference-request  svd  matrix-factorization 

Saya sedang mengerjakan perpustakaan matriks header-only untuk memberikan tingkat kemampuan aljabar linier yang masuk akal dalam paket sesederhana mungkin, dan saya mencoba mensurvei seperti apa keadaan seni saat ini: menghitung SVD dari sebuah matriks kompleks. Saya sedang melakukan dekomposisi dua fase, bidiagonisasi diikuti oleh perhitungan nilai singular. Saat ini saya …

12 linear-algebra  matrix  svd 

Adakah yang bisa merekomendasikan metode untuk masalah kuadrat-terkecil berikut: R ∈ R3 × 3R∈R3×3R \in \mathbb{R}^{3 \times 3}∑i = 0N( R xsaya- bsaya)2→ min∑saya=0N(Rxsaya-bsaya)2→min\sum\limits_{i=0}^N (Rx_i - b_i)^2 \rightarrow \minRRR Saya bisa mendapatkan solusi perkiraan dengan meminimalkan (arbitrer ), dengan mengambil matriks dan:∑i = 0N( A xsaya- bsaya)2→ min∑saya=0N(SEBUAHxsaya-bsaya)2→min\sum\limits_{i=0}^N (Ax_i - …

11 linear-algebra  least-squares  svd 

Apakah ada algoritma SVD terpotong yang menghitung nilai singular satu per satu? Masalah saya: Saya ingin menghitung nilai singular pertama (dan vektor singular) dari matriks padat , tapi saya tidak tahu apa nilai tepat. itu besar, jadi untuk alasan efisiensi, saya lebih suka tidak mengevaluasi SVD lengkap hanya untuk memotong …

11 linear-algebra  algorithms  svd 

Misalnya pustaka matriks jarang C ++ yang saya gunakan - Eigen dan SuiteSparse, mereka tampaknya tidak memiliki funcitionality SVD untuk matriks jarang. Jadi hanya ingin tahu, apakah SVD lebih sulit daripada QR / LU untuk matriks jarang?

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

Untuk reduksi model, saya ingin menghitung vektor singular kiri terkait dengan - katakanlah 20 - nilai singular terbesar dari matriks , di mana N ≈ 10 6 dan k ≈ 10 3 . Sayangnya, matriks A saya akan padat tanpa struktur apa pun.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx 10^3AAA Jika saya …

10 linear-algebra  python  reference-request  numpy  svd 

Saya telah menggunakan Intel MKL SVD ( dgesvdmelalui SciPy) dan memperhatikan bahwa hasilnya sangat berbeda ketika saya mengubah presisi di antara float32dan float64ketika matriks saya dikondisikan secara buruk / bukan peringkat penuh. Apakah ada panduan tentang jumlah minimum regularisasi yang harus saya tambahkan untuk membuat hasil tidak sensitif terhadap float32-> …

10 numerical-analysis  stability  svd  intel-mkl  numerical-limitations