Discussion :

[Globoxx]

Bonjour,

Je fais un script qui permet d'empêcher que 2 points aléatoires ne se trouvent aux mêmes coordonnées.

Ce code devrait donc placer 9 points bien alignés qui remplissent une surface de 3 fois 3 avec un plot.

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
22
 
N = 9;
 
x = round(1 + (3-1) .* rand(N,1));
y = round(1 + (3-1) .* rand(N,1));
 
for j = 1:N
    for i = 1:N
 
    if i == j
        continue
 
    else
        while  x(j) == x(i) && y(j) == y(i)
                x(j) = round(1 + (3-1) * rand(1));
                y(j) = round(1 + (3-1) * rand(1));
        end
 
    end
 
    end
end

Mais comme vous vous en doutez, ça ne marche pas.

Si vous pouviez m'aider, merci !

[Dombrai]

je crois deja que tu t'es peut être emmelé un peu dans tes indices i et j.

Ensuite, le problème est que au fur et à mesure que tu avances tu ne prends plus en compte les lignes précédemment validées, par exemple : si tu trouves que la ligne 4 est identique à la ligne 3, tu retires la ligne 4 jusqu'a ce qu'elle soit différente de la ligne 3 mais tu ne vas pas tester si cette nouvelle ligne 4 est différente de la ligne 2.

une solution :

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
22
23
24
25
26
 
clear all
close all
clc
tic
N = 9;
 
x = round(1 + (3-1) .* rand(N,1));
y = round(1 + (3-1) .* rand(N,1));
 
XY = [x,y];
ligne = 2;
while ligne<=N
    ligne_verifiee = 0;
    for i = 1 :ligne-1
        if any(XY(ligne,:) ~= XY(i,:))
            ligne_verifiee = ligne_verifiee + 1;
        end
    end
    if ligne_verifiee == ligne-1
        ligne = ligne+1;
    else
        XY(ligne,:) = round(1 + (3-1) * rand(1,2));
    end
end
toc

il y a probablement plus élégant, mais ca marche

[Globoxx]

Merci beaucoup

[Dombrai]

ou sinon, beaucoup plus simplement, avec la fonction randperm :

Code :

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

1
2
3
4
5
6
 
x = [1;1;1;2;2;2;3;3;3];
y = [1;2;3;1;2;3;1;2;3];
ind = randperm(9);
x = x(ind);
y = y(ind);

[Jerome Briot]

ou encore :

Code :

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

1
2
3
ind = randperm(9);
 
[x,y] = ind2sub([3 3],ind);

[Globoxx]

Merci, mais le problème avec les deux codes que vous me proposez, c'est que si l'écart n'est pas de 3, mais prenons de 10, cela donnera quelque chose comme ça en résultat :

x = [1,4,7,9,6,8,3,5,2]

y = [1,1,1,1,1,1,1,1,1]

Alors premièrement, cela empêche d'avoir deux fois la même position en x ce qui n'est absolument pas un problème et que je ne veux donc pas empêcher.
Et deuxièmement, la coordonnée de y est toujours de 1.
Peut-on y remédier facilement ?

[Jerome Briot]

La solution que j'ai donnée fonctionne parfaitement :

Code :

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

1
2
3
4
5
N = 10;
 
ind = randperm(N*N);
 
[x,y] = ind2sub([N N],ind);

Ou alors je n'ai pas bien compris ton problème...

[Globoxx]

En fait mon but n'est pas de construire un joli bloc de N*N, j'ai juste donné cet exemple afin que l'on comprenne bien ce que mon script devrait faire si N valait 9 et le maximum du Random valait 3.
Ce n'est surement pas clair.
Ce que je fais précisément, c'est une simulation de colonisation spatiale, N étant le nombre de système à coloniser. Mon script doit simplement servir à empêcher que deux systèmes ne se trouvent aux mêmes coordonnées.
D'ailleurs celui de Dombrai fonctionne parfaitement, même si compliqué.

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
N = 9;
ecart = 10;
 
ind = randperm(ecart*ecart);
 
[x,y] = ind2sub([ecart ecart],ind(1:N));

[Globoxx]

C'est nickel merci !