The way to solve a singular matrix
66 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Hi
There is any one know how the method to decompose the singular square matrix using Matlab. Someone told me the Matlab have something like a ready Forthran subroutine. Does anyone know how to use it in Matlab?
0 Commenti
Risposte (3)
John D'Errico
il 23 Ago 2014
Modificato: John D'Errico
il 23 Ago 2014
help pinv
Not much more to say, since you give very little info to help you on. Note that computing the inverse of a matrix is almost never recommended. The backslash operator is a better choice always than inv. But pinv is a good tool for this purpose, when backslash (and surely also inv) will fail.
A = ones(2);
A\[1;1]
Warning: Matrix is singular to working precision.
ans =
NaN
NaN
inv(A)*[1;1]
Warning: Matrix is singular to working precision.
ans =
Inf
Inf
pinv(A)*[1;1]
ans =
0.5
0.5
2 Commenti
Einat Shoval
il 24 Dic 2017
Thank you so much for this!! Was stuck on this for two days now until I found your answer :)
Jess
il 22 Mar 2016
% Goal: solve A*x == b for x
% Set up some matrix A (I used a sparse matrix) -- do yourself
% Set up the vector b -- do yourself
% Perform SVD on A
[U,S,V] = svd(A);
% A == U*S*V' % Not needed, but you can check it yourself to confirm
% Calc number of singular values
s = diag(S); % vector of singular values
tolerance = max(size(A))*eps(max(s));
p = sum(s>tolerance);
% Define spaces
Up = U(:,1:p);
%U0 = U(:,p+1:Nx);
Vp = V(:,1:p);
%V0 = V(:,p+1:Nx);
%Sp = spdiags( s(1:p), 0, p, p );
SpInv = spdiags( 1.0./s(1:p), 0, p, p );
% Calc AInv such that x = AInv * b
AInv = Vp * SpInv * Up';
x = AInv * b; % DONE!
1 Commento
Xin Shen
il 10 Lug 2019
When I tried your idea to solve my problem, I got an error "SVD does not support sparse matrices. Use SVDS to compute a subset of the singular values and vectors of a sparse matrix". Does SVDS work for your idea?
Vedere anche
Categorie
Scopri di più su Operating on Diagonal Matrices in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!