Discussion :

[matretro]

Bonjour,

J'ai un soucis sur lequel je bute depuis plusieurs jours sous MATLAB. J'ai vu qu'il y avait des cracks par ici alors je vous soumets mon problème. Si quelqu'un trouve une solution, je lui en serai très reconnaissant.

Voici mon problème : j'ai 2 matrices comme si dessus. Elles ont le même nombre de colonne mais pas le même nombre de ligne. Sur chaque colonne, il y a le même nombre de valeurs non nulles.
La première matrice

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
A=[
     1  2  9
    58  6  5
     0 65 98
     0  0 65
     6  6  0
     1  2  0
     0  5  8
     0  0  0
     2  3  0];
 
B=[
     0   6    0
     0   0    0 
     0  65.1  0
     1   0    0
     0   0    9.2
     0   0    0
     0   6    0
     0   2    0
    58   0    0 
     0   0    0
     0   0    5
     0   5    0 
     6   0   98
     0   0   65
     1   3    0 
     2   0    8 ];

J'aimerai comparer ces valeurs colonne par colonne mais sans faire de boucle.

L'objectif, c'est d'avoir des matrices ou il n'y aurait plus de zéros.
Pour faire la différence des deux matrices et connaitre les colonnes ou la différence est nulle.

exemple de ce que je veux :

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
A1=[
     1   2  9
    58   6  5
     6  65  98
     1   6  65
     2   2  8
     0   5  0
     0   3  0];
 
B1=[
     1   2     9.2
    58   6     5
     6  65.1  98
     1   6    65
     2   2    8
     0   5    0
     0   3    0];

Code :

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

sum(abs(A1-B1),1)

Il va sans dire que mes matrices sont en réalités bien plus grandes (200x8000) environ.

J'ai essayer des codes comme :

Code :

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

1
2
3
4
 
A1=A(A==0);
B1=B(B==0);
sum(abs(A1-B1),1);

mais je perds la forme de la matrice et je ne plus savoir quel colonnes correspond à une somme de différences de valeurs absolues nulle.

Merci d'avance pour votre aide.

[Jerome Briot]

Je ne suis pas sûr de bien tout comprendre.

Essaie néanmoins ceci :

Code :

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

1
2
3
4
5
6
7
A1(1:nnz(A(:,1)),1) = nonzeros(A(:,1));
A1(1:nnz(A(:,2)),2) = nonzeros(A(:,2));
A1(1:nnz(A(:,3)),3) = nonzeros(A(:,3));
 
B1(1:nnz(B(:,1)),1) = nonzeros(B(:,1));
B1(1:nnz(B(:,2)),2) = nonzeros(B(:,2));
B1(1:nnz(B(:,3)),3) = nonzeros(B(:,3));

ce qui donne

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
>> A1
 
A1 =
 
     1     2     9
    58     6     5
     6    65    98
     1     6    65
     2     2     8
     0     5     0
     0     3     0
 
>> B1
 
B1 =
 
    1.0000    6.0000    9.2000
   58.0000   65.1000    5.0000
    6.0000    6.0000   98.0000
    1.0000    2.0000   65.0000
    2.0000    5.0000    8.0000
         0    3.0000         0

[matretro]

Merci pour cette réponse mais si je ne me trompe pas, cela revient à faire une boucle sur les colonnes. Telle que tu présentes le code, tu as fait la boucle à la main.

Dans mon cas réel, j'ai des matrices 200x8000, ce qui me fait faire une boucle sur 8000 colonnes et que j'aimerai éviter.

Sans faire de boucle, je ne sais pas implémenter ta solution. As-tu une idée ?

[Jerome Briot]

L'ordre des valeurs dans les colonnes a-t-il de l'importance ?

[matretro]

L'ordre n'a pas vraiment d'importance du moment que c'est le même entre la matrice A et B afin que je puisse comparer les sommes de différences colonnes par colonne.

[Jerome Briot]

Je n'ai pas une machine très puissante, pourtant même avec une boucle for-end, la solution est très rapide :

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
A = [1  2  9
    58  6  5
    0 65 98
    0  0 65
    6  6  0
    1  2  0
    0  5  8
    0  0  0
    2  3  0];
 
A = repmat(A,20,3000);
 
for k = 1:5
    tic
    for n = 1:size(A,2)
        A1(1:nnz(A(:,n)),n) = nonzeros(A(:,n));
    end
    t(k) = toc;
    clear A1
end

Code :

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

1
2
3
4
5
>> t
 
t =
 
    0.0824    0.0798    0.0796    0.0763    0.0776

[matretro]

J'ai effectivement peut être chercher compliqué là ou il fallait faire simple. Je n'ai pas penser à regarder le temps d'exécution car je pensais à priori qu'en faisant des boucles ce serait trop long.
J'essaye d'implémenter tout ça et je vous tiens au courant.