24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sapphirechina

铜虫 (正式写手)

应助: 13 (小学生)
金币: 624.8
散金: 178
红花: 1
帖子: 385
在线: 332.9小时
虫号: 431292
注册: 2007-08-11
专业: 环境与生态统计学

[求助] 如何存储和计算超大规模稀疏矩阵已有1人参与

本人在R语言环境中生成了一个18万行*18万列的单位矩阵，然后用循环给这个单位矩阵赋值，结果循环到8万多次就报错无法分配内存了。
代码如下：
#生成单位矩阵，用普通矩阵无法创建，所以用sparseMatrix来创建。
num_samples=180000
rownumber=c(1:num_samples)
colnumber=c(1:num_samples)
value=c(rep(1,each=num_samples))
M=sparseMatrix(rownumber,colnumber,x=value)

#给单位矩阵赋值，全部赋值完成后是一个4*50*180000个数值的对称稀疏矩阵。
for(i in 1:num_samples){
w=W[,i]                            #有50个浮点型数值。
n=neighbors[,i]                #有50个整型数值。
M[i,n]=M[i,n]-t(w)
M[n,i]=M[n,i]-w
M[n,n]=M[n,n]+w%*%t(w)
}

请教上述的情况该怎么处理呢？我的机器内存是4G，32位操作系统。
如果这个矩阵能顺利赋值完毕，我随后还需要计算这个矩阵的特征值和特征向量。
请教可能实现么？

回复此楼

» 猜你喜欢

Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有6人回复
2025冷门绝学什么时候出结果已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
论文投稿，期刊推荐已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复

1楼2016-11-26 09:18:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sapphirechina

铜虫 (正式写手)

应助: 13 (小学生)
金币: 624.8
散金: 178
红花: 1
帖子: 385
在线: 332.9小时
虫号: 431292
注册: 2007-08-11
专业: 环境与生态统计学

引用回帖:

9楼: Originally posted by FMStation at 2016-12-11 21:13:41
https://www.mathworks.com/help/matlab/math/constructing-sparse-matrices.html

S = sparse(i,j,s,m,n)

S = sparse(,,,4,4)

D = sparse(1:n,1:n,-2*ones(1,n),n,n);
E = sparse(2:n,1:n-1,ones(1,n-1 ...

太好了，这是matlab代码吧？我试试看。