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

Calcul matriciel avec la multiplication

C

Vue hybride

Discussions similaires

Partager

Partager