UNIVERSITE MOHAMMED PREMIER

ECOLE NATIONALE DES SCIENCES APPLIQUEES

OUJDA-MAROC
Filière : GIndus3 Année universitaire : 2020-2021
Module : Méthodes Numériques Enseignant : M. Derouich

Fiche TD N ◦ : 2
Exercice 1 (Interpolation de Newton). :
On interpole f ( x) = ln( x) par un polynôme aux points d’interpolation

x0 = 1, x1 = 2, x2 = 3, x3 = 4, x4 = 5.

1. Trouver une expression algébrique de ce polynôme en utilisant la méthode de

Newton.
2. Estimer la valeur de f (6, 32) avec le polynôme trouvé en 1.
3. Calculer l’erreur absolue.
4. Combien de points d’interpolation à intervalle régulier de 0,5 faudrait-il ajouter,
en partant de x5 = 5, 5, afin que l’erreur absolue de l’estimé de f (6, 32) obtenu en
2. diminue d’un facteur 100.
5. Sur l’intervalle [3, 4], le graphe du polynôme trouvé en 1. est-il au-dessus de celui
de f ( x), en dessous, ou se croisent-ils ?

Exercice 2 (Interpolation de Lagrange).

Soit les trois points (0, 1), (π /16, cos(π /16) et (π /8, cos(π /8)) de la fonction
f ( x) = cos( x).
1. Obtenir à l’aide de l’interpolation de Lagrange, le polynôme de degré 2 qui passe
par les 3 points et en déduire une approximation de cos(π /32).
2. Calculer le développement de Taylor de degré 2 de la fonction f(x)= cos(x) autour
de x0 = 0 et en déduire une approximation de cos(π /32).
3. Sachant que f 0 (0) = 0 , calculer le polynôme de degré 2, passant par les points
(0, 1) et (π /8, cos(π /8) et dont la dérivée en x = 0 est égale à 0 et en déduire une
approximation de cos(π /32).
4. Des trois approximations cos(π /32). que vous avez obtenues, qu’elle est la plus
précise ? Pourquoi ?

Exercice 3. :
En relevant toutes les 10 secondes la vitesse d’écoulement de l’eau dans une conduite
cylindrique, on a obtenu

t 0 10 20 30
v 2 1.89 1.72 1.44
1. Trouver une approximation de la vitesse en t = 15 via un polynôme interpolant
de degré 2 ;
2. Répéter l’opération avec un polynôme de degré 3.

1
 Corrigé

Exercice 1. :
On interpole f ( x) = ln x par un polynôme, aux points d’interpolation 1, 2, 3, 4, 5 .
1. Il y a 5 points d’interpolation, donc le degré du polynôme est 4. Le polynôme de
Newton est donné par :

Pn ( x) = f ( x0 ) + f [ x0 , x1 ] N1 ( x) + .... + f [ x0 , ..., xn ] Nn ( x)

Avec N0 ( x) = 1, N1 ( x) = ( x − x0 ), · · · , Nn ( x) = ( x − x0 ) · · · ( x − xn−1 )
Donc on construit la table des différences divisées comme suit :

i xi f [ xi ] f [ xi −1 , xi ] f [ xi −2 , xi −1 , xi ] f [ xi −3 , xi −2 , xi −1 , xi ] f [ xi −4 , xi −3 , · · · , xi ]
0 1 0
1 2 0,693 147 1806 0,693 147 1806
2 3 1,098 612 289 0,405 465 1084 -0,143 841 0361
3 4 1,386 294 361 0,287 682 072 -0,058 891 5182 0,028 316 505 97
4 5 1,609 437 912 0,223 143 551 -0,032 269 2605 0,008 874 085 90 -0,004 860 605 018
Le polynôme de Newton de degré 4 est donc :
P4 ( x) = D0 + D1 N1 ( x) + D2 N2 ( x) + D3 N3 ( x) + D4 N4 ( x)
D0 = f [ x0 ] = 0,
D1 = f [ x0 , x1 ] = 0, 6931471806,
D2 = f [ x0 , x1 , x2 ] = −0, 1438410361,
D3 = f [ x0 , x1 , x2 , x3 ] = 0, 02831650597,
D4 = f [ x0 , x1 , x2 , x3 , x4 ] = −0, 004860605018.
=⇒ P4 ( x) = 0, 6931471806( x − 1) − 0, 1438410361( x − 1)( x − 2)
+0, 02831650597( x − 1)( x − 2)( x − 3)
−0, 004860605018( x − 1)( x − 2)( x − 3)( x − 4)
P4 ( x) = −0.004860605018x4 + 0.07692255615x3 − .4838612475x2
+1.679182105x − 1.267382809
2. Pour l’estimation, il suffit d’évaluer le polynôme P4 en x = 6, 32.
On obtient alors P4 (6, 32) = 1, 68190203 . Or f (6, 32) = ln(6, 32) ' 1, 843719208.
L’erreur absolue est donc E = |1, 68190203 − 1, 843719208| ' 0.161817178.
3. On veut maintenant diminuer cette erreur d’un facteur 100 et donc obtenir une
erreur absolue de 0, 001618 .

2
 Le polynôme de Newton de degré 5 obtenu en ajoutant le nœud x5 = 5, 5 est :

P5 ( x) = P4 ( x) + f [ x0 , · · · , x5 ]( x − 1)( x − 2) · · · ( x − 5)
= P4 ( x) + 0, 78558 × 103 ( x − 1)( x − 2)( x − 3)( x − 4)( x − 5)

de sorte que P5 (6, 32) = 1, 681902 + 0, 183563 = 1, 865465.

L’erreur absolue est alors E = |ln(6, 32) − 1, 865465| = 0, 021746
et il faut encore ajouter un nœud ( x6 = 6, 0) pour obtenir :

P6 ( x) = P5 ( x) + f [ x0 , · · · , x6 ]( x − 1)( x − 2) · · · ( x − 5)( x − 5, 5)
= P5 ( x) − 0, 11905 × 10−3 ( x − 1)( x − 2) · · · ( x − 5, 5)

On a alors P6 (6, 32) = 1, 865465 − 0, 2281 = 1, 842654 et l’erreur absolue est don-
née par E = |1, 842654 − ln(6, 32)| = 0, 001065,
ce qui est mieux que la précision requise.
4. On sait que l’erreur exacte peut s’écrire
1 (5)
E= f (ξ ( x))( x − 1)( x − 2)( x − 3)( x − 4)( x − 5); ξ ∈ [1; 5]
5!
et nous sommes intéressés au signe de l’erreur.
(−1)n−1 (n − 1)!
 
( 5 ) 24 (n)
Or la fonction f (t) = 5 car ln ( x) =
t xn
et donc le signe de l’erreur ne dépend que de celui de

( x − 1)( x − 2)( x − 3)( x − 4)( x − 5)

Si x ∈ [3, 4], alors trois des facteurs sont positifs et deux sont négatifs. Par consé-
quent, l’erreur est positive et le graphe de ln( x) est au dessus de celui de P4 ( x).

Exercice 2. :

1. Rappelons que le polynôme de Lagrange basé sur les points d’interpolation d’abs-
cisses x0 , x1 , ...., xn est de degré n et s’écrit :
Pn ( x) = L0 ( x) f ( x0 ) + L1 ( x) f ( x1 ) + .... + Ln ( x) f ( xn )
j=n
(x − x j )
avec Li ( x) = ∏
j=0 ( xi − x j )
j 6 =i
On a les points d’interpolation donnés par :
x0 = 0, x1 = π /16 et x2 = π /8 ,
f ( x0 ) = 1 , f ( x1 ) = cos(π /16) , f ( x2 ) = cos(π /8) , détermineront donc un
polynôme de Lagrange de degré 2, celui-ci s’écrit :

P2 ( x) = L0 ( x) f ( x0 ) + L1 ( x) f ( x1 ) + L2 ( x) f ( x2 )

3
 ( x − x1 )( x − x2 ) ( x − π /16)( x − π /8) ( x − π /16)( x − π /8)
Avec L0 ( x) = = =
( x0 − x1 )( x0 − x2 ) (−π /16) × (−π /8) π 2 /128
( x − x0 )( x − x2 ) x( x − π /8) − x( x − π /8)
L1 ( x) = = =
( x1 − x0 )( x1 − x2 ) π /16(π /16 − π /8) π 2 /256
( x − x0 )( x − x1 ) x( x − π /16) x( x − π /16)
L2 ( x) = = =
( x2 − x0 )( x2 − x1 ) π /8(π /8 − π /16) π 2 /128
finalement

P2 ( x) = L0 ( x) f ( x0 ) + L1 ( x) f ( x1 ) + L2 ( x) f ( x2 )
( x − π /16)( x − π /8) − x( x − π /8) x( x − π /16)
= 2
+ 2
cos(π /16) + cos(π /8)
π /128 π /256 π 2 /128
( x − π /16)( x − π /8) − x( x − π /8)
 
= 128 +2 cos(π /16)
π2 π2
x( x − π /16)
+128 cos(π /8)
π2
cos(π /32) ' P2 (π /32)
(π /32 − π /16)(π /32 − π /8) −π /32(π /32 − π /8)
 
' 128 +2 cos(π /16)
π2 π2
π /32(π /32 − π /16)
+128 2
cos(π /8)
π
128π 2 (1 − 2)(1 − 4) −(1 − 4) (1 − 2)
 
' + 2 cos ( π / 16 ) + cos ( π / 8 )
322 π2 π2 π2
1
' (3 + 6 cos(π /16) − cos(π /8))
8
3 3 1
' + cos(π /16) − cos(π /8)
8 4 8
' 0.995104018738512
1
2. On a cos( x) = 1 − x2 + o( x2 ) ce qui donne
2
1
cos(π /32) ' 1 − (π /32)2 ' 0.995180857226031
2
2
3. on a P2 ( x) = αx + βx + γ donc 

 P2 (0) = 1 
 γ = 1
  π 2 π π
P2 (π /8) = cos(π /8) ⇐⇒ α ( ) + β + γ = cos( )
8 8 8
P20 (0) = 0

 

 β = 0


 γ = 1
64  π 

⇐⇒ α = − 1 + cos ( )

 π2 8
β = 0


64  π 
=⇒ P2 ( x) = 2 −1 + cos( ) x2 + 1
π 8
π 64  π  π2 15 1 π
=⇒ P2 ( ) = 2 −1 + cos( ) 2
+1 = + cos( )
32 π 8 32 16 16 8
π
=⇒ cos(π /32) ' P2 ( ) = 0.995242470781955
32

4
 4. on a cos(π /32) = 0.995184726672197 Le développement de Taylor de dégré 2
donne la meilleure approximation. Cela s’explique par le fait que π /32 soit proche
de 0 et que le polynôme de Taylor utilise les valeurs exactes de f (0) , f 0 (0) et f ”(0).

Exercice 3. :

i xi f [ xi ] f [ xi −1 , xi ] f [ xi −2 , xi −1 , xi ] f [ xi −3 , xi −2 , xi −1 , xi ]
0 0 2
1 10 1.89 -0.011
2 20 1.72 -0.017 -0.0003
3 30 1.44 -0.028 -0.00055 -0.00000833333

1. Pour avoir un polynôme de degré 2 on va prendre trois points d’interpolation, le

polynôme de Newton de degré 2 est donné par :

P2 ( x) = f ( x0 ) + f [ x0 , x1 ]( x − x0 ) + f [ x0 , x1 , x2 ]( x − x0 )( x − x1 )

=⇒ P2 ( x) = 2 − 0.011x − 0.0003x( x − 10)

=⇒ P2 (15) = 1.8125

2. Pour avoir un polynôme de degré 3 on va prendre quatre points d’interpolation,

le polynôme de Newton de degré 2 est donné par :

P3 ( x) = P2 ( x) + f [ x0 , x1 , x2 , x3 ]( x − x0 )( x − x1 )( x − x2 )

=⇒ P3 (15) = 1.8125 + −0.833333 × 10−5 x( x − 10)( x − 10)

= 1.8125 + 0.0031125
=⇒ P3 (15) = 1.815625