原文首发于blog.scloudyy.com,转载请注明出处
记一次Matlab处理大型系数矩阵过程中的探索。
创建一个系数矩阵有很多方法,最简单的如
A = sparse(1000, 1000)
这会创建一个1000*1000的系数矩阵,虽然很大,其实没有占什么内存,因为稀疏矩阵只会存储非零的值,刚创建的稀疏矩阵没有存储任何值,默认全为零。
赋值时可以按普通矩阵的方法赋值
A(1,1) = 2
讲稀疏矩阵的第一行第一列赋值为2
但如果要往稀疏矩阵里放很多值,不推荐以上做法,因为稀疏矩阵内部是类似链表的实现,一次赋值太多速度会非常慢,我之前忘稀疏矩阵里添加大约十万个值发现一个下午都不一定能跑完
如果要往系数矩阵里赋大量的值,正确的做法应该是这样的。
将元素的行坐标,列坐标,还有元素的值分别存入三个向量中,三个向量中各个分量一一对应,然后使用spar