向量化“for”循环

Question

这段代码按照我想要的方式工作，但是本着良好的 MATLAB 代码的精神，有没有一种方法可以对其进行向量化（之前是 akx 1 向量）：

start = zeros(k,1);
for i = 2:length(previous)
    if (previous(i-1) == -1)
        start(previous(i))= start(previous(i))+1;
    end    
end

一般来说，什么是对代码进行向量化的直观方法MATLAB？

Answer 1

使用 MATLAB 中的 find 命令，该命令返回向量为 TRUE 的索引 (i)。所以：

% precache indices i where previous(i-1) == 1
idx = find(previous==-1)+1;
start(previous(idx)) = start(previous(idx))+1;

我预缓存 idx 的原因是，如果 previous 是一个大向量，并且执行 find 需要一段时间。否则你可以这样做

start( find(previous==-1)+1 ) = start( find(previous==-1) ) + 1;

Answer 2

您可以在不使用 find 的情况下执行此操作，以获得最佳性能：

I = [false; previous(1:end-1) == -1];
idx = previous(I);
start(idx) = start(idx) + 1;

这还可以避免 previous(end) == -1 的风险，这会导致索引超出范围替代方案中的错误。

请注意，如果 idx 包含重复索引，则其工作方式与原始索引不同。

Answer 3

我会做类似的事情：

prev_is_minus1 = [false; previous(1:end-1)==-1]
start = accumarray(previous(prev_is_minus1), 1, k)

我相信这与您发布的代码循环具有相同的效果。