匿名用户
1级
2018-10-19 回答
矩阵太大了,我很怀疑你的K矩阵能表示出来吗?更不要说后面的求逆及相乘的运算了。
按照你贴出来的数据,K矩阵的维数是 94507*94507,每个元素按照double类型考虑,占8个字节,则K矩阵占用的内存就是
>> 94507*94507*8/1024^3
ans =
66.545也就是66.5G,你的机器内存的确挺多了,但离这个需求还差很远吧?
这么大的矩阵,其元素都有什么样的规律?如果大多数为0,只有少量的非零元素,可以考虑使用稀疏矩阵,例如:
>> K=sparse(94507,94507)
K =
All zero sparse: 94507-by-94507
>> whos
Name Size Bytes Class
K 94507x94507 378044 double array (sparse)
ans 1x1 8 double array
Grand total is 2 elements using 378052 bytes
追问:
我就是用的稀疏矩阵啊。。。。。已经表示好了,就是在求解的时候出现了inv(K)*p的时候出现了问题
追问:
我就是用的稀疏矩阵啊。。。。。已经表示好了,就是在求解的时候出现了inv(K)*