Sort Index - Bug

Question

Rainer Ng 2019-9-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481239-sort-index-bug

评论： Walter Roberson 2021-4-14

[a, sort_index]=sort([100;20;30;40])

The sort_index should return

But it does not.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2019-9-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481239-sort-index-bug#answer_392692

>> [a, sort_index]=sort([100;20;30;40])
a =
    20
    30
    40
   100
sort_index =
     2
     3
     4
     1

This output is correct. It tells you that a(1) == 20 came from element #2 of the input, and that a(2) == 30 came from element 3 in the input, and that a(3) == 40 came from element 4 of the input, and that a(4) == 100 came from element 1 of the input.

To get the output you want, use

b(sort_index) = 1:length(sort_index);

4 个评论
显示 2更早的评论隐藏 2更早的评论

Walter Roberson 2021-4-13

test=[1 5;6 2;1 2]
test = 3×2
     1     5
     6     2
     1     2
 [testgeo,testidx]=sort(test,1)
testgeo = 3×2
     1     2
     1     2
     6     5
testidx = 3×2
     1     2
     3     3
     2     1
 test(testidx(:,2),2)
ans = 3×1
     2
     2
     5
 test([3 2 1],2)
ans = 3×1
     2
     2
     5

This tells us that for the second column, [2 3 1] and [3 2 1] lead to the same output for sorting purposes, and so both are correct.

When there are multiple items with the same value, there is the question of which order sort algorithms should output. Some algorithms output the items in the order they happen to encounter them during sorting -- and especially with parallel sorts, or sort/merge algorithms, that can lead to later copies of the identical values happening to get sorted before the earlier copies of the value. Other algorithms take the time to order identical values according to the original order; those are generally referred to as "stable" algorithms. Algorithms that are not "stable" have the issue that if you sort an already sorted array, then the sort order might not be just the indices, but that never happens for "stable" sort algorithms.

(You can see "stable" as an option for unique() where it has a slightly different meaning https://www.mathworks.com/help/matlab/ref/double.unique.html#bs_6vpd-1-setOrder )

Even in algorithms that are not "stable", there is no reason why [3 2 1] would necessarily have to be output instead of [2 3 1] for the indices: it could just happen to return [2 3 1] instead.

Franziska Höppe 2021-4-13

Thank you for you detailed answer! :-)

The issue is if you let run the loop, you get a different result. You get:

2

5

2

Do you know why?

Walter Roberson 2021-4-14

I never got that output; which code were you using, and which MATLAB version?

请先登录，再进行评论。

Answer 2

James Tursa 2019-9-19

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481239-sort-index-bug#answer_392690

The sort index gives the location in the original array of the sorted values. I.e., the sort results "a" are "original_array(sort_index)"

>> x = [100;20;30;40];
>> [a,sort_index] = sort(x)
a =
    20
    30
    40
   100
sort_index =
     2
     3
     4
     1
>> x(sort_index)
ans =
    20
    30
    40
   100
>> isequal(x(sort_index),a)
ans =
     1