Calcul matriciel avec la multiplication

**argon** · 02/05/2007, 11h14

Bonjour tout le monde

Je un probleme, j'arrive à résoudre sur le papier les calculs sur les matrices.

Mais je n'arrive pas à le faire en programme.

Voici le 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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
 
#include <stdio.h>
#include <stdlib.h>
 
#define a 4
#define b 3
#define c 3
#define d 4
 
int multipli_matrice(int m1[a][b], int m2[c][d])
{
 int i,j,k,l;      
int m3[0][0];
 
puts("Multiplication de M1 par M2\n");
 
 
    for(i=0;i<4;i++)
    {
 
    puts("\n");
         for(j=0;j<4;j++)
         {
 
         m3[i][j]=m1[i][j]+ m2[j][i];
         }
         }
 
 
 
puts("Resultat M3\n");
 
 
         for(k=0;k<4;k++)
         {
 
         puts("\n");                
              for(l=0;l<4;l++)
              {
               printf(" %d",m3[k][l]);    
 
              }            
 
              }
 return 0;
}
 
void entre_donnee_mat()
{
 
    int mu[a][b];
 
 
    int md[c][d];
 
 
    int n,i,j,l,m,v,w;
    int partSize   = (n == RAND_MAX) ?        1 : 1 + (RAND_MAX - n)/(n+1);
   int maxUsefull = partSize * n + (partSize-1);
   int draw;
   srand(time(NULL)); 
 
 
     for(i=0;i<a;i++)
     {
 
         for(j=0;j<b;j++)
         {
           draw = draw = (int)(100.0*(double)rand()/(double)RAND_MAX);  
           mu[i][j]= draw; 
 
         }               
     }
 
     for(l=0;l<c;l++)
     {
          for(m=0;m<d;m++)
          {
               draw = draw = (int)(100.0*(double)rand()/(double)RAND_MAX);
               md[l][m]= draw;
 
               }
 
               }
 
     puts("Voici la matrice M1\n");
 
     for(v=0;v<4;v++)
     {
 
     puts("\n");
 
          for(w=0;w<3;w++)
               {
               printf(" %d",mu[v][w]);
 
     }
     }
 
     puts("\n");
 
     puts("Voici la matrice M2\n");
 
     for(v=0;v<3;v++)
     {
 
     puts("\n");
 
          for(w=0;w<4;w++)
               {
               printf(" %d",md[v][w]);
 
     }
     }
 
     puts("\n");
 
     multipli_matrice( mu, md);
}
 
int main(int argc, char *argv[])
{
 
 
    puts("calcul matriciel\n");
 
    entre_donnee_mat();
 
 
 
 
  system("PAUSE");	
  return 0;
}

Le calcul dois ce faire dans la fonction multipli_matrice.

Je n'est pas d'erreur de compile mais, le résultat est faux.

Quel est mon erreur svp?

Merci.

Cordialement

**souviron34** · 02/05/2007, 11h41

Envoyé par argon

....
Je n'est pas d'erreur de compile mais, le résultat est faux.
...

euhhhh........

c -Wall -pedantic ./test5.c -lm -o test5
./test5.c: In function `multipli_matrice':
./test5.c:13: warning: ISO C forbids zero-size array `m3'
./test5.c:13: warning: ISO C forbids zero-size array `m3'
./test5.c: In function `entre_donnee_mat':
./test5.c:61: warning: implicit declaration of function `time'
./test5.c:69: warning: operation on `draw' may be undefined
./test5.c:79: warning: operation on `draw' may be undefined
./test5.c:59: warning: unused variable `maxUsefull'

**femtosa** · 02/05/2007, 12h04

Juste une petite remarque sur le code : Il y a-t-il une raison au fait que tu utilise des variables de contrôles de boucle chaque fois différentes ? Tu utilise par exemple 'i', 'j' et 'k', 'l' dans ta fonction . Pourquoi n'utilise-tu pas juste 'i' et 'j'.

Même remarque pour l'autre fonction !

**Darwiz** · 02/05/2007, 12h42

Rappel (corrige moi si je me trompe) sur le produit de matrice :
(
x1, y1,
x2, y2
)
et
(
x1', y1',
x2', y2'
)

la matrice de multiplication =
(
x1*x1' + y1*x2', x1*y1' + y1*y2',
x2*x1' + y2*x2', x2*y1' + y2*y2'
)

En terme de tableau C ça devrait donner
m1[0][0] m1[1][0]
m1[0][1] m1[1][1]

m2[0][0] m2[1][0]
m2[0][1] m2[1][1]

Calcul :
m1[0][0]*m2[0][0] + m1[1][0]*m2[0][1], m1[0][0]*m2[1][0] + m1[1][0]*m2[1][1],
m1[0][1]*m2[0][0] + m1[1][1]*m2[0][1], m1[0][1]*m2[1][0] + m1[1][1]*m2[1][1]

Ton code ne fait pas du tout ce calcul, essaie un truc comme ça (j'ai pas testé) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
int i,j, k;
for (i = 0; i < 4; i++)
{
	for (j=0; j < 4; j++)
	{
		m3[i][j] = 0;
		for (k=0; k < 4: k++)
		{
			m3[i][j] = m3[i][j] + m1[k][j]*m2[i][k];
		}
	}
}

**argon** · 02/05/2007, 20h17

Bonsoir tout le monde

souviron34 je n'est pas c'est erreur lors de la compil, je signal que je travail avec devcpp, car la fac l'utilise.

Je n 'ais pas fais de beaucoup de modification.

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
37
38
39
40
 
 
int multipli_matrice(int m1[a][b], int m2[c][d])
{
 int i,j,k;      
int m3[0][0];
puts("Multiplication de M1 par M2\n");
 
 
    for(i=0;i<4;i++)
    {
 
    puts("\n");
         for(j=0;j<4;j++)
         {
 
         for(k=0;k<4;k++)
         {
         m3[i][j]+=m1[i][k]* m2[k][j];
         }
         }
         }
 
 
puts("Resultat M3\n");
 
 
         for(i=0;i<4;i++)
         {
 
         puts("\n");                
              for(j=0;j<4;j++)
              {
               printf(" %d",m3[i][j]);    
 
              }            
 
              }
 return 0;
}

Voici le résultat.

Resultat M3

-1731784084 151198941 -1340018197 -1608011244

-1731784084 151198941 -1340018197 -1608011244

-1731784084 151198941 -1340018197 -1608011244

-1731784084 151198941 -1340018197 -1608011244

Je ne comprends pas pourquoi il fait ca?

A +

**Emmanuel Delahaye** · 02/05/2007, 21h06

Envoyé par argon

souviron34 je n'est pas c'est erreur lors de la compil, je signal que je travail avec devcpp, car la fac l'utilise.

Problème de réglage :

http://emmanuel-delahaye.developpez....tm#cfg_compilo

Sinon :

Code :

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

int m3[0][0];

Tu n'as pas l'impression de faire n'importe quoi là ?

Si la taille est 'en dur' (connue à la compilation), définit un tableau statique de la bonne taille. Si elle est connue à l'exécution, crée une matrice dynamique :

http://emmanuel-delahaye.developpez.....htm#tabdyn_2d

**argon** · 02/05/2007, 21h12

Re

le probleme vient de moi Emmanuel Delahaye?

A+

**Emmanuel Delahaye** · 02/05/2007, 21h16

Envoyé par argon

le probleme vient de moi Emmanuel Delahaye?

Quel problème ? Apprends à

, ça ira plus vite...

**argon** · 02/05/2007, 21h24

Re

Est ce que l'erreur de souviron34 vient de la configuration de son compilateur.

J'ai moddifié:

Code :

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

1
2
 
int m3[4][4];

Voici le résultat m'aintenant.

Voici la matrice M1

99 14 34

88 98 70

52 34 83

74 74 29

Voici la matrice M2

31 87 97 94

94 63 69 99

57 95 13 10

Multiplication de M1 par M2

Resultat M3

2013342529 13958 4157531 18778

2013291503 4159456 2311376 2013350678

16894 2013426765 10986 2013520796

-515502658 -1876034641 -261137819 1092628012

Ce résultat est faux.

A+

**Emmanuel Delahaye** · 02/05/2007, 21h37

Envoyé par argon

Ce résultat est faux.

Ton application, je ne la connais pas ('multiplication de matrice', pour moi, ça n'a aucun sens, mais étant quasiment suis nul en math, c'est normal). Donc, j'ai soulevé un problème de C. Le reste, je ne n'ai pas regardé. Il y a rarement une seule erreur dans un code...

Une vérification du domaine des indices permet en général de lever les derniers problèmes de C...

assert() est ton ami.

http://emmanuel-delahaye.developpez....tes.htm#assert

**souviron34** · 03/05/2007, 01h38

Envoyé par argon

Est ce que l'erreur de souviron34 vient de la configuration de son compilateur.

disons plutôt que l'absence d'erreur de ton côté vient de la mauvaise utilisation des flags du compilateur de TON côté

Dans devcpp il suffit que tu ailles dans le panneau Options Makefile ou Flags de compilations et tu rajoutes (au minimum) -Wall -pedantic

Et comme c'est noté dans les diagnostics et comme Emmanuel l'a re-fait remarquer, tu définis un tableau comme étant de taille (0,0), ce qui ne PEUT PAS passer...

Ensuite tu multiplies des "int", ça risque fort de dépasser la capacité d"un "int", sauf si tu sais que les valeurs sont faibles...

Et enfin, n'as-tu pas oublié d'initialiser m3 à zéro, maintenant que tu lui a mis des dimensions non nulles ?

**Darwiz** · 03/05/2007, 10h34

Bon j'ai testé le bout de programme que j'ai donné et ça m'a donné des résultats cohérents par contre étant donné qu'on travail sur des matrices 4x4 tu dois mettre des 0 dans les cases qui ne te servent pas sur les matrices m1 et m2. Si tu ne mets pas 0 dans les cases non utilisées tu vas obtenir des résultats incohérents car tu vas accéder à de la mémoire non initialisée.

Voici le résultat que j'ai obtenu :

Matrice M1 :
99 14 34 0
88 98 70 0
52 34 83 0
74 74 29 0

Matrice M2 :
31 87 97 94
94 63 69 99
57 95 13 10
0 0 0 0

Matrice M3 :
6323 12725 11011 11032
15930 20480 16208 18674
9539 14551 8469 9084
10903 13855 12661 14572

Ton code avec les remarques de chacuns :

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
int multipli_matrice(int m1[a][b], int m2[c][d])
{
	int i,j,k;      
	int m3[4][4]; // Initialiser la taille du tableau pour le résultat
	puts("Multiplication de M1 par M2\n");
 
 
	for (i = 0; i < 4; i++)
	{  
		for (j=0; j < 4; j++)
		{  
			m3[i][j] = 0; // Initialisation à 0
			for (k=0; k < 4: k++)  
			{  
				m3[i][j] += + m1[k][j]*m2[i][k];
			}  
		} 
	}
 
	puts("Resultat M3\n");
 
	for(i=0;i<4;i++)
	{
 
		puts("\n");                
		for(j=0;j<4;j++)
		{
			printf(" %d",m3[i][j]);    
		}            
 
	}
	return 0;
}