How can I approximate this large matrix?
5 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I have a sparse symmetric matrix M, roughly 10^6 x 10^6, but only about 100 values in each column (or row) are non-zero, so they'll all fit memory / disk space. Values are complex, and all magnitudes are < 10^-1 I need an approximation of
M * inv(eye() - M) = M + M^2 + M^3 + M^4 + ...
The inverse in there won't be sparse, but the produc hopefully will be.
So, if I have the nonzero values and their locations tabled up, what's the best way to go about this calculation in MATLAB?
Are there library tools avaiable for large sparse matrices?
And should I try to approximate the LHS or the RHS? I know the RHS can be approximated by calculating M^2, M^4, M^8,... and then multiplying M * (I+M) * (I+M^2) * (I+M^4) * (I+M^8) = M + M^2 + ... + M^16. But I'm hoping there might be some iterative way to approximate the LHS. And in either case I'd still need to know how to set M up as a sparse matrix and take advantage of M's sparseness.
0 commentaires
Réponses (1)
Bruno Luong
le 4 Nov 2022
Modifié(e) : Bruno Luong
le 4 Nov 2022
What about
T = M;
niter = 5:
for k=1:niter
T = M*(speye(size(M)) + T);
end
This returns
T = M + M^2 + .... + M^6
6 commentaires
Bruno Luong
le 4 Nov 2022
Modifié(e) : Bruno Luong
le 4 Nov 2022
But MATLAB sparse stores only non-zeros elements.
If your matrix are big and not sparse enough to be store as sparse in the memory of your PC, then you are simply stuck.
I don't know what you meant by tool, MATLAB matrix library comprises a comprehensive array of functions that can deal with sparse matrix.
Voir également
Catégories
En savoir plus sur Sparse Matrices dans Help Center et File Exchange
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!