Algo

Inverse d'une matrice

Recherche :

Mot : Pseudo : Filtrer
Page : 1 2 Page Suivante Page Précédente Bas de page
Auteur	Sujet : Inverse d'une matrice

sanglier04

boujour ,

je dois faire un algo qui calcule l'inverse d'une matrice. pour les matrices carrées c simple mais pour des matrices non carrées je n'y arrive pas. le seul indice que l'on met donné c qu'il ya suivant la matrice soit une inverse a guauche, soit une inverse a droite :??: . jai beau chercher mais je ne trouve pas

Publicité

Ace17

Je vais peut etre dire une connerie ( arretez moi au cas! ) mais je vais quand meme essayer d'aider

L'inverse d'une matrice A c'est bien la matrice qui donne l'identité quand on la multiplie par A c'est ca? Dans ce cas retrouver un inverse de matrice ne se résume-t-il pas a résoudre un systeme tout simplement?

Ace17

http://www.gel.ulaval.ca/~fortier/ [...] rices.html

"L'inverse d'une matrice ne peut être obtenue qu'à partir d'une matrice carrée"

sanglier04

surement mais dans ce cas je ne sais po le resoudre quand la matrice n'est pas carré

Ace17

Ace17 a écrit a écrit :

http://www.gel.ulaval.ca/~fortier/ [...] rices.html

"L'inverse d'une matrice ne peut être obtenue qu'à partir d'une matrice carrée"

Désolé c'est une connerie

edit : finalement ca doit etre possible quand les deux matrices ont des dimensions échangées

Message édité par Ace17 le 23-10-2002 à 16:41:42

Toucouch

L'inverse d'une matrice ce n'est pas la transposée de la comatrice sur le déterminant de A?

Problème: je ne me souviens plus comment on calcule le déterminant d'une matrice de dimension supérieure à 3x3. :sweat: :honte: :cry:

sanglier04

daprès ce que jai trouvé en cherchant, lorsque ce n'est pas carré, il y a une infinité de solutions. Il est probable que ce n'est pas une matrice "complètement inversse" mais "pseudo inverse" ou "demi inverse".

l'inverse a droite d'une matrice A est du type X A = B et a guauche c du typs A X = B. mais après ca je n'en sais po plus

Pitounet

A^-1 = 1/det(A) * transposée(Comatrice(A)).

Il faut que A soit carrée.

Toucouch

Pitounet a écrit a écrit :

A^-1 = 1/det(A) * transposée(Comatrice(A)).

Il faut que A soit carrée.

Et voilà Joe La Science qui vient faire son quetquet avec ses formules de la mort qui tue...

Ca geht's?

darkoli

Le Petit Dinosaure Bleu

Toucouch a écrit a écrit :

Et voilà Joe La Science qui vient faire son quetquet avec ses formules de la mort qui tue...

Ca geht's?

Tout à fait thierry ...
Mais jamais deux sans trois alors me voilà !!!
Sinon j'ai oublié comment on faisait alors ... Mais je vais chercher ça dont mes fichier C ce soir (on sait jamais).

---------------
Le site de l'année :D (XHTML 1.0 strict) : http://darkoli.free.fr/index.html

Publicité

Toucouch

DarkOli a écrit a écrit :

Tout à fait thierry ...
Mais jamais deux sans trois alors me voilà !!!
Sinon j'ai oublié comment on faisait alors ... Mais je vais chercher ça dont mes fichier C ce soir (on sait jamais).

:lol: Excellent! Décidément, ça bosse dur chez les anciens de MT43... N'empêche, quelqu'un peut-il me réexpliquer vite fait comment on calcule un déterminant de matrice > 3x3?

tomlameche

Et pourquoi pas ?

sanglier04 a écrit a écrit :

boujour ,

je dois faire un algo qui calcule l'inverse d'une matrice. pour les matrices carrées c simple mais pour des matrices non carrées je n'y arrive pas. le seul indice que l'on met donné c qu'il ya suivant la matrice soit une inverse a guauche, soit une inverse a droite :??:

. jai beau chercher mais je ne trouve pas

Ben déjà, toute matrice, même carré, n'est pas inversible. Ca depend si les vecteurs formés par les colonnes sont indépendant ou pas.
Sinon, le meilleurs axe de recherche pour trouver de bons algo de calcul d'inverse de matrice, c'est de chercher des algos de recherche opérationnel. Tu devrai trouver ton bonheur la dedans.

bjone

Insert booze to continue

ce qui a été dit est correct.

par contre, petite astuce, pour un matrice 3x3 qui effectue une rotation, transposer la matrice donne la matrice inverse de la rotation (valable uniquement dans le cas d'une rotation)

farib

tomlameche a écrit a écrit :

Ben déjà, toute matrice, même carré, n'est pas inversible. Ca depend si les vecteurs formés par les colonnes sont indépendant ou pas.
Sinon, le meilleurs axe de recherche pour trouver de bons algo de calcul d'inverse de matrice, c'est de chercher des algos de recherche opérationnel. Tu devrai trouver ton bonheur la dedans.

fo écrire la matrice a inverser et une matrice identité cote a cote, puis ensuite il faut réduire la premiere matrice en une matrice identité, et reporter les opératiosn qu'on fait sur les lignes sur la matrice identité qu'on a écrit a dote de la matrice a inverser.

cela dit, c'est des maths pures ca.

sanglier04

farib a écrit a écrit :

fo écrire la matrice a inverser et une matrice identité cote a cote, puis ensuite il faut réduire la premiere matrice en une matrice identité, et reporter les opératiosn qu'on fait sur les lignes sur la matrice identité qu'on a écrit a dote de la matrice a inverser.

cela dit, c'est des maths pures ca.

J utilise ce principe pour les matrices carré mais pour les matrices qui ne le sont pas je suis pas sur que cela fonctionne.

darkoli

Le Petit Dinosaure Bleu

Toucouch a écrit a écrit :

:lol:

Excellent! Décidément, ça bosse dur chez les anciens de MT43... N'empêche, quelqu'un peut-il me réexpliquer vite fait comment on calcule un déterminant de matrice > 3x3?

Je l'ai noté sur mon poly mais je ne sais plus où il est ce con !!!

En attendant j'ai retrouvé ça (j'ai le même en assembleur ). Par contre ce n'est pas très clair mais c'est une version 'optimisée' et spécifique aux matrices 4x4.

Euh la notation de la matrice est (je crois) :

Code :

00 01 02 03
04 05 06 07
08 09 10 11
12 13 14 15

Code :

//---------------------------------------------------------------------------
// inversion de matrices : matrice de taille 4x4 exclusivement
void Inversion_de_matrice(/*double* matrice,double* resultat_*/) {
// variables temporaires
double t[6]={0,0,0,0,0,0};
// calcul des cofacteurs de la ligne 1
t[0]=matrice[10]*matrice[15] - matrice[14]*matrice[11];
t[1]=matrice[9]*matrice[15] - matrice[13]*matrice[11];
t[2]=matrice[9]*matrice[14] - matrice[13]*matrice[10];
t[3]=matrice[8]*matrice[15] - matrice[12]*matrice[11];
t[4]=matrice[8]*matrice[14] - matrice[12]*matrice[10];
t[5]=matrice[8]*matrice[13] - matrice[12]*matrice[9];
resultat[0]=matrice[5]*t[0] - matrice[6]*t[1] + matrice[7]*t[2];
resultat[1]=matrice[4]*t[0] - matrice[6]*t[3] + matrice[7]*t[4];
resultat[2]=matrice[4]*t[1] - matrice[5]*t[3] + matrice[7]*t[5];
resultat[3]=matrice[4]*t[2] - matrice[5]*t[4] + matrice[6]*t[5];
// calcul du determinant
td=matrice[0]*resultat[0] - matrice[1]*resultat[1] + matrice[2]*resultat[2] - matrice[3]*resultat[3];
// Form1->Edit1->Text=FloatToStr(td);
// test du determiant : si non nul on continue
if (td==0) {
resultat_[ 0]=1;
resultat_[ 1]=0;
resultat_[ 2]=0;
resultat_[ 3]=0;
resultat_[ 4]=0;
resultat_[ 5]=1;
resultat_[ 6]=0;
resultat_[ 7]=0;
resultat_[ 8]=0;
resultat_[ 9]=0;
resultat_[10]=1;
resultat_[11]=0;
resultat_[12]=0;
resultat_[13]=0;
resultat_[14]=0;
resultat_[15]=1;
return;
}
// calcul des cofacteurs de la ligne 2
resultat[4]=matrice[1]*t[0] - matrice[2]*t[1] + matrice[3]*t[2];
resultat[5]=matrice[0]*t[0] - matrice[2]*t[3] + matrice[3]*t[4];
resultat[6]=matrice[0]*t[1] - matrice[1]*t[3] + matrice[3]*t[5];
resultat[7]=matrice[0]*t[2] - matrice[1]*t[4] + matrice[2]*t[5];
// calcul des cofacteurs de la ligne 3
t[0]=matrice[2]*matrice[7] - matrice[6]*matrice[3];
t[1]=matrice[1]*matrice[7] - matrice[5]*matrice[3];
t[2]=matrice[1]*matrice[6] - matrice[5]*matrice[2];
t[3]=matrice[0]*matrice[7] - matrice[4]*matrice[3];
t[4]=matrice[0]*matrice[6] - matrice[4]*matrice[2];
t[5]=matrice[0]*matrice[5] - matrice[4]*matrice[1];
resultat[8] =matrice[13]*t[0] - matrice[14]*t[1] + matrice[15]*t[2];
resultat[9] =matrice[12]*t[0] - matrice[14]*t[3] + matrice[15]*t[4];
resultat[10]=matrice[12]*t[1] - matrice[13]*t[3] + matrice[15]*t[5];
resultat[11]=matrice[12]*t[2] - matrice[13]*t[4] + matrice[14]*t[5];
// calcul des cofacteurs de la ligne 4
resultat[12]=matrice[9]*t[0] - matrice[10]*t[1] + matrice[11]*t[2];
resultat[13]=matrice[8]*t[0] - matrice[10]*t[3] + matrice[11]*t[4];
resultat[14]=matrice[8]*t[1] - matrice[9]*t[3] + matrice[11]*t[5];
resultat[15]=matrice[8]*t[2] - matrice[9]*t[4] + matrice[10]*t[5];
// inversion
td=1/td;
resultat_[0]=resultat[0]*td;
resultat_[1]=-resultat[4]*td;
resultat_[2]=resultat[8]*td;
resultat_[3]=-resultat[12]*td;
resultat_[4]=-resultat[1]*td;
resultat_[5]=resultat[5]*td;
resultat_[6]=-resultat[9]*td;
resultat_[7]=resultat[13]*td;
resultat_[8]=resultat[2]*td;
resultat_[9]=-resultat[6]*td;
resultat_[10]=resultat[10]*td;
resultat_[11]=-resultat[14]*td;
resultat_[12]=-resultat[3]*td;
resultat_[13]=resultat[7]*td;
resultat_[14]=-resultat[11]*td;
resultat_[15]=resultat[15]*td;
// fin
}
//---------------------------------------------------------------------------

Message édité par darkoli le 23-10-2002 à 20:15:35

---------------
Le site de l'année :D (XHTML 1.0 strict) : http://darkoli.free.fr/index.html

farib

tu fais marcher ton imagination et tu cherches sur google !

sanglier04

farib a écrit a écrit :

tu fais marcher ton imagination et tu cherches sur google !

si je poste ici c que jai deja chercher ...

Carbon_14

Pour 3x3, y a http://www.geocities.com/SiliconVa [...] ach03.html

Je me souviens très bien de "A^-1 = 1/det(A) * transposée(Comatrice(A))" bien que très ancien, mais suis guère plus capable de l'implémenter :pt1cable: .

sanglier04

carbon_14 a écrit a écrit :

Pour 3x3, y a http://www.geocities.com/SiliconVa [...] ach03.html

Je me souviens très bien de "A^-1 = 1/det(A) * transposée(Comatrice(A))" bien que très ancien, mais suis guère plus capable de l'implémenter :pt1cable:

merci, mais c tjs pour des matrices carrés, ce qui me pose problème c surtout pour celles qui ne sont pas carrés.

Kristoph

Pas carré => pas inversible.

Ca a deja ete dit ici même.

B est l'inverse de A ssi AB = I et BA = I. Or c'est pas possible avec des matrices pas carrées. Point.

chicha

Kristoph a écrit a écrit :

Pas carré => pas inversible.

Ca a deja ete dit ici même.

B est l'inverse de A ssi AB = I et BA = I. Or c'est pas possible avec des matrices pas carrées. Point.

+1
evidemment qu'une matrice n'est inversible que si elle est carré ! condition NECESSAIRE mais pas suffisante bien sûr...

Ace17

Est-ce que le produit de deux matrices non carrées donne/peut donner une matrice carrée?

sanglier04

comme je l'ai dit plus haut, il semble que ce ne soit po une inverse complète mais une sorte de demi inverse (jaurai plus de précision vendredi) merci a tout ce qui ont essayer de m'aider.

Carbon_14

Ca permettra de savoir quelle est la transformation envisagée. Il doit y avoir un "abus de langage". La vraie inversion ne pouvant QUE concerner des matrices carrées.

darkoli

Le Petit Dinosaure Bleu

Ace17 a écrit a écrit :

Est-ce que le produit de deux matrices non carrées donne/peut donner une matrice carrée?

Pour multiplier A par B il faut que le nombre de colonnes de A soit identiques au nombres de lignes de B.

La matrice obtenue aura autant de lignes que la matrice A et autant de colonne que la matrice B.

Donc avec une matrice A de 7 colonnes et 4 lignes et une matrice B de 4 colonnes et 7 lignes tu obtiens bien une matrice carrée de 4 colonnes et 4 lignes.

Conclusion pour obtenir une matrice carrée avec une multiplication de deux matrices A et B (AxB). Les matrices doivent avoir les mêmes dimensions mais inversées ... (lignes A = colonnes B et colonnes A = lignes B).

Message édité par darkoli le 24-10-2002 à 08:52:43

---------------
Le site de l'année :D (XHTML 1.0 strict) : http://darkoli.free.fr/index.html

Ace17

Ok , dans ce cas il est possible d'avoir deux matrices non carrées, aux dimensions inversées l'une par rapport a l'autre, dont le produit donne la matrice identité; n'est-ce pas? Pourquoi ne pourrait-on pas dire qu'elles sont inverses l'une de l'autre?

Kristoph

Parceque A est inversible ssi il existe B telle que AB = I et BA = I et ca c'est impossible si A n'est pas carrée. Faut vraiment se repeter ici

Ace17

Ca y est j'ai compris
En fait on peut trouver B pour que AB = I mais a ce moment la il y n'y a pas forcément BA = I; d'ou les termes d'inverse a droite ou a gauche mentionnés plus haut

Pitounet

Toucouch a écrit a écrit :

Et voilà Joe La Science qui vient faire son quetquet avec ses formules de la mort qui tue...

Ca geht's?

Héhé, il faut bien resortir le peu qui me reste.
D'ailleurs, je ne sais vraimentpas pourquoi cette formule est restée alors que tout le reste est parti !

Je peux t'appeller ce soir au tel ?

Pitounet

pour le determinant de la matrice 3*3, je crois que c'et un truc comme ça :

a b c
A= d e f
i j k

det(A) =aek + bfi + cdj - cei -afj -dbk

Attention, c'est à prendre avec des pincettes comme formule, c'est suuuuuuuper vieux !

Message édité par Pitounet le 24-10-2002 à 11:15:32

darkoli

Le Petit Dinosaure Bleu

pitounet a écrit a écrit :

pour le determinant de la matrice 3*3, je crois que c'et un truc comme ça :

a b c
A= d e f
i j k

det(A) aek + bfi + cdj - cei -afj -dbk

Attention, c'est à prendre avec des pincettes comme formule, c'est suuuuuuuper vieux !

Sinon tu dois pouvoir retrouver la formule dans mon bout de code (qui fonctionne correctement !!!) mais comme j'avais tout décomposer les calculs à l'époque, il faut faire l'opération inverse

---------------
Le site de l'année :D (XHTML 1.0 strict) : http://darkoli.free.fr/index.html

farib

pour le déterminant, revenez a la définition avec les formes p-linéaires alternées

Toucouch

pitounet a écrit a écrit :

Héhé, il faut bien resortir le peu qui me reste.
D'ailleurs, je ne sais vraimentpas pourquoi cette formule est restée alors que tout le reste est parti !

Ben c'était pas ce jour-là que tu suivais le cours?

Citation :

Je peux t'appeller ce soir au tel ?

Caramba! Ce soir c'est notre 5° anniversaire de fiançailles! (Ca fait un peu rétro, mais bon, on est sentimental ou on l'est pas ). Donc ce soir, ch'tit resto sympa, etc. donc cramé...

Mais demain soir c'est avec grand plaisir que j'écouterai de nouveau le son mélodieux et suave de ta jolie voix... (J't'avais dit que j'étais sentimental... )

PS: Oli, t'as une belle voix aussi, hein! Si tu veux appeler aussi...

Message édité par Toucouch le 24-10-2002 à 15:18:37

bobricard

Ah bah ouais mais bon !

Soit la matrice

On peut calculer son déterminant par 2 manières:

1)
Sachant que , on en déduit:

det(A)=a*(e*k-j*f)-d*(b*k-j*c)+i*(b*f-e*c)
det(A)=a*e*k+d*j*c+i*b*f-a*j*f-d*b*k-i*e*c

2)On rajoute les deux premières lignes en dessous de la matrice:

On multiplie entre eux les termes des diagonales en rouge, on fait la somme, on obtient M=a*e*k+d*j*c+i*b*f
On fait la même opération pour les diagonales en bleu, on obtient N=i*e*c+a*j*f+d*b*k

On a alors det(A)=M-N=a*e*k+d*j*c+i*b*f-i*e*c-a*j*f-d*b*k

Pitounet

Toucouch a écrit a écrit :

Ben c'était pas ce jour-là que tu suivais le cours?

Citation :

Je peux t'appeller ce soir au tel ?

Bon anniversaire à tous les deux alors
J'essaye de t'appeler demain soir, sans pb.

darkoli

Le Petit Dinosaure Bleu

Toucouch a écrit a écrit :

Ben c'était pas ce jour-là que tu suivais le cours?

Citation :

Je peux t'appeller ce soir au tel ?

Caramba! Ce soir c'est notre 5° anniversaire de fiançailles! (Ca fait un peu rétro, mais bon, on est sentimental ou on l'est pas ). Donc ce soir, ch'tit resto sympa, etc. donc cramé...

Mais demain soir c'est avec grand plaisir que j'écouterai de nouveau le son mélodieux et suave de ta jolie voix... (J't'avais dit que j'étais sentimental... )

PS: Oli, t'as une belle voix aussi, hein! Si tu veux appeler aussi...

Sauf que je ne suis pas devin ... [:cupra] (Le numéro de téléphone ???)

---------------
Le site de l'année :D (XHTML 1.0 strict) : http://darkoli.free.fr/index.html

Toucouch

bobricard a écrit a écrit :

Soit la matrice http://www.ifrance.com/manupap/matrice.jpg

On peut calculer son déterminant par 2 manières:

1) http://www.ifrance.com/manupap/formule.jpg
Sachant que http://www.ifrance.com/manupap/det.jpg, on en déduit:

det(A)=a*(e*k-j*f)-d*(b*k-j*c)+i*(b*f-e*c)
det(A)=a*e*k+d*j*c+i*b*f-a*j*f-d*b*k-i*e*c

2)On rajoute les deux premières lignes en dessous de la matrice:
http://www.ifrance.com/manupap/astuce.jpg
On multiplie entre eux les termes des diagonales en rouge, on fait la somme, on obtient M=a*e*k+d*j*c+i*b*f
On fait la même opération pour les diagonales en bleu, on obtient N=i*e*c+a*j*f+d*b*k

On a alors det(A)=M-N=a*e*k+d*j*c+i*b*f-i*e*c-a*j*f-d*b*k

Oui, mais pour une matrice plus grande?

Toucouch

Pitounet > Ca roule! A demain, donc...

Darkoli > Demande à Pitounet...

farib

Toucouch a écrit a écrit :

Oui, mais pour une matrice plus grande?

ce sotn des cours de math qu'il faut chercher

Publicité

Page : 1 2

Page Suivante

Page Précédente

Haut de page

FORUM HardWare.fr

Programmation

Algo

Inverse d'une matrice

Sujets relatifs
Mettre une image de type BMP ou PCX dans une matrice?? (Newbies)	utilisation de var et boucle for en sens inverse
[C++Builder] Matrice de structures	Determinant d'une matrice[n][n]???
[C++] Comment accèder à la fonction inverse de cosinus	[c] probleme initialisation matrice
[C] Matrice Dynamique URGENT svp ...	[MPEG2] PTS inversé ???
[openGL] Que contient une matrice de visualisation-modelisation ?	Comment charger et recuperer une image sous forme de matrice ?
Plus de sujets relatifs à : Inverse d'une matrice

Page générée en 0.106 secondes