清空筛选条件
清空筛选条件

I find a problem about accumarray​(subs,val,​[],@(x)sum​(diff(x)))​, maybe a bug

1 次查看（过去 30 天）
jacky chen
jacky chen 2014-1-7
评论： Roger Stafford 2014-1-12
in the example of function accumarray, when use the diff ,it's so weird that the index of the val ,sometime the input value of x is Positive order ,but sometime is Reverse order,(u can run the code ,then pay attention to the value of y(the input value of x)),in the B3，u can get y(the input of diff ) like this
y =
103
y =
104
106
y =
102
101
it's so weird that 104 106 and 102 101，why it's not 101 102 ,it will lead to a wrong diff .
the code followed :
clc;clear all;
val = [101 102 103 104 105 106];
subs=[1 2; 1 2; 3 1; 4 1; 4 4; 4 1];
B1 = accumarray(subs,val,[],@(x)sum(my(x)))
B2 = accumarray(subs,val,[],@(x)sum(my(diff(x))))
subs= [1 2; 3 1; 1 2; 4 4; 4 1; 4 1;];
B3 = accumarray(subs,val,[],@(x)sum(my(x)))
B4 = accumarray(subs,val,[],@(x)sum(my(diff(x))))
function [ y] = my( x )
y=x
end

请先登录，再回答此问题。

采纳的回答

Roger Stafford
Roger Stafford 2014-1-8
When the 'subs' values are not in what Mathworks considers as "sorted order", the 'accumarray' algorithm calls for some kind of sorting operation to place 'subs' in sorted order, and it subsequently reorders the 'val' vector accordingly. However, apparently the sorting operation which is used does not necessarily behave in the same way Mathworks' own 'sort' function does with respect to equal entries. For equal entries, matlab's 'sort' operation guarantees that the corresponding sort indices are a sequence of successively increasing integers. That is, if one does
[B,IX] = sort([7,7,4]);
the value for IX is guaranteed to be [3,1,2] and never [3,2,1]. Such is apparently not true for whatever sorting method Mathworks is using with 'accumarray'. Sometimes the index order for equal entries is reversed as apparently happens in their Example 2 in which the indexing is reversed for the two [1,2] equal entries, leading to a sequence of 'val' values of 102,101 instead of the expected 101,102. It would be interesting to find out just what sorting algorithm Mathworks is using in 'accumarray' that it produces such reversals for equal entries.

更多回答（2 个）

Sean de Wolski
Sean de Wolski 2014-1-8
编辑：Sean de Wolski 2014-1-8
Explicitly sort either subs or vals (inside of the function) when you need sorted values in accumarray()
Example two and the note discuss this briefly.
http://www.mathworks.com/help/releases/R2013b/matlab/ref/accumarray.html
  1 个评论
jacky chen
jacky chen 2014-1-8
编辑：jacky chen 2014-1-8
Thanks,,when I read the example two(in the help window) that I find the question ,in the example two I can understand B(4,1)=2,it comes from val(6)-val(4)=106-104=2, but why the B(1,2)= -1, why it's not val(2)-val(1)=102-101=1?

请先登录，再进行评论。

jacky chen
jacky chen 2014-1-12
Hi,,everyone ,finally ,I find the answer ,it's about the definition of the differential, when the number is the first one of a line(use L for a single line,the size is (n,1)),we can use it(Forward difference) as diff: L(1)=L(1)-L(2); and when the number is the last one of a line(L),we can use it(Backward difference) as diff: L(n)=L(n)-L(n-1); where other number of the line, the diff: L(m)=L(m)-L(m-1),(1<m<n). Jsut enjoying.
  1 个评论
Roger Stafford
Roger Stafford 2014-1-12
The problem is not with the 'diff' function, Jacky. Accumarray exhibits the same anomalous behavior with other functions. For example, someone on the internet wrote
r = accumarray(subs, vals', [], @(x){x'})
with the non-sorted subs = [2,2,1,1,1,1]' and val = [10,13,11,14,12,10]' with the result
r{1} = 12 11 14 10
r{2} = 13 10
which shows directly the strange orderings the 'accumarray' sorting procedure performs with equal values of 'subs'.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Matrices and Arrays 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by