C++

inversion matricielle

Recherche :

Mot : Pseudo : Filtrer
Page : 1 2 Page Suivante Page Précédente Bas de page
Auteur	Sujet : inversion matricielle

Elmoricq

Modérateur

Reprise du message précédent :
De toute façon, même sans %e manquant, c'est horrible de faire comme ça.
Horrible question élégance, mais aussi (et surtout) parce que c'est casse-gueule.
Aucun contrôle d'erreur, possibilité d'obtenir absolument ce qu'on veut...

Mieux vaut un bon vieux fgets() et on travaille sur la chaîne lue à partir de là.

Publicité

skelter

ben ouai, fgets + strtod

el muchacho

Comfortably Numb

nicosmos a écrit :

bon effectivement je n ai pas suivi les instructions, donc honte a moi!!
effectivement ce n est pas du c++ mais bien du C, enfin, si je comprends bien que le C++ est de la prog orienté objet??
le prog tourne quand meme mais avec une initialisation par "random"de la matrice.
je dois avoir un problem sur la lecture du fichier contenant les valeurs numeriques. j utilise "fscanf".j'affiche bien les valeurs de la matrice a inverser mais lorsque je l inverse elle me sort des nan nan...
PS : j ai corrigé le %e qui manquait a l appel.
voila.

merci

C'est normal, tu passes tout sur la pile (de plus, ta syntaxe me parait louche, ça ne peut pas compiler comme ça), elle explose. Le C est chouillas plus sophistiqué que le Pascal au niveau gestion mémoire. Il y a déjà eu une explication à ce sujet sur le forum. Fais une recherche sur "pile tas".

Message édité par el muchacho le 21-05-2005 à 07:55:32

chicotruss

nicosmos a écrit :

Bonjour,
j aai essayé Gauss jordan mais sans réel succès!!

La methode de Gauss necessite de prendre quelques précautions sinon elle renvoid des résultats completements faux.
De mémoire, il faut que le pivot soit une valeur la plus grande possible, ce qui implique des inversions de lignes et|ou colonnes au préalable.

Perso, j'utilise la methode de Cramer (seule restriction, le determinant ne doit pas etre nul), par contre c'est beaucoup plus long pour des matrices de taille a partir de 9x9 (sur mon implémentation)

pour les formules par exemple: http://membres.lycos.fr/emauvais/idm/MatOpe.htm

nicosmos

Elmoricq a écrit :

De toute façon, même sans %e manquant, c'est horrible de faire comme ça.
Horrible question élégance, mais aussi (et surtout) parce que c'est casse-gueule.
Aucun contrôle d'erreur, possibilité d'obtenir absolument ce qu'on veut...

Mieux vaut un bon vieux fgets() et on travaille sur la chaîne lue à partir de là.

oui mais un fgets est pour une chaine de caractere et là, j'ai des nombres positif ou negatifs, donc pas le meme nombre d'octets.
quel controle est ce que je peux mettre? je veux bien en mettre!!

nicosmos

chicotruss a écrit :

je ne suis pas avard en temps, je veux bien essayé autre chose??
est ce que tu as un source avec la methode de cramer en C ??
merci

chicotruss

nicosmos a écrit :

je ne suis pas avard en temps, je veux bien essayé autre chose??

Je precise quand meme que dans la cas d'une matrice 10x10 de float ca prend environ 1min sur un 3GHZ l'inversion par Cramer.

Citation :

est ce que tu as un source avec la methode de cramer en C ??

En C non mais je peu te mettre les formules de math, c'est recursif et ca passe tout seul

Code :

pour une matriice carrée uniquement
Si determinant(matrice originale) != 0
alors c'est inversible
matrice inverse = transposee(comatrice(matrice originale)) / determinant(matrice originale) //division terme à terme
transposee:
boucle sur tous les elemnents de la matrice
element(i, j) devient element(j, i) dans la matrice transposee
comatrice:
pour chaque element(i, j) de la matrice originale correspond un cofacteur(i, j) = (-1)^(i+j) * determinant(mineur(i, j))
avec i et j commencant à 1 donc si i et j commencent à 0 tu remplace (-1)^(i+j) par (-1)^(i+j+2) (mais bon c'est juste au niveau logique, ca change rien au resultat)
mineur:
mineur(i, j) = matrice sans les lignes i et j
determinant:
somme pour j fixé quelconque (1 par exemple) de((-1)^(i+j) * determinant(mineur(i, j)) * element(i, j))
avec i et j commencant à 1
donc recursif -> condition d'arret: quand nb_lignes (= nb_colonnes) = 2 -> alors determinant = element(1, 1)*element(2, 2) - element(2, 1)*element(1, 2)
Attention (-1)^(i+j) signifie (-1) puissance (i+j) et pas (-1) * 10 puissance (i+j)
-> cela dit, pour optimiser tu dois pouvoir remplacer par (((i+j)&0x01)?(-1):(1))
-> voila, ca passe tout seul si tu fais une classe avec ces methodes

Message édité par chicotruss le 26-05-2005 à 09:46:28

skelter

nicosmos a écrit :

oui mais un fgets est pour une chaine de caractere et là, j'ai des nombres positif ou negatifs, donc pas le meme nombre d'octets.

on a dis fgets + strtol (ou strtod)

cat stevens

dnner moi le programme en c du calcul de l'inverce d'une matrice
merci

push

/dev/random

cat stevens a écrit :

dnner moi le programme en c du calcul de l'inverce d'une matrice
merci

et un mars dans le cul avec ça ?

Publicité

cat stevens

connassez vous la methode gausse seideel

nicosmos

cat stevens a écrit :

connassez vous la methode gausse seideel

non je ne connais pas, dis m en plus...
je suis toujours la recherche d un code d 'inversion de matrice n*n qui fonctionne.
voil la fonction que j utilise mais qui ne fonctionne pas;je l ai faite a partir d un source en pascal mais ici il y un soucis de division par zero.

Code :

void
invers(double a[nmax][nmax])
{
int k,j,i;
double b[nmax][nmax],tot;
for (k=0;k<nmax;k++)
{
for(j=0;j<nmax;j++)
b[k][j]=0;
b[k][k]=1;
}
// b est initialisé a la matrice identité
for (k=0;k<nmax;k++)
{
tot=a[k][k];
for(j=0;j<nmax;j++)
{ a[k][j]=a[k][j]/tot;
b[k][j]=b[k][j]/tot;
}
for (i=0;i<nmax;i++)
{
if (i!=k)
{
tot=a[i][k];
for (int j=0;j<nmax;j++)
{
a[i][j]=a[i][j]-tot*a[k][j];
b[i][j]=b[i][j]-tot*b[k][j];
}
}
}
}
//on inverse la matrice
for(int z=0;z<nmax;z++)
{
for (int s=0;s<nmax;s++)
a[z][s]=b[z][s];
}
}

nicosmos

push a écrit :

et un mars dans le cul avec ça ?

je crois bien que si c est pour faire ce genre de reflexion tu peux ne rien faire ce sera mieux.
par contre si tu deviens constructif et que tu me trouves un prog d'inversion matricielle en C pour des matrice de rang n>100 ... je te serais reconnaissant...
si je suis sur ce forum c est quand meme pour que l on m'aide.
d'ailleurs je suis sur l inversion matricielle des numerical recipices et je n'arrive pas a compiler leur inversion gauss jordan... si...
bref
merci
nico

blastman

just me !

nicosmos a écrit :

je crois bien que si c est pour faire ce genre de reflexion tu peux ne rien faire ce sera mieux.
par contre si tu deviens constructif et que tu me trouves un prog d'inversion matricielle en C pour des matrice de rang n>100 ... je te serais reconnaissant...
si je suis sur ce forum c est quand meme pour que l on m'aide.
d'ailleurs je suis sur l inversion matricielle des numerical recipices et je n'arrive pas a compiler leur inversion gauss jordan... si...
bref
merci
nico

tu acceptera son mars

---------------
http://www.blastmanu.info

nicosmos

blastman a écrit :

tu acceptera son mars

...

blastman

just me !

ça veut dire oui ?

vas y push trouve lui un prog. d'inversion matricielle en C pour des matrice de rang n>100 comme ça tu aura son autorisation pour le sodomiser

fouette le de ma part

:lol:

Message édité par blastman le 02-06-2005 à 15:38:55

---------------
http://www.blastmanu.info

nicosmos

...

FrenchFrogger

http://antony-ibanez.eg2.fr/

nicosmos a écrit :

Code :

.....
for (k=0;k<nmax;k++)
{
for(j=0;j<nmax;j++)
b[k][j]=0;
b[k][k]=1; // Aucune erreur ici, mais je pense que tu peux optimiser ces instructions...
}
.....
for (k=0;k<nmax;k++)
{
tot=a[k][k]; // tot peut donc être égal à 0?
for(j=0;j<nmax;j++)
{ a[k][j]=a[k][j]/tot; // Si tot = 0, ça coince...
b[k][j]=b[k][j]/tot;
}
......

Code :

Pour inverser une matrice, si l'on parle de la même chose, j'ai vu cette méthode:
A = (e1 e2) matrice dont ont connait la valeur des n*n membres.
(e3 e4)
B = inverse de A =
(x1) si A a n*n éléments, cette matrice doit avoir 1*n variables (x1,...,xn)
(x2) (en fait elle en a n*n éléments mais on verra plus tard)
Puisque B est l'inverse de A,
A*B = B*A = (y1) on s'enfout des valeurs de y1 et y2...
(y2)
Donc A*B = (x1*e1 + x2*e2) = (y1) => x1*e1 + x2*e2 = y1
(x1*e3 + x2*e4) = (y2) => x1*e3 + x2*e4 = y2
Puis, par exemple, si e1 = 2, e2 = 3, e3 = -2, e4 = 9, on a:
2 * x1 + 3 * x2 = y1 => y1 + y2 = 12 * x2
-2 * x1 + 9 * x2 = y2 => 3 * y1 - y2 = 4 * x1
=> x1 = 3/4 * y1 - 1/4 * y2 => B = ( 3/4 -1/4 )
x2 = 1 / 12 * y1 + 1 / 12 * y2 ( 1/12 1/12 )

Je sais pas si tu as compris...

Message édité par FrenchFrogger le 07-06-2005 à 13:27:31

---------------
http://antony-ibanez.eg2.fr/

Publicité

Page : 1 2

Page Suivante

Page Précédente

Haut de page

FORUM HardWare.fr

Programmation

C++

inversion matricielle

Sujets relatifs
inversion de liste chaine	Inversion de couleurs [Résolu]
[algo] inversion d'une matrice, cas "particulier"	Produit de 2 matrices de format différent et inversion de matrice
Plus de sujets relatifs à : inversion matricielle

Page générée en 0.070 secondes