Mldivide gives rank deficency when solving for large full system
8 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
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?
0 comentarios
Respuestas (1)
Star Strider
el 10 de Dic. de 2019
2 comentarios
Star Strider
el 11 de Dic. de 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.
Ver también
Categorías
Más información sobre Linear Algebra en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!