Средство ссылок на разреженную матрицу

Question

Итак, я работаю с разреженной матрицей, и мне нужно узнать различную информацию о очень большом (размер 10^6), и мне нужно узнать среднее значение исходящих ссылок. Просто, чтобы быть уверенным, что я имею в виду то, что вы получаете от 3 + 4 + 5/3 = 4, 4 - это среднее значение.

Я думал, что-то вроде этого:

[row,col] = find(A(:,2),1,'first') 

, а затем я хотел бы сделать 1/numberInThatIndex или что-то подобное, так как это S-матрица (уверен, что это называется, что).

И я буду перебирать столбцы по столбцам, но по какой-то причине он не дает мне первое число в каждом столбце, если я нахожу (A (:, 1), 1, 'first'), это дает мне первый в первом столбце, но не во втором, если я изменю его на A (:, 2).

Мне также нужно что-то, чтобы сохранить этот индекс, чтобы получить доступ к значению, я подумал о векторе 2xN, но я думаю, что это не лучшая идея. Я имею в виду, что поиск будет давать мне индекс, но мне нужно значение в этом индексе, а затем сохранить это или показать его. Не уверен, что я объясняю себя правильно, но я пытаюсь, извините.

Просто, чтобы быть понятным как при вводе A (:, 1), так и A (:, 2), он дает мне индекс из первого столбца, и я не хочу этого, я хочу, чтобы первый элемент был найден из каждого столбца, поэтому я могу вычислить среднее из числа в этом индексе.

редактирование: все кажется, что это действительно работает, но когда я проверял результаты, я помещал 3817 вместо 3871, который был данным ответом, и поэтому я нашел 0, когда мне захотелось чего-то, что не является нулем. Не уверен, что я должен удалить все это.

Answer 1

Чтобы решить вашу проблему, вы можете сделать следующее:

numberNonZerosPerColumn = sum(S~=0,1); 
meanValue = nanmean(1./numberNonZerosPerColumn); 

• Подсчитайте число ненулевых элементов в каждом столбце n(i)
• Compute значение v(i), которые хранятся там, которые определяются v(i) := 1/n(i)
• Возьмите среднее значение тех значений, где n(i) не равно нулю (т.е. суммируя все эти значения, где v(i) не является NaN и разделить на число столбцов, которые содержат по крайней мере один нуль)

Если вы хотите, чтобы рассматривать столбцы без любого ненулевого входа в v(i):= 0, но по-прежнему использовать их в среднем, вы можете использовать:

numberNonZerosPerColumn = sum(S~=0,1); 
meanValue = nansum(1./numberNonZerosPerColumn)/size(S,2);