Discussion :

[zicos]

Bonjour,
Je voudrais comparer toutes les valeurs d'une matrices par toutes les valeurs d'une autre matrice.Abituer du C, j'ai realiser 4 boucles for avec des matrices de 30 par 5 mais cela me fait 30x30x5x5=22500 execution. Je me suis rendu compte que cette solution etait trop lente a l'execution. Il a t'il des solutions pour palier a ce probleme, comme par exemple:

Code :

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

Diff=R(1:30,1:5)-C(1:30,1:5)

Bon, cela ne marche pas ainsi que toutes les solutions auquels j'ai penser. J'ai pas l'impression que ce soit realisable sans boucle mais bon pas sur non plus.
Si qq'un a une idee a ce probleme, je suis preneur.
MERCI d'avance pour vos reponses.

[Pingu1]

si c'est juste la différence qui t'intéresse, alors tu fais directement

Code :

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

Diff = R-C

[Jerome Briot]

Essaie ceci

Code :

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

all(all(R(1:30,1:5)==C(1:30,1:5)))

[zicos]

Ok je me suis sans doute un peu mal exprimer. Je ne souhaite pas comparer deux matrices mais bien calculer la difference entre chaque valeur de toutes les deux tables pour ensuite en faire des statistiques et des tracer de courbe.

Voici ce que je veux faire en MATLAB mais avec une philosophie de C, vous alez mieux comprendre.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
e=1;
for c=1:Nb_C
   for r=1:Nb_R
      for i=1:Nb_ngh
         for j=1:Nb_ngh
	    DiffVal(e)=TOC_C(c,j)-TOC_R(r,i);
            e=e+1;
         end
      end
   end
end

[rostomus]

Salut,

Essaie:

Code :

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

1
2
 
Dif=kron(A(:),ones(size(B(:))))-kron(ones(size(A(:))),B(:))

c'est ça?

[zicos]

Non ca pas l'air d'etre ca.
Avec mes boucles for, j'obtiens DiffVal qui fait 1x24025 pour deux matrices de 31x5.
Avec ta commande,j'obtiens la valeur Diff qui fait 384400x1.

[rostomus]

tu es sûr?

refais le test.
un exemple:

Code :

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

1
2
3
4
5
6
7
8
9
 
>> A=randn(31,5);
>> B=randn(31,5);
>> Dif=kron(A(:),ones(size(B(:))))-kron(ones(size(A(:))),B(:));
>> size(Dif)
 
ans =
 
       24025           1

[rostomus]

Un petit exemple:

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
A =
 
     1     2
     3     4
 
 
B =
 
     5     3
 
>> Dif'
 
ans =
 
    -4    -2    -2     0    -3    -1    -1     1

[zicos]

ca ma l'air d'etre bon, autant pour moi je m'ettais tromper de matrice.
Plus qu'a verifier l'ordonnancement et le tour est jouer.
Je te remercie et te souhaite un bon Week End.
MERCI