Permutation sans répétition

Version imprimable

13/01/2009, 17h31
vincent_LSCP

Permutation sans répétition
Bonjour,

J'ai un petit problème sous matlab. J'ai ce vecteur:

Code:

x = {'-1' '-0.5' '0' '0.5' '1' '-1' '-0.5' '0' '0.5' '1' '-1' '-0.5' '0' '0.5' '1'};

Je fais une permutation aléatoire des éléments:

Code:

y = x(randperm(length(x)))

Le problème, c'est que dans ce vecteur y, il y a souvent deux valeurs identiques qui se suivent. Par exemple:
Code:

1 2 3 4 5 6 7 8 9 y = Columns 1 through 9 '1' '0.5' '1' '-1' '-1' '-0.5' '-0.5' '0.5' '-0.5' Columns 10 through 15 '0' '0' '1' '0' '-1' '0.5'
Les valeurs -1, -0.5, et 0 se suivent apparaissent deux fois consécutivement.

Avez-vous une idée d'une méthode pour éviter cela?

Merci,
Vincent
13/01/2009, 17h42
Jerome Briot

Tu peux nous en dire plus sur le vecteur x ?

Quel phénomène physique représente-t-il ?
Comment est-il généré ?
Pourquoi utiliser un tableau de cellules pour stocker un vecteur ?
Est-ce normale que ce soit toujours les mêmes valeurs qui soient stockée dedans ?
13/01/2009, 18h08
vincent_LSCP

Donc le vecteur x représente des valeurs d'abscisses, il y a 5 valeurs possibles (-1 -0.5 0 0.5 1). Chaque valeur est présente 3 fois dans le vecteur.

J'utilise un format cellule pour pouvoir faire la permutation: obtenir un ordre aléatoire des 15 valeurs. L'ordre est ici primordial car l'ordre d'apparition des valeurs dans le vecteur correspond à l'ordre temporel. Ma question vient du fait que je ne veux pas que 2 objets consécutifs apparaissent à la même position.

Une première méthode un peu brutale ( et pas forcément très robuste) :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
x = {'-1' '-0.5' '0' '0.5' '1' '-1' '-0.5' '0' '0.5' '1' '-1' '-0.5' '0' '0.5' '1'};
 
nx = numel(x);
 
xnum = str2double(x)
 
idx = true;
k = 0;
while (idx == 1) && k<=10
    idxperm = randperm(nx);
    y = xnum(idxperm);
    d = diff(y);
    idx = any(d==0);
    k = k+1;
end
 
if k>10
    disp('Mauvaise permutation')
else
    X = x(idxperm)
end

13/01/2009, 19h21
vincent_LSCP

Super, merci. J'ai trouvé une solution de mon côté aussi, pas très esthétique (je la poste dès que c'est amélioré).
-Vincent