Эффективное создание логической разреженной матрицы

Question

Моя цель - создать прямоугольную сетку и присвоить значение да/нет каждой точке. Я занимаюсь довольно большой сеткой, поэтому решил сохранить ее как разреженную матрицу. Мой подход рассматривает N x N сетку, только с отношением r точек, чтобы быть истинными (так что Nact = r*N^2 баллов):

N = 200; 
r = 0.25; 
NAct = r*N^2; 
ss = spalloc(ceil(N),ceil(N),ceil(NAct)); 
for j = 1:N, for i = 1:N, ss(i,j) = rand < r; end; end; 
ssL = logical(ss); 

Однако это, кажется, не эффективен для больших N.

Я также попытался это:

N = 200; 
r = 0.25; 
NAct = r*N^2; 
iAct = randi(N,1,N); 
jAct = randi(N,1,N); 
sssL = sparse(iAct,jAct,true,ceil(N),ceil(N),ceil(NAct)); 
nnz(sssL) 

в качестве альтернативы, но у меня есть две проблемы: true точки N вместо NAct, и возникает ошибка, если две пары элементов iAct и jAct совпадают ,

Answer 1

N = 200; 
Grid(N,N) = false; %// preallocate grid 
Grid = sparse(Grid); %// add to make sparse 
r = 213; 
TruePoints = randi(N^2,r,1); %// get random indices of points to be true 
Grid(TruePoints) = true; 
Grid = logical(Grid); %// switch to logical 

Answer 2

Я наконец использовал sprand

N = 200; 
r = 0.25; 
NAct = r*N^2; 
ss = sprand(N,N,r); 
ssL = logical(ss);