Universite de Nice Sophia-Antipolis
Licence L3 Mathematiques Annee 2008/2009
Analyse Numerique
TD 6
EXERCICE 1
Matrices diagonales, triangulaires
1.1 Matrices diagonales
Soit D = (d
ii
)
i=1,...,n
une matrice diagonale dordre n > 0. Donner une condition necessaire
et susante pour que D soit inversible.
1.2 Matrices triangulaires inferieures
Soit L une matrice triangulaire inferieure dordre n > 0.
a. Sous quelle condition necessaire et susante L est-elle inversible ?
b. On suppose que la matrice triangulaire inferieure L est inversible. Soit b un vecteur
colonne ayant n composantes.
Donner un algorithme qui permet de resoudre lequation dinconnue y :
Ly = b .
Quel est le co ut de cet algorithme en termes doperations elementaires (additions, multi-
plications, divisions) ?
1.3 Matrices triangulaires superieures
On consid`ere une matrice triangulaire superieure U dordre n > 0.
a. Donner une condition necessaire et susante pour que U soit inversible.
b. On suppose que la matrice triangulaire superieure U est inversible. Soit y un vecteur
colonne donne ayant n composantes.
Ecrire un algorithme qui permet de resoudre lequation dinconnue x :
U x = y .
Donner la complexite de cet algorithme.
EXERCICE 2
Methode delimination de Gauss
2.1 Des exemples
Eectuer une elimination de Gauss sur les syst`emes lineaires suivants
2 4 4
1 3 1
1 5 6
x
1
x
2
x
3
2
1
6
,
1
Universite de Nice Sophia-Antipolis
Licence L3 Mathematiques Annee 2008/2009
1 0 6 2
8 0 2 2
2 9 1 3
2 1 3 10
x
1
x
2
x
3
x
4
6
2
8
4
.
2.2 Cas general
On consid`ere maintenant le cas general dun syst`eme lineaire Ax = b.
a.
Ecrire un algorithme de resolution de ce syst`eme par la methode de Gauss.
b. Donner la complexite de cet algorithme.
EXERCICE 3
Factorisation LU
3.1 Un exemple
On revient sur la premi`ere matrice donnee dans lexercice 2 :
2 4 4
1 3 1
1 5 6
.
Eectuer une factorisation LU de cette matrice o` u L est une matrice triangulaire inferieure
ayant des 1 sur sa diagonale et U est une matrice triangulaire superieure.
3.2 Cas general
a. Montrer que le produit de deux matrices triangulaires inferieures de meme ordre est
une matrice triangulaire inferieure.
b. Soit L une matrice triangulaire inferieure et inversible.
Montrer que son inverse L
1
est egalement une matrice triangulaire inferieure.
c. Soit A une matrice carree reguli`ere possedant une decomposition LU, avec L une
matrice triangulaire inferieure ayant des 1 sur sa diagonale et U une matrice triangulaire
superieure.
Montrer que cette factorisation A = LU est unique.
3.3 Factorisation LU dune matrice tridiagonale
Soit A une matrice tridiagonale inversible
a
i1 i
= a
i
, a
i i
= b
i
, a
i i+1
= c
i
Ecrire un algorithme de factorisation LU de A.
Quelle est la complexite de cet algorithme ?
2
Universite de Nice Sophia-Antipolis
Licence L3 Mathematiques Annee 2008/2009
Application : eectuer la decomposition A = LU de la matrice
2 1 0 0 0
4 5 2 0 0
0 3 1 1 0
0 0 2 4 1
0 0 0 2 2
.
EXERCICE 4
Localisation des valeurs propres dune matrice
On introduit la denition suivante.
Denition 4.1. Soit A = (a
kj
)
k,j=1,...,n
une matrice carree dordre n.
On appelle disque de Gerschgorin centre en a
kk
lensemble
D
k
= {z C/ |z a
kk
|
n
j=1
j=k
|a
kj
|} .
On donne le theor`eme suivant qui sera demontre dans la premi`ere partie de cet exercice.
Theor`eme (theor`eme de Gerschgorin)
Soit A une matrice carree dordre n.
Les valeurs propres de A appartiennent `a lunion
des n disques de Gerschgorin du plan complexe :
valeur propre de A k {1, ..., n} , D
k
.
4.1 Demonstration du theor`eme
a. Soit une valeur propre de A et u un vecteur propre associe `a cette valeur propre.
Soit k tel que |u
k
| = max
1jn
|u
j
|.
Montrer que | a
kk
|
n
j=1
j=k
|a
kj
|.
Conclure.
4.2
Etude dun exemple
On consid`ere la matrice suivante
A =
1 + i i 2
3 2 + i 1
1 i 6
.
a. Dessiner les 3 disques de Gerschgorin et localiser les valeurs propres de A.
3
Universite de Nice Sophia-Antipolis
Licence L3 Mathematiques Annee 2008/2009
b. En remarquant que A et
t
A ont les memes valeurs propres, representer les disques
de Gerschgorin associes aux valeurs propres de
t
A.
c. Donner une majoration des valeurs absolues des valeurs propres de A.
4.3 Matrice `a diagonale dominante
Denition 4.2. Soit A = (a
kj
)
k,j=1,...,n
une matrice carree dordre n.
On dit que A est `a diagonale dominante si k , 1 k n, |a
kk
|
n
j=1
j=k
|a
kj
| .
La matrice A est dite `a diagonale strictement dominante si
k , 1 k n, |a
kk
| >
n
j=1
j=k
|a
kj
| .
a. Montrer, en utilisant le theor`eme de Gerschgorin, que si A est `a diagonale strictement
dominante alors elle est inversible.
b. Montrer que la matrice suivante est `a diagonale dominante, mais nest pas inversible :
2 1 1
1 3 2
1 0 1
.
4
