我有以下代码,但我必须运行它 n>10^6。对于 n=10^5,内存结束。我有 16GB。有没有有效的方法来做到这一点?
n = 10^6;
m = n/2;
H = sparse(m,n);
parfor jj=1:n
pos = sample(cumulative);
H(:,jj) = (rand(m,1) > 1 - degrees(pos)/m);
end