Mldivide gives rank deficency when solving for large full system

13 vues (au cours des 30 derniers jours)
Giulio Tagliaferro
Giulio Tagliaferro le 10 Déc 2019
Commenté : Star Strider le 11 Déc 2019
I have tried to solve a full linear system of approx 3 millions equation and 3 uknown.
If I try to solve the system (A\y) using mldivide i get:
Warning: Rank deficient, rank = 1, tol = 4.893629e+05.
And the results are completely wrong
However when forming the normal equation system (A'*A) x = A'*y the solver runs perfectly fine which is not surprising. (I know that in general is less numerically stable but in this case works fine)
From the documentation I read mldivide should use the QR solver on rectangular matrices. If i try to call QR on the A matrix I get:
Error using qr
Requested 3107227x3107227 (71934.3GB) array exceeds maximum array size preference.
My question are: what is mldivive trying to do? why it says that a rank deficecy is present?

Réponses (1)

Star Strider
Star Strider le 10 Déc 2019
See if the lsqr function will do what you want.
Another option is to use pinv to do the matrix inversion (not normally an approach I would recommend, however this appears to be an unusual situation).
  2 commentaires
Giulio Tagliaferro
Giulio Tagliaferro le 11 Déc 2019
Playing with tolernace and number of iteration i was able to get the correct result using lsqr.Also solving uisng pinv seem to give a resonable result. This still does not explain the behaviour of mldivide.
Star Strider
Star Strider le 11 Déc 2019
This still does not explain the behaviour of mldivide.
It does, actually. See the Algorithms section of the documentation. The mldivide function can deal with some sparse matrices, although not matrices that are badly scaled or nearly singular. That is likely the reason that lsqr and the similar functions are included in the core MATLAB functions.

Connectez-vous pour commenter.

Catégories

En savoir plus sur Linear Algebra dans Help Center et File Exchange

Produits


Version

R2016a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by