　ここで、k番目までは昇順に並んでいるとします。s=kとして、kをsにも格納します。k番目の要素を最小値minとして、k+1番目から最後の要素まで、minより小さいかを比較します。もし、n番目の要素がminよりも小さい場合、min=a(n)として、minを置き換えます。また、s=nとして、k番目の要素位置sをnと置き換えます。最後に、t=a(k); a(k)=a(s); a(s)=t;として、k番目の要素と、最小であったn番目の要素を置き換えます。kを最初の要素番号1から最後の要素の手前まで変えながら、以上の操作を繰り返していくと、要素は昇順に並びます。

　また、コマンドticとtocをforループの前後に置いて、計算に要した時間を計測しています。

clear;
a=randperm(1000);b=a;
tic;
for k=1:length(a)-1
    min=a(k);
    s=k;
    for n=k+1:length(a)
        if(a(n)<min)
            min=a(n);s=n;
        end
    end
    t=a(k);a(k)=a(s);a(s)=t;
end
time1=toc;
disp(['The calculation needs ',num2str(time1),'second']);

ex317.m

＞＞「ex317.m」ダウンロード

　FreeMatでは「sort」コマンドで、昇順に並べ替えることができます。sortコマンドを使って、ex317.mのforループを置き換えたものが下記です。とても簡単です。

tic;
sort(b);
time2=toc;
disp(['The calculation needs ',num2str(time2),'second']);

　双方の計算効率を調べてみましょう。ex317.mに上記のプログラムを付け加えて動作させると、例えば下記のような結果が得られます。

The calculation needs 5.238second
The calculation needs 0.001second

　forループによる処理では、5.238秒で、sortによる処理では、0.001秒であることが分かります。先ほどの探索アルゴリズムでの比較よりも大きな差が生じているのは、2重のループを使っているためです。先のアルゴリズムは、初歩的なアルゴリズムで効率があまり良くないとはいえ、sortコマンドを使った方が、簡単で計算が速いことが分かると思います。なお、これは数値をランダムに並べたものから該当の数値を見つけ出すプログラムであり、ハードウェアの処理能力の違いなども考えられるため、同様の計算を行っても、同じ結果にはなりませんのであらかじめご了承ください。

　ちなみに、sortコマンドは、文字列をソートすることも可能です。例えば、sort('How much wood a woodstock chuck?')とすると、?Haccccddhhkkmoooooostuuwwwと、特殊記号、大文字アルファベット、小文字アルファベットの順に並び変わります。

　セル配列のソートも同様です。例えば、sort({'abc','def','ghi','jkl'},2,'descend')とすると、[jkl] [ghi] [def] [abc]と、文字列がアルファベット順に降順に並び変わっています。

　マージは配列の結合とソートの組み合わせで実現できます。例えば、ab=randperm(1000); a=ab(1:500); b=ab(501:1000);とすると、1～1000までの整数がランダムに並んだ2つの配列aとbが出来上がります。plot(a,'o',b,'*')とすると、図1のように、○で示すaの要素と＊で示すbの要素とがランダムに並んでいることが分かります。