Scribd
Importer
25 vues1 page

Correction Exercice1

Le document présente un système d'équations linéaires et sa représentation matricielle. Il décrit les méthodes de résolution, y compris la méthode directe de Gauss et la méthode itérative de Jacobi, en utilisant MATLAB. Enfin, il compare la convergence des méthodes de Jacobi et de Gauss-Seidel, notant que cette dernière converge généralement plus rapidement.

Transféré par

mohamedaliammar035
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd
25 vues1 page

Correction Exercice1

Le document présente un système d'équations linéaires et sa représentation matricielle. Il décrit les méthodes de résolution, y compris la méthode directe de Gauss et la méthode itérative de Jacobi, en utilisant MATLAB. Enfin, il compare la convergence des méthodes de Jacobi et de Gauss-Seidel, notant que cette dernière converge généralement plus rapidement.

Transféré par

mohamedaliammar035
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd

Correction - Exercice 1

Soit le système d'équations linéaires suivant :

10x₁ - x₂ + 2x₃ = 6
-x₁ + 11x₂ - x₃ + 3x₄ = 25
2x₁ - x₂ + 10x₃ - x₄ = -11
3x₂ - x₃ + 8x₄ = 15

1. Forme matricielle

Ce système peut s’écrire sous la forme AX = b, où :


A = [10 -1 2 0; -1 11 -1 3; 2 -1 10 -1; 0 3 -1 8]
X = [x₁; x₂; x₃; x₄]
b = [6; 25; -11; 15]

2. Méthode directe de Gauss avec MATLAB

On utilise la fonction gauss vue en TP pour résoudre le système :


A = [10 -1 2 0; -1 11 -1 3; 2 -1 10 -1; 0 3 -1 8];
b = [6; 25; -11; 15];
X = gauss(A, b);

3. Méthode itérative de Jacobi

Utilisation de la fonction Jacobi :


X0 = [0; 0; 0; 0]; % Vecteur initial
nmax = 40; % Nombre max d'itérations
tol = 1e-4; % Tolérance

[X, niter] = Jacobi(A, b, X0, nmax, tol);

4. Convergence : Jacobi vs Gauss-Seidel

En général, la méthode de Gauss-Seidel converge plus rapidement que la méthode de Jacobi,


car elle utilise les valeurs les plus récentes dans les calculs, ce qui accélère la convergence.

Vous aimerez peut-être aussi