Discussion :

[hidies]

Bonjour,
J aimerais comparer toutes les lignes au sein d une meme matrice, et si plusieurs lignes sont identiques, j aimerais recuperer le numero de la premiere ligne. Mon idee etait de prendre chaque ligne de la matrice et de voir si elle appartient a la matrice sans cette ligne. Voici les matrices A et B identiques :

Code :

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

1
2
3
4
5
6
7
8
A et B =
 
     1     1     1     1     1
     1     1     1     1     2
     1     3     1     1     1
     1     3     1     1     2
     1     3     1     1     1
     1     3     1     1     2

et mon code

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
d=length(B);
for i=1:d-1
    % On prend chaque ligne de B
    Ligne=B(i,:);
    % On reduit la dimension de la matrice
    A(1,:)=[];
    % On compare la ligne a la nouvelle matrice, si oui on prend l indice
    if (ismember(Ligne,A,'rows'))
        result(i)=i;
    end
end
% resultat vertical
resultat1 = result';
% resultat sans 0
result2 = resultat1(logical(resultat1));

Resultat :

Code :

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

1
2
3
4
t =
 
     3
     4

Cela fonctionne. Toutefois mon probleme est le temps d execution pour des matrices de grandes dimensions. Y a t il un moyen d ameliorer ce code ou meme supprimer la boucle ? merci bien

[salseropom]

Salut avec ton code, en faisant un tic toc j'obtiens un temps de calcul de 0.068825 seconds. Voici mon code :

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
 
A=[
1     1     1     1     1
     1     1     1     1     2
     1     3     1     1     1
     1     3     1     1     2
     1     3     1     1     1
     1     3     1     1     2];
 
 
tic
d=size(A,1);
dm1=d-1;
count=1;
 
for i=1:dm1
  test=1;
  j=i+1;
  x=A(i,:);
  while(test && j<=d)
    if(isequal(x,A(j,:)))
      test=0
    else
      j=j+1;
    end
  end
 
  if(~test)
    out(count)=i;
    count=count+1;
  end
end
toc
 
out

temps de calcul : 0.000191 seconds.

[Jerome Briot]

Voila une solution vectorisée :

Code :

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

1
2
3
4
5
6
7
8
9
M =[  1     1     1     1     1
     1     1     1     1     2
     1     3     1     1     1
     1     3     1     1     2
     1     3     1     1     1
     1     3     1     1     2]
 
[u,a,b]=unique(M,'rows');
idx=find(histc(b,unique(b))>1)

[salseropom]

Voila une solution vectorisée :

Code :

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

1
2
3
4
5
6
7
8
9
M =[  1     1     1     1     1
     1     1     1     1     2
     1     3     1     1     1
     1     3     1     1     2
     1     3     1     1     1
     1     3     1     1     2]
 
[u,a,b]=unique(M,'rows');
idx=find(histc(b,unique(b))>1)

Salut, je ne connaissais pas la fonction unique. Elle porte bien son nom

Mais en temps de calcul, elle est un peu longue (hidies veut une méthode très rapide). Avec ton algo j'ai un temps de calcul de 0.074878 seconds (toujours en faisant un tic toc).

[Jerome Briot]

en temps de calcul, elle est un peu longue

As-tu essayé avec des "matrices de grandes dimensions" comme le spécifie hidies dans son premier message ?
Les solutions vectorisées ne sont souvent efficaces qu'avec des problèmes de taille suffisante (ce qui peut devenir pénalisant en mémoire)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
M =[  1     1     1     1     1
     1     1     1     1     2
     1     3     1     1     1
     1     3     1     1     2
     1     3     1     1     1
     1     3     1     1     2];
 
n=300;
M=repmat(M,n,1);
 
<Test des solutions>

Regardes aussi la robustesse de ta solution...

[salseropom]

Salut, voici les tests pour une matrice 100*100
(1e ligne : algo de hidies, 2e ligne : mon algo, 3e ligne : ton algo)

Elapsed time is 0.100790 seconds.
Elapsed time is 0.014229 seconds.
Elapsed time is 0.074729 seconds.


pour une matrice 500*500

Elapsed time is 8.197833 seconds.
Elapsed time is 1.035057 seconds.
Elapsed time is 0.117466 seconds.

pour une matrice 1000*1000

Elapsed time is 73.572814 seconds.
Elapsed time is 7.952947 seconds.
Elapsed time is 0.268672 seconds.

j'ai arrêté le massacre à 1000*1000. Je m'incline devant ton aglo (qui plus est, ne prend que 2 lignes de code)

[Caro-Line]

Et puis la solution de Dutmatlab est la seule qui fonctionne quand on a une ligne qui apparait plus de 2 fois...Enfin si j'ai bien compris l'intitulé...

[salseropom]

Et puis la solution de Dutmatlab est la seule qui fonctionne quand on a une ligne qui apparait plus de 2 fois...Enfin si j'ai bien compris l'intitulé...

exact...
la prochaine fois je vérifirais par deux fois mon algo...

[hidies]

Merci a vous trois pour vos idees et commentaires . Etant donne les resultats je me suis attardee sur la methode de Dutmatlab. Toutefois cela ne fonctionne pas avec cette matrice :

Code :

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

1
2
3
4
5
6
7
8
9
M =
     1    13     3     0     1
     1    13     3     0     2
     1    13     3     0     1
     1    13     3     0     2
     1    13     5     0     1
     1    13     5     0     2
     1    13     5     0     1
     1    13     5     0     2

Normalement on devrait trouver comme resultat res et on obtient idx

Code :

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

1
2
3
4
5
res =       idx =
     1            1
     2            2
     5            3
     6            4

Cela vient du parametre de sortie b de la fonction unique qui donne

Code :

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

b = 1     2     1     2     3     4     3     4

Mais le calcul de ce parametre est interne a Matlab, donc je ne vois pas comment le modifier...

[Jerome Briot]

Code :

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

1
2
[u,a,b]=unique(M,'rows');
idx=setdiff(1:size(M,1),a)

Vu que la fonction UNIQUE renvoie les indices des lignes uniques... les autres indices ne sont donc pas uniques...

[hidies]

Merci c etait exactement ca : trouver les indices des lignes qui ne sont pas uniques ! Rapide et efficace