1. 程式人生 > >隨機打亂和逆過程

隨機打亂和逆過程

實驗中,我們經常需要對樣本進行隨機打亂,一般情況下,首先生成一個隨機序列,然後通過下標索引的方式得到打亂後的序列。

A 表示原樣本,B表示打亂的樣本,通過 rand_ind,我們可以知道:B ==> A (B到A的對映)。
而現在我們想知道 A ==> B 的對映。


clear all;

A = [3 4 5 6 1 2 ];

disp('A:');
disp(A);

if ~exist('tst_rand.mat','file')
    ind = randperm(length(A));
    save('tst_rand.mat','ind');
else
    load tst_rand.
mat; end B = A(ind); disp('B: '); disp(B); % create map map = zeros(length(A),2); reverse_map = zeros(length(A),2); for i = 1:length(ind) map(i,1) = i; map(i,2) = ind(i); % B ==> A end for i = 1:length(ind) reverse_map(i,1) = map(i,2); reverse_map(i,2) = map(i,1); % A == >B
end ss = sortrows(reverse_map,1); re_ind = ss(:,2); re_A = B(re_ind); disp('re_A:'); disp(re_A);

結果:

A:
     3     4     5     6     1     2

B: 
     4     1     2     6     5     3

re_A:
     3     4     5     6     1     2