矩阵操作以提取某些子列

Question

M = [1007  1007  4044  1007  4044  1007  5002 5002 5002 622 622;
      552   552   300   552   300   552   431  431  431 124 124 ; 
     2010  2010  1113  2010  1113  2010  1100 1100 1100  88  88;
        7    12    25    15    12    30     2   10   55  32  12]

X = {[2 5 68 44],[2 10 55 9 17],[1 55 6 7 8 9],[32 12]}

A = [1007  4044  5002  622
      552   300   431  124
     2010  1113  1100   88
        7    25     2   32
       12    12    10   12
       15          55
       30                 ]

A 是一个解释我想要什么的实体。

A 包含M(1:3,:) 的唯一列向量，以及M(4,:) 中的对应值

A(1:3,:) = unique(M(1:3,:)','rows')'

我希望找到A(1:3,:) 的列向量，其在M(4,:) 中的对应值不属于单元格X 的向量之一（并且显然不等于这些向量之一）。

对于我的示例，所需的结果是矩阵：

 [1007  4044; 
   552   300; 
  2010  1113;]

列向量[5002;431;1100] 被删除，因为[2;10;55] 包含在X{2} = [2 10 55 9 17] 中

列向量[622;124;88] 被消除，因为[32 12] = X{4}

Answer 1

输入：

M = [1007  1007  4044  1007  4044  1007  5002 5002 5002 622 622;
      552   552   300   552   300   552   431  431  431 124 124; 
     2010  2010  1113  2010  1113  2010  1100 1100 1100  88  88;
        7    12    25    15    12    30     2   10   55  32  12];

X = {[2 5 68 44],[2 10 55 9 17],[1 55 6 7 8 9],[32 12]};

这样做（你做了什么）

A(1:3,:) = unique(M(1:3,:).','rows').';

给予：

>> A

A =

     622        1007        4044        5002
     124         552         300         431
      88        2010        1113        1100

然后使用unique 和accumarray

[~, ~, subs] = unique(M(1:3,:)','rows');

A4 = accumarray(subs(:),M(4,:).',[],@(x) {x});

现在我们将A4 作为单元阵列

>> A4

A4 = 

[2x1 double]    [4x1 double]    [2x1 double]    [3x1 double]

然后使用cellfun、ismember、all 和any

%// getting a mask of which columns we want
idxC(length(A4)) = false;
for ii = 1:length(A4)
    idxC(ii) = ~any(cellfun(@(x) all(ismember(A4{ii},x)), X));
end

显示我们想要的列

out = A(:,idxC)

结果：

>> out

out =

    1007        4044
     552         300
    2010        1113

我建议您按照@Dan 的建议自己尝试一下。如果你被困在某个地方，你可以参考这个。如果您有任何澄清/修改，请告诉我:)