La résolution par la méthode de gauss se déroule en deux étapes.
Tout d'abord dans la transformation de la matrice en une matrice triangulaire supérieure. Ax = b => A'x=b'
Ensuite dans la résolution du systeme linéaire A'x = b' ou A' est une matrice triangulaire supérieure.
La seconde partie est très simple.
L'algorithme est appelé algorithme de la remontée. Il consiste sur un couple de matrice de ce genre :
A' b'
1 2 3 4|| 5
0 6 8 7|| 6
0 0 2 1|| 3
0 0 0 1|| 1
La derniere composante du vecteur x est facile a obtenir ( 1*x4 = 1).
Tu remonte d'un cran : (2*x3+1*x4=3) <=> x3 = (3-1*x4)/2
Et ainsi de suite.
xn = b'n-sum(n+1=>Taille(A))[A2j*xj]
En francais, ta valeur xn c'est la valeur de b' de la meme ligne moins tout les x qui lui sont supérieurs affectées par leur multiplicateur.
Pour la triangulation de la matrice de gauss, tu peux la trouver sur le net facilement normalement, je te laisse chercher, j'ai la flemme
Message édité par henker92 le 09-02-2011 à 17:30:13