USTHB 2020-2021 Semestre 1 Analyse numérique 1

Faculté de Mathématiques 2ème année Lic Maths

Examen final - 30 mars 2021. Durée : 1 heure

Nom et Prénom :...........................................................................................................................................................

Matricule :............................................................................................................................................................................. 20
Exercice 1 (7 pts.) : Soient α ∈ R, h > 0 et f ∈ C 3 ([α − 1, α + 1] , R).

a) Calculer lim f [α − h, α] et montrer qu’on a lim f [α − h, α, α + h] = 12 f 00 (α) .

h→0 h→0

b) Soit p2 le polynôme d’interpolation associé à f aux nœuds α − h, α, α + h écrit sous forme de Newton.

Déterminer alors q (x) où q (x) = lim p2 (x).

h→0

c) Vérifier que q est le seul polynôme, de degré inférieur ou égal à 2, tel que

q (α) = f (α) , q 0 (α) = f 0 (α) et q 00 (α) = f 00 (α) .

Réponse.

a) La différence divisée f [α − h, α] est définie par

f (α) − f (α − h) f (α) − f (α − h)
f [α − h, α] = = .
α − (α − h) h

Alors
f (α) − f (α − h)
lim f [α − h, α] = lim = f 0 (α) .
h→0 h→0 h
La différence divisée f [α − h, α, α + h] est définie par
f (α+h)−f (α)
f [α, α + h] − f [α − h, α] h
− f (α)−fh(α−h)
f [α − h, α, α + h] = =
α + h − (α − h) 2h
f (α + h) − 2f (α) + f (α − h)
= .
2h2

Effectuons un développement de Taylor d’ordre 2 de f autour de α,

h2 00 h3
f (α + h) = f (α) + hf 0 (α) + f (α) + f 000 (ξ1 ) , ξ1 ∈ ]α, α + h[ ,
2 6
2 3
h h
f (α − h) = f (α) − hf 0 (α) + f 00 (α) − f 000 (ξ2 ) , ξ2 ∈ ]α − h, α[ .
2 6

Alors f [α − h, α, α + h] s’écrit comme

1 h 000
f [α − h, α, α + h] = f 00 (α) + (f (ξ1 ) − f 000 (ξ2 )) .
2 12

D’où
1
lim f [α − h, α, α + h] = f 00 (α) .
h→0 2
1/6
b) Le polynôme d’interpolation p2 associé à f aux nœuds α − h, α, α + h sous forme de Newton est

donné par

p2 (x) = f (α − h) + f [α − h, α] (x − (α − h)) + f [α − h, α, α + h] (x − (α − h)) (x − α) .

Alors
1
q (x) = lim p2 (x) = f (α) + f 0 (α) (x − α) + f 00 (α) (x − α)2 ,
h→0 2
qui est un polynôme de degré 2.

c) On a

q 0 (x) = f 0 (α) + f 00 (α) (x − α) et q 00 (x) = f 00 (α) .

On remarque le polynôme q vérifie

q (α) = f (α) , q 0 (α) = f 0 (α) et q 00 (α) = f 00 (α) .

Le polynôme q est unique car il s’agit du polynôme de Taylor du second ordre de f autour de α.

2/6
 R1 2
Exercice 2 (7 pts.) : On recherche une approximation de l’intégrale I = 0
e−x dx, dont on ne connaı̂t

pas la primitive.

a) Donner une approximation de cette intégrale en appliquant la méthode du trapèze sur 4 sous-

intervalles (n = 4) et en utilisant 9 chiffres significatifs avec arrondi.

b) Indiquer le terme d’erreur que l’on fait par cette méthode de trapèze.

c) Déterminer alors le nombre d’intervalles minimal nécessaire pour déterminer l’approximation de

l’intégrale I par la méthode du trapèze avec une erreur de 10−4 .

d4 −x2
d) Refaire les mêmes questions pour la méthode de Simpson. (Indication sup dx4
e = 12).
x∈[0,1]
Réponse.
b−a 2
On a a = 0, b = 1, n = 4, h = = 14 , xi = 4i , i = 0, ..., 4 et f (x) = e−x .
n
On obtient alors le tableau de données suivantes

1 1 3
xi 0 4 2 4
1
2 1 1 9
yi = f (xi ) = e−xi 1 e− 16 e− 4 e− 16 e−1

a) L’approximation de l’intégrale I par la méthode des trapèzes est

n−1
!
h X
Tn (f ) = f (a) + 2 f (xi ) + f (b) .
2 i=1

Pour n = 4, on obtient
 3
 1 
h  1 1 9
 
f (a) + 2 f (xi ) + f (b) = 4 1 + 2 e− 16 + e− 4 + e− 16 + e−1 ' 0.742 984 098 .
P
T4 (f ) =
2 i=1 2

b) L’erreur théorique commise par la méthode des trapèzes est donnée par

b−a 2 00 h2 2
sup 2 2x2 − 1 e−x .


|ET (h)| ≤ h sup |f (x)| =
12 x∈[a,b] 12 x∈[0,1]

000 d  −x2

2
On a f (x) = 2
2 (2x − 1) e = −4x (2x2 − 3) e−x = 0 sur [0, 1] implique x = 0. Alors
dx
 
00 2

−x2 00 00 2
sup |f (x)| = sup 2 2x − 1 e = max (|f (0)| , |f (1)|) = max |−2| , = 2.
x∈[0,1] x∈[0,1] e

D’où
h2
|ET (h)| ≤ .
6
c) Pour avoir une erreur inférieure à ε > 0, il suffit donc que

h2
≤ε
6

ce qui est équivalent à

1 1
n= ≥√ .
h 6ε
3/6
 Pour ε = 10−4 , on obtient
1
n≥ √ ' 40.82 .
6 × 10−4
Puisque n est entier, on choisit n = 41.

d) L’approximation de l’intégrale I par la méthode de Simpson est

!
h X X
Sn (f ) = f (a) + f (b) + 2 f (xi ) + 4 f (xi ) .
3 i pair i impair

Pour n = 4, on obtient

h

f (0) + f (1) + 2f 12 + 4f 41 + 4f 3



S4 (f ) = 4
3
1  
1 1 9
= 4 1 + e−1 + 2e− 4 + 4e− 16 + 4e− 16
3
' 0.746 855 380 .

L’erreur théorique commise par la méthode de Simpson est donnée par

b−a 4 h4 h4 h4
|ES (h)| ≤ h sup f (4) (x) = sup f (4) (x) = × 12 = .
2880 x∈[a,b] 2880 x∈[0,1] 2880 240

Pour avoir une erreur inférieure à 10−4 , il suffit que

h4
≤ 10−4
240

ce qui est équivalent à

1 1
n= ≥ √
4
' 2.54.
h 240 × 10−4
Puisque n est un entier pair, on choisit n = 4.

4/6
Exercice 3 (6 pts.) : Soit f : [0, 1] → R une fonction continue strictement décroissante telle que f (0) = 1

et f (1) = −3.

a) 1. Sachant que f (0.3131421438) = 0, déterminer la suite des premiers quatre étirés de la méthode

de la dichotomie dans l’intervalle [0, 1] pour l’approximation du zéro de f .

On pourra utiliser le tableau ci-dessous :

n an cn bn signe f (an ) signe f (cn ) signe f (bn )

2. Combien d’itérations faut-il effectuer pour approcher le zéro de f à 2−5 près ?

b) Considérons l’équation ex = x (3 + ex ).

1. Montrer que cette équation admet une unique solution réelle r dans [0, 1].

2. Écrire la méthode de Newton pour approcher la solution r.

3. Effectuer quatre itérations avec la méthode de Newton en démarrant de x0 = 0.5 et en utilisant 9

chiffres significatifs avec arrondi.

Réponse.

a)

1. On a

n an cn bn signe f (an ) signe f (cn ) signe f (bn )

0 0 0.5 1 + − −

1 0 0.25 0.5 + + −

2 0.25 0.375 0.5 + − −

3 0.25 0.3125 0.375 + + −

4 0.3125 0.34375 0.375 + − −

donc, après quatre itérations, le zéro de f est approché par 0.34375.

 1−0 
Log
2. Il faut effectuer au moins N = Log2−5 2
= 5 itérations.

5/6
b)

1. Soit f : R → R définie par f (x) = ex − x (3 + ex ). Alors f (x) = 0 si et seulement si x est solution

de l’équation donnée. La fonction f est de classe C ∞ (R) , f (0) = 1 > 0 et f (1) = −3 < 0 donc,

d’après le théorème des valeurs intermédiaire, la fonction f admet au moins une racine r ∈ [0, 1]. De

plus, f est strictement monotone sur [0, 1] car f 0 (x) = −xex − 3 < 0 pour tout x ∈ [0, 1], donc cette

racine est unique.

2. La méthode de Newton est


 x0 ∈ [0, 1] ,

f (x ) exn − xn (3 + exn ) exn
 xn+1 = xn − 0 n = xn −
 = (x 2
n − x n + 1) .
f (xn ) −xn exn − 3 xn exn + 3

3. En utilisant l’algorithme de Newton on obtient le tableau suivant

x0 x1 x2 x3 x4

0.5 0.323332727 0.313169525 0.313142144 0.313142144

6/6