Why matlab returns inaccurate response for rank of my matrix?
14 views (last 30 days)
This is my matrix:
6.07787530038946e+15 0 0 0 0 0 0 0
0 6.07787530038946e+15 0 0 0 0 0 0
0 0 6.07787530038946e+15 0 0 0 0 0
0 0 0 1215578336877.89 0 0 0 0
0 0 0 0 1215578336877.89 0 0 0
0 0 0 0 0 0 0 0
0 0 0 0 0 0 -1.47264000000000e-12 0
0 0 0 0 0 0 0 -1.96352000000000e-16
When ask matlab to say the rank of A it returns 5, while it is 7 as you see. why it returns 5, and how to fix it?
Jan on 28 May 2017
If you have a 2000x2000 matrix with elements on the diagonal, which vary by a factor of 10^30 on the diagioal, and do not want to process the data symbolically, you are simply lost. This is the wrong method to solve your problem. It can neither be done symbolically nor numerically. The only solution will be to do it mathematically, which means using your brain.
You have to determine the meaning of the values during the creation already. You need a method to decide, if 1e-15 is a rounding error or a meaningful value. Note that this is not trivial, because then you have to consider 0 also as potential meaningful value, when it was created e.g. by 1e17 + 1 - 1e17.
A standard approach in numerics is scaling the values: Together with e.g. a vector of positions, another vector is created carrying scaling factors such that the numbers get about 1. Then e.g. a variation of the values is much easier, or the calculation of the Hessian matrix. Scaling is used in optimization and simulation of dynamic systems, example: The integration of the equations of motion of a lunar orbiter. To control vibrations of the engines, millimeters matter. Right at the start the position is known in the magnitud of centimeters, but when the lunar orbit is reached, such a precision would be nonsense. A scaling of the different components of the trajectory allows the integrator to consider them with a physically meaningful precision.
Summary: Think twice. Your problem cannot be solved by standard methods.
Star Strider on 27 May 2017
The result of the rank function is correct. With respect to the first 5 values of diag(A), the last three are vanishingly small and evaluate to zero. For all intents and purposes, rank(A) is 5.
John D'Errico on 28 May 2017
Edited: John D'Errico on 28 May 2017
As others have said, just wanting magic to happen will not suffice. In MATLAB, you need either to use symbolic tools for that computation, or numerical tools. The numerical tools in MATLAB CANNOT be used to give you higher than double precision accuracy, because they are written in double precision. Your problem is wildly beyond the domain of double precision arithmetic.
Of course, if your problem is too large, then symbolic tools will also not suffice, as they will take far too long for any possible solution.
There is no law that says that ALL problems are solvable. If you choose to create a problem that is too big and with too wide of a dynamic range, then expect your problem to fall into the set that do not have a solution, at least in MATLAB. And tomorrow, we would probably find out that you have now decided to make your problem even more difficult, with a wider dynamic range yet, because if it was easy then anybody can do it.
Nothing stops you from writing or finding linear algebra code in some other language that supports a highly extended precision arithmetic. Even an implementation of quad precision would probably be insufficient here. But but be careful, as infinite precision tools (like the JAVA tools) tend to do poorly on problems like this, as those computations generate numbers with more and more digits at each step of the process.
A solution could be to write a version of rank that used my HPF tool, it would probably be fine on your problem if you decided to work in 50 or 100 decimal digits of precision. You could write a simple version of rank then, based on a column pivoted Q-less QR algorithm. Or use a rref variation, that will work on HPF inputs. (SVD would be slower, and take a bit more effort. rref would be easy enough to hack to allow HPF. I might even write one to put in HPF if I have some free time.) But at least using a variable precision tool like HPF, you can choose the precision to work in.
(Actually, rref is REALLY easy to hack to work with HPF inputs).