Как я могу посчитать количество элементов данного значения в матрице?

Question

Как я могу посчитать количество элементов данного значения в матрице?

66

Кто-нибудь знает, как подсчитать количество раз, когда значение появляется в матрице?

Например, если у меня есть матрица размером 1500 x 1 M (вектор), которая хранит значения будних дней (1-7), как я могу подсчитать, сколько воскресений (1), понедельников (2),..., Суббота (7) хранятся в M?

Niko Gamulin 21 май 2010, в 08:55

Источник

Теги:

matlab

matrix

count

7 ответов

73

Здесь перечислены все способы, с помощью которых можно было бы подсчитать уникальные элементы:

M = randi([1 7], [1500 1]);

Вариант 1: таблица

t = tabulate(M);
counts1 = t(t(:,2)~=0, 2);

Вариант 2: hist/histc

counts2_1 = hist( M, numel(unique(M)) );
counts2_2 = histc( M, unique(M) );

Вариант 3: накопитель

counts3 = accumarray(M, ones(size(M)), [], @sum);
%# or simply: accumarray(M, 1);

Вариант 4: sort/diff

[MM idx] = unique( sort(M) );
counts4 = diff([0;idx]);

Вариант 5: arrayfun

counts5 = arrayfun( @(x)sum(M==x), unique(M) );

Вариант 6: bsxfun

counts6 = sum( bsxfun(@eq, M, unique(M)') )';

Вариант 7: разреженный

counts7 = full(sparse(M,1,1));

Amro 21 май 2010, в 20:08

2

Второй вариант - это то, что я искал, спасибо!
Martin Vseticka 15 март 2012, в 11:21
4

Я предпочитаю этот ответ тому, который связан с математическими работами, поскольку он очень ясен и самодостаточен.
gakera 28 апр. 2014, в 14:26
0

Я сделал немного времени, и Вариант 3 обычно самый быстрый.
Roun 21 июль 2014, в 23:11
0

@Amro При M = [1 2 100000] , вариант 2, counts2_1 = hist( M, numel(unique(M)) ) дает [2 0 1] что недостаточно хорошо, да? Это должно быть просто counts2_1 = hist( M, unique(M) ) .
Evgeni Sergeev 09 авг. 2014, в 04:29
0

@EvgeniSergeev: Ты прав. Второй аргумент HIST может иметь различные значения: 1) это может быть скаляр, обозначающий количество бинов hist(data,nbins) (в этом случае диапазон [min(data),max(data)] делится на равные размер N бинов) 2) он также может указывать центры hist(data,xvalues) виде hist(data,xvalues) , предполагая, что xvalues является вектором равномерно распределенных значений 3) или, если вектор xvalues не является равномерно распределенным, то средние точки между последовательные значения используются в качестве ребер бина.
Amro 09 авг. 2014, в 14:46
0

@EvgeniSergeev: во всех трех случаях код в hist внутренне заканчивает тем, что histc функцию histc (проверьте edit hist.m исходного кода edit hist.m , тогда как histc - встроенная функция). Поэтому я думаю, что здесь лучше избегать hist и напрямую вызывать histc с явными ребрами бина, как я показал выше.
Amro 09 авг. 2014, в 14:47
1

Примеры, которые я привел выше, были в основном предназначены для данных, показанных OP, а некоторые методы предполагают, что значения данных являются целыми числами и предпочтительно начинаются с 1 . См. Следующие ответы о том, как обрабатывать больше случаев в целом: stackoverflow.com/a/4093228/97160 , stackoverflow.com/a/6933863/97160
Amro 09 авг. 2014, в 15:14

Показать ещё 5 комментариев

9

Один из способов выполнения этой операции для всех значений с 1 по 7 - это использовать функцию ACCUMARRAY:

>> M = randi(7,1500,1);  %# Some random sample data with the values 1 through 7
>> dayCounts = accumarray(M,1)  %# Will return a 7-by-1 vector

dayCounts =

   218       %# Number of Sundays
   200       %# Number of Mondays
   213       %# Number of Tuesdays
   220       %# Number of Wednesdays
   234       %# Number of Thursdays
   219       %# Number of Fridays
   196       %# Number of Saturdays

gnovice 21 май 2010, в 15:16

5

Предположим, что w содержит номера недель ([1: 7])

n = histc(M,w)

если вы не знаете диапазон чисел в M:

n = histc(M,unique(M))

Он является таким, как SQL Group по команде!

Babak 15 апр. 2012, в 22:56

4

Это было бы идеальной причиной того, что мы делаем операцию над матрицей, и ответ должен быть единственным числом

sum(sum(matrix==value))

radhe_shyam 27 июнь 2013, в 15:42

3

Это очень хороший файл функций, доступный на центральном файловом сервере Matlab.

ссылка countmember.m

Этот файл функции полностью векторизован и, следовательно, очень быстрый. Кроме того, по сравнению с функцией, которая упоминается в ответе на aioobe, эта функция не использует функцию аккумулирования, поэтому она даже совместима со старыми версиями Matlab. Кроме того, он работает как с массивами ячеек, так и с числовыми массивами.

РЕШЕНИЕ: Вы можете использовать эту функцию в сочетании со встроенной функцией matlab, "unique".

occurance_count = countmember (unique (M), M)

occurance_count будет представлять собой числовой массив с тем же размером, что и у уникального (M), а различные значения массива occurance_count будут соответствовать счету соответствующих значений (одного и того же индекса) в уникальном (M).

Abhinav 29 сен. 2011, в 10:45

2

Используйте nnz вместо суммы. Нет необходимости в двойном вызове, чтобы свернуть матрицы к векторам и, скорее всего, быстрее, чем сумма.

nnz(your_matrix == 5)

Doc

Andy Campbell 31 авг. 2015, в 21:23

Ещё вопросы

Второй вариант - это то, что я искал, спасибо!
Я предпочитаю этот ответ тому, который связан с математическими работами, поскольку он очень ясен и самодостаточен.
Я сделал немного времени, и Вариант 3 обычно самый быстрый.
@Amro При M = [1 2 100000] , вариант 2, counts2_1 = hist( M, numel(unique(M)) ) дает [2 0 1] что недостаточно хорошо, да? Это должно быть просто counts2_1 = hist( M, unique(M) ) .
@EvgeniSergeev: Ты прав. Второй аргумент HIST может иметь различные значения: 1) это может быть скаляр, обозначающий количество бинов hist(data,nbins) (в этом случае диапазон [min(data),max(data)] делится на равные размер N бинов) 2) он также может указывать центры hist(data,xvalues) виде hist(data,xvalues) , предполагая, что xvalues является вектором равномерно распределенных значений 3) или, если вектор xvalues не является равномерно распределенным, то средние точки между последовательные значения используются в качестве ребер бина.
@EvgeniSergeev: во всех трех случаях код в hist внутренне заканчивает тем, что histc функцию histc (проверьте edit hist.m исходного кода edit hist.m , тогда как histc - встроенная функция). Поэтому я думаю, что здесь лучше избегать hist и напрямую вызывать histc с явными ребрами бина, как я показал выше.
Примеры, которые я привел выше, были в основном предназначены для данных, показанных OP, а некоторые методы предполагают, что значения данных являются целыми числами и предпочтительно начинаются с 1 . См. Следующие ответы о том, как обрабатывать больше случаев в целом: stackoverflow.com/a/4093228/97160 , stackoverflow.com/a/6933863/97160

aioobe · Accepted Answer · 2010-05-21T11-13-00.000Z

Посмотрите Определите и подсчитайте уникальные значения массива.

Или, чтобы подсчитать количество вхождений 5, просто выполните

sum(your_matrix == 5)

Хороший и простой ответ. А чтобы расширить код для работы с несколькими матрицами строк, просто используйте: sum (sum (your_matrix == 5))
@AlexB. Правильное расширение, которое работает для матриц произвольных измерений, будет sum(your_matrix(:) == 5)
также для подсчета вхождений нескольких значений выполните sum(your_matrix == [ 1, 2, 3 ])