Co to jest prosty algorytm obliczania SVD macierzy ?2 × 22×22 \times 2 Idealnie chciałbym mieć solidny algorytm liczbowy, ale chciałbym zobaczyć zarówno proste, jak i nie tak proste implementacje. Kod C został zaakceptowany. Wszelkie odniesienia do artykułów lub
Pracuję nad biblioteką macierzy zawierającą tylko nagłówki, aby zapewnić pewien rozsądny stopień możliwości algebry liniowej w tak prostym pakiecie, jak to możliwe, i próbuję zbadać, jaki jest obecny stan techniki: obliczanie SVD złożona macierz. Robię dwufazowy rozkład, dwukieragonalizację, a...
Czy istnieje skrócony algorytm SVD, który oblicza wartości osobliwe pojedynczo? Mój problem: chciałbym obliczyć pierwsze wartości pojedynczych (i wektorów pojedynczych) dużej gęstej macierzy , ale nie wiem, jaka byłaby odpowiednia wartość . jest duży, więc ze względu na wydajność wolałbym nie...
Czy ktoś mógłby polecić metodę dla następującego problemu najmniejszych kwadratów: znajdź który minimalizuje: , gdzie jest jednostką (obrót) matryca.N ∑ i = 0 ( R x i - b i ) 2 → minR ∈ R3 × 3R∈R3×3R \in \mathbb{R}^{3 \times 3}∑i = 0N.( R xja- bja)2)→ min∑i=0N(Rxi−bi)2→min\sum\limits_{i=0}^N (Rx_i...
W celu zmniejszenia modelu chcę obliczyć lewe wektory osobliwe związane z - powiedzmy 20 - największymi wartościami osobliwymi macierzy , gdzie N ≈ 10 6 i k ≈ 10 3 . Niestety moja matryca A będzie gęsta bez żadnej struktury.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx...
Na przykład biblioteki rzadkich macierzy C ++, których użyłem - Eigen i SuiteSparse, wydają się nie mieć żadnej funkcji SVD dla rzadkiej macierzy. Ciekawe, czy SVD jest trudniejsze niż rzadka matryca QR /
Korzystałem z SVD Intela MKL ( dgesvdprzez SciPy) i zauważyłem, że wyniki są znacząco różne, kiedy zmieniam precyzję pomiędzy float32i float64kiedy moja matryca jest źle uwarunkowana / nie w pełni ustawiona. Czy istnieje przewodnik dotyczący minimalnej ilości regularyzacji, którą powinienem dodać,...