sparse matrix problem in solve systems of linear equations
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
NA
am 24 Jan. 2019
Bearbeitet: John D'Errico
am 24 Jan. 2019
I have F= 1.0e+05 *[-6.2281; -2.6137;0.1177; -0.4726]
J= 1.0e+05 *
(1,1) 1.2527
(2,1) -0.2980
(3,1) -0.1361
(4,1) 0.0240
(1,2) -0.2980
(2,2) 0.3484
(3,2) -0.0156
(4,2) -0.0416
(1,3) -0.1361
(2,3) -0.0156
(3,3) 1.2470
(4,3) -0.2760
(1,4) 0.0240
(2,4) -0.0416
(3,4) -0.2760
(4,4) 0.4438
when I want to calculate
dx = (J \ F);
gives Warning: Matrix is singular to working precision.
I change it to dx = pinv(J)*F
gives this error (SVD does not support sparse matrices. Use SVDS to compute a subset of the
singular values and vectors of a sparse matrix). so I change it to dx = pinv(full(J))*F. Does it correct ?
0 Kommentare
Akzeptierte Antwort
John D'Errico
am 24 Jan. 2019
Bearbeitet: John D'Errico
am 24 Jan. 2019
J is not sparse. In fact, J is essentially a full matrix. It has no zeros in it at all! So why are you using sparse matrix storage? And especially on such a small, tiny matrix? I think you are using sparse storage for no good reason. In fact, sparse is costing you capability, time, and even memory. Not a lot. But it is just a bit silly here.
As far as using pinv(full(J))*F, yes, that is a way to resolve the singularity issues, but it may be a poor choice. A singular matrix is usually a sign that something is wrong. Yes, you can close your eyes, pushing forth, full speed ahead. But far too often that just results in a train wreck.
As a wild guess, the use of the variable names J & F here indicates to me that you might be trying to do some sort of Newton iteration. A singular matrix there suggests that you are getting into a bad place in that search, that something is wrong. But then, what do I know, since you have told us nothing at all, except that you seem to be very confused in why you are doing what you are doing. So maybe, just maybe, it is time to slow down. Look at the warnings, and take that curve ahead with the brakes applied. (In other words, think very carefully about what you are doing.)
0 Kommentare
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Creating and Concatenating Matrices finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!