matlab矩阵查找,在matlab矩阵中查找子矩阵的一般方法

这是低性能,但(据称)任意维度函数.它使用find来创建总计潜在匹配位置的(线性)索引列表,然后只检查总大小的匹配子块是否匹配.

function loc = matrixFind(total, sub)

%matrixFind find position of array in another array

% initialize result

loc = [];

% pre-check: do all elements of sub exist in total?

elements_in_both = intersect(sub(:), total(:));

if numel(elements_in_both) < numel(unique(sub))

% if not, return nothing

return

end

% select a pivot element

% Improvement: use least common element in total for less iterations

pivot_element = sub(1);

% determine linear index of all occurences of pivot_elemnent in total

starting_positions = find(total == pivot_element);

% prepare cell arrays for variable length subscript vectors

[subscripts, subscript_ranges] = deal(cell([1, ndims(total)]));

for k = 1:length(starting_positi