U.I.

T Faculté des sciences Filière SMP : S3

Contrôle Final du Module Analyse Numérique

Session de Janvier 2016. De : 10h30 à 12h

Exercice 1 (4 pt.) Donner la décomposition LU de la matrice A suivante

 
2 1 −1
A =  4 6 −1  .
2 −7 3
 
3
Résoudre le système Ax = b où b= 5 
−1

Exercice 2 On donne la fonction f (x) = x + x − 1, on veut approcher son zéro positif

(6 pt.) 2

≈ 0.61803.
√
5−1
xexacte= 2
1. En partant de a = 0.5 et b = 1. Donner les deux premiers intervalles obtenus en utilisant
la méthode, de dichotomie et de Lagrange et préciser l'erreur dans chaque itération.
0 0

2. En prenant x = 1, donner les deux premières itérations obtenues en utilisant la méthode

de Newton et préciser l'erreur dans chaque itération.
0

Exercice 3 (4 pt.) Donner le polynôme interpolant la fonction f aux points (−1, −1), (0, −1),
(1, −1) et (2, 5) en utilisant
1. la base de Lagrange
2. la méthode des fraction divisées de Newton

Soit I = 1 +1 x dx.
Z 1
Exercice 4 (6 pt.) 2

1. Rappeler les formules composites du point milieu, des trapèzes et de simpson ainsi que
0

l'erreur commise dans chacune des méthodes,

2. Combien faut-il de subdivisions de [0, 1] pour évaluer I à 10 près en utilisant la méthode
−3

de Simpson
3. Donner la valeur approchée de I , par la méthode de Simpson composite en utilisant la
subdivision x = 0, x = et x = 1.
0 1
1
2 2

4. Quelle est alors l'ordre de l'erreur commise?

M. Moussa Page 1/1

U.I.T Faculté des sciences Filière SMP : S3

Corrigé du contrôle Final du Module Analyse Numérique

Session de Janvier 2016. De : 10h30 à 12h

Corrigé de l'exercice 1 (4 pt.) On suivant l'algorithme on obtient

   
1 0 0 2 1 −1
4 (1)
L1 =  − 2 1 0 , A = L1 A =
  0 4 1 ,
− 22 0 1 0 −8 4
   
1 0 0 2 1 −1
L2 =  0 1 0  , A(2) = L2 A(1) = L2 L1 A =  0 4 1 .
0 2 1 0 0 6
Donc,    
1 0 0 2 1 −1
( 1 pt.) L =  2 1 0 , et U = 0 4 1  1 pt.) (
1 −2 1 0 0 6
On a bien,     
2 1 −1 1 0 0 2 1 −1
A= 4 6 −1  =  2 1 0  0 4 1  = LU
2 −7 3 1 −2 1 0 0 6
Pour résoudre le système Ax = b , on résout le système échelonné Ly = b et U x = y
successivement. Donc,
      
1 0 0 y1 3  y1 = 3  y1 = 3
 2 1 0   y2  =  5  ⇒ 2y1 + y2 = 5 ⇒ y2 = −1 (1 pt.)
1 −2 1 y3 −1 y1 − 2y2 + y3 = −1 y3 = −6
 

Puis,
      
2 1 −1 x1 3  2x1 + x2 − x3 = 3  x1 = 1
 0 4 1   x2  =  −1  ⇒ 4x2 + x3 = −1 ⇒ x2 = 0 (1 pt.)
0 0 6 x3 −6 6x3 = −6 x3 = −1
 

Finalement on a bien     
2 1 −1 1 3
 4 6 −1   0  =  5 
2 −7 3 −1 −1

Corrigé de l'exercice 2 On pose f (x) = x + x − 1 2

1. Pour la première itération de la dichotomie, on pose a = 0.5 et b = 1, on a f (a ) =

f (0.5) = −0.25 < 0 et f (b ) = f (1) = 1 > 0. On calcule f (
0 0 0
a0 +b0 2
) = f (0.75) = (0.75) +
0.75 − 1 = 0.3125 > 0. On pose a = a = 0.5 et b = = 0.75. ( 0.5 pt.)
0 2
a0 +b0
1 0 1
On a alors 0.5 < x < 0.75 l'erreur est donc inférieure à 0.75 − 0.5 = 0.25 ( 0.5 pt.)
2
exacte
Pour la deuxième itération de la dichotomie, on calcule f ( ) = f (0.625) = (0.625) +
a1 +b1 2

0.625 − 1 = 0.015625 > 0. On pose a = a = 0.5 et b = = 0.625. ( 0.5 pt.)

2
a0 +b0
2 1 2 2

M. Moussa Page 2/2

U.I.T Faculté des sciences Filière SMP : S3

On a alors 0.5 < x exacte < 0.625 l'erreur est donc inférieure à 0.625 − 0.5 = 0.125 ( 0.5
pt.)
Pour la première itération de Lagrange, on a a = 0.5 et b = 1, on pose c = a −
f (a ), donc c = 0.5 − (−0.25) = or f (c ) =
0 0 0 0
b0 −a0 1−0.5 0.5 2
0 0f (0.5) = 0.5 − 0
f ( ) = > 0 . On pose alors a = a = 0.5 et b = c = . ( 0.5 pt.)
f (b0 )−f (a0 ) f (1)−f (0.5) 1−0.25 3
2 1 2
1 0 1 0
On a alors < x < l'erreur est donc inférieure à − = ≈ 0.166666 ( 0.5 pt.)
3 9 3
1 2 2 1 1
2 exacte 3 3 2 6

Pour la deuxième itération de, on pose c = a −

2 1
− b1 −a1 1 3 2 1
f (a ) = − 1 f( ) =
1 f (b1 )−f (a1 ) 1 2 2 1 2
f( ) − f( ) 3 2

or f ( on pose alors a et b . ( 0.5

1
1 6 1 8 8 1 8 2
2
− 13 (− 4 ) = 13 13
) = − 169 <0 2 = c1 = 13 2 = b1 = 3
36
pt.)
On a alors 8
13
< xexacte < 2
3
l'erreur est donc inférieure à 2
3
− 8
13
= 2
39
≈ 0.051282 ( 0.5
pt.)
2. Pour la méthode de Newton, la suite est donnée par x = x n+1 n − f (xn )
f 0 (xn )
,x 0 =1donné.
Donc la première itération donne x = x − = 1− = 1 0
f (x0 )
f 0 (x0 )
1
3
2
3
x1 = 2
3
( 0.5 pt.) et
e1 = |x1 − xexacte | = 0.0486326 ( 0.5 pt.)
la deuxième itération x2 x1 = x1 − ff0(x 1)
= 2
− 1/9
= 2
− 1
= 13
≈ 0.619047 ( 0.5 pt.)
et
(x1 ) 3 7/3 3 21 21
e2 = |x2 − xexacte | = 0.00101 ( 0.5 pt.)

Corrigé de l'exercice 3 1. Le polynôme d'interpolation de Lagrange de degré n sur l'en-

semble des n + 1 points (x , y ) s'écrit i
i=n
i i=0

j=n
i=n
!
X Y x − xj
Pn (x) = yi (0.5 pt.)
i=0 j=0,j6=i
xi − xj

Ici n = 3 on a alors,
x(x − 1)(x − 2) x3 − 3x2 + 2x
L1 (x) = =− (0.25 pt.)
(−1 − 0)(−1 − 1)(−1 − 2) 6
(x + 1)(x − 1)(x − 2) x3 − 2x2 − x + 2
L2 (x) = = (0.25 pt.)
(0 + 1)(0 − 1)(0 − 2) 2
x(x + 1)(x − 2) x3 − x2 − 2x
L3 (x) = =− (0.25 pt.)
(1 + 1)(1 − 0)(1 − 2)) 2
x(x + 1)(x − 1) x3 − x
L4 (x) = = (0.25 pt.)
(2 + 1)(2 − 0)(2 − 1) 6
le polynôme d'interpolation est donné par
P3 (x) = −L1 (x) − L2 (x) − L3 (x) + 5L2 (x) = x3 − x − 1 (0.5 pt.)
2. Pour la méthode des fractions divisées de Newton on a
P3 (x) = a0 + a1 (x + 1) + a2 (x + 1)(x − 0) + a3 (x + 1)(x − 0)(x − 1) on dresse le tableau
suivant,
C'est à dire que P (x) = −1 + x(x + 1)(x − 1) = x
3
3
− x − 1 (2 pt.)

M. Moussa Page 3/3

U.I.T Faculté des sciences Filière SMP : S3

x f (x)
−1 -1 = a0
0 −1
−1 − (−1)
0 − (−1)
= = a1 0
1 −1
−1 − (−1)
1−0
=0 =
0−0
1 − (−1)
= 0 = a2

2 5
5 − (−1)
2−1
=6
6−0
2−0
=3
3−0
2 − (−1)
= 1 = a3

Corrigé de l'exercice 4 1. Soit f : [a, b] → R continue, n ∈ N, n ≥ 2 et h = , on pose b−a

n

x = a, x = b, x = x +kh et x = , pour 1 ≤ k ≤ n. Pour évaluer f (x) dx,

x +x k−1 k R b
0 n k 0 k
on utilise l'une des méthodes numériques2usuelles suivantes la formule composite du
a

point milieu est donnée par

n
X
IPc M (f ) =h f (xk ), (0.5 pt.)
k=1

Si en plus f est de classe C l'erreur commise est estimée par

2

Z b
h2

max |f 00 (x)| (0.5 pt.)
 c

 f (x) dx − IPM (f ) ≤ (b − a)

a
 24 x∈[a,b]

la formule composite des trapèzes est donnée par

n
hX
ITc (f ) = f (xk−1 + f (xk )), (0.5 pt.)
2 k=1

Si en plus f est de classe C l'erreur commise est estimée par

2

Z b  2
 ≤ (b − a) h max |f 00 (x)| (0.5 pt.)
 c

 f (x) dx − IT (f )

a
 12 x∈[a,b]

et la formule composite de Simpson est donnée par

n
hX
ISc (f ) = [f (xk−1 + 4f (xk ) + f (xk )] , (0.5 pt.)
6 k=1

Si en plus f est de classe C l'erreur commise est estimée par

4

Z b
h4

 

 f (x) dx − IS
c
(f )  ≤ max |f (4) (x)| (0.5 pt.)
 2880 x∈[a,b]
a

2. Pour évaluer I à 10 près par la méthode de Simpson composite il sut que l'erreur
−3

commise soit inférieure à 10 . −3

M. Moussa Page 4/4

U.I.T Faculté des sciences Filière SMP : S3

Or dans ce cas (b−a) = 1, et f (x) = − (4) 24

+ 144x+48 96
− (1+x donc max |f (4) (x)| ≤
312.
(1+x2 )3 (1+x2 )4 2 )5 x∈[0,1]

Pour avoir l'erreur désirée on devrait avoir

2880
h4 ≤ × 10−3 ⇒ h ≤ 0.31
312
Il sut de prendre h = pour obtenir la précision voulue. On aura donc besoin de 4
1

points de subdivisions pour cette méthode. (1 pt.)

4

3. La valeur approchée de I , par la méthode de Simpson composite en utilisant la subdivision

x = 0, x = et x = 1 est donnée par h =
0 1
1
2 2
1
2
h
f (0) + f (1) + 2f ( 12 ) + 4 f ( 41 ) + f ( 34 )



ISc (f ) = 6

C'est à dire que,

  
1 1 5 16 16
c
IS (f ) = 1+ +2× +4 + = 0.86039 (1 pt.)
12 2 4 17 25
4. la valeur exacte de est , l'erreur commise donc est
I
π
4
π
|I − ISc (f )| = | − 0.86039| ≤ 0.075
4
(1 pt.)

M. Moussa Page 5/5