Discussion :

[clara1815]

Bonjour a tous,
J'ai fait un programme matlab qui permet de faire le graphe de gabriel sur des figures rentrées par un nuages de points ou par des coordonnées. J'aurais voulu savoir si ce programme vous parraissait juste et si c'est le cas si on ne peut pas l'ameliorer.
Merci par avance pour votre aide.
Je mets juste la fonction, si vous avez besoin du programme en entier demandez le, je vous le donnerais :

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
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
 
function GG=GrapheGabriel(x,y,Triang)
 
NbTr=size(Triang,1);
MC=[];
for l=1:NbTr
    Tr=Triang(l,: );
    i=Tr(1,1);j=Tr(1,2);k=Tr(1,3);
    MC=[MC;[i,j];[j,k];[k,i] ];   
end;
SML=sort(MC,2); % 
SML=sortrows(SML);
N=MC(1,: );
n=size(MC,1);
for i=2:n 
    if MC(i-1,: )==MC(i,: ) ;
    else N=[N;MC(i,: )]; end;
end;  
 
GG=[];
n=size(MC,1);
for k=1:n
    Cote=MC(k,: );i=Cote(1);j=Cote(2);
    n=numel(x);
    A=[x(i) y(i)]; B=[x(j) y(j)]; 
    M=1/2*(A+B);
    R=(1/2)*norm(B-A);
    Go=1;k=0;Bool=0;
    while Go
        k=k+1;
        if ismember(k,[i j]); 
        else
          C=[x(k) y(k)];d=norm(C-M);Bool=d<R; 
        end;    
        Go=not(Bool) & (k<n);
    end;    
    if not(Bool) GG=[GG;Cote];end;
end;  
clf;
plot(x,y,'.r');
 
n=size(GG,1);
for k=1:n
    i=GG(k,1);j=GG(k,2);
    plot([x(i) x(j)],[y(i) y(j)],'-b');hold on;
end;

[tug83]

Pour avoir une idée comment MATLAB te suggère de l'améliorer , tu peux utiliser le mlint (Check M-files for possible problems)

Code :

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

mlint('GrapheGabriel.m')

Ce qui te donnera déjà des indications comme:

L 8 (C 5-6): 'MC' might be growing inside a loop. Consider preallocating for speed.
L 11 (C 1-3): The value assigned here to variable 'SML' might be unused.
L 16 (C 10): 'N' might be growing inside a loop. Consider preallocating for speed.
L 27 (C 10): Loop index 'k' is changed inside a FOR loop.
L 29 (C 9): Loop index 'k' is changed inside a FOR loop.
L 36 (C 18-19): 'GG' might be growing inside a loop. Consider preallocating for speed.
L 36 (C 18): For better readability, use newline, semicolon, or comma before this statement.

Tapes

Code :

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

doc mlint

pour avoir plus de détails sur cette fonctionnalité.

[clara1815]

Merci beaucoup pour votre aide, j'ai fait quelques modifications, mais je me demandais pourquoi matlab ne m'a pas donné les mêmes indications que vous?
Encore merci je ne connaissais pas cette fonction.

[tug83]

je me demandais pourquoi matlab ne m'a pas donné les mêmes indications que vous?

Les messages indiqués par le M-lint peuvent varier d'une version MATLAB à une autre. Moi j'ai utilisé la version actuelle de MATLAB la R2008b. De plus il faut savoir qu'on peut controler les parametres du Mlint dans les Preferences.
En allant dans File-> Preferences -> M-Lint (depuis toutes les versions égales ou supérieures à la R2006b) tu verras qu'il y a une cinquante de cases que tu peux cocher ou décocher à ta guise.