MATLAB accumarray的使用

先看看subs和val的具体内容

subs = [1 1 1; 2 1 2; 2 3 2; 2 1 2; 2 3 2];

subs =

     1     1     1

     2     1     2

     2     3     2

     2     1     2

     2     3     2

val = (101:105)'

val =

   101

   102

   103

   104

   105

val的元素个数与subs的行数是一致的。

A = accumarray(subs, val)的实现过程分成2步。

第一步是把val中的元素，按照subs对应行所给出的下标放到一个新的cell矩阵B中（cell是为了方便解释，也就是说B矩阵中的每个位置可以放入多个数值），注意，subs的值是B的下标，不是val的。举例来说，subs第一行[ 1     1     1],意思就是把val中第一个元素（val（1））放入到B（1，1，1）的位置，依次类推，val（2）放入到B(2     1     2)，val（3）放入到B(2     3    2)，val（4）放入到B(2     1     2)，val（5）放入到B(2    3   2)。此时，可以看到B（1，1，1）中有1个数（val（1））；B(2     1    2)有2个数（val（2），val（4））；B(2    3    2)也有2个数（val（3），val（5））。

第二步，把B中每个单元中的数分别累加，并放入到A的对应位置。

注：accumarray默认的是把每个单元中的数累加，因为对每个单元中的数的默认处理函数是sum。可以通过A = accumarray(subs,val,[],@fun)的调用格式来指定其他的处理函数，比如说mean。对指定的fun函数的要求是，接受列向量输入，输出单个的数值型，、字符型或逻辑型变量。

A的维数与B相同，A中的元素默认为零。A的大小为max(subs(1))×max(subs(2))×max(subs(3))...

因此，

A（1，1，1）=val（1）=101；

A (2 , 1 , 2)=val（2）+val（4）=102+104=206;

A (2, 3 ,  2)=val（3）+val（5）=1-03+105=208,

A的其他位置都为零。输出A就得到

A(:,:,1) =

   101     0     0

     0     0     0

A(:,:,2) =

     0     0     0

   206     0   208

很显然，A的维数与subs的列数相等。

A = accumarray(subs, val)

A = accumarray(subs,val,sz)

sz 可以用来指定A大小，但是不能小于A = accumarray(subs, val)得到的A的大小。比如A = accumarray(subs, val)的到A是一个3×4的二维矩阵，那么sz应当为一个包含2个元素的向量sz=[m1,m2] （sz向量的长度和A的维数相等），其中，m1大于等于3，m2大于等于4. 但是，当得到的A是一个p×1的一维向量时，sz=[m，1]，m大于等于p。另外，sz可以赋值为空，表示由函数自动决定A的大小。

A = accumarray(subs,val,sz,fun)

fun可以指定专门的处理函数，默认的处理函数为sum

A = accumarray(subs,val,sz,fun,fillval)

fillval指定A中元素的默认值。可以等于NaN

A = accumarray(subs,val,sz,fun,fillval,issparse)

isspares选择A是否使用稀疏矩阵的格式

A = accumarray({subs1, subs2, ...}, val,...)

{subs1, subs2, ...}，等同于A = accumarray(subs, val,...)，此时，subs=[subs1, subs2, ...]或者=[subs1；subs2； ...]