UNIVERSITÉ ALIOUNE DIOP DE BAMBEY UFR SATIC

LICENCE 2 DE MATHS-PHYSIQUE-INFORMATIQUE

FICHE 1 DE CALCUL NUMÉRIQUE

Les méthodes Directes pour les systèmes linéaires et la théorie du conditionnement

Exercice 1. Détermination de LU
Soit A ∈ Mn (R) une matrice dont les sous-matrices principales sont inversibles. D’après
le théorème du cours, il existe une unique matrice L ∈ Mn (R) triangulaire inférieure avec
lii = 1, ∀i = 1, ..., n et une unique matrice U ∈ Mn (R) triangulaire supérieure telle que
A = LU .
1. Montrer que pour i = 1, ∀j = 1, ..., n u1j = aij (première ligne de U).

2. Vérifier que u11 6= 0.

aij
3. Montrer que pour j = 1, ∀i = 2, ..., n li1 = u11
(première colonne de L).

4. Montrer par récurrence que la i − eme ligne de U et la j − eme colonne de L peuvent

être exprimées selon :

i−1
X

aij − lik ukj , ∀j = i, ..., n

pour i allant de 1 a n : uij =
 k=1
0, ∀j = 1, ..., i − 1




 0, ∀i = 1, ..., j − 1
1, ! i=j


pour j allant de 1 a n : lij = j−1
X
1
aij − lis usj , ∀i = j + 1, . . . , n


ujj


s=1

Exercice 2. Soient α ∈ R et A ∈ M3 (R), b ∈ R3 définis par

   
1 −α 1 1
A =  −α 2α2 −α − α2  , b =  −α  . (0.1)
1 −α − α 1 + α2 + α4
2
1 + α2

1. Pour α 6= 0, résoudre Ax = b par la méthode d’élimination de Gauss.

2. En déduire le déterminant de A.

3. Déterminer la factorisation A = LU , lii = 1, i = 1, 2, 3.

4. Résoudre Ax = b par la factorisation A = LU .

   
1 1 1 1 0
 1 2 3 4 
 et b ∈ R4 tel que b =  0 .
 
Exercice 3. Soit A =   1 3 6 10   1 
1 4 10 20 0
1. Triangulariser le système Ax = b et en déduire le déterminant de A et le vecteur x ∈ R4 .
2. Montrer que A est symétrique définie positive.

TRAVAUX DIRIGÉS DE CALCUL NUMÉRIQUE 1 L2 MPI M. S. DIALLO, A. DIOP et S. SHAMS c 2022-2023

En déduire qu’il existe un unique couple (L, U ) telle que A = LU où L est une matrice
triangulaire inférieure et U une matrice triangulaire supérieure.
3. Nous rappelons que la méthode de Doolittle (lii = 1) fournit d’abord la i-ième ligne de U,
puis la i-ième colonne de L selon
i−1
X
pour i = 1, . . . , n uij = aij − lik ukj j = i, . . . , n (i ≤ j)
k=1
j−1
!
X
et pour j = 1, . . . , n lij = u1jj aij − lis usj i = j + 1, . . . , n (i > j).
s=1
a) Utiliser cet algorithme pour déterminer les matrices L et U .
b) Reprendre la résolution de Ax = b en utilisant la décomposition LU .
c) Sans calculer explicitement l’inverse de la matrice, dire que représente la solution x ∈ R4
pour A−1 .

Exercice 4. On considère la matrice A et le vecteur b suivants :

   
2 −1 0 1
A =  −1 3 −1  , b =  2  .
0 −1 2 1

1. Déterminer la matrice triangulaire supérieure en utilisant la méthode d’élimination de

Gauss sans imposer des 1 sur la diagonale.

2. En déduire le déterminant de la matrice A.

3. Déterminer la factorisation de LU de la matrice A par la méthode de Croût.

On rappelle que pour i = 1, ..., n
i−1
X
lki = aki − lkr uri k = i, ..., n (0.2)
r=1

i−1
!
1 X
uij = aij − lir urj j = i + 1, ..., n. (0.3)
lii r=1

4. Calculer le déterminant en utilisant la décomposition LU .

5. Résoudre Ax = b.

6. Montrer que la matrice A est symétrique définie positive.

Rappeler le théorème de la décomposition de Choleski.

7. Soit D une matrice diagonale inversible définie par D = diag(U ) (Dii = uii 6= 0, ∀i) et
R = D−1 U (matrice triangulaire supérieure).

a) Montrer que la matrice R est à diagonale unité.

b) Montrer que A admet une factorisation LDLT , Lii = 1, la déterminer.

c) En déduire la factorisation de Choleski A = CC T , Cii > 0, où C est une matrice

triangulaire inférieure.

TRAVAUX DIRIGÉS DE CALCUL NUMÉRIQUE 2 L2 MPI M. S. DIALLO, A. DIOP et S. SHAMS c 2022-2023

Exercice 5. Soit  > 0 et A ∈ M4,4 (R) définie par
 
1+ 1 1 1
 1 1 +  1 1 
A=  1
. (0.4)
1 1+ 1 
1 1 1 1+

1. Calculer k A k∞ .

2. Calculer A−1 en utilisant les co-matrices.

3. Calculer Cond∞ (A).

4. Pour quelles valeurs de , A est-elle mal conditionnée ? bien conditionnée ?

Exercice 6. Après avoir modélisé un problème de diffusion stationnaire, nous avons obtenu
le système suivant Au = b où u ∈ Rn le vecteur des inconnues de ce problème , b ∈ Rn , et A
une matrice carrée d’ordre n à coefficients dans R définie par
 
2 −1 0 . . . . . . 0
 .. 
 −1 2 −1 0 . 0 
 .. 
1  0 −1 2 −1 0

. 

A= 2 . . .
h  .. . . . . . . . . . . . ... 
 
 .. . . 
 . ... . −1 2 −1 
0 . . . . . . 0 −1 2
1
Avec h = n+1 , soit k . k une norme matricielle induite.
1. Montrer que la matrice A × h2 est symétrique définie positive.
2. Pour u ∈ Rn , k u k= max{| ui | i ∈ {1, ..., n}}.
En admettant que k u k≤ (1/8) k b k, donner une majoration de k A−1 k∞ .
3. Montrer que le conditionnement de A est exprimée par
1
K∞ (A) ≤ .
2h2
4. Soient b, δb ∈ Rn et soient u, δu ∈ Rn tels que Au = b et A(u + δu) = b + δb.
Montrer que
k δu k k δb k
≤ K(A) .
kuk kbk
5. En perturbant cette fois ci la matrice A, nous obtenons le système suivant (A + δA)(u +
δu) = b où δA ∈ Mn (R) et u, δu ∈ Rn .
a) Montrer que δu satisfait la relation (I + A−1 δA)δu = −A−1 δAu où I est une matrice
identité.
b) On suppose que k A−1 δA k≤ 1, et (I + A−1 δA) est inversible.
Montrer que k I k= 1, puis établir que
1
k (I + A−1 δA)−1 k≤ .
1− k A−1 δA k
kA−1 k
c) En déduire que k δu k≤ 1−kA−1 δAk
k δA kk u k.

BON COURAGE!

TRAVAUX DIRIGÉS DE CALCUL NUMÉRIQUE 3 L2 MPI M. S. DIALLO, A. DIOP et S. SHAMS c 2022-2023