What will be a suitable option to circumvent a non-invertible matrix?
6 views (last 30 days)
Show older comments
Line 157 of this code-file (relevant part is reproduced below) used to calculate the optical flow between two images occassionally returns an error:
% .........
% Solve this system of linear equations, adding a small value along the
% diagonal to avoid potentially having a singular matrix
diag_small_value = sparse(1:2*N, 1:2*N, 1e-10);
A = A + diag_small_value;
xexact = A\b;
% .........
>>> Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate.
So, any downsides to using
pinv(A)*B
under a try-except block? Or, maybe we avoid the iterative computation since we already have the eigen vector though that will be weird .....
1 Comment
Walter Roberson
on 27 Sep 2023
Note that "Matrix is singular" is a warning not an error, so try ; catch would not catch it.
I do not know anything about the mathematics of optical flow, so I do not know what equations are being solved or why they might give problems.
Answers (2)
Matt J
on 27 Sep 2023
The approach is valid if A is symmetric, positive definite, but 1e-10 may not be large enough to add substantial conditioning.
2 Comments
Bruno Luong
on 27 Sep 2023
When you encount warning message like this, the first thing to ask is not "how do I change linear inversion algorithm?", but "what condition I forgot to take into account in order for my problem to be wellposed?".
In the optical flow, the answer could be that you forgot to impose some boundary condition so that the flow is uniquely determine.
If you won't fix that, any algorithm will fail to give you the correct answer.
See Also
Categories
Find more on Linear Algebra in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!