我希望找到
A(1:3,:)
的列向量,其在M(4,:)
中对应的值不是单元格X
的其中一个向量的一部分(显然也不等于这些向量之一)。如果X
非常大,我希望找到一种快速的解决方案。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];
这里直接取 A
:
A = [1007 4044 5002 622;
552 300 431 124;
2010 1113 1100 88];
A
包含 M(1:3,:)
的唯一列向量
X = {[2 5 68 44],[2 10 55 9 17],[1 55 6 7 8 9],[32 12]};
[~, ~, subs] = unique(M(1:3,:)','rows');
A4 = accumarray(subs(:),M(4,:).',[],@(x) {x});
%// 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
由于[2;10;55]
包含在X{2} = [2 10 55 9 17]
中,因此列向量[5002;431;1100]
被排除。
由于[32 12] = X{4}
,因此列向量[622;124;88]
被排除。
另一个例子:使用相同的X
M = [1007 4044 1007 4044 1007 5002 5002 5002 622 622 1007 1007 1007;
552 300 552 300 552 431 431 431 124 124 552 11 11;
2010 1113 2010 1113 2010 1100 1100 1100 88 88 2010 20 20;
12 25 15 12 30 2 10 55 32 12 7 12 7];
X = {[2 5 68 44],[2 10 55 9 17],[1 55 6 7 8 9],[32 12]};
A = [1007 4044 5002 622 1077;
552 300 431 124 11;
2010 1113 1100 88 20];
结果: (带有scmg答案)
如果按照第一行对 A
进行排序,我得到了正确的结果:
out =
1007 1007 4044
11 552 300
20 2010 1113
如果我不对矩阵
A
进行排序,结果会是错误的:(错误的结果)out =
4044 5002 622
300 431 124
1113 1100 88
应该消除列向量A(:,4) = [622;124;88]
,因为[32 12] = X{4}
。
应该消除列向量[5002;431;1100]
,因为[2;10;55]
包含在X{2} = [2 10 55 9 17]
中。