Question

我有矩阵X，mX2，我想得到一个大小为

的矩阵S.

size(unique(X(:,2),1) X size(unique(X(:,2),1)

对于每个S（i，j）我想计算i，j出现在一起的次数。例如：

X = [1 11 ; 
     2 11; 
     3 11; 
     5 23; 
     6 23;
     1 23; 
     9 24;
     9 25;
     3 23;
    10 23]
unique(X(:,2)) 
11
23
24
25

应该是：

 0     2     0     0
 0     0     0     0
 0     0     0     1
 0     0     0     0

（我不关心对角线，它可能有或没有它们，S也可以是对称的。）

S(1,2) = 2

因为11和23（位于1,2位）一起出现两次（即在X（：，1）中具有相同的值）。

由于

Answer 1

这是一种方法：

[~, ~, n1] = unique(X(:,1));
[~, ~, n2] = unique(X(:,2));
B = accumarray([n2 n1],1);
S = B*B';

这给出了完整的矩阵：

>> S
S =
     3     2     0     0
     2     5     0     0
     0     0     1     1
     0     0     1     1

要删除对角线和下三角形，您可以使用

S = triu(B*B',1);

产生

>> S
S =
     0     2     0     0
     0     0     0     0
     0     0     0     1
     0     0     0     0

Answer 2

尝试以下方法：

% convert each columns to indices starting from 1
[a,~,aa] = unique(X(:,1));
[b,~,bb] = unique(X(:,2));

% group occurences of col2 according to values of col1
C = accumarray(aa, bb, [], @(x){x});

% keep only occurences of two or more values
C = C(cellfun(@numel,C) > 1);

% in case of three or more values co-occured, generate all pairs
C = cellfun(@(v) nchoosek(v,2), C, 'UniformOutput',false);

% concatenate all pairs
C = cell2mat(C);

% build count matrix
C = sparse(C(:,[1 2]), C(:,[2 1]), 1);
C = full(C);

这种情况下的结果（显然是对称矩阵）：

>> C
C =
     0     2     0     0
     2     0     0     0
     0     0     0     1
     0     0     1     0

或漂亮打印的行/列标题：

>> [{[]} num2cell(b'); num2cell(b) num2cell(C)]
ans = 
      []    [11]    [23]    [24]    [25]
    [11]    [ 0]    [ 2]    [ 0]    [ 0]
    [23]    [ 2]    [ 0]    [ 0]    [ 0]
    [24]    [ 0]    [ 0]    [ 0]    [ 1]
    [25]    [ 0]    [ 0]    [ 1]    [ 0]

计算矩阵matlab中对的频率

2 个答案: