0

我有一个 ND 数组，对于数组中的每个元素，我需要找到它下面的向量中最大元素的索引。我每次都在做这个，所以我真的对它尽可能快地超级感兴趣。

我写了一个函数locate，我用一些有代表性的示例数据调用它。（我使用arrayfun增加timeit运行该函数的次数来最小化随机波动。）

Xmin = 5;
Xmax = 300;
Xn = 40;
X = linspace(Xmin, Xmax, Xn)';
% iters = 1000;
% timeit(@() arrayfun(@(iter) locate(randi(Xmax + 10, 5, 6, 6), X), 1:iters, 'UniformOutput', false))
timeit(@() locate(randi(Xmax + 10, 5, 6, 6), X))

我的原始版本locate如下所示：

function indices = locate(x, X)
    % Preallocate
    indices = ones(size(x));

    % Find indices
    for ix = 1:numel(x)
        if x(ix) <= X(1)
            indices(ix) = 1;
        elseif x(ix) >= X(end)
            indices(ix) = length(X) - 1;
        else
            indices(ix) = find(X <= x(ix), 1, 'last');
        end
    end
end

我可以召集的最快版本如下所示：

function indices = locate(x, X)
    % Preallocate
    indices = ones(size(x));

    % Find indices
    % indices(X(1) > x) = 1;  % No need as indices are initialized to 1
    for iX = 1:length(X) - 1
        indices(X(iX) <= x & X(iX + 1) > x) = iX;
    end
    indices(X(iX) <= x) = length(X) - 1;
end

你能想到其他更快的方法吗？

4

1 回答 1

1

你需要的基本上是第二个输出histc

[pos,bin]=histc(magic(5),X);
bin(bin==0)=1;
于 2020-04-04T13:51:02.477 回答