1

我有矩阵 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

S 应该是：

 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) 中具有相同的值）。

谢谢

4

2 回答 2

4

这是一种方法：

[~, ~, 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
于 2013-09-10T06:52:10.787 回答
0

尝试以下操作：

% 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]
于 2013-09-10T07:21:00.910 回答