_{碰到函数accumarray是在rocwood的关于寻找完全数的程序中。没找到M文件，Document中有语法分析，没看懂。结合例子，作简单说明。}

_Syntax

_{A = accumarray(sub, val)}

_{A = accumarray(sub, val, sz)}

_{A = accumarray(sub, val, sz, fun)}

_{A = accumarray(sub, val, sz, fun, fillvalue)}

  sub：提供累计信息的指示向量

  val：提供累计数值的向量

  sz：控制输出向量A的size

  fun：用于计算累计后向量的函数，默认为@sum，即累加

  fillvalues: 填补A中的空缺值，默认为0

_{Matlab的Document 里这么说：}

_{accumarray groups elements from a data set and applies a function to each group. A = accumarray(subs,val) creates an array A by accumulating elements of the vector val using the elements of subs as indices. The position of an element in subs determines which value of vals it selects for the accumulated vector; the value of an element in subs determines the position of the accumulated vector in the output.}

_{对于A = accumarray(subs,val)这么一个调用，有这么几个问题，理解清楚，就理解了这个函数。}

_例子：

_{val    =     [     1     2     3     4     5    6]}

_{subs =     [     1     2     4     2     4    3]'  % subs是列向量}

_{Q: accumarray总体是干嘛的？}

_{A: 笼统的说，是用subs向量中的信息从val中提取数值做累加，累加完的结果放到A中。}

_{Q: subs是干嘛的？}

_{A: subs是一个累加指示向量。}

_{subs提供的信息由两个：}

_{(a). subs向量中的每个位置对应val的每个位置；}

_{(b). subs中元素值相同的，val中的对应元素累加，元素值是累加完后放到A的什么地方。}

       _{如：上面的例子中，subs(2),subs(4)都是2，所以，val(2)和val(4)累加起来，放到A(2)这个位置上。}

_{Q: val是干嘛的？}

_{A: val是提供累加数值的，谁累加呢？就是A中的数值累加。选哪些数进行累加呢？subs向量中数值相同的对应位置的数。累加完后放到哪里呢？放到subs中指示的位置。}

_{Q: A是怎么出来的？A的维度是什么？A的内容如何确定？}

_{A: A的维度是subs中表示维度的数值最大的那个，如例子中size(A,1)==4，因为max(subs)==4。当然，这只是一维的情况。}

_{最后A的结果就是：}

_A =  

     _{1   % subs(1)==1，所以，A(1) = val(1)。}  

     _{6   % subs(2)==subs(4)==2，所以，A(2)=val(2)+val(4)}  

     _{6   % subs(6)==3的值， A(3) = Val (6) = 6} 

    _{8   % subs(3)==subs(5)==4，所以，A(4)=val(3)+val(5)}  

_{产生一个2*3*4的矩阵}

ind = [1 1 1; 2 1 2; 2 3 4; 2 3 4];

A = accumarray(ind,11:14)

A(:,:,1) =

     11     0     0

     0     0     0

A(:,:,2) =

      0     0     0

    12     0     0

A(:,:,3) =

     0     0     0

     0     0     0

A(:,:,4) =

      0     0     0

     0     0    27

_{另外，关于二维和高维矩阵，只是大概知道A矩阵的维度，数值的变化不详细得知。}

_References:

_{[1] accumarray - The MathWorks http://www.mathworks.com/help/techdoc/ref/accumarray .html}

_{[2] 向量化操作的又一重要函数accumarray的用法总结http://www.simwe.com/forum/thread-811616-1-3.html}

_{[3] rocwood关于寻找完全数的程序}

_第一种：

_{clear;clc;close all}

_a=1:10000;

_{a(2*a'==accumarray(a',a',[],@(x)sum(a(abs(x./(1:x)-round(x./(1:x)))<eps))));}