Introduction to Numerical Analysis
TP 3
Question 1
Soit le système Ax = b où :
! " ! "
0, 780 0, 563 0, 217
A= ; b=
0, 913 0, 659 0, 254
dont la solution est ! "
1
x=
→1
On propose deux approximations de la solution :
! " ! "
0, 999 0, 341
x1 = ; x2 =
→1, 001 →0, 087
Calculer les résidus et interpréter le résultat.
Question 2 (Janvier 2012 )
4 1 →2 0.5 →0.5 0
Soit A = 2 1 →2 dont l’inverse est A→1 = →228 207 →2 .
→249 →11 21 →113.5 102.5 →1
(a) Calculer le nombre de conditionnement de A en norme ↑.
(b) On souhaite résoudre le système Ax = b par élimination de Gauss, où b = (1.02 4.75 53.2) est
connu avec trois chi!res significatifs. Quelle est la précision sur le résultat ?
(c) On équilibre les équations. Quelle sera la nouvelle précision sur le résultat ?
(d) On utilise à présent un préconditionneur basé sur la somme des valeurs absolues des éléments de
chaque ligne, c’est-à-dire que l’on prémultiplie A par D = diag( 17 , 15 , 281
1
). Quelle est la nouvelle
précision sur le résultat ?
Question 3
Soit le système linéaire
! "! " ! "
0.1036 0.2122 x1 0.7381
= .
0.2081 0.4247 x2 0.9327
Calculez la solution de ce système par élimination de Gauss en arrondissant chaque résultat in-
termédiaire à quatre chi!res significatifs. Comparez le résultat obtenu à la résolution avec une
double précision.
Questions supplémentaires
Question 4
Soient
1 1 1 1.1223 →0.1235 0.0011
A = 1 10 100 ; dont l’inverse est A→1 = →0.1235 0.1247 →0.0012 ;
1 100 10000 0.0011 →0.0012 0.0001
1, 0000
1 1 1
b = 0, 3642 ; D = diag( ; ; )
3 111 10101
0, 8525
Les éléments de b sont correctement arrondis. Ceux de A et D sont exacts.
Calculer la précision relative que l’on a sur x dans les cas suivants :
(a) Ax = b
(b) DAx = Db
Question 5 On considère la matrice
1 0.5 3
A = 0.5 25 75
3 75 2500
dont l’inverse vaut
56875 →1025 →37.5
1
A→1 = →1025 2491 →73.5
56250
→37.5 →73.5 24.75
(a) Calculer le nombre de conditionnement, en norme infinie, de A
1 1 1
(b) Calculer le nombre de conditionnement de DA si D = diag( 4.5 ; 100.5 ; 2578 )
3.25
(c) Si b = 6.91 (arrondi avec 2 chi!res exacts après la virgule) quelle précision relative
8.32
peut-on atteindre sur la solution des systèmes Ax = b et DAx = Db ?