«svd» etiketlenmiş sorular

matrisin SVD'sini hesaplamak için basit bir algoritma nedir ?2 × 22×22 \times 2 İdeal olarak, sayısal olarak sağlam bir algoritma istiyorum, ancak hem basit hem de basit olmayan uygulamaları görmek istiyorum. C kodu kabul edildi. Makalelere veya kodlara referans var mı?

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

Mümkün olduğunca basit bir pakette makul derecede doğrusal cebir yeteneği sağlamak için sadece başlık matris kütüphanesi üzerinde çalışıyorum ve sanatın şu andaki durumunun ne olduğunu araştırmaya çalışıyorum: SVD'nin hesaplanması karmaşık matris. İki fazlı ayrışma, tekdüze değer hesaplaması ve ardından tek köşeli değer hesaplaması yapıyorum. Şu anda bidiagonalization için hanehalkı yöntemini …

12 linear-algebra  matrix  svd 

Herkes aşağıdaki en küçük kareler sorunu için bir yöntem önerebilir: bulmak o en aza indirir: , R, bir bütün olan (rotasyon) matris.R ∈ R3 × 3R∈R3×3R \in \mathbb{R}^{3 \times 3}Σi = 0N-( R xben- bben)2→ dk∑i=0N(Rxi−bi)2→min\sum\limits_{i=0}^N (Rx_i - b_i)^2 \rightarrow \minR,RR Ben en aza indirerek, bir yaklaşık çözüm alabilir Σi …

12 linear-algebra  least-squares  svd 

Tekil değerleri birer birer hesaplayan kesilmiş bir SVD algoritması var mı? Benim sorun: Ben ilk hesaplamak istiyorum kkk büyük yoğun matris tekil değerler (ve tekil vektörleri) MMM , ama uygun bir değer bilmiyorum kkk olurdu. MMM büyüktür, bu nedenle verimlilik nedenleriyle, tam SVD'yi yalnızca daha sonra en küçük SV'leri kesmek …

11 linear-algebra  algorithms  svd 

Örneğin, kullandığım C ++ seyrek matris kütüphaneleri - Eigen ve SuiteSparse, seyrek matris için herhangi bir SVD işlevine sahip görünmüyorlar. Çok meraklı, seyrek matris için SVD QR / LU'dan daha mı zor?

10 linear-algebra  sparse-matrix  matrix  svd  eigen 

Modelin indirgenmesi için, matrisi olan bir matrisin en büyük tekil değerlerini - say 20 ile ilişkili olan sol tekil vektörleri hesaplamak istiyorum , burada ve . Ne yazık ki, matrisim herhangi bir yapı olmadan yoğun olacak. N ≈ 10 6 k ≈ 10 3 AA ∈ RN-, kA∈RN,kA \in \mathbb …

10 linear-algebra  python  reference-request  numpy  svd 

Intel MKL'nin SVD'sini ( dgesvdSciPy aracılığıyla) kullanıyorum ve matrisim kötü koşullu / tam sıralı değilken hassasiyeti değiştirdiğimde float32ve sonuçların önemli ölçüde farklı olduğunu fark ettim float64. Ben sonuçları için duyarsız hale getirmek için eklemek gerekir regülarizasyon minimum miktarda bir rehber var mı float32> - float64değişiklik? Özel olarak, kullanıcının durumu , …

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