Discussion :

[sonialand]

j'espere que c'est le bon endroit pour poser ma question

je dois ecrire un code pour calculer toutes les distances entre les vecteurs d'une matrice et apres trouver la distance minimale.

j'ai ecris un code mais le probleme j'arrive pas a stocker toutes les valeurs ( cad toutes les distances) et du coup la valeur minimale est toujours la derniere valeur car je perds tout les resultats!

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
%cette partie est pour la creation de la matrice
Sym = [-1 -1/3 1/3 1];
X = zeros(256,4);
n = 1;
for i=1:4
    for j=1:4
        for l=1:4
            for m=1:4
                X(n,: ) = [Sym(i) Sym(j) Sym(l) Sym(m)];
                n = n+1;
            end
        end
    end
end
%calcul de la distance
for f=1:255
    for k=(f+1):256
        M_dist(k,1) = sqrt( (X(f,1)-X(k,1))^2  + (X(f,2)-X(k,2))^2  + (X(f,3)-X(k,3))^2 + (X(f,4)-X(k,4))^2 );
    end
        dmin = min(M_dist)
end

quelqu'un pourrait m'aider pour corriger mon code
merci
sonia

[Jean Dumoncel]

Bonjour,

Tu peux faire par exemple ceci :

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
%cette partie est pour la creation de la matrice
Sym = [-1 -1/3 1/3 1];
X = zeros(256,4);
n = 1;
for i=1:4
    for j=1:4
        for l=1:4
            for m=1:4
                X(n,: ) = [Sym(i) Sym(j) Sym(l) Sym(m)];
                n = n+1;
            end
        end
    end
end
%calcul de la distance
dmin = Inf;
for f=1:255
    M_dist = [];
    for k=(f+1):256
        M_dist(k-f,1) = sqrt( (X(f,1)-X(k,1))^2  + (X(f,2)-X(k,2))^2  + (X(f,3)-X(k,3))^2 + (X(f,4)-X(k,4))^2 );
    end
        dmin = min([M_dist ; dmin]);
end

[sonialand]

Merci bcp pour la reponse
encore une chose, le dmin doit etre different de 0!
comment est ce que je peux coder ca?
merci d'avance

[Jean Dumoncel]

Si M_dist prends précisément la valeur 0, pour les enlever, tu peux faire :

Code :

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

dmin = min([M_dist(M_dist~=0) ; dmin]);

Après attention, les nombres dans matlab ont une précision donnée (Pourquoi 0.3-0.2-0.1 est-il différent de 0 ?) et il sera peut-être nécessaire de faire un test en utilisant une tolérance (comme c'est expliqué dans le lien de la faq).