3

我有一个像这样的矩阵:

A=
    10 31 32 22
    32 35 52 77
    68 42 84 32

我需要一个类似mode但有范围的函数,例如返回 30 的 mymode(A,10),在 0-10、10-20、20-30 范围内找到最频繁的数字,并返回范围内的大多数数字。

4

4 回答 4

4

您可以使用histc将数据分箱到您想要的范围内,然后max在输出中找到成员最多的箱histc

ranges = 0:10:50;                 % your desired ranges
[n, bins] = histc(A(:), ranges);  % bin the data
[v,i] = max(n);                   % find the bin with most occurrences

[ranges(i) ranges(i+1)]           % edges of the most frequent bin

对于您的具体示例,这将返回

ans =

    30    40

这与您所需的输出相匹配,因为 A 中的大多数值介于 30 到 40 之间。

于 2012-09-14T07:24:18.943 回答
1
[M,F] = mode(  A((A>=2) & (A<=5))  ) %//only interested in range 2 to 5

...其中 M 将为您提供模式,F 将为您提供出现频率

于 2012-09-14T07:15:37.223 回答
1
> A = [10 31 32 22; 32 35 52 77; 68 42 84 32]
A =

   10   31   32   22
   32   35   52   77
   68   42   84   32

> min = 10
min =  10
> max = 40
max =  40
> mode(A(A >= min & A <= max))
ans =  32
> 
于 2012-09-14T07:17:15.537 回答
1

我猜根据不同答案的数量,我们可能会错过您的目标。这是我的解释。

如果您想要有多个范围并且希望为每个范围输出最频繁的数字,请创建一个包含所有所需范围的单元格(它们可能重叠)并使用 cellfun 为每个范围运行 mode() 。您还可以使用 arrayfun 以类似的方式创建具有所需范围的单元格:

A = [10 31 32 22; 32 35 52 77; 68 42 84 32];

% create ranges
range_step = 10;
range_start=[0:range_step:40];
range=arrayfun(@(r)([r r+range_step]), range_start, 'UniformOutput', false)

% analyze ranges
o = cellfun(@(r)(mode(A(A>=r(1) & A<=r(2)))), range, 'UniformOutput', false)

o = 
[10]    [10]    [22]    [32]    [42]
于 2012-09-14T08:04:07.927 回答