INSA TOULOUSE Analyse Numérique
Méthodes itératives de résolution de systèmes linéaires
Dans cette feuille, on liste les questions de cours/exercices types relatifs au chapitre sur la résolution
de systèmes linéaires à l’aide de méthodes indirectes que vous devez connaı̂tre/savoir faire.
Questions de cours
1. Décrire le principe des méthodes itératives de résolution de systèmes linéaires.
2. A quelle condition une méthode itérative de résolution du système Ax = b basée sur la
décomposition A = M − N converge-t-elle ?
3. Ecrire l’algorithme associé à la méthode de Jacobi.
4. Ecrire l’algorithme associé à la méthode de Gauss-Seidel
5. Ecrire l’algorithme associé à la méthode de relaxation.
6. Donner le critère nécessaire de convergence de la méthode de relaxation.
Exercice.
1 a a
Soit a ∈ R et A = a 1 a .
a a 1
1. Pour quelles valeurs de a, la matrice A est-elle définie positive ?
2. Pour quelles valeurs de a la méthode de Gauss-Seidel est-elle convergente ?
3. Ecrire la matrice J de l’itération de Jacobi.
4. Pour quelles valeurs de a la méthode de Jacobi converge-t-elle ?
5. Ecrire la matrice L1 de l’itération de Gauss Seidel. Calculer ρ(L1 ).
6. Pour quelles valeurs de a la méthode de Gauss-Seidel converge-t-elle plus vite que celle de
Jacobi.
7. Ecrire la matrice Lω liée à la méthode de relaxation et donner des conditions sur a et ω pour
avoir convergence.
Exercice.
Etant donné une matrice A ∈ Mn (R), on rappelle que la norme de A subordonnée à la norme infinie
est donnée par
n
X
kAk∞ = max |Aij |, et |Ax|∞ ≤ kAk∞ |x|∞ , ∀x ∈ Rn .
i=1,...,n
j=1
1
On considère le système Ax = b avec
3 1 0 0 0
1 2 1 0 0
A=
0 2 3 1 0 .
0 0 1 4 3
0 0 0 1 1
1. Donner la décomposition LU de A et en déduire que Ax = b possède une solution unique x∗ .
2. Ecrire une itération du schéma de Gauss-Seidel associé à la résolution de Ax = b.
3. Soit (Xn )n∈N ∈ R5 les itérées du schéma de Gauss-Seidel. On note en = Xn − x∗ . Montrer qu’il
existe a ∈ [0, 1[ tel que
|en+1 |∞ ≤ a |en |∞ , ∀n ∈ N.
En déduire la convergence de la suite.
4. Déterminer la matrice L1 de Gauss-Seidel associée à A et calculer kL1 k∞ . En déduire la conver-
gence de (Xn )n∈N vers x∗ .