22

行列がAあり、この行列の行を並べ替えるとします。マトリックス上で同じ順序を複製するにはどうすればよいですBか（もちろん同じサイズ）？

例えば

A = rand(3,4);
[val ind] = sort(A,2);
B = rand(3,4);
%// Reorder the elements of B according to the reordering of A

これは私が思いついた最高のものです

m = size(A,1);
B = B(bsxfun(@plus,(ind-1)*m,(1:m)'));

好奇心から、他に選択肢はありますか？

更新： 2008a（XP）でプロファイルされたJonasの優れたソリューション：

n = n

0.048524       1.4632       1.4791        1.195       1.0662        1.108       1.0082      0.96335      0.93155      0.90532      0.88976

n = 2m

0.63202       1.3029       1.1112       1.0501      0.94703      0.92847      0.90411       0.8849       0.8667      0.92098      0.85569

これは、 JITAのおかげで（おそらく）ループがMATLABプログラマーにとってもはや嫌悪感ではないことを示しています。

4

3 に答える 3

17

これを行うためのやや明確な方法は、ループを使用することです

A = rand(3,4);
B = rand(3,4);
[sortedA,ind] = sort(A,2);

for r = 1:size(A,1)
   B(r,:) = B(r,ind(r,:));
end

興味深いことに、ループバージョンは、小さい（<12行）および大きい（>〜700行）正方形アレイ（r2010a、OS X）の方が高速です。行に比べて列が多いほど、ループのパフォーマンスは向上します。

これが私がテストのためにすぐにハックしたコードです：

siz = 10:100:1010;
tt = zeros(100,2,length(siz));

for s = siz
    for k = 1:100

        A = rand(s,1*s);
        B = rand(s,1*s);
        [sortedA,ind] = sort(A,2);

        tic;
        for r = 1:size(A,1)
            B(r,:) = B(r,ind(r,:));
        end,tt(k,1,s==siz) = toc;

        tic;
        m = size(A,1);
        B = B(bsxfun(@plus,(ind-1)*m,(1:m).'));
        tt(k,2,s==siz) = toc;

    end
end

m = squeeze(mean(tt,1));

m(1,:)./m(2,:)

正方形配列の場合

ans =

    0.7149    2.1508    1.2203    1.4684    1.2339    1.1855    1.0212    1.0201    0.8770       0.8584    0.8405

行の2倍の列（同じ行数）の場合

ans =

    0.8431    1.2874    1.3550    1.1311    0.9979    0.9921    0.8263    0.7697    0.6856    0.7004    0.7314
于 2010-04-21T00:04:30.250 に答える
6

Sort（）は、並べ替えたディメンションに沿ったインデックスを返します。行を安定させたままにする他のディメンションのインデックスを明示的に作成してから、線形インデックスを使用して配列全体を再配置できます。

A = rand(3,4);
B = A; %// Start with same values so we can programmatically check result

[A2 ix2] = sort(A,2);
%// ix2 is the index along dimension 2, and we want dimension 1 to remain unchanged
ix1 = repmat([1:size(A,1)]', [1 size(A,2)]); %//'
%// Convert to linear index equivalent of the reordering of the sort() call
ix = sub2ind(size(A), ix1, ix2) 
%// And apply it
B2 = B(ix)
ok = isequal(A2, B2) %// confirm reordering
于 2010-04-21T13:22:50.130 に答える
0

あなたはこれをすることができませんか？

[val ind]=sort(A);
B=B(ind);

私があなたの問題を間違って理解していない限り、それは私のために働いた。

于 2012-06-05T17:03:36.667 に答える