MATLAB的答案

我怎么能近似这个大型矩阵?

3视图(30天)
杰瑞Guern
杰瑞Guern 2022年11月4日
评论道: 杰瑞Guern2022年11月4日
我有一个稀疏对称矩阵M,大约10 ^ 6 x 10 ^ 6,但只有100的值在每个列(或行)是零,所以他们会适合所有内存/磁盘空间。价值观是复杂的,大小是我需要一个< 10 ^ 1 近似 的
M *发票(眼()- M) = M + M ^ 2 + M M ^ ^ 3 + 4 +…
不会有稀疏的逆,但名单的希望。
所以,如果我有非零值及其位置提出,最好的方法是什么在MATLAB计算吗?
有图书馆大型稀疏矩阵的用具吗?
和我应该试图近似lh和RHS吗?我知道的皇家可以近似计算M ^ 2, M ^ 4 M ^ 8,……然后乘以M *(我+ M) * (I + M ^ 2) * (I + M ^ 4) * (I + M ^ 8) = M + M ^ 2 +…+ M ^ 16。但是我希望会有一些迭代方法来近似lh。和在这两种情况下我仍然需要知道如何设置米作为一个稀疏矩阵,利用M的稀疏。

在回答这个问题。

答案(1)

布鲁诺陈德良
布鲁诺陈德良 2022年11月4日
编辑:布鲁诺陈德良 2022年11月4日
是什么
T = M;
硝石= 5:
为k = 1:硝石
T = M * (speye(大小(M)) + T);
结束
这将返回
T = M + M ^ 2 + ....+ M ^ 6
6个评论
显示5年长的评论
杰瑞Guern
杰瑞Guern 2022年11月4日
啊,谢谢你,这是有益的稀疏矩阵处理部分。我现在意识到我可能应该张贴这是两个独立的问题,一个在稀疏矩阵的工具和一个在近似。

登录置评。

类别

找到更多的在稀疏矩阵在帮助中心和文件交换

标签

下载188bet金宝搏

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!