隨機打亂和逆過程
阿新 • • 發佈:2019-01-06
實驗中,我們經常需要對樣本進行隨機打亂,一般情況下,首先生成一個隨機序列,然後通過下標索引的方式得到打亂後的序列。
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