Замените ненулевые элементы разреженной матрицы на единицы
R = spones(S)
R = spones(S) генерирует матричный R с той же структурой разреженности как S, но с 1 в ненулевых положениях.
c = sum(spones(S)) является количеством ненулей в каждом столбце.
r = sum(spones(S'))' является количеством ненулей в каждой строке.
sum(c) и sum(r) равны, и равны nnz(S).