向量化“for”循环

Question

这段代码就像我想要的那样工作，但是本着好的MATLAB代码的精神，有没有一种方法来向量化它(之前的是一个k x 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

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

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

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

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

Answer 2

使用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 3

我会这样做：

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

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