Chapitre 3 : Interpolation polynômiale
0.1 Position du problème
Soit f une fonction dénie et continue sur un intervalle [a,b]. On suppose que l'on
connait seulement sa valeur en n+1 points xi de [a,b] que l'on suppose ordonnés :
a ≤ x0 < x2 < x3 < . . . < xn < xn ≤ b
Interpoler la fonction f par un polynôme de degré n aux points x0 , x1 , . . . , xn
consiste à trouver un polynôme p de degré ≤ n tel que
f (xi ) = pn (xi ), 0≤i≤n (1)
Les points (xi , f (xi )) sont appelés points d'interpolation ou d'appui.
Théorème 1 Une condition nécessaire et susante pour qu'il existe un polynôme
Pn unique interpolant f est que les points d'interpolation xi soient tous distincts.
Théorème 2 Il existe un seul polynôme Pn de degré n interpolant f aux points
d'interpolation xi , i = 0, . . . , n.
0.2 Le polynôme d'interpolation de Lagrange
0.2.1 Construction
Soient x0 , x1 , . . . , xn n+1 réels donnés distincts. On dénit n+1 polynômes li
pour i = 0, ..., n par
(x − x0 )(x − x1 )...(x − xn ) x − xj
li (x) = = Πnj=0,j6=i
(xi − x0 )(xi − x1 )...(xi − xn ) xi − x j
Le numérateur de chacun de ces polynômes est un produit de n termes (x−xk ),
il est donc un polynôme de degré n. Le dénominateur est une constante. Donc
chaque polynôme li vérie :
1. li est un polynôme de degré n.
2. li (xk ) = 0 si i 6= k pour k = 0, ..., n.
3. li (xk ) = 1 si i = k (li (xi ) = 1).
1
Proposition
Les polynômes l0 (x), l1 (x), ..., ln (x) forment une base de Rn [X]
Démonstration :
Il sut de montrer que la famille l0 (x), l1 (x), ..., ln (x) est libre et génératrice.
Donc le polynôme pn vériant le problème (1.1) peut s'écrire sous la forme
n
X
pn (x) = αi li (x)
i=0
puisque les li forment une base de Rn [X]. En prenant x = xk pour k = 0, ..., n et
en utilisant le fait que li (xk ) = δi,k , on obtient
αk = pn (xk ) = f (xk )
Proposition :
L'unique solution du problème (1.1) est donc
n
X
pn (x) = f (xi )li (x)
i=0
Dénition :
Le polynôme pn s'appelle l'interpolant de Lagrange de la fonction f de degré n aux
points x0 , x1 , ..., xn . Il s'appelle aussi le polynôme d'interpolation de Lagrange aux
points x0 , x1 , ..., xn .
Exercice :
Déterminer le polynôme d'interpolation de Lagrange satisfaisant au tableau ci-
dessous :
x 0 2 3 5
f(x) -1 2 9 87
0.3 Le polynôme d'interpolation de Newton
0.3.1 Base de Newton
Les polynômes ek de la base de Newton sont dénis comme suit :
k−1
ek = Πi=0 (x − xi ) = (x − x0 )(x − x1 )...(x − xk−1 ), k = 1, ...n
avec pour convention
e0 = 1
en outre
e1 = (x − x0 )
2
e2 = (x − x0 )(x − x1 )
...
...
...
en = (x − x0 )(x − x1 )...(x − xn−1 )
L'ensemble des polynomes (ek )0≤k≤n (base de Newton) forment une base de
l'espace Pn des polynomes de degré au plus n, puisqu'il s'agit d'une famille de
(n+1) polynomes de degré échelonné.
Donc le polynome d'interpolation de Newton de degré n relatif à la subdivision
(x0 , y0 ), (x1 , y1 ), ..., (xn , yn ) s'écrit :
n
X
Pn (x) = αk ek (x)
k=0
avec
Pn (xi ) = f (xi ), i = 0, ..., n
Il faut alors déterminer les coecients (αk )0≤k≤n , qui sont basés sur les dié-
rences divisées.
0.3.2 Diérences divisées
Dénition 1 Soit f une fonction dont on connait les valeurs en n+1 points dis-
tincts x0 , x2 , . . . , xn . On appelle diérence divisée d'ordre 0, 1, . . . , n de la fonction
f les expressions suivantes :
Ordre Notation Dénition
0 f [xi ] f (xi )
f (xi ) − f (xj )
1 f [xi , xj ], i 6= j
xi − xj
f [xi , xj ] − f [xj , xk ]
2 f [xi , xj , xk ] i 6= j, j 6= k, i 6= k
xi − xk
... ... ...
f [x1 , . . . , xm−1 ] − f [x2 , . . . , xm ]
m-1 f [x1 , . . . , xm ]
x1 − xm
... ... ...
f [x1 , . . . , xn ] − f [x2 , . . . , xn+1 ]
n f [x1 , . . . , xn+1 ]
x1 − xn+1
3
4
5
Exercice :
1. Ecrire le polynôme d'interpolation de f, noté P4 , construit sur les données
suivantes : (-1,-1), (0,1), (1,0) et (2,0) en utilisant la formule de Newton et
les diérences divisées.
2. Calculer Pi (20) pour i = 1, . . . , 4 et comparer à f (20).