Université Mohamed BOUDIAF de M’SILA Année universitaire 2023/2024
ème
Faculté des Mathématiques et de l'Informatique 2 Année Licence Informatique
Département de l’informatique Module : Méthodes numériques
Corrigé type de l’examen
Exercice 1 : (7 points)
1. La forme matricielle AX=b. (1 points)
1 0 −2 7 𝑥1 11
2 −1 3 4 𝑥2 9
( ) (𝑥 ) = ( )
3 −3 1 5 3 8
2 1 4 4 𝑥4 10
︸ ︸ ︸
𝐴 𝑥 𝑏
2. La résolution du système (S) par la méthode de Gauss. (4 points)
1 0 −2 7 11 1 0 −2 7 11
4 9 𝐿 ← 𝐿 − 2𝐿
Etape 1 : (2 −1 3 | ) 2 2 1
; Etape 2 : (0 −1 7 −10| −13) 𝐿 ← 𝐿 − 3𝐿
3 −3 1 5 8 𝐿3 ← 𝐿3 − 3𝐿1 0 −3 7 −16 −25 3 3 2
2 1 4 4 10 𝐿4 ← 𝐿4 − 2𝐿1 0 1 8 −10 −12 𝐿4 ← 𝐿4 + 𝐿2
1 0 −2 7 11 1 0 −2 7 11
0 −1 7 −10 −13 0 −1 7 −10 −13
Etape 3 : ( | ) ; On obtient : ( | )
0 0 −14 14 14 15
0 0 −14 14 14
0 0 15 −20 −25 𝐿4 ← 𝐿4 + 14 𝐿3 0 0 0 −5 −10
𝑥1 − 2𝑥3 + 7𝑥4 = 11 𝑥1 = −1
−𝑥 + 7𝑥3 − 10𝑥4 = −13 𝑥 =0
(𝑆) ⇔ { 2 Donc la solution de (S) est: { 2
−14𝑥3 + 14𝑥4 = 14 𝑥3 = 1
−5𝑥4 = −10 𝑥4 = 2
3. La décomposition de la matrice A, (A=LU). (1 points)
La matrice 𝑈 = 𝐴̂ obtenue dans (Q2), alors,
1 0 −2 7 1 0 0 0
0 −1 7 −10 2 1 0 0
𝑈=( ) ; et L= 3 3 1 0
0 0 −14 14
−15⁄14
0 0 0 −5 (2 −1 1)
4. Le déterminant de la matrice A. (1 points)
det(𝐴) = det(𝐿)det(𝑈) = 1(1)(−1)(−14)(−5) = −70.
Exercice 2 : (6 points)
1. Vérification que les processus itératifs de Jacobi et Gauss-Seidel convergent, quelque soit le vecteur initial
𝑿(𝟎) ∈ ℝ𝟑 . (1 points)
|0,6| > |0,2| + |−0,1|
On a : { |0,4| > |0,1| + |0,1| . Donc la matrice A est à diagonale strictement dominante par ligne. Alors, les processus
|0,5| > |−0,2| + |0,1|
itératifs de Gauss-Seidel associés au système linéaire AX=b convergent, quelque soit le vecteur initial 𝑋 (0) ∈ ℝ3 .
2. les processus itératifs de Gauss-Seidel associés au système linéaire AX=b: (1 points)
(𝑘+1) (𝑘) (𝑘)
𝑥1 = (0,7 − 0,2𝑥2 + 0,1𝑥3 )/0,6
(𝑘+1) (𝑘+1) (𝑘)
𝑥2 = (1,2 − 0,1𝑥1 − 0,1𝑥3 )/0,4
(𝑘+1) (𝑘+1) (𝑘+1)
{𝑥3 = (1,5 + 0,2𝑥1 − 0,1𝑥2 )/0,5
Page 1 sur 2
3. Les quatre (4) premiers itérés de la méthode de Gauss-Seidel : (4 points)
K 0 1 2 3 4
x1 0.9 1.0167 0.9935 0.9992 0.9996
x2 1.9 2.0208 2.0010 2.0009 2.0002
X3 2.9 3.0025 2.9972 2.9995 2.9998
L'erreur 0.1208 0.0232 0.0057 0.0007
Exercice 3 : (7 points)
1. L’algorithme qui permet de résoudre un système linéaire par la méthode de Gauss. (4 points)
Algorithme SL_Methode_Gauss
Variables :
A : Tableau(N×N) des réels
X,b : Tableaux(N×1) des réels
N,k,i,j : Entiers
S,α : Réels
Début
Ecrire (‘Entrer la matrice A, et le vecteur b’)
Lire (A, b)
Pour k=1 jusqu’à N-1 (pas +1) faire
Pour i=k+1 jusqu’à N (pas +1) faire
α A(i,k)/A(k,k)
b(i) b(i)-α*b(k)
Pour j=k+1 jusqu’à N (pas +1) faire
A(i,j) A(i,j)-α*A(k,j)
Fin pour(j)
Fin pour(i)
Fin pour(k)
X(N) b(N)/A(N,N)
Pour i=N-1 jusqu’à 1 (pas -1) faire
S 0
Pour j=i+1 jusqu’à N (pas +1) faire
S S+A(i,j)*X(j)
Fin pour(j)
X(i) (b(i)-S)/A(i,i)
Fin pour(i)
Ecrire (‘La solution est X= ‘, X)
Fin
2. Calcul du nombre d'opérations arithmétiques (Calcul de complexité arithmétique). (3 points)
……………………………………………………………………………
Page 2 sur 2