先看看subs和val的具体内容
subs = [1 1 1; 2 1 2; 2 3 2; 2 1 2; 2 3 2];
subs =
val = (101:105)'
val =
val的元素个数与subs的行数是一致的。
A = accumarray(subs, val)的实现过程分成2步。
第二步,把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的其他位置都为零。输出A就得到
A(:,:,1) =
A(:,:,2) =
很显然,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; ...]