替换为这个MATLAB脚本中的循环

Question

有没有什么办法可以避免在下面的代码中使用循环？ 我正在做一个减少总和，即分组a和b的所有相同指数的值。

a = [ 1 3 2 3 1 2] 
b = [0.1 0.2 0.3 0.4 0.5 0.6] 
for i = unique(a) 
    c(i) = sum(b(ismember(a, i))) 
end 

Answer 1

这正是*意味着什么accumarray为：

>> accumarray(a.',b.') 
ans = 

    0.60000 
    0.90000 
    0.60000 

* _{其实，这是最简单的件事accumarray一样。它默认为@sum，但您可以指定自己的函数来告诉它如何处理累计索引。}

Answer 2

这是一个非常昂贵避免环路的：

rsum = b*(bsxfun(@eq,a.',unique(a))); 

的总和在a的元件的增加的顺序给出。此外，a未在b索引的情况下适用。