Discussion :

[vincent_LSCP]

Bonjour,

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[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 ?

[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.

[Jerome Briot]

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[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