在MATLAB中找到互补向量的快捷方法

时间:2017-12-03 16:22:02

标签: algorithm matlab combinations

我有一排N行二进制向量,即

mymatrix = [ 1 0 0 1 0;
             1 1 0 0 1;
             0 1 1 0 1;
             0 1 0 0 1;
             0 0 1 0 0;
             0 0 1 1 0;
             ....       ]

我想在哪里找到行的组合,当它们加在一起时,会让我完全理解:

[1 1 1 1 1]

因此,在上面的示例中,可行的组合是1/31/4/52/6

我现在的代码是:

i = 1;
for j = 1:5
    C = combnk([1:N],j); % Get every possible combination of rows
    for c = 1:size(C,1)
        if isequal(ones(1,5),sum(mymatrix(C(c,:),:)))
            combis{i} = C(c,:);
            i = i+1;
        end
    end
end

但正如您想象的那样,这需要一段时间,特别是因为那里的combnk

什么可能是一个有用的算法/功能,可以帮助我加快速度?

1 个答案:

答案 0 :(得分:1)

M = [
 1 0 0 1 0;
 1 1 0 0 1;
 0 1 1 0 1;
 0 1 0 0 1;
 0 0 1 0 0;
 0 0 1 1 0;
 1 1 1 1 1
];

% Find all the unique combinations of rows...
S = (dec2bin(1:2^size(M,1)-1) == '1');

% Find the matching combinations...
matches = cell(0,1);

for i = 1:size(S,1)
    S_curr = S(i,:);

    rows = M(S_curr,:);
    rows_sum = sum(rows,1);

    if (all(rows_sum == 1))
        matches = [matches; {find(S_curr)}];
    end
end

以良好的程式化方式显示您的比赛:

for i = 1:numel(matches)
    match = matches{i};

    if (numel(match) == 1)
        disp(['Match found for row: ' mat2str(match) '.']);
    else
        disp(['Match found for rows: ' mat2str(match) '.']);
    end
end

这将产生:

  

找到匹配的行:7。

     

找到行匹配:[2 6]。

     

找到行的匹配项:[1 4 5]。

     

找到行匹配:[1 3]。

就效率而言,在我的机器中,这个算法正在完成大约2 milliseconds的匹配检测。