Singular Value Decomposition (SVD) ist eine Zerlegung (Faktorisierung) einer rechteckigen reellen oder komplexen Matrix in das Produkt einer einheitlichen Rotationsmatrix, einer diagonalen Skalierungsmatrix und einer zweiten einheitlichen Rotationsmatrix.
Was ist ein einfacher Algorithmus zur Berechnung der SVD von 2×22×22 \times 2 Matrizen? Idealerweise hätte ich gerne einen numerisch robusten Algorithmus, aber ich würde gerne sowohl einfache als auch nicht so einfache Implementierungen sehen. C-Code akzeptiert. Verweise auf Papiere oder...
Ich arbeite an einer Nur-Header-Matrixbibliothek, um einen angemessenen Grad an linearer Algebra in einem möglichst einfachen Paket bereitzustellen, und ich versuche, einen Überblick über den aktuellen Stand der Technik zu erhalten: die SVD von a komplexe Matrix. Ich mache eine zweiphasige...
Gibt es einen abgeschnittenen SVD-Algorithmus, der die Singularwerte einzeln berechnet? Mein Problem: Ich möchte die ersten kkk Singularwerte (und Singularvektoren) einer großen dichten Matrix berechnen MMM, weiß aber nicht, was ein geeigneter Wert von kkk wäre. MMM ist groß, daher würde ich aus...
Könnte jemand eine Methode für das folgende Problem der kleinsten Quadrate empfehlen: finde R∈R3×3R∈R3×3R \in \mathbb{R}^{3 \times 3} , das minimiert: ∑i=0N(Rxi−bi)2→min∑i=0N(Rxi−bi)2→min\sum\limits_{i=0}^N (Rx_i - b_i)^2 \rightarrow \min , wobei RRR eine einheitliche (Rotations-) Matrix ist. Ich...
Zur Modellreduktion möchte ich die linken Singularvektoren berechnen, die den - sagen wir 20 - größten Singularwerten einer Matrix , wobei und . Leider wird meine Matrix ohne Struktur dicht sein. N ≤ 10 6 k ≤ 10 3 A.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx 10^3AAA...
Zum Beispiel scheinen die von mir verwendeten C ++ - Sparse-Matrix-Bibliotheken - Eigen und SuiteSparse - keine SVD-Funktionalität für Sparse-Matrix zu haben. Also nur neugierig, ist SVD für eine spärliche Matrix schwieriger als QR /
Ich habe die SVD von Intel MKL ( dgesvdüber SciPy) verwendet und festgestellt, dass die Ergebnisse erheblich voneinander abweichen, wenn ich die Genauigkeit zwischen float32und float64wenn meine Matrix schlecht konditioniert ist / nicht den vollen Rang hat. Gibt es einen Leitfaden zum Mindestmaß an...
Ich löse ein System mit Singularwertzerlegung. Die Singularwerte (vor der Skalierung) sind: 1.82277e+29 1.95011e+27 1.15033e+23 1.45291e+21 4.79336e+17 7.48116e+15 8.31087e+12 1.71838e+11 5.63232e+08 2.17863e+08 9.02783e+07 1.72345e+07 1.73889e+05 8.09382e+02 2.16644e+00 Ich habe festgestellt, dass...
Gibt es eine Bibliothek / Toolbox, in der inkrementelle SVD in MATLAB implementiert ist? Ich habe dieses Papier implementiert , es ist schnell, funktioniert aber nicht gut. Ich habe dies versucht, aber auch hier breitet sich der Fehler schnell aus (beim Aktualisieren von 5-10 Punkten ist der Fehler...
Bei einer großen Matrix EINAA mit Eigenwerten σ1≥ σ2≥ …σ1≥σ2≥…\sigma_1\ge \sigma_2 \ge \dotsc möchte ich nur eine Teilmenge dieser Werte bestimmen, z. B. σ5, σ8σ5,σ8\sigma_5,\sigma_8 und σ19σ19\sigma_{19} . Gibt es einen Algorithmus, der dies kann, oder ist es am besten, die Top-19-Eigenwerte zu...