Discussion :

[pedrosanchau]

Bonjour,
J'aimerais savoir comment effectuer un tirage sans remise.
Par exemple:

Code :

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

1
2
 a=   [1   2   3   4   5
6   7   8   9   10]

Dans cette matrice, je voudrais pouvoir selectionner un des éléments au hasard qui lors d'un deuxième tirage au hasard ne pourrait pas être sélectionner.
Je pensais que la commande 1:a marcherait mais je viens de découvrir que non.

Merci.

[Caro-Line]

Pour la sélection "au hasard" je pense qu'il faut que tu regardes du côté de la fonction RAND, à adapter à ton problème évidemment.
Après pour le 2ème tirage, il faut enlever l'élément déjà pris au 1er tirage.
As-tu vraiment besoin que ce soit une matrice ou peux-tu la transformer en vecteur (car sinon on ne peut enlever simplement un élément) ?

[mr_samurai]

il y a la fonction randperm

Code :

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

1
2
3
4
>> randperm(6)
ans =
     2     4     3     6     5     1 
>>

C'est une permutation des 1:6 .

Pour ton cas, un truc du style :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
A =  [1   2   3   4   5 ; 6   7   8   9   10];
 
% Generer les indices
r = randperm(size(A,1)*size(A,2));
 
% Prendre le i-eme
 
tempo = A(r(i)); % indexage linéaire ici 
 
% Plus encore : 
Aperm = A(r);  % Aperm est un vecteur qui contiens tous les elements de la matrice A permuté.
 
%  ... utilise la fonction  'reshape' si tu veux retrouver un matric a partir de  Aperm

[pedrosanchau]

Non, pas du tout, ça peut être un vecteur mais je ne sais pas vraiment comment les gérer, j'ai toujours utiliser des matrices.

[Jerome Briot]

Code :

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

1
2
3
 
% Generer les indices
r = randperm(size(A,1)*size(A,2));

Il est plus judicieux d'utiliser la fonction NUMEL ici :

Code :

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

1
2
% Generer les indices
r = randperm(numel(A));

Note : NUMEL évite aussi l'emploi de PROD(SIZE(A))

[mr_samurai]

j'avais changé juste apres ... je teste mes codes avant de les poster.

Code :

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

numel(A)

bien pratique celle la .

[pedrosanchau]

Utiliser le randperm est une bonne idée pour générer le hasard. Il suffirait ensuite de choisir le premier chiffre lors du tirage, le deuxième lors du deuxième tirage et ainsi de suite. Est ce possible d'intégrer cela dans une boucle car j'ai 60 tirages?
Cela m'éviterais de devoir procéder ainsi: "premier tirage, essai, deuxième tirage, essai..." ce qui ferait un programme fastidieux a ècrire.

[mr_samurai]

60 tirages ?

60 experiences de tirage (1) , ou tu tire 60 element de la matrice (2)?

pour (1) :

Code :

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

1
2
3
4
5
6
7
8
tirage = 0;
 
while tirage < 60 
     r = randperm(numel(A));
     % corps de l'expérience
     % extraction resultats
     tirage = tirage + 1;
end

pour (2) :

Code :

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

1
2
3
4
5
6
7
8
% ....
r = randperm(numel(A));
 
for i=1:60  % /!\ ta matrice A doit contenir au moins 60 elements
     tirage = A(r(i));
     % corps de l'expérience
     % extraction resultats
end

[vinc-mai]

Code :

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

1
2
i=randperm(numel(a));
result=a(i)

result est un vecteur contenant tous tes tirages!
a(1) le premier tirage, a(2) le deuxième...

Ce n'est pas ce que tu attends?

[pedrosanchau]

Je vais tenter ta deuxième proposition, elle m'a l'air pas mal du tout.