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