MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

4. DÉVELOPPEMENT DE TAYLOR EN UNE VARIABLE

ET NOTION D’ORDRE (§1.6.1; 1.5.3)
4.1 Rappels

Définition : Soit f une fonction dérivable n fois en un point x0 . Le

polynôme de Taylor de dégré n de f autour de x0 est donné par :
n (k )
f ( x0 )
Pn ( x ) = ∑ ( x − x0 ) k
k =0 k!
f ′′( x 0 )
= f ( x 0 ) + f ′( x 0 )( x − x 0 ) + ( x − x0 ) 2
2!
(n)
f ( x0 )
+ ... + ( x − x0 ) n .
n!

Remarques
• f (k )
( x0 ) = Pn( k ) ( x0 ), ∀k ≤ n :
C’est-à-dire que f et Pn ont les mêmes
dérivées jusqu’à l’ordre n en x0 .
• Pn est le «meilleur» polynôme d’approximation de f en x0 .

Quelques développements de Taylor x0 = 0.

n
xk
f (x) = ex ⇒ Pn (x) = ∑
k =0 k!

x2k +1
n
f (x) = sin(x) ⇒ Pn (x) = ∑(−1) k

k =0 (2k +1)!
n
x2k
f (x) = cos(x) ⇒ Pn (x) = ∑(−1)k
k =0 (2k)!
n
x2k +1
f (x) = arctan(x) ⇒ Pn (x) = ∑(−1)k , x ∈[−1;1]
k =0 2k +1
n
x ∈]−1;1[
1
f (x) = ⇒ Pn (x) = ∑xk ,
1− x k =0

xk n n
xk +1
f (x) = ln(1+ x) ⇒ Pn (x) = ∑(−1)k −1
= ∑(−1)k , x ∈]−1;1]
k =1 k k =0 k +1
α α(α −1) 2 α(α −1)(α − 2) 3 α(α −1)(α − 2)...(α − n +1) n
f (x) = (1+ x) = 1+ x + x ∈]−1;1[
α
x + x +...+ x,
1! 2! 3! n!

Partie I : Fondements – MTH2201A (P. M.) Page 1

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

Théorème 1 :

Soit f une fonction dérivable n + 1 fois en un point x0 . Alors :

• f ( x) = Pn ( x) + Rn ( x)
avec :
f ( n +1) (ξ ( x) )
• Terme d ' erreur : Rn ( x) = ( x − x0 ) n+1 ,
(n + 1)!
• ξ ( x) : compris entre x0 et x.

Remarques

• f ( x) = Pn ( x) + Rn ( x) ⇒ Égalité
• f ( x) ≈ Pn ( x) ⇒ Approximation.
• ξ ( x ) : dépend de x.
• Le théorème garantit l’existence de ξ ( x ) mais en général il n’est
pas possible de déterminer explicitement ce point.

Définition: Erreur de troncature

L’erreur commise lors de l’approximation de f ≈ Pn est appelée

erreur de troncature.

Forme pratique du développement de Taylor autour de x0

En posant h = x − x0 , ⇒ x = x0 + h on a :

• f ( x ) = f ( x 0 + h ) = Pn ( h ) + R n ( h )
où :
f ′′( x 0 ) 2 f (n ) ( x0 ) n
• Pn ( h ) = f ( x 0 ) + f ′( x 0 ) h + h + ... + h .
2! n!
f ( n +1) (ξ ( h ) ) n +1
• Terme d ' erreur : R n ( h ) = h .
(n + 1)!
• ξ ( h ) : compris entre x 0 et x 0 + h.

Partie I : Fondements – MTH2201A (P. M.) Page 2

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

Exemple 1: Polynômes de Taylor e x : x0 = 0, n = 0,...,5

h2 h3 hn e ε ( h ) h n +1
• f ( x) ≈ Pn (h) = 1 + h + + + ... + et Rn ( h) = , ε (h) ∈ [0, h ].
2 3! n! (n + 1)!
h
h n +1
• Borne sup . Rn (h) : e ε ( h ) ≤ e h dans [0, h] ⇒ Rn (h) ≤ e .
(n + 1)!
• Estimer : e 0,1 avec h = 0,1 et P3 (0,1) = 1,1051667 ⇒ e 0,1 − P3 (0,1) = 0,420 × 10 −5 .
⇒ R3 (0,1) ≤ 0,460 × 10 −5 et P3 (0,1) = 1,1051667 possède 6chiffres significatifs.
• Estimer : e 0,05 avec h = 0,05 et P3 (0,05) = 1,051270833 ⇒ e 0, 05 − P3 (0,05) = 0,263 × 10 −6 .
e 0,1 − P3 (0,1) 0,420 × 10 −5
• Rapport des erreurs absolues : = −6
= 16,14 ≈ 2 4. Pourquoi ?
e 0,05 − P3 (0,05) 0,263 × 10

Partie I : Fondements – MTH2201A (P. M.) Page 3

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

Exemple 2:
x ln(1 + t )
Soit la fonction g ( x) = ∫0 dt.
t

a) Obtenir le développement de Taylor de la fonction g (x) autour

de x = 0.
b) Donner une borne supérieure de l’expression analytique du
terme d’erreur pour le développement de Taylor de degré n de
1
g (x) pour x ≤ .
2
c) Déterminer n pour que l’approximation donnée par ce
développement de Taylor ait une erreur qui soit au plus
 1 1
10 −7 sur − , .
 2 2

Solution

ln(1 + t ) x 1 ∞ tk  x ∞ t k −1 
dt = ∫ ∑ (−1) k −1  dt = ∫ ∑ (−1) k −1
x
a) g ( x) = ∫ dt
0 t 0 t
 k =1 k 0
 k =1 k 
t=x
∞
tk  x2 x3 x4 n −1 x
n
= ∑ (−1) k −1
 = x − + − + ... + ( − 1) + ...
k =1 k 2  t =0 2 2 32 4 2 n2

x n +1
b) Série altérnée : Rn ( x) ≤ (−1) n −1
(n + 1) 2

Pour x≤
1
⇒ Rn ( x) ≤ (−1) n −1
(1 2) = n +1
1
.
2 (n + 1) 2
(n + 1) 2 2 n+1

1
c ) Rn ( x ) ≤ 2 n +1
≤ 10 − 7.
(n + 1) 2
Solveur : n = 14.37 ⇒ R14 ( x) = 1,355 × 10 −7 et R15 ( x) = 0,961 × 10 −7 ⇒ n = 15.

Partie I : Fondements – MTH2201A (P. M.) Page 4

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

4.2 Notion d’ordre

Définition

Une fonction f (h) est un grand ordre de h n , noté f (h) = O(h n ) , au

voisinage de 0, s’il existe une constante C > 0 telle que

f ( h)
≤ C ⇒ f (h) ≤ C × h n .
hn

Caractéristique de f (h) = O(h n ) : proche de 0, f décroit à la même

vitesse que h n .

Idée du comportement de f (h) = O(h n )

n
h h 1
Si h → ⇒ f (h) ≤ C × h n devient C ×   = n × C × h n .
2 2 2

Théorème (Remarque 1.12)

L’erreur d’approximation d’une fonction f par son polynôme de

Taylor de degré n est d’ordre n + 1. Autrement dit,

f ( x0 + h) = Pn (h) + O (h n +1 )

Remarque

Pn (h) est généralement (mais pas toujours) une approximation

d’ordre n + 1 de f (x).

L’ordre est également le degré du premier terme non nul du terme

d’erreur.

Partie I : Fondements – MTH2201A (P. M.) Page 5

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

Exemple 3:
2 x
Soit la fonction f ( x) = ∫ e −t dt. Pour obtenir le développement de
2

π 0

Taylor, on peut suivre les étapes suivantes :

a) Obtenir le développement de Taylor de la fonction e x .
b) Déduire de (a) le développement de Taylor de e −t .
2

c) Déduire de (b) le développement de Taylor de f (x).

d) Donner une approximation de f (1) en utilisant les 4 premiers
termes de son développemt de Taylor.
e) Quel est l’ordre de précision de l’approximation obtenue en (d).
f) Donner le nombre de chiffres significatifs de l’approximation
obtenue en (d) en la comparant avec la valeur exacte
f (1) = 0,842701.

Solution à corriger
x 2 x3 x 4 x5 xn
a) ex = 1 + x + + + + + .. + + ...
2! 3! 4! 5! n!

t 4 t 6 t 8 t 10 t 2n
e −t = 1 − t 2 + − + − + .. + (−1) n + ....
2
b)
2! 3! 4! 5! n!

2 x
−t 2 2  x t 4 t 6 t 8 t 10 n t
2n

∫ dt = ∫  1− t + − + − + .. + (−1) + ... dt
2
c) e
0 
π 0
π  2! 3! 4! 5! n! 
2  t3 t5 t7 t9 t 2 n+1 
=.  t − + − + + ... + (−1) n + .....  t = x
π  3 × 1! 5 × 2! 7 × 3! 9 × 4! (2n + 1) × n!  t =0
t=x
2  t3 t5 t7 t9  2  x3 x5 x7 x9 
=. t − + − + + ... = 
 x − + − + + ... 
π  3 10 42 216  t =0 π  3 10 42 216 

2  1 1 1 
d ) f (1) ≈ 1 − + −  = 0,838224524.
π  3 10 42 

e) P7 ( h) et R 7 ( h) = O( h 9 ) ⇒ Approximation d ' ordre 9.

f) x − x * = 0,842701 − 0,838224524 = 0.004476476 = 0.45 × 10 − 2 ≤ 0.5 × 10 − 2
∴ 2chiffres significat ifs

Partie I : Fondements – MTH2201A (P. M.) Page 6

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

Exemple 4:

Soit la fonction f ( x) = 1 + x .

a) Obtenir le développement de Taylor d’ordre 3 de f (x) autour

de x0 = 0.
b) Donner l’expression analytique du terme d’erreur pour le
développemt de Taylor obtenu en (a).
c) À l’aide du polynôme trouvé en (a), donner des
approximations de 1,1 et 1,025 et calculer les erreurs
absolues e1 et e2 commises en comparant les valeurs exactes
correspondantes.
e1
d) Effectuer le rapport et expliquer le résultat.
e2
e) Calculer une borne supérieure de l’expression analytique du
terme d’erreur calculé en (b) et utiliser cettte expression pour
estimer l’erreur sur l’évaluation de 1,1 et 1,025 . Comparer
avec les erreurs exactes e1 et e2 .

Solution

x x2
a ) Développement de Taylor : dégré = 2 et ordre = 3 : p 2 ( x) = 1 + −
2 8
b) Puisque f ′′′( x) = (1 + x) −5 2 ⇒ R2 ( x) = (1 + ξ ( x)) −5 2 x 3 , pour ξ ( x) ∈ ]0, x[.
3 1
8 16
c) e1 = 1,1 − p 2 (0.1) ) ≤ 0,588 × 10 − 4 ≤ 0,5 × 10 −3 , avec : p 2 ( x) = 1,041875
⇒ p 2 ( x) possède 4chiffres significatifs.
De même : e2 = 1,025 − p 2 (0.025) ) ≤ 0,962 × 10 − 6 ≤ 0,5 × 10 −5 , avec : p 2 ( x) = 1,012421875
⇒ p 2 ( x) possède 6 chiffres significatifs.

e1 1,1 − p 2 (0.1) ) h
d) = = 61,17 ≈ 64 = 4 3 ⇒ appr. ordre = 3 et h→ .
e2 1,025 − p 2 (0.025) ) 4

(1 + ξ ( x)) −5 2 x 3 , pour ξ ∈ ]0, x[ ⇒ Borne

1 1
e) Puisque R2 ( x) = Sup. : e( x) = x 3 .
16 16
1 1
⇒ x = h = 0,1 :⇒ e(0,1) = (0,1) 3 = 0,625 × 10 − 4 et e(0,025) = (0,025) 3 = 0,97656 × 10 −6 .
16 16

Partie I : Fondements – MTH2201A (P. M.) Page 7

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

4.3 Propagation de l’erreur en général

Soit la quantité x approximée par x * avec erreur ∆x, c’est-à-dire que
x = x * ± ∆x. Soit f une fonction de la variable x.

L’erreur ∆f = f ( x) − f ( x * ) commise lors de l’approximation de f (x) est

estimée par

∆f ≈ f ′( x * ) ∆x. ⇔ f ( x) ≈ f ( x * ) ± f ′( x * ) ∆x

Exemple 5:

Estimer l’erreur commise dans l’évaluation de la fonction

suivante f ( x) = arctan( x) avec x * = 1,0100. Tous les chiffres fournis sont
significatifs. Indiquer le nombre de chiffres significatifs du
résultat.

Solution

En pratique ∆x ≤ 0,5 × 10 −4 , car tous les chiffres fournis sont

significatifs.

∆f ≈ f ′( x * ) ∆x =
1
(0,5 × 10 ) =
−4 1
(0,5 × 10 ) = 0,495029988× 0,5 × 10
−4 −4

( )
1+ x ∗ 2
1 + (1,0100)
2

∆f = 0,2475× 10−4 ≤ 0,5 × 10−4 et f ( x ∗ ) = arctan(1,0100) = 0,79037

Donc f ( x∗) = 0,79037 possède 4 chiffres significatifs (0,7903).

Généralisation :

Si f est une fonction de plusieurs variables, par exemple deux

variables x et y , on aura alors

∂f ( x * , y * ) ∂f ( x * , y * )
∆f = ∆x + ∆ y.
∂x ∂y

Partie I : Fondements – MTH2201A (P. M.) Page 8

 MTH2210A - Calcul scientifique pour ingénieurs

ANALYSE D’ERREURS

Exemple 6:

Le volume d’un cylindre de hauteur H = 5m et de rayon R = 1,01m est

donné par V = πR 2 H . Si H est connue à 2% près et que R possède 2
chiffres significatifs, calculer la valeur de V correspondante et
donner le nombre de chiffres significatifs.

Solution

V = πR 2 H
∂V ∂V
∆V = ∆R + ∆H = 2πRH ∆R + πR 2 ∆H
∂R ∂H
∆H 2
Er = = ⇒ ∆H = E r × H = 0,02 × 5 = 0,1 et ∆ R ≤ 0,5 × 10 −1
H 100
∆V ≤ 2 × π × 1,01 × 5 × 0,5 × 10 −1 + (1,01) 2 × π × 0,1 ≤ 0,19069 × 101
⇒ ∆V ≤ 0,19069 × 101 ≤ 0,5 × 101 ⇒ V = 16,023693 : 1 seul chiffre significatif (1).

Propagation de l’erreur pour les opérations élémentaires

Si ∆x et ∆y sont des erreurs sur x et y respectivement, alors

Opération Erreur absolue Erreur relative
x± y ∆x + ∆y ∆x + ∆y
x± y
x× y y ∆x + x ∆y ∆x ∆y
+
x y
x÷ y y ∆x + x ∆y ∆x ∆y
+
y
2 x y

Exemple 7:

Soit x = 0,12345 × 10 −4 un nombre qui possède 3 chiffres significatifs et le

nombre y = 0,67890 × 10 2 est tel ∆y ≤ 0,7. En évaluant l’expression
( x + 1) y 2 , combien de chiffres significatifs obtiendrez-vous?
Réponse: ∆f = 0,95 ×102 < 0,5 ×103 ⇒ f ( x * , y * ) = 0,460910×104 (un seul chiffire significatif «4»).

Partie I : Fondements – MTH2201A (P. M.) Page 9