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.