创建0s和1s的大型稀疏矩阵

Question

我有以下代码，但我必须运行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