matlab的val,Matlab的accumarray(subs, val) 解释

先看看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； ...]