Discussion :

[dreamdam]

bonjour tt le monde voila dit fois que j'essaye de programmer une procedure pour multiplier deux matrice si qqun arrive a faire ca ou l'a deja pouvez vous m'aidez svp ??sachant qu'une matrice de taille 2lignes mais nombres de colonne indefini... merci pr votre aide ...

[Eric Sigoillot]

Il suffit d'utiliser la forumule de multiplication...

Si (a) et (b) sont deux matrices de tailles m,n et n,p alors si on pose (c) = (a)*(b), on a :

(c)i,j = Somme(k=1, n, (a)i,k * (b)k,j)

Bon, transposons tout ça en Pascal à présent. On prend A: array[1..m, 1..n] of Real; (m lignes, n colonnes) et B: array[1..n, 1..p] of Real; (n lignes, p colonnes) et on pose C: array[1..m, 1..p] of Real; (m lignes, p colonnes), ça nous donne donc :

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
 
var
  A: array[1..m, 1..n] of Real;
  B: array[1..n, 1..p] of Real;
  C: array[1..m, 1..p] of Real;
  i, j, k: Integer;
  R: Real;
 
begin
  for i := 1 to m do
    for j := 1 to p do
    begin
      R := 0;
      for k := 1 to n do R := R + A[i, k] * B[k, j];
      C[i, j] := R;
    end;
end;

Avec un peu de méthode, on arrive à tout

A+

[charly]

surtous avec de bon niveau en math

[Axiome]

Normalement, ça ne devra poser aucun problème.
Tu déclare trois matrices (tableau à deux dimensions)
les deux premières seront les matrices données et la troisième la matrice résultat.

il suffit de parcourir les deux matrices par une boucle et de faire les multiplications et les sommes nécessaires et d'affecter ça dans l'emplacement correct de la troisième matrice.

[dreamdam]

merci bien c'est avec joie que je vais faire tourner ca mais j'ai pas mal de mal avec les matrices en faite lol j'ai poste une autres question sur comment inverse une matrice avec le pivot de gauss mais les gens on ete assez sec alors je vais m'arrache sur la tache si je trouve je posterais pour a la fin vous poste un programme de crypto

[Bibipulse]

Bonjour

J'ai fait une procédure qui multiplie 2 matrices de taille identique dans une 3e, mais au moment de l'éxecution il y a un RANGE CHECK ERROR. J'ai beau chercher je ne trouve pas d'où ça vient.
Le type matrice a été crée auparavant de la manière suivante:

Code :

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

1
2
3
4
5
 
type matrice = record
        L : integer;
        C : integer;
       data : array [1..Lmax, 1..Cmax] of Zm {Zm est un type crée également, qui au final est un entier}

Ma procédure est alors :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
procdedure matmult (var a,b: matrice; var c : matrice); 
var i,j,k,r : integer;
 
begin
     for i:= 1 to Cmax do
          begin
               for j:=1 to Lmax do
                   begin
                      r:=0;
                      for k:=1 to Cmax do r:= r + (a.data[i,k]) * (b.data[k,j]);
                      c.data[i,j]:= r;
                   end;
           end;
end;

Le Range Check Error intervient apparamment juste après le premier passage dans la boucle k.

Si vous voyez où peut se trouver l'erreur, ca serait sympa.

[Giovanny Temgoua]

Salut,

Il y'a une technique que j'utilise souvent pour le débogage, elle consiste à mettre des writeln() partout, histoire de savoir exactement quelle sont les valeurs de mes variables (dans ton cas, il s'agira notemment des indices de boucles).

Bref, si LMax et CMax sont différents alors tu obtiens l'erreur à cette ligne

Code :

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

1
2
 
r:= r + (a.data[i,k]) * (b.data[k,j]);

car comme k et i sont compris entre 1 et CMax alors tu essaie d'atteindre l'élément [CMax, CMax] du tableau, opération impossible si CMax > LMax

Ainsi, successivement tu accedes aux éléments [CMax,CMax], [CMax,LMax] opération impossible si CMax > LMax (impossible d'atteindre l'élément [CMax,CMax]) et impossible aussi si LMax > CMax (impossible d'atteindre [CMax,LMax]) puisque le tableau est déclaré comme array [1..Lmax, 1..Cmax]

==> RANGE CHECK ERROR

[Eric Sigoillot]

Y'a quand même le programme complet juste au-dessus avec des matrices de taille quelconque. Alors après, c'est juste une question du type "je sais adapter du code" ou pas.

A+

[Binga]

Il y aune chose que je veux ajouter c'est qu'il ya une condition importante avant de faire la multiplication de deux matrices:
Le nombre de colone de la 1ére matrice doit etre égale à celui de lignes de la 2ème matrice.

[forbin]

Salut à tous,

J'ai lu qu'un certain Volker Strassen(*) avait mis au point un algorithme de calcul du produit de matrices plus rapide que l'algorithme algébrique classique.

Quelqu'un l'aurait-il testé ?

A+

4BIN

______________

(*) "Gaussian Elimination is not Optimal", Numerische Mathematik, N°14(3), pp.354-56, 1969.

[Binga]

J'ai lu qu'un certain Volker Strassen(*) avait mis au point un algorithme de calcul du produit de matrices plus rapide que l'algorithme algébrique classique.

Le produit de deux matrices est une définition donc on ne peut pas la changer (voir un peu le cours d'algebre). A mon avis Je pense pas qu'il existe une autre methode de calcule du produit différente de la methode classique (la définition).
a+

[Eric Sigoillot]

Il est possible qu'il existe d'autres algorithmes permettant d'améliorer les performances de calcul.
Quoiqu'il en soit c'est une discussion du forum Algorithmes si tel est le cas.

Et merci de ne pas déterrer les sujets résolus depuis plus d'un mois

A+