Chapitre 2 Interpolation polynomiale

Interpolation polynomiale

2.1 Introduction
Dans ce chapitre, on dispose d'une fonction f connue uniquement par ses valeurs en
certains points, et on cherche à approcher f par une fonction polynôme. Le problème de
l’interpolation polynomiale consiste à trouver un polynôme de degré inférieur ou égal à n dont
la courbe passe par les n + 1 points donnés.
Supposons par exemple que nous connaissons les valeurs d’une fonction 𝑓 en un nombre
fini de points distincts x0, x1,…, xn de l'intervalle [a, b] selon le tableau suivant :
xi x0 x1 … xn
f(xi) f(x0) f( x1) … f( xn)

Pour estimer la valeur de f en un point quelconque x  IR , on peut construire un polynôme

P de degré inférieur ou égal à 𝑛 tel que :

P( xi )  f ( xi ) pour i  0,1, 2, ....n et utiliser l’approximation P ( x)  f ( x) .

C’est ce qu’on appelle Interpolation polynomiale : interpolation de la fonction f par le
polynôme Pn(x) aux points x0, x1,…, xn.

Les points x0, x1,…, xn sont appelés points d'interpolation (nœuds).

2.2 Méthodes Numériques

2.2.1 Interpolation de Lagrange
[Link] Théorème et définition
Soient n+1 points distincts x0, x1,…, xn de l'intervalle [a, b] et et f une fonction dont les
valeurs sont f(x0), f(x1), … … , f(xn), il existe un unique polynôme Pn de degré inférieur ou
égal à n tel que :
f ( xi )  Pn ( xi ) pour i  0, 1, 2, .... n
Le polynôme s'écrit :
n
Pn ( x)   f ( xi ) Li ( x)  f ( x0 ) L0 ( x)  f ( x1 ) L1 ( x)  ....  f ( xn ) Ln ( x)
i 0

avec yi  f ( xi )

Où
20
 Chapitre 2 Interpolation polynomiale

n (x  x j ) ( x  x0 )( x  x1 )....(x  xi 1 )( x  xi 1 )....(x  xn )
Li ( x)  
j 0, j i ( xi  x j )

( xi  x0 )( xi  x1 )....(xi  xi 1 )( xi  xi 1 )....(x i  xn )

Le polynôme Pn est appelé polynôme d'interpolation de Lagrange de la fonction f aux points

x0, x1,…, xn.
Les polynômes Li(x) sont appelés coefficients polynômes de Lagrange.

[Link] Quelques exemples simples

- Polynôme de degré 1 (n = 1) :
Le polynôme de Lagrange, passant par 2 points, est une droite.

P1 ( x)  f 0 L0 ( x)  f1 L1 ( x) P1(x)
f1
( x  x1 ) ( x  x0 )
 f0  f1
( x0  x1 ) ( x1  x0 )
f0
Ce qui s’écrit encore :

( y1  y0 )
P1 ( x)  ( x  x0 )  y 0 x0 x1
( x1  x0 )

- Polynôme de degré 2 (n = 2) :
Le polynôme de Lagrange, passant par 3 points, est une f
2

parabole. P2(x)
f1
P2 ( x)  f 0 L0 ( x)  f1 L1 ( x)  f 2 L2 ( x) f0

( x  x1 )( x  x2 ) ( x  x0 )( x  x2 )
 f0  f1 x0 x1 x2
( x0  x1 )( x0  x2 ) ( x1  x0 )( x1  x2 )
( x  x0 )( x  x1 )
 f2
( x2  x0 )( x2  x1 )

Exemple :
Déterminer le polynôme d’interpolation de Lagrange satisfaisant au tableau ci-dessous :

xi x0 =0 x1=1 x2=2 x3=3

f(xi ) f0=-4 f1=-2 f2=2 f3=14

21
 Chapitre 2 Interpolation polynomiale

Solution :
Dans ce cas on a 4 points, donc le polynôme de degré inférieur ou égal à 3.

Le polynôme d’interpolation de Lagrange est :

3
P3 ( x)   f ( xi ) Li ( x)  f ( x0 ) L0 ( x)  f ( x1 ) L1 ( x)  f ( x2 ) L2 ( x)  f ( x3 ) L3 ( x)
i 0
On calcul les coefficients polynômes de Lagrange :
( x  x1 )( x  x2 )( x  x3 ) ( x  1)( x  2)( x  3) 1
L0 ( x)     ( x 3  6 x 2  11x  6)
( x0  x1 )( x0  x2 )( x0  x3 ) (0  1)(0  2)(0  3) 6
( x  x0 )( x  x2 )( x  x3 ) ( x  0)( x  2)( x  3) 1 3
L1 ( x)    ( x  5 x 2  6 x)
( x1  x0 )( x1  x2 )( x1  x3 ) (1  0)(1  2)(1  3) 2
( x  x0 )( x  x1 )( x  x3 ) ( x  0)( x  1)( x  3) 1
L2 ( x)     ( x3  4 x 2  3x)
( x2  x0 )( x2  x1 )( x2  x3 ) (2  0)(2  1)(2  3) 2
( x  x0 )( x  x1 )( x  x2 ) ( x  0)(x  1)( x  2) 1 3
L3 ( x)    ( x  3x 2  2 x )
( x3  x0 )( x3  x1 )( x3  x2 ) (3  0)(3  1)(3  2) 6
Finalement on remplace les coefficients polynômes et on obtient :
 1  1   1 
P3 ( x)  ( 4)  ( x 3  6 x 2  11x  6)   (2) ( x 3  5 x 2  6 x)   2  ( x 3  4 x 2  3 x) 
 6  2   2 
1 
 14 ( x 3  3 x 2  2 x) 
6 

P3 ( x)  x 3  2 x 2  3x  4

2.2.2 Interpolation de Newton

[Link] Différences divisées d’une fonction
Forme générale :
La forme de Newton du polynôme d’interpolation est :
Pn ( x)  a0  a1 ( x  x0 )  a2 ( x  x0 )( x  x1 )  ....  an ( x  x0 )( x  x1 )....( x  xn1 )
Où Pn est le polynôme de degré inférieur ou égal à n qui interpole f aux points x0, x1,…, xn.
Les coefficients an seront trouvés en utilisant des différences divisées d’ordre n de la fonction
f.

Pn ( x0 )  a0  f ( x0 )
f ( x1 )  f ( x0 )
Pn ( x1 )  a0  a1 ( x1  x0 )  f ( x0 )  a1 ( x1  x0 )  f ( x1 )  a1 
x1  x0

On procède de la même manière jusqu’à l’obtention de an.

22
Chapitre 2 Interpolation polynomiale

Définition 2.1:
Soit une fonction définie sur [a, b] et x0, x1,…, xn (n+1) points de [a, b] distincts. On
appelle différences divisées d’ordre k de f les relations de récurrences suivantes :

- d’ordre 0 : f xi   f ( xi )
f ( xi 1 )  f ( xi )
- d’ordre 1 : f xi , xi 1  
xi 1  xi
f xi 1 , xi  2   f xi , xi 1 
- d’ordre 2 : f xi , xi 1 , xi  2  
xi  2  x i
…
f xi 1 , xi  k   f xi , xi  k 1 
- d’ordre k : f xi , xi 1 ,..., xi  k  
xi  k  xi
D’après cette définition on a :
a 0  f x 0   f ( x 0 )

a1  f x 0 , x1   f ( x1 )  f ( x0 )
 x1  x0

 f ( x 2 )  f ( x1 ) f ( x1 )  f ( x 0 )

 f x1 , x 2   f x0 , x1  x 2  x1 x1  x 0
a 2  f x 0 , x1 , x 2   
 x2  x0 x2  x0
 

a n  f x 0 , x1 ,..., x n 




Ces coefficients an sont les différences divisées d’ordre n de la fonction f.

Théorème 2.1:
Le polynôme d’interpolation Pn(x) de degré inférieur ou égal à n passant par les points (xi ,
f(xi)), i=0,1,.., n peut s’écrire:

Pn ( x)  f ( x0 )  f x0 , x1 ( x  x0 )  f x0 , x1 , x 2 ( x  x0 )( x  x1 )  .... 

 f x0 , x1 ,..., xn ( x  x0 )( x  x1 )....( x  x n1 )
Ce polynôme d’interpolation est appelé forme de Newton du polynôme d’interpolation par les
différences divisées.

Calcul des différences divisées :

Le tableau des différences divisées permet de calculer les différences divisées d'une fonction
selon le schéma suivant :

23
Chapitre 2 Interpolation polynomiale

xi f xi  f xi , xi 1  f xi , xi 1 , xi  2  … f xi ,..., xi  n 

x0 f x0 
f x0 , x1 
x1 f x1  f x0 , x1 , x2 
f x1 , x2  …
f x0 , x1 ,..., xn 
x2 f x2  …
...
…
… … f xn  2 , xn 1 , xn 
f xn 1 , xn 
xn f xn 

Exemple :
Trouver le polynôme d’interpolation de Newton qui passe par les points suivants (0,1), (2,5)
et (4,17).
Solution:
On a 3 points, donc le degré du polynôme est 2.

Le polynôme d’interpolation de Newton est :

P2 ( x)  f ( x0 )  f x0 , x1 ( x  x0 )  f x0 , x1 , x2 ( x  x0 )( x  x1 )

On construit le tableau des différences divisées de f :

xi f xi  f xi , xi 1  f xi , xi 1 , xi  2 

x0  0 f x0   1
5 1
f x0 , x1   2
20
x1  2 f x1   5 62
f x0 , x1 , x2   1
17  5 40
f x1 , x2   6
42

f x2   17
x2  4

On obtient :

P2 ( x )  1  2 x  x ( x  2)

P2 ( x)  1  x 2

24
Chapitre 2 Interpolation polynomiale

[Link] Différences finies progressives (cas des points équidistants)

Dans ce cas les points d’interpolation sont en progression arithmétiques, i.e. :
x0 , x1  x0  h, x2  x0  2h, ...., xn  x0  nh, h0

Définition 2.2:
Soient f ( xi )  yi pour i=0, 1, …, n des nombres réels. On appelle différences finies
d’ordre k de f les relations de récurrences suivantes :

- d’ordre 1 :  f ( xi )  f ( xi 1 )  f ( xi ) pour i  0,1, ....., n  1

d’ordre 2 :  f ( xi )   f ( xi 1 )   f ( xi ) pour i  0,1, ....., n  2
2
-
…
k 1
- d’ordre k :  f ( xi )  
k
f ( xi 1 )  k 1 f ( xi ) pour i  0,1, ....., n  k

Convention  f ( xi )  f ( xi ) , pour i  0,1, ....., n

0

Table des différences finies progressives :

xi f ( xi )  f ( xi ) 2 f ( xi ) 3 f ( xi ) …
n f ( xi )

x0 f ( x0 )
 f ( x0 )
x1 f ( x1 ) 2 f ( x0 )
 f ( x1 ) 3 f ( x0 )
…
x2 f ( x2 ) 2 f ( x1 )
… n f ( x0 )
 f ( x2 )
…
x3 f ( x3 ) …
... 3 f ( xn  3 )

… … 2 f ( xn2 )
 f ( xn 1 )
xn f ( xn )

Relation entre les différences finies progressives et les différences

divisées :
Théorème 2.2:
Soit f une fonction dont on connait les valeurs f ( xi )  yi , i=0, 1, …, n, avec
xi  xi 1  ih, h  0 . Alors :
25
Chapitre 2 Interpolation polynomiale

k ( f ( xi ))
f xi , xi 1 , ...., xi  k   pour 0  i  i  k  n
h k k!
Où f xi , xi 1 , ...., xi  k  est la différence divisée d’ordre k de f aux points xi , xi 1 , ...., xi  k et
k ( f ( xi )) est la différence finie d’ordre k au point f ( xi ) .

Polynôme d’interpolation de Newton par les différences finies :

Théorème 2.3:
Soient x0, x1,…, xn points équidistants. Le polynôme d’interpolation Pn(x) de degré
inférieur ou égal à n passant par ces points peut s’écrire:
 f ( x0 ) 2 f ( x0 )
Pn ( x)  f ( x0 )  ( x  x0 )  ( x  x0 )( x  x1 )  .... 
1!h 2 !h 2
n f ( x0 )
 ( x  x0 )( x  x1 )....( x  x n1 )
n !h n
Ce polynôme d’interpolation est appelé forme de Newton du polynôme d’interpolation par les
différences finies.
Exemple :
Trouver le polynôme d’interpolation de Newton passant par les points suivants :
xi 0 2 4
f(xi) 1 5 17
Solution :
On a 3 points, donc le degré de polynôme est n=2.

On remarque que les points donnés sont équidistants avec le pas d'interpolation h =2, alors
sous la forme de Newton par les différences finies, le polynôme d'interpolation de f est donné
par :
 f ( x0 ) 2 f ( x0 )
P2 ( x)  f ( x0 )  ( x  x0 )  ( x  x0 )( x  x1 )
1!h 2 !h 2
Où : h  x1  x0  x2  x1  2

Le calcul des différences finies se fait comme suit :

xi f ( xi )  f ( xi ) 2 f ( xi )
x0  0 f ( x0 )  1
 f ( x0 )  5  1  4
x1  2 f ( x1 )  5 2 f ( x0 )  12  4  8
 f ( x1 )  17  5  12
x2  4 f ( x2 )  17

26
Chapitre 2 Interpolation polynomiale

D'où,

P2 ( x )  1  2 x  x ( x  2)

P2 ( x)  1  x 2

2.3 Estimation de l’erreur d’interpolation

Soit f : a , b   IR de classe C . Notons Pn son polynôme d’interpolation aux nœuds x0,
n 1

x1,…, xn dans a, b  . Alors pour tout x  a , b  , il existe c  a, b  tel que :

f n 1 (c) n M n
E n ( x)  f ( x)  Pn ( x)   ( x  xi )  (n  1)! 
( n  1)! i 0 i 0
( x  xi )

n 1
Où M  sup f ( x)
xa ,b 

2.4 Exercices
Exercice 1:
Soit les points d’interpolation suivants : (0,1), (1,3), (2,7).

1/ Quel est le degré n du polynôme d’interpolation P ?

2/ Trouver le polynôme d’interpolation passant par les trois points :
a/ Par une méthode d’identification.
b/ à l’aide des polynômes de Lagrange.
Exercice 2 :
1
Soit la fonction f ( x) 
x 1
1/ Déterminer le polynôme de degré 2 qui interpole la fonction f (x) en x0  0, x1  1, x2  3

2/ Estimer la valeur de f(1.5) en utilisant le polynôme trouvé en (a).

3/ Donner une borne supérieure de l’erreur commise.

Exercice 3 :
Les valeurs d’une fonction f (x) pour trois valeurs de x sont données par :
X 3.5 4.0 4.5

f (x ) 0.9086 1.0000 1.0772

27