L'inverse d'une matrice en langage C
mpiaya
Messages postés
3
Date d'inscription
Statut
Membre
Dernière intervention
-
mpiaya Messages postés 3 Date d'inscription Statut Membre Dernière intervention -
mpiaya Messages postés 3 Date d'inscription Statut Membre Dernière intervention -
Bonjour les amis, SVP y'a il quelqu'un qui peut me traduire cette fonction en C car je cherche la fonction qui calcul l'inverse d'une matrice en C mais je l'ai trouvé en Matlab alors que moi je ne métrise pas ce dernier, merci beaucoup d'avance.
*************************************************************************
determinant=0;
if size(M,1)==1
determinant = M;
else
for p=1:size(M,1)
Mm=M;
Mm(:,1)=[];
Mm(p,:)=[];
determinant = determinant + (-1)^(p+1)*M(p,1)*calcul_determinant(Mm);
end
end
***************************************************************************
function inverse = calcul_inverse(M)
inverse=zeros(size(M,1));
for p=1:size(M,1)
for q=1:size(M,1)
Mm=M;
Mm(p,:)=[];
Mm(:,q)=[];
inverse(p,q) = (-1)^(p+q)*calcul_determinant(Mm);
end
end
inverse = inverse.'/calcul_determinant(M);
****************************************************************************
*************************************************************************
determinant=0;
if size(M,1)==1
determinant = M;
else
for p=1:size(M,1)
Mm=M;
Mm(:,1)=[];
Mm(p,:)=[];
determinant = determinant + (-1)^(p+1)*M(p,1)*calcul_determinant(Mm);
end
end
***************************************************************************
function inverse = calcul_inverse(M)
inverse=zeros(size(M,1));
for p=1:size(M,1)
for q=1:size(M,1)
Mm=M;
Mm(p,:)=[];
Mm(:,q)=[];
inverse(p,q) = (-1)^(p+q)*calcul_determinant(Mm);
end
end
inverse = inverse.'/calcul_determinant(M);
****************************************************************************
A voir également:
- Inverse d'une matrice en c
- Annuaire inversé - Guide
- Ecran inversé - Guide
- Annuaire inversé sfr - Forum SFR
- Slash inversé ✓ - Forum Clavier
- Mac comment taper l'inverse d'un slash (/) ? ✓ - Forum MacOS
2 réponses
Matlab travaille directement avec des variables matrices, ce qu'il faut implémenter à la main avec des pointeurs. C'est long à faire, en revanche t'expliquer ce que fait le code Matlab est plus rapide :
Explication de la sélection des cases des matrices :
Mm(_,_)
avant la virgule on sélectionne les lignes, après la virgule on sélectionne les colonnes
exemple:
2:4 -> de la 2e à la 4e (ligne ou colonne selon si c'est à gauche ou à droite de la virgule)
: -> toutes (les lignes ou les colonnes selon si c'est à gauche ou à droite de la virgule)
Le code situé avant la fonction inverse dans ton post présente le calcul de déterminant
function inverse = calcul_inverse(M) // M est la matrice passée en paramètre
// et c'est la variable inverse qui est renvoyée par la fonction Matlab à la fin de la fonction
inverse=zeros(size(M,1));
// size(M,1) renvoie le nombre de lignes de la matrice M
// size(M,2) renverrait le nombre de colonnes de la matrice M
// zeros(N) renvoie la matrice nulle carrée NxN
for p=1:size(M,1) // p va de 1 au nombre de lignes de M
for q=1:size(M,1) // p va de 1 au nombre de lignes de M
// attention ce code ne marche que pour une matrice carrée sinon on ferait for q=1:size(M,2)
Mm=M; // copie de la matrice M (toute une procédure en C)
Mm(p,:)=[]; // suppression de la p-ième ligne de la matrice Mm
Mm(:,q)=[]; // suppression de la q-ième colonne de la matrice Mm
inverse(p,q) = (-1)^(p+q)*calcul_determinant(Mm);
// dans la case à la p-ième ligne et la q-ième colonne on calcule et on range le coefficient
end
end
inverse = inverse.'/calcul_determinant(M);
Explication de la sélection des cases des matrices :
Mm(_,_)
avant la virgule on sélectionne les lignes, après la virgule on sélectionne les colonnes
exemple:
2:4 -> de la 2e à la 4e (ligne ou colonne selon si c'est à gauche ou à droite de la virgule)
: -> toutes (les lignes ou les colonnes selon si c'est à gauche ou à droite de la virgule)
Le code situé avant la fonction inverse dans ton post présente le calcul de déterminant
