Qu'est-ce qu'un algorithme simple pour calculer la SVD de matrices?2×22×22 \times 2 Idéalement, j'aimerais un algorithme numériquement robuste, mais j'aimerais voir à la fois des implémentations simples et pas si simples. Code C accepté. Des références à des articles ou à du...
Je travaille sur une bibliothèque de matrices uniquement en-tête pour fournir un certain degré raisonnable de capacité d'algèbre linéaire dans un package aussi simple que possible, et j'essaie de faire le point sur l'état actuel de la technique: le calcul de la SVD d'un matrice complexe. Je fais...
Existe-t-il un algorithme SVD tronqué qui calcule les valeurs singulières une par une? Mon problème: je voudrais calculer les premières valeurs singulières (et vecteurs singuliers) d'une grande matrice dense M , mais je ne sais pas ce que serait une valeur appropriée de k . M est grand, donc pour...
Quelqu'un pourrait-il recommander une méthode pour le problème des moindres carrés suivant: trouver qui minimise: , où R est unitaire (rotation) matrice.R∈R3×3R∈R3×3R \in \mathbb{R}^{3 \times 3}∑i=0N(Rxi−bi)2→min∑i=0N(Rxi−bi)2→min\sum\limits_{i=0}^N (Rx_i - b_i)^2 \rightarrow \minRRR Je pourrais...
Pour la réduction du modèle, je veux calculer les vecteurs singuliers de gauche associés aux - disons 20 - plus grandes valeurs singulières d'une matrice , où N ≈ 10 6 et k ≈ 10 3 . Malheureusement, ma matrice A sera dense sans aucune structure.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx...
Par exemple, les bibliothèques de matrice clairsemée C ++ que j'ai utilisées - Eigen et SuiteSparse, elles ne semblent pas avoir de fonctionnalité SVD pour la matrice clairsemée. Donc, juste curieux, SVD est-il plus difficile que QR / LU pour une matrice
J'ai utilisé le SVD d'Intel MKL ( dgesvdvia SciPy) et j'ai remarqué que les résultats sont considérablement différents lorsque je change de précision entre float32et float64lorsque ma matrice est mal conditionnée / n'est pas pleine. Existe-t-il un guide sur le montant minimum de régularisation à...