Discussion :

[dwight]

Bonjour,

Voici ce que je cherche à faire. Je dispose d'une série de couples, disons (5,7), (3,4) et (10,1), et je souhaite associer à chacun de ces points un identifiant, simplement représenté par un entier croissant :
1 => (5,7)
2 => (3,4)
3 => (10,1)
4 => ...

De plus, j'aimerais pouvoir accéder facilement à l'id à partir du couple, mais aussi au couple à partir de l'id. Comment faire ça le plus simplement et le plus proprement en matlab ?

J'ai pensé à former une matrice à partir des points :
m =
(5 7
3 4
10 1)

De cette façon, on peut récupérer le couple correspondant à l'id 2 en faisant m(2,:). Mais comment faire le contraire ? Comment savoir que (3, 4) correspond à la seconde ligne de la matrice ?

Merci d'avance !

[nahouto]

De cette façon, on peut récupérer le couple correspondant à l'id 2 en faisant m(2,. Mais comment faire le contraire ? Comment savoir que (3, 4) correspond à la seconde ligne de la matrice ?

essaie :

Code :

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

find(m == [3,4])

[dwight]

Merci nahouto, mais ça ne fonctionne pas.

Code :

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

1
2
??? Error using ==> eq
Matrix dimensions must agree.

Ce qui est somme toute assez logique.

Au cas où, j'ai essayé toutes les combinaisons possibles (virgule, pas virgule, parenthèses), rien à faire. Si tu vois une autre façon d'aborder le problème, n'hésite pas à m'en faire part.


Edit : je précise quand même que les couples sont nécessairement uniques.

[Jerome Briot]

Jette un oeil aux fonctions SORTROWS ou UNIQUE

[dwight]

Merci pour les pistes, mais je coince.

La fonction unique n'a pas l'air de vouloir faire ce que je veux. Je l'utilise forcément comme ceci :

Code :

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

1
2
 
[b, i, j] = unique(m, 'rows');

Etant donné que les lignes de m sont toutes uniques, b a le même contenu que m, et i et j sont des vecteurs-colonne contenant [1 2 3 4]'. Ce qui ne m'avance pas du tout.

Résultat à peu près similaire avec sortrows.

Où est-ce que j'ai faux ?

[Caro-Line]

Il y a peut-être plus joli mais essayes ça :

Code :

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

find(ismember(m,[3 4],'rows'))

[dwight]

Super Caro-Line, ça fonctionne ! Connaissais pas ismember...

Grand merci, ainsi qu'à Dut et à nahouto, et bonne soirée.

[Jerome Briot]

On pourra se passer de FIND ici au besoin pour accélérer le code si on doit accéder plusieurs fois aux lignes de m :

Code :

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

1
2
3
4
5
6
m = [5 7 ; 3 4 ; 10 1]
 
n = 1:size(m,1);
idx = ismember(m,[3 4],'rows');
 
n(idx)

Il y a aussi la fonction INTERSECT :

Code :

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

1
2
3
4
 
m = [5 7 ; 3 4 ; 10 1]
 
[x,u] = intersect(m,[3 4],'rows')

Mais je doute que cette dernière solution soit plus rapide

[dwight]

Ah ! Oui, ce fameux indexage logique, j'aurais certainement fini par y penser (mais quand ? ). Merci pour les précisions, Dut.

[dwight]

Hello. Je re-poste parce que j'ai de sérieux problèmes de performances avec mon code, utilisant les méthodes préconisées ci-dessus.

Voilà, j'ai une boucle qui contient 4 segments de ce genre :

Code :

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

1
2
3
4
5
 
        if m(i(1), i(2)) > 0
            r = ismember(t,[i(1) i(2)-1],'rows');
            a(k, n(r)) = 1;
        end

Où :
- m est une matrice creuse 30 000 x 30 000
- t est une matrice à 2 colonnes et à 10 000 lignes
- n est un vecteur de 10 000 éléments
- a est une matrice 10 000 x 10 000

La boucle fait 10 000 itérations. Problème : un temps d'exécution avoisinant les 7 minutes. Sachant que matlab est conçu pour traiter de grandes matrices, je me pose de sérieuses questions... C'est l'appel à ismember aussi souvent qui pose problème ?

En simulant l'exécution d'une itération, c'est-à-dire en tapant une à une les instructions de mon programme dans l'invite de commandes matlab, tout a l'air d'être instantané...

Merci bien

[nahouto]

tu peux montrer les fameux 4 segments ?

[dwight]

Ce sont presque les mêmes, seuls les indices de m changent. Tant que j'y suis, voilà la boucle complète. Le contexte, c'est la résolution de l'équation de Laplace.

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
 
for k = 1 : sz_t
        i = t(k, :);
 
        if m(i(1), i(2)) == 4
            a(k, k) = 1; 
            b(k) = 30;  
            continue;
        end     
 
        if m(i(1), i(2)) == 3 || m(i(1), i(2)) == 2
            b(k) = 15; 
        end
 
        if m(i(1)-1, i(2)) > 0
            r = ismember(t,[i(1)-1 i(2)],'rows'); 
            a(k, n(r)) = 1;
         end
 
        if m(i(1)+1, i(2)) > 0
            r = ismember(t,[i(1)+1 i(2)],'rows');            
            a(k, n(r)) = 1;
        end
 
        if m(i(1), i(2)-1) > 0
            r = ismember(t,[i(1) i(2)-1],'rows');
            a(k, n(r)) = 1;
        end
 
        if m(i(1), i(2)+1) > 0
            r = ismember(t,[i(1) i(2)+1],'rows');
            a(k, n(r)) = 1;
        end
 
        a(k, k) = -sum(a(k,:));
end

[Caro-Line]

ISMEMBER n'est effectivement pas une des fonctions les plus rapides qui soient, mais il peut y avoir d'autres paramètres comme par exemple ton IF.
Il y a surement moyen d'améliorer ton code mais il faudrait au moins toute la boucle.[EDIT]Ah bah voilà

As-tu essayé d'utiliser le PROFILER pour bien cerner ce qui prenait le plus de temps ?

[Jerome Briot]

Sachant que matlab est conçu pour traiter de grandes matrices, je me pose de sérieuses questions...

MATLAB est initialement fait pour traiter des problèmes de résolution de système linéaire et de problèmes de moindres carrés... pas pour trier des matrices

[dwight]

Bingo. D'après le profiler, c'est clairement l'appel à ismember qui est problématique. Au passage, c'est assez étrange, parce qu'il m'annonce 78 000 appels à cette fonction, alors que je n'en fais au maximum que 40 000 (cf. ma boucle donnée ci-dessus).

Et j'ai aussi découvert que Dut parlait par énigmes :

Jette un oeil aux fonctions SORTROWS ou UNIQUE

... or le profiler m'apprend que ismember appelle en interne ces deux fonctions, et je présume que c'est ce qui consomme pas mal.

Bon, j'aurai appris pas mal de choses sur matlab avec un seul topic. Merci à vous, je m'oriente vers un autre algo.

[Jerome Briot]

C'est la le grand avantage de MATLAB... mais aussi son principale inconvénient

La plupart des fonctions sont faites pour fonctionner avec des entrées de différents type. Ce qui nécessite généralement de multiplier les tests en début de fonction pour sélectionner la bonne méthode à appliquer. D'où pas mal de condition parfois bien inutile lorsqu'on applique la fonction sur des cas spécifiques.