U.S.T.H.B.

2018=2019
2éme année- Licence Genie Civil - Section C - Méthodes Numériques - Epreuve Finale
Nom & Prénoms : ..................................Matricule: ......................Gr:..
——————————————————————————————————
Exercice 1 On se propose de calculer des approximations des racines de la
fonction:

f (x) = ex x 2
1 ) Montrer que l’équation f (x) = 0 admet une unique racine r dans l’intervalle
[1 ; 2]
2 ) Montrer qu’on a f (x) = 0 , '(x) = x; avec '(x) = log(x + 2)
3 ) L’algorithme du point …xe xn+1 = '(xn ) est-il convergent vers r quelque
soit x0 2 [1 ; 2] : Justi…er votre réponse.
4 ) Trouver le nombre d’itération n nécessaire pour approcher r avec une erreur
inférieure à 10 6 en utilisant les itérations du point …xe précédent en partant
de x0 = 1:
——————————————————————————————————

1 ) f (x) = ex x 2 ) f 0 (x) = ex 1 > 0 sur l’intervalle [1; 2] ; f (1) =

e 3 = 0:281 72; f (2) = e2 4 = 3:389 1

f (x) continue, strictement croissante sur l’intervalle [1; 2] ; et f (1) f (2) <
0 ) 9! une racine unique dans l’intervalle [1; 2]- d’aprés le théorème des valeurs
intermédiaires.
2 f (x) = 0 , ex x 2 = 0 , ex = x + 2 , log(ex ) = log(x + 2) , x =

log(x + 2) = '(x)
1 1
3 ) '0 (x) = ) max j'0 (x)j = < 1 ; et '([1; 2]) = [log(3); log(4)] =
x+2 1 x 2 3
[1:098 ; 1:386 ] [1; 2] )L’algorithme du point …xe xn+1 = '(xn ) est conver-
gent vers r quelque soit x0 2 [1 ; 2] et ceci d’aprés le théorème du point …xe -
voir par exemple le cours de T. Takahashi de l’Ecole des Mines de Nancy page
22-24 dont une copie est dans ma page perso:usthb:dz=~rbellout:
4 ) Le nombre d’itération n nécessaire pour approcher r avec une erreur in-
férieure à 10 6 s’obtient d’aprés l’estimation d’erreur:

jxn xexacte j Kn jx0 xexacte j

1 1
ou K = max j'0 (x)j = < 1 est la constante de contraction) ( )n
1 x 2 3 3
6 log(10)
10 ,n 6 = 12:575 ) n 13:
log(1=3)
———————————————————————————— R 1— —2 — — — —
Exercice 2 1 ) Calculer des valeurs approchées de l’intégrale 0 ex dx à l’aide
de :

1
a) la méthode du rectangle b) la méthode du trapéze c) la méthode de Simpson.
2) Quel est le nombre de points équidistants N qu’il faut utiliser dans la quadra-
ture des trapézes pour obtenir une approximation de cette intégrale 10 3 ?
——————————————————————————————————
f (a) = f (0) = e0 = 1:0 ; f (b) = f (1) = e1 = 2:718 ; f ( a+b 2 ) = e
1=4
=
R 1 x2
1:284 ; valeur exacte = 0 e dx = 1:462 7
a) R = (b a) f (a) = 1:0; ; erreur = 1:462 7 1:0 = 0:462 7
1 1
b) T = (b a) [f (a) + f (b)] = (1:0 + 2:718) = 1:859 ; erreur = 1:462 7
2 2
1:859 = 0:396 3
1 1
c) S = (b a) f (a) + 4f ( a+b 2 ) + f (b) = 6 (1:0 + 4 1:284 + 2:718) =
6
1:475 7 ; erreur = 1:462 7 1:475 7 = 0:013
2)
f "( )h2 (b a)
ETN =
12
2 2 2 2
f (x) = ex ) f 0 (x) = 2xex ) f "(x) = 2ex 2x2 + 1 ) f 000 (x) = 4ex x 2x2 + 3
h2 (b a) 1
M2 = max jf "(x)j = 6e ; ) M2 10 3 ) N 2 = 2
0 x 1
r 12 h
1 3 2 6e 3 e
M2 10 ) N 10 ) N 103 = 36:867 ) N 37:
12 12 2
——————————————————————————————————
1
Exercice 3 1 ) Interpoler la fonction f (x) = sur l’intervalle [ 1 ; 1] par
1 + x2
un polynôme P1 (x) aux points x0 = 1 ; x1 = 1:Donner une majoration de
l’erreur E1 (x) = jf (x) P1 (x)j sur [ 1 ; 1]
2 ) Faire de même pour un polynôme P2 (x) aux points x0 = 1 ; x1 = 1 ; x2 =
0 et une majoration de l’erreur E2 (x) = jf (x) P2 (x)j sur [ 1 ; 1] :
——————————————————————————————————
1 1
f (x0 ) = f ( 1) = ; f (x1 ) = f (1) = ;
f (x2 ) = f (0) = 1:0
2 2
1 1
(f (x1 ) f (x0 ) 1 (2 2) 1
1) P1 (x) = f (x0 )+ (x x0 ) = + (x+1) ) P1 (x) =
(x1 x0 ) 2 (1 + 1) 2
2) P2 (x) = P1 (x)+f [x0 ; x1 ; x2 ] (x x0 )(x x1 ) = P1 (x)+f [ 1; 1; 0] (x+1)(x 1)
1 2 1 1 2 1 2
P2 (x) = P1 (x) (x 1) = (x 1) = 1 x
2 2 2 2
1 2
P2 (x) = 1 x
2
1 x 2
f (x) = ) f 0 (x) = 2 00
2 ) f (x) = 3 3x
2
1 )
1 + x2 2
(x + 1) 2
(x + 1)
M2 =max jf "(x)j = 2
1 x 1

000 x2 1
f (x) = 24x 4 ) M3 = max jf 000 (x)j = 4:668
(x2 + 1) 1 x 1

2
Erreur d’interpolation :
i=n
f (n+1) ( ) Y
En (x) = (x xi )
(n + 1)! i=0

max jf "(x)j 2
E1 (x) = jf (x) P1 (x)j max j(x x0 )(x x1 )j = 2 max(x2
(2)! 1 x 1
1) ) E1 (x) = 1
max jf 000 (x)j
E2 (x) = jf (x) P2 (x)j max j(x x0 )(x x1 )(x x2 )j =
(3)!
p 1 x 1
4:66
6 max x(x2 1) = 4:66
6
2
9 3 ) E2 (x) = 0:298 94
1
f (x) =
1 + x2
1
P1 (x) =
2
1 2
P2 (x) = 1 x
2

1.0
y

0.9

0.8

0.7

0.6

-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0
x