Algorithme de jarvis

Version imprimable

bonosir ,
ces lignes d'instruction sont l'implémentation de l'algorithme de jarvis mais ça ne donne rien , qui peut me dire où est l'erreur ; mer6
je vous laisse avec le code :

Code:

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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
clear all ;
axis([0 10 0 10]);
%hold on ;
 
[y,x]=getpts ;
n=length(x);   %la longueur = le nbr de pt dessiné
[ymin i] = min(y);  %le pt le plu min + l'indice de ymin
 
 
pery=y(1);
y(1)=ymin ;  %%%de permutation de y pour débuter avec le pt le plus bas
y(i)=pery;
 
 
 
perx=x(1) ;
x(1)=x(i) ;  %%%  de permutation de x
x(i)=perx;
 
firstx=x(1) ;
firsty=y(1);
 
for k=2:n
 
    t=zeros(1,length(x)-1);   %declaration d'un vecteur dé 0
 
 
    for j=1:length(x)-1
        %%%% calculeeeer la tangente
 
        tangente =y(j+1)-firsty/x(j+1)-firstx;
        t(j)=tangente;  %%%% mettre lé tg ds un tableau t
 
    end
 
    [tmin j] = min(t); % la tg min + l'indice de la tag min
 
    %line([firstx, x(j+1)],[firsty, y(j+1)],'LineWidth',1,'Color','b')
    plot(firstx,firsty,x(j+1),y(j+1)) ;  % j car j c l'indice du tgmin
 
    firstx=x(k);  %aprés finir le 1er pt
    firsty=y(k);  %aller vers lé suivants
 
end

22/04/2011, 01h02
Franck Dernoncourt

Citation:

Envoyé par marocco.mimi

mais ça ne donne rien

= ?
22/04/2011, 01h57
Riham Darine

Oui , justement je veux connaitre pourquoi il ne donne pas de résultat
22/04/2011, 01h59
Franck Dernoncourt

Sur la figure ou dans la console ?
22/04/2011, 02h18
Riham Darine

je sais pa moi ,
pourquoi il ne dessine rien ???

Essaie :

Code:

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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
clear all ;
figure
hold on ;
 
axis([0 10 0 10]);
[y,x]=getpts ;
n=length(x);   %la longueur = le nbr de pt dessiné
[ymin i] = min(y);  %le pt le plu min + l'indice de ymin
 
plot(x,y)
pery=y(1);
y(1)=ymin ;  %%%de permutation de y pour débuter avec le pt le plus bas
y(i)=pery;
 
 
 
perx=x(1) ;
x(1)=x(i) ;  %%%  de permutation de x
x(i)=perx;
 
firstx=x(1) ;
firsty=y(1);
 
for k=2:n
 
    t=zeros(1,length(x)-1);   %declaration d'un vecteur dé 0
 
 
    for j=1:length(x)-1
        %%%% calculeeeer la tangente
 
        tangente =y(j+1)-firsty/x(j+1)-firstx;
        t(j)=tangente;  %%%% mettre lé tg ds un tableau t
 
    end
 
    [tmin j] = min(t); % la tg min + l'indice de la tag min
 
    %line([firstx, x(j+1)],[firsty, y(j+1)],'LineWidth',1,'Color','b')
    plot(firstx,firsty,x(j+1),y(j+1)) ;  % j car j c l'indice du tgmin
 
    firstx=x(k);  %aprés finir le 1er pt
    firsty=y(k);  %aller vers lé suivants
 
end

22/04/2011, 02h34
Riham Darine

nn ,ça dessine les lignes avec l'ordre de clique , et ne respecte pas le parcourt des indices dans le code
22/04/2011, 02h35
Invité

Il se dessinait bien quelque chose, mais seulement quelques points à peine visibles

Citation:

plot(firstx,firsty,x(j+1),y(j+1)) ; % j car j c l'indice du tgmin

à remplacer par

Code:

plot(firstx,firsty,'r*',x(j+1),y(j+1),'cp') ; % j car j c l'indice du tgmin

Pour mieux voir ces points
22/04/2011, 02h55
Riham Darine

1-pourquoi les positions des points se chengent ?
2-il ne dessine pas un enveloppe convexe , il juste relie entre des points (n'existe pas une fonction "polygon qui fait ça ?")

Fait un diff avec ce code, ca te donnera peut-être des idées de debug :

Code:

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
27
28
29
30
31
32
33
clear all ;
figure
hold on ;
 
axis([0 10 0 10]);
[x,y]=getpts ;
n=length(x);   %la longueur = le nbr de pt dessiné
%[ymin i] = min(y);  %le pt le plu min + l'indice de ymin
 
 
 
for k=1:n
 
    firstx=x(k);  %aprés finir le 1er pt
    firsty=y(k);  %aller vers lé suivants
 
    t=zeros(1,length(x)-1);   %declaration d'un vecteur dé 0
 
    for j=1:length(x)-1
        %%%% calculeeeer la tangente
 
        tangente =y(j+1)-firsty/x(j+1)-firstx;
        t(j)=tangente;  %%%% mettre lé tg ds un tableau t
 
    end
 
    [tmin j] = min(t); % la tg min + l'indice de la tag min
 
    %line([firstx, x(j+1)],[firsty, y(j+1)],'LineWidth',1,'Color','b')
    plot(firstx,firsty,'r*',x(j+1),y(j+1),'cp') ;  % j car j c l'indice du tgmin
 
 
end

22/04/2011, 03h14
Riham Darine

vraiment je suis désolé , mais ça aussi ne marche pas :cry::cry::cry::cry:
22/04/2011, 10h17
Invité

Tu parles d'enveloppe convexe, mais quels sont les variables qui définissent cette enveloppe? Je n'en vois aucune.
Tu ne fait qu'afficher les points que tu obtients avec getpts.
Il faudrait peut-être que tu revoies l'algorithme.
23/04/2011, 01h13
Riham Darine

revoies l'algorithme ? comment ;
mon code normalement :
1-l'user dessine des points
2-le code doit :
_récupérer les coordonnées des points dessinés et les mettre dans [x, y].
_choisir le point le plus en bas pour débuter (c pour ça j'ai fais ymin et la permutation)
_calculer l'angle entre l'horizon et le premier point et les autres points(x2,x3...xn) par la tangente
_mettre les tangentes dans un tableau et choisir la tangente min .
_ relier le point x(i) qui nous a donner la tangente min avec le x(1)
et reboucler x(i)=x(1) et ainsi de suite
24/04/2011, 21h16
Invité

Je ne suis pas tout à fait d'accord avec ces lignes

Citation:

_calculer l'angle entre l'horizon et le premier point et les autres points(x2,x3...xn) par la tangente
_mettre les tangentes dans un tableau et choisir la tangente min .

Je te laisse regarder ceci: [ame="http://www.youtube.com/watch?v=rsbxHqhyKD0"]Jarvis Algorithm [/ame] afin de bien comprendre le principe.
Bon il a commencé par un xmin, mais c'est pareil.
On fait un parcours 'double' de tous les points différents du précédent, et dès que la ligne rouge passe à l'extérieur de la verte, le point à l'extrémité verte devient celui à l'extrémité rouge.
Lorsqu'on a fini ce parcours, on ajoute au contour, l'extrémité verte.
Et ainsi de suite jusqu'à ce que le contour revienne à son point de départ.