A N A L Y S E S NUMERIQUE

´ 1- 2eme CP ST
ENP/ALGER 2024/2025

Solution TD 01 Analyse d’erreurs

Exercice 1 enfin, on obtient :

Si pour une approximation, l’erreur ∆x 6 0.5 × 10m 1 √ 1

b
alors le premier chiffre significatifs a u s ens étroit corres- ∆S = √ ∆a + 2 √ ∆b
a+ b a+ b
pond à la mieme puissance de dix. Si ∆x 6 10m on parle de
chiffre significatif au sens large. Dans tout ce qui suit on
adopte la première définition. mais ∆a = ∆b = 0.5 × 10−1 . L’application numérique nous
donne : ∆S = 0.5 × 10−2 ce qui donne S ' 2.52.
1 x1 = 0.1342 ⇒ ∆x1 6 0.5 × 10−4
0.5 × 10−4 Exercice 4
⇒ δ1 = ' 4 × 10−2 %
0.1342
2 x2 = 7.480 ⇒ ∆x1 6 0.5 × 10−3 1 Le développement d’ordre quatre prend la forme :
0.5 × 10−3 f ( x ) = f ( x0 + h) = f ( x0 ) + f 0 ( x0 )h + f 00 ( x0 )
h2 h3
+ f (3) ( x 0 ) + R 3 ( h )
⇒ δ2 = ' 7 × 10−3 % 2! 3!
7.480
remarquant que bien que ∆x1 < ∆x2 mais δ1 > δ2 ce qui les dérivées d’ordre n et leurs valeurs au point x0 = 0 de la
fait apparaı̂tre l’interet du concept de l’erreur relative. fonction f ( x ) = ln(2x + 1) sont regroupées dans le tableau
3 x3 = 3.1415 ⇒ ∆x3 6 0.5 × 10−4 suivant :
0.5 × 10−4 n f (n) ( x ) f ( n ) (0)
⇒ δ3 = ' 1.6 × 10−3 %
3.1415
4 x4 = 0.34235 × 10−3 ⇒ ∆x4 6 0.5 × 10−8 0 ln(2x + 1) 0
0.5 × 10−8
⇒ δ4 = ' 1.5 × 10−3 % 2
0.34235 × 10−3 1 2
2x + 1
5 x5 = 7.000 ⇒ ∆x5 6 0.5 × 10−3
−4
0.5 × 10−3 2 -4
⇒ δ5 = ' 7 × 10−3 % (2x + 1)2
7.000
16
3 16
6 x6 = 0.22754 × 109 ⇒ ∆x6 6 0.5 × 104 (2x + 1)3
0.5 × 104
⇒ δ6 = ' 2.2 × 10−3 % −96
0.22754 × 109 4
(2x + 1)4
Exercice 2 en remplaçant dans l’équation précédente, on obtient :
∆x 8
a) x = 3.14156 mais δ = 0.1% = 10−3 = ⇒ f ( x ) = f (h) = 2h − 2h2 + h3 + R3 (h)
x 3
∆x ' 0.3 × 10−2 6 0.5 × 10−2 ⇒ x ' 3.14. Si
∆x 2 La formule du reste de Taylor-Lagrange permet
δ = 1% = 10−2 = ⇒ ∆x ' 0.3 × 10−1 6
x d’écrire :
0.5 × 10−1 ⇒ x ' 3.1.
f (4) ( c ) 4
∆x ∃ c ∈ [0, h] tel que : R3 (h) = h
b) x = 392.12 mais δ = 0.01% = 10−4 = ⇒ 4!
x
−1 −1
∆x ' 0.4 × 10 6 0.5 × 10 ⇒ x ' 392.1. Si l’expression de la quatrième dérivée étant calculée dans le
∆x tableau précédent, on obtient :
δ = 1% = 10−2 = ⇒ ∆x ' 0.4 × 101 6
x
0.5 × 101 ⇒ x ' 3.9 × 102 . On utilisant la notation −96 h4 −4h4
R3 ( h ) = 4
=
scientifique on évite l’ambiguı̈té d’écrire x ' 390 car (2c + 1) 4! (2c + 1)4
le zéro n’est pas significatif1 .
3 On peut calculer l’approximation de f ( x ) au point
x = 0.1 en utilisant le polynôme de Taylor obtenu.
Exercice 3 Puisque x = x0 + h et x0 = 0 on a x = h = 0.1 donc :
 √   √ 
S = ln a + b = ln 10.3 + 4.4 = 2.51750433. 8
y∗ = f (0.1) ' 2(0.1) − 2(0.1)2 + (0.1)3 = 0.1826666667
Mais a et b étant entachés d’erreurs, ces erreurs se pro- 3
pagent dans le calcul et on doit déterminer les chiffre si-
4 Pour déterminer les chiffres significatifs de l’approxi-
gnificatif du résultat final. Pour celà, on doit évaluer ∆S.
mation précédente on peut calculer l’erreur par comparai-
On utilisant la notion de la différentielle :
son avec la valeur exacte y = f (0.1) = 0.182321557. On
∂S ∂S arrive a :
dS = da + db
∂a ∂b ∆y = |y∗ − y| ' 0.35 × 10−3 6 0.5 × 10−3
√ 1
1 b signifiant que le chiffre des millième est significatif. Donc :
dS = √ da + 2 √ db f (0.1) ' 0.182 avec trois chiffres significatifs.
a+ b a+ b
1 On note que certains manuels utilise cette notation x ' 390 on indiquant que le zéro n’est pas significatif.

Page 1/4 [Link]

 ANALYSE NUMERIQUE 1- 2eme CP ST

5 Si on suppose que la valeur exacte est inconnue, on Exercice 5

peut quand même obtenir une estimation de l’erreur en
1 Le développement d’ordre trois prend la forme :
partant de l’expression du reste R3 (h). Le problème est
qu’on trouve le nombre c dont la valeur est inconnue dans h2
f ( x ) = f ( x0 + h) = f ( x0 ) + f 0 ( x0 )h + f 00 ( x0 ) + R2 ( h )
l’expression de ce reste. Cependant, on peut se débarasser 2!
de ce nombre pour obtenir une majoration du reste. Pour les dérivées d’ordre n et leurs valeurs au point x0 = 0 de
√ 1
celà on note qu’on peut majorer la valeur absolue de l’er- la fonction f ( x ) = 3 2x + 1 = (2x + 1) 3 sont regroupées
reur plutôt que l’erreur elle même puisque le signe de cette dans le tableau suivant :
dernière ne nous interesse pas : n f (n) ( x ) f ( n ) (0)

−4h4 4h4 1
R3 ( h ) = ⇒ | R 3 ( h )| = 0 (2x + 1) 3 1
(2c + 1)4 (2c + 1)4
2 2
1 (2x + 1)− 3 2
3
en effet : 3
8 5

c ∈ [0, h] ⇒ 06c6h 2 − (2x + 1)− 3 − 89

9
⇒ 1 6 (2c + 1)4 6 (2h + 1)4
80 8
1 1 3 (2x + 1)− 3
⇒ 1> > 27
(2c + 1)4 (2h + 1)4
4h4 4h4 en remplaçant, on obtient :
⇒ 4h4 > 4
>
(2c + 1) (2h + 1)4 2 4
f ( x ) = f ( h ) = 1 + h − h2 + R2 ( h )
3 9
on reconnait l’expression de la valeur absolue du reste au
milieu de la dernière inégalité. Donc : 2 La forme analytique du reste :

f (3) ( c ) 3
4 4h4 ∃ c ∈ [0, h] tel que : R2 (h) = h
4h > | R3 (h)| > 3!
(2h + 1)4
l’expression de la troisième dérivée étant calculée dans le
c’est la partie gauche de l’inégalité qui nous interesse. On tableau précédent, on obtient :
déduit que pour x = h = 0.1 on a : h380 40 h3
R2 ( h ) = = 8 8
| R3 (0.1)| 6 4(0.1)4 = 0.4 × 10−3 6 0.5 × 10−3 27(2c + 1)− 3
3! 81(2c + 1)− 3
√
3 Pour calculer l’approximation de 3 1.4 on remarque
on retrouve le résultat précédent f (0.1) ' 0.182 avec trois √
chiffres significatifs. On note que l’erreur réelle, calculée que 3 1.4 = f (0.2) et puisque x = x0 + h = h on déduit
précédement 0.35 × 10−3 , est effectivement inférieur à la que h = 0.2. Donc :
majoration obtenue 0.4 × 10−3 . √ 2 4
1.4 = f (0.2) ' 1 + 0.2 − 0.22 = 1.115555556
3

6 L’approximation pour x = h = −0.1 : 3 9

4 On peut estimer l’erreur en majorant le reste :
8
f (−0.1) ' 2(−0.1) − 2(−0.1)2 + (−0.1)3 = 0.2226666667
3 c ∈ [0, h] ⇒ 06c6h
8 8
on peut estimer l’erreur commise dans cette approxima- ⇒ 1 6 (2c + 1) 3 6 (2h + 1) 3
tion par le même procédé de majoration précédent, mais 1 1
⇒ 1> 8 > 8
en notant que h < 0 donc : (2c + 1) 3 (2h + 1) 3
40h3 40h3 40h3
⇒ > 8 > 8
c ∈ [h, 0] ⇒ 0>c>h 81 81(2c + 1) 3 81(2c + 1) 3
⇒ 1 > (2c + 1)4 > (2h + 1)4 on reconnait l’expression de la valeur absolue du reste au
1 1 milieu de la dernière inégalité. Donc :
⇒ 16 6
(2c + 1)4 (2h + 1)4
40h3 40h3
4h4 4h4 > | R2 (h)| >
4h4
8
⇒ 6 6 81 81(2c + 1) 3
(2c + 1) 4 (2h + 1)4
On déduit que pour x = h = 0.2 on a :
Donc :
40(0.2)3
4 4h4 | R2 (0.2)| 6 ' 0.4 × 10−2 6 0.5 × 10−2
4h 6 | R3 (h)| 6 81
(2h + 1)4 √
donc 3 1.4 = f (0.2) ' 1.12 avec trois chiffres significatifs.
à présent, c’est la partie droite de l’inégalité qui nous in-
teresse. On déduit que pour x = h = −0.1 on a : Exercice 6
4(−0.1)4 1 Le développement d’ordre trois prend la forme :
| R3 (−0.1)| 6 = 0.98 × 10−3 6 0.5 × 10−2
(2(−0.1) + 1)4 h2
f ( x ) = f ( x0 + h) = f ( x0 ) + f 0 ( x0 )h + f 00 ( x0 ) + R2 ( h )
2!
enfin, f (−0.1) ' 0.22 avec seulement deux chiffres signi- les dérivées d’ordre n et leurs valeurs au point x0 = 0 de
ficatifs. la fonction f ( x ) = sin x :

Page 2/4 [Link]

ENP/ALGER

n f (n) ( x ) f ( n ) (0) Exercice 7

0 sin x 0 Prenant le développement d’ordre deux :

1 cos x 1 y = f ( x ) = f ( h ) = 1 + R1 ( h )
2 − sin x 0
notant que c’est bien un polynôme de degré 1 dont le
coéfficient de h est nul. On obtient un reste :
3 − cos x
f (2) ( c ) 2 h2 h2
en remplaçant, on obtient : R1 ( h ) = h = −(cos c) ⇒ | R1 (h)| 6
2! 2! 2!
f ( x ) = f ( h ) = h + R2 ( h )
avec une erreur relative pour x = 0.3 :
on remarque que c’est bien un développement d’ordre
trois mais le coéfficient de h2 est nul. h 2
2 La forme analytique du reste : ∆y 0.32
δ= = 2! = ' 4.5 % > 0.1 %
|y| 1 2!
f (3) ( c ) 3 h3
∃ c ∈ [0, h] : R2 (h) = h = −(cos c)
3! 3! on essaye le développement suivant qui est d’ordre
quatre :
3 puisque x = x0 + h = h on déduit que h = 0.1 :
f (0.1) = sin 0.1 ' 0.1 h2
y = f ( x ) = f (h) = 1 − + R3 ( h )
on note que h doit être en radian. 2!
4 Pour estimer le nombre de chiffres significatifs de cette
on note aussi que c’est un polynôme de degré 3 dont le
approximation on peut majorer le reste :
coéfficient de h3 est nul. On obtient un reste :
c ∈ [0, h] ⇒ 06c6h
⇒ 1 > cos c > cos h f (4) ( c ) 4 h4 h4
R3 ( h ) = h = (cos c) ⇒ | R3 (h)| 6
h 3 h3 h3 4! 4! 4!
⇒ > cos c > cos h
3! 3! 3!
avec une erreur relative :
h3
⇒ > | R2 (h)| > cos h
3! ∆y h 4
0.3 4

δ= = 4! h2 = 4!
2 ' 0.04% < 0.1 %
On déduit que pour x = h = 0.1 on a : |y| 1 − 2! 1 − 0.3
2!

0.13
| R2 (0.1)| 6 ' 0.17 × 10−3 6 0.5 × 10−3 Donc le développement d’ordre quatre cos x = 1 − x2!
2
6
donne une précision inférieure à 0.1 % pour x = 0.3 rad.
donc y = f (0.1) = sin 0.1 ' 0.100 avec trois chiffres si-
∆y On remarque qu’on peut même aller encore au delà de
×10−3 =
gnificatifs et une erreur relative δ = |y| = 0.170.1 x = 0.3 tout en gardant cette précision.
0.17 × 10−2 = 0.17 %
4 Si on calcul l’erreur pour h = 0.2 on obtient : Exercice 8
0.23 2 2
| R2 (0.2)| 6 ' 0.13 × 10−2 6 0.5 × 10−3 En utilisant l’approximation cos x ' 1 − x2 = 1 − h2
6 on commet une erreur calculée dans l’exercice précédent :
avec une erreur relative :
∆y 0.13 × 10−2 f (4) ( c ) 4 h4 h4 x4
δ= = = 0.65 × 10−2 = 0.65 % R3 ( h ) = h = (cos c) ⇒ | R3 (h)| 6 =
|y| 0.2 4! 4! 4! 4!
on voit que l’erreur et l’erreur relative augmentent à
pour une précision au centième prés :
mesure que h augmente ( on s’éloigne du centre du
développement x0 ). Si on souhaite que l’erreur relative
reste inférieur à 1% on doit avoir : x4 1
6 0.5 × 10−2 ⇒ x 6 (4!(0.005)) 4 ' 0.59 rad
∆y ∆y 4!
δ= 6 1% = 10−2 ⇒ 6 10−2
|y| y
pour une précision au millième prés :
h3 x3
mais ∆y = = parce que x = h et y = sin x ' x.
3! 3! x4 1
Donc : 6 0.5 × 10−3 ⇒ x 6 (4!(0.0005)) 4 ' 0.33 rad
4!
x3
∆y √ pour une précision au dix-millième prés :
= 3! 6 10−2 ⇒ x 6 0.06 ' 0.25 rad ' 14˚
|y| x
on peut donc utiliser l’approximation sin x ' x jusqu’à x4 1
6 0.5 × 10−4 ⇒ x 6 (4!(0.00005)) 4 ' 0.19 rad
l’angle x ' 14˚ en étant certain que l’erreur relative ne 4!
dépassera pas 1%.

Page 3/4 [Link]

ENP/ALGER ANALYSE NUM ERIQUE 1- 2eme CP ST

Exercice 9 f (3) ( c ) 3 x3 x3
R2 ( x ) = x = sin c ⇒ | R2 ( x )| 6
3! 3! 3!
Z 1/2
sin x
I= dx f (4) ( c ) 4 x4 x4
0 x R3 ( x ) = x = cos c ⇒ | R3 ( x )| 6
4! 4! 4!
on peut utiliser une développement de Taylor de f ( x ) =
sin x autour de x0 = 0. Cependant, on doit déterminer f (5) ( c ) 5 x5 x5
R4 ( x ) = x = − sin c ⇒ | R4 ( x )| 6
le degré du polynôme necessaire pour atteindre la 5! 5! 5!
précision demandée de 0.5 × 10−4 . Lorsqu’on utilise un
f (6) ( c ) 6 x6 x6
développement de degré n tel que sin x = Pn ( x ) + Rn ( x ) R5 ( x ) = x = − sin c ⇒ | R5 ( x )| 6
(on remarque que h = x) on obtient une approximation de 6! 6! 6!
l’intégrale In : f (7) ( c ) 7 x7 x7
R6 ( x ) = x = − sin c ⇒ | R6 ( x )| 6
Z 1/2
Pn ( x )
Z 1/2
Rn ( x ) 7! 7! 7!
I= dx + dx
x x on déduit les erreurs correspondantes ∆In :
|0 {z } |0 {z }
approximation In l’erreur ∆In Z 1 2
x 1 3
x /3 dx = = 0.15 > 0.5 × 10−4
1
donc : ∆I1 6
Z 1/2 0 2! 2! 10
Rn ( x )
∆In 6 dx
0 x Z 1 3
x 1 3
x /3 dx = = 0.04 > 0.5 × 10−4
1
on calcul cette erreur pour n = 1, 2, 3, . . . jusqu’à obtenir ∆I2 6
0 3! 3! 13
∆In 6 0.5 × 10−4 .
Z 1 4
x 1 3
x /3 dx = =' 0.8 × 10−2 > 0.5 × 10−4
1
∆I3 6
f (2) ( c ) x2 x2 0 4! 4! 16
R1 ( x ) = x2 = − sin c ⇒ | R1 ( x )| 6 Z 1 5
2! 2! 2! x 1 3
x /3 dx = =' 0.13 × 10−2 > 0.5 × 10−4
1
∆I4 6
f (3) ( c ) x3 x3 0 5! 5! 19
R2 ( x ) = x3 = − cos c ⇒ | R2 ( x )| 6
3! 3! 3! Z 1 6
x 1 3
x /3 dx = =' 0.2 × 10−3 > 0.5 × 10−4
1
f (4) ( c ) 4 x4 x4 ∆I5 6
R3 ( x ) = x = sin c ⇒ | R3 ( x )| 6 0 6! 6! 22
4! 4! 4! Z 1 7
x 1 3
f (5) ( c ) 5 x /3 dx = =' 0.2 × 10−4 6 0.5 × 10−4
1
x5 x5 ∆I6 6
R4 ( x ) = x = cos c ⇒ | R4 ( x )| 6 0 7! 7! 25
5! 5! 5!
on déduit les erreurs correspondantes ∆In : on doit donc utiliser le polynôme de degré six :
Z 1/2 x2
1 1 1 1 x2 x4 x6
∆I1 6 2!
dx = = ' 0.6 × 10−1 > 0.5 × 10−4 cos x ' P6 ( x ) = 1 − + −
0 x 2! 2 22 16 2! 4! 6!
Z 1/2 x3
1 1 1 1 on arrive à l’approximation :
∆I2 6 3!
dx = 3
= ' 0.7 × 10−2 > 0.5 × 10−4
0 x 3! 3 2 144 !
x2 x4 x6
Z 1/2 x4 Z 1√ Z 1
1 1 1 −3 −4 I ' I = 3
x P ( x ) dx = x /3
1
1− + − dx
∆I3 6 4!
dx = 4
' 0.7 × 10 > 0.5 × 10 6
0
6
0 2! 4! 6!
0 x 4! 4 2
Z 1/2 x5
1 1 1
∆I4 6 dx =5!
' 0.5 × 10−4 6 0.5 × 10−4
0 x 5! 5 25 !
x /3 x /3 x /3
Z 1 7 13 19
1/3
on doit donc utiliser le polynôme de degré quatre : ⇒ I6 = x − + − dx
0 2! 4! 6!
h3 x3 #1
sin x ' P4 ( x ) = h − = x−
"
3 x /3 3 x /3 3 x /3
10 16 22
3! 3! 3 4/3
⇒ I6 = x − + −
4 10 2! 16 4! 22 6!
on arrive à l’approximation : 0
3 3 1 3 1 3 1
Z 1/2
P4 ( x )
Z 1/2
x − x3!
3
⇒ I6 = − + − = 0.6076
I ' I4 = dx = dx 4 10 2! 16 4! 22 6!
0 x 0 x
avec quatre chiffres significatifs. On remarque, cependant,
x
 3 1/2
que le polynôme de degré 6 est le même que le polynôme
I4 = x − = 0.4931
18 0 de degré 7 dont lequel le coéfficient de h7 est nul. Celà
avec quatre chiffres significatifs. nous permet d’obtenir une meilleure estimation de l’er-
reur de l’approximation précédente :
Exercice 10 Z 1 8
x 1 3
x /3 dx = =' 0.3 × 10−5 6 0.5 × 10−5
1
Comme dans l’exercice précédent, on doit d’abord ∆I7 6
0 8! 8! 28
déterminer le polynôme à utiliser pour approcher f ( x ) =
cos x au voisinage du zéro : ce qui permet de rajouter un autre chiffre significatif :
f (2) ( c ) x2 x2
R1 ( x ) = x2 = − cos c ⇒ | R1 ( x )| 6 I ' 0.60762
2! 2! 2!

Page 4/4 [Link]