找到相交多阵列

Question

的。如果我有几个数组：

A = [7 1 7 7 4]; 

B = [7 0 4 4 0]; 

C = [0 0 4 1 5]; 

D = [5 7 2 4 0]; 

我知道在Matlab“相交”可以找到两个矩阵之间的共享内容与他们的索引。如果我想将它用于四个矩阵，我该怎么做？

注：这可以用于两个矩阵为： [K，IA，IB] =相交（A，B）

http://uk.mathworks.com/help/matlab/ref/intersect.html

Answer 1

你可以串联所有输入阵列（矢量）转换成2D数组，然后尝试查找所有输入数组中存在的唯一元素。接下来介绍的bsxfun基于代码试图达到相同的 -

%// Concatenate all vector arrays into a 2D array 
M = cat(1,A,B,C,D) 

%// Find unique values for all elements in all arrays 
unqvals = unique(M(:),'stable')' %//' 

%// Find which unqiue elements are common across all arrays, which is the 
%// desired output 
out = unqvals(all(any(bsxfun(@eq,M,permute(unqvals,[1 3 2])),2),1)) 

代码输出 -

M = 
    7  1  7  7  4 
    7  0  4  4  0 
    0  0  4  1  5 
    5  7  2  4  0 
unqvals = 
    7  1  4  0  5  2 
out = 
    4 

为了验证针对intersect基于代码，它的一种形式是这样的 -

out1 = intersect(intersect(intersect(A,B,'stable'),C,'stable'),D,'stable') 

对于给定的输入，它会给 -

out1 = 
    4 

为了进一步验证，假设你介绍一个7为C即C = [0 7 4 1 5]，使得在所有输入数组可用7，你将有输出[7 4]。

---

如果你想使bsxfun工作与2D阵列，可以更高效存储，这里是一个另类 -

%// Concatenate all vector arrays into a 2D array 
M = cat(1,A,B,C,D) 

%// Find unique values for all elements in all arrays 
unqvals = unique(M(:),'stable')' %//' 

[m1,m2] = size(M) %// Get size of concatenated 2D array 

%// Matches for all elements in all arrays against the unique elements 
matches = bsxfun(@eq,reshape(M',[],1),unqvals) %//' 

%// Desired output 
out = unqvals(all(any(permute(reshape(matches,m2,m1,[]),[1 3 2]),1),3))