ESIM

TD 2
Interpolation polynômiale
A.U. 2023-2024

Exercice 1

Déterminer le polynôme d’interpolation de Lagrange passant par les points (0, 6), (2, 0) et (6, 12).

Corrigé
Soit P2 ∈ P2 [X] l’interpolant de Lagrange à déterminer. Pour tout x ∈ R, on a

X
2
P2 (x) = yi Li (x) = f(x0 )L0 (x) + f(x1 )L1 (x) + f(x2 )L2 (x),
i=0

avec
(x − x1 )(x − x2 )
L0 (x) =
(x0 − x1 )(x0 − x2 )
(x − 2)(x − 6)
=
−2 × −6
1
= (x − 2)(x − 6)
12
(x − x0 )(x − x2 )
L1 (x) =
(x1 − x0 )(x1 − x2 )
1
= − x(x − 6)
8
(x − x0 )(x − x1 )
L2 (x) =
(x2 − x0 )(x2 − x1 )
1
= x(x − 2).
24
Le polynôme de Lagrange est :
6 12
P2 (x) = (x − 2)(x − 6) + x(x − 2)
12 24
= x2 − 5x + 6.

Exercice 2

Soit les trois points d’appuis (0, 1); (1, 0.5) et (3, 0.25) de la fonction f(x). Déterminer:

1. Le polynôme de Lagrange passant par ces points.

2. Une approximation de f(1.5).

Corrigé
ESIM

page 1/5
 1. Soit P2 ∈ P2 [X] l’interpolant de Lagrange à déterminer. Pour tout x ∈ R, on a
X
2
P2 (x) = yi Li (x) = y0 L0 (x) + y1 L1 (x) + y2 L2 (x),
i=0
avec
(x − x1 )(x − x2 )
L0 (x) =
(x0 − x1 )(x0 − x2 )
(x − 1)(x − 3)
=
−1 × −3
1
= (x − 1)(x − 3)
3
(x − x0 )(x − x2 )
L1 (x) =
(x1 − x0 )(x1 − x2 )
1
= − x(x − 3)
2
(x − x0 )(x − x1 )
L2 (x) =
(x2 − x0 )(x2 − x1 )
1
= x(x − 1).
6
Le polynôme de Lagrange est :
1 0.5 0.25
P2 (x) = (x − 1)(x − 3) − x(x − 3) + x(x − 1)
3 2 6
= 0.125x2 − 0.625x + 1

2. Approximation de P2 (1.5) = 0.125(1.5)2 − 0.625(1.5) + 1 = 0.344

Exercice 3

1. Déterminer P3 le polynôme d’interpolation de Lagrange satisfaisant au tableau ci-dessous

xi 0 2 3 5
f(xi ) -1 2 9 87

2. Donner une valeur approchée de f(4.2) par cette interpolation.

Corrigé

1. Soit P3 ∈ P3 [X] l’interpolant de Lagrange à déterminer. Pour tout x ∈ R, on a

X
3
P3 (x) = yi Li (x) = y0 L0 (x) + y1 L1 (x) + y2 L2 (x) + y3 L3 (x),
i=0
avec
x−2 x−3 x−5 (x − 2)(x − 3)(x − 5)
L0 (x) = × × =−
−2 −3 −5 30
x x−3 x−5 x(x − 3)(x − 5)
L1 (x) = × × =
2 −1 −3 6
x x−2 x−5 x(x − 2)(x − 5)
L2 (x) = × × =−
ESIM

3 1 −2 6
x x−2 x−3 x(x − 2)(x − 3)
L3 (x) = × × = .
5 3 2 30
page 2/5
 Ainsi
(x − 2)(x − 3)(x − 5) x(x − 3)(x − 5) 3 87
P3 (x) = + − x(x − 2)(x − 5) + x(x − 2)(x − 3)
30 3 2 30
53 3 2 253
= x − 7x + x − 1.
30 30

2. Une valeur approchée de f(4.2) : f(4.2) ≃ P3 (4.2) = 41, 8288.

Exercice 4

Le tableau suivant donnant les valeurs v de la viscosité cinématique (m2 .s−1 ) de l’eau en fonction
de la température T (◦ C).
T 15 16 17 18 19 20 21 22 23 24 25
v 1.14 1.11 1.08 1.06 1.03 1.01 0.983 0.960 0.938 0.917 0.894

1. Par interpolation linéaire, à quelle température la viscosité vaut-elle v = 1 m2 .s−1 ?

2. Par interpolation quadratique, quelle est la viscosité à 24.5 ◦ C ?

Corrigé

1. On considère le tableau suivant:

T 20 21
v 1.01 0.983

P1 (T ) = L0 (t)f(T0 ) + L1 (t)f(T1 ) = 1.55 − 0.027T.

Si v = 1, alors 1 = 1.55 − 0.027T , donc T = 20.37◦ C.

2. On considère le tableau suivant:

T 23 24 25
v 0.938 0.917 0.894

P2 (T ) = L0 (t)f(T0 ) + L1 (t)f(T1 ) + L2 (T )f(T2 ) = 0.863 + 0.026T − 0.001T 2 .

Si T = 24.5◦ C, alors P2 (24.5) = 0.863 + 0.026(24.5) − 0.001(24.5)2 = 0.8997 m2 .s−1 .

Exercice 5

1. Déterminer le polynôme de Lagrange qui passe par les points donnés par le tableau

xi 0 2 4 6
f(xi ) 0 4 0 4
ESIM

2. Retrouver ce polynôme en utilisant la formule des différences divisées.

3. Retrouver ce polynôme en utilisant la formule des différences finies progressives.

page 3/5
 Corrigé

1. Le polynôme de Lagrange
X
3
P3 (x) = f(xi )Li (x)
i=0
= f(x0 )L0 (x) + f(x1 )L1 (x) + f(x2 )L2 (x) + f(x3 )L3 (x),

avec
(x − x1 )(x − x2 )(x − x3 ) (x − 2)(x − 4)(x − 6)
L0 (x) = =−
(x0 − x1 )(x0 − x2 )(x0 − x3 ) 48
(x − x0 )(x − x2 )(x − x3 ) x(x − 4)(x − 6)
L1 (x) = =
(x1 − x0 )(x1 − x2 )(x1 − x3 ) 16
(x − x0 )(x − x1 )(x − x3 ) x(x − 2)(x − 6)
L2 (x) = =−
(x2 − x0 )(x2 − x1 )(x2 − x3 ) 16
(x − x0 )(x − x1 )(x − x2 ) x(x − 2)(x − 4)
L3 (x) = =
(x3 − x0 )(x3 − x1 )(x3 − x2 ) 48
Ainsi
4 4
P3 (x) = x(x − 4)(x − 6) + x(x − 2)(x − 4)
16 48
x3 20x
= − 3x2 + .
3 3

2. Le polynôme de Newton :

P3 (x) = f (x0 ) + ∆ [x0 , x1 ] (x − x0 ) + ∆2 [x0 , x1 , x2 ] (x − x0 ) (x − x1 )

+ ∆3 [x0 , x1 , x2 , x3 ] (x − x0 ) (x − x1 ) (x − x2 )


Les différences divisées DDi , i = 1, 2, 3 sont calculées dans le tableau suivant :

xi f(xi ) DD1 DD2 DD3

x0 = 0 f(x0 ) = 0
x1 = 2 f(x1 ) = 4 ∆ [x0 , x1 ] = 2
x2 = 4 f(x2 ) = 0 ∆ [x1 , x2 ] = −2 ∆2 [x0 , x1 , x2 ] = −1
1
x3 = 6 f(x3 ) = 4 ∆ [x2 , x3 ] = 2 ∆2 [x1 , x2 , x3 ] = 1 ∆3 [x0 , x1 , x2 , x3 ] = 3

On déduit que :
1
P3 (x) = 0 + 2 (x − x0 ) − 1 (x − x0 ) (x − x1 ) + (x − x0 ) (x − x1 ) (x − x2 )
3
1
= 0 + 2(x − 0) − 1(x − 0)(x − 2) + (x − 0)(x − 2)(x − 4)
3
3
x 20x
= − 3x2 + .
3 3

3. Remarquons que les points donnés sont équidistants avec le pas d’interpolation h = 2, alors
sous la forme de Newton progressive, le polynôme d’interpolation de f est donné par :
∇f (x0 ) ∇2 f (x0 ) ∇3 f (x0 )
ESIM

P3 (x) = f (x0 )+ (x − x0 )+ (x − x 0 ) (x − x 1 )+ (x − x0 ) (x − x1 ) (x − x2 ) .
h 2!h2 3!h3
Le calcul des différences finies progressives se fait comme suit:

page 4/5
 xi ∇f(xi ) ∇1 f(xi ) ∇2 f(xi ) ∇3 f(xi )
x0 = 0 f(x0 ) = 0
x1 = 2 f(x1 ) = 4 ∇f(x0 ) = 4
x2 = 4 f(x2 ) = 0 ∇f(x1 ) = −4 ∇2 f(x0 ) = −8
x3 = 6 f(x3 ) = 4 ∇f(x2 ) = 4 ∇2 f(x1 ) = 8 ∇3 f(x0 ) = 16

D’où,
4 −8 16
P3 (x) = 0 + (x − x0 ) + (x − x0 ) (x − x1 ) + (x − x0 ) (x − x1 ) (x − x2 )
2 2! × 4 3! × 8
x3 20x
= − 3x2 + .
3 3
ESIM

page 5/5