How can I approximate this large matrix?

1 次查看（过去 30 天）

Jerry Guern 2022-11-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1843443-how-can-i-approximate-this-large-matrix

评论： Jerry Guern 2022-11-4

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 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Bruno Luong 2022-11-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1843443-how-can-i-approximate-this-large-matrix#answer_1091278

编辑：Bruno Luong 2022-11-4

在 MATLAB Online 中打开

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 个评论
显示 4更早的评论隐藏 4更早的评论

Bruno Luong 2022-11-4

编辑：Bruno Luong 2022-11-4

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.

Jerry Guern 2022-11-4

Ah, thank you, this was helpful on the sparse matrix handling part. I realize now I probably should have posted this as two separate questions, one on the sparse matrix tools and one on the approximation.

请先登录，再进行评论。

请先登录，再回答此问题。