Faculté des sciences et ingénierie (Toulouse III) Année universitaire

Département de mathématiques 2019-2020

L2 Maths, UE d’Analyse numérique

Feuille de TD 1 - Correction : Interpolation de Lagrange

Exercice 1. (Identification)
On considère x, y ∈ R4 donnés par : x = [−2, 0, 1, 2] et y = [4, 0, 0, 4]. Parmi les poly-
nômes suivants, lequel est le polynôme d’interpolation P aux points x, y (justifiez votre
réponse) ?
1. P1 (X) = X 4 − 23 X 3 − 3 X 2 + 38 X
2. P2 (X) = 43 X 2 − 4
3
3. P3 (X) = 13 X 3 + X 2 − 34 X.

Correction : On ne demande pas ici de calculer le polynôme mais de l’identifier. On va

donc utiliser la caractérisation équivalente (liée à l’unicité) du polynôme d’interpolation
de Lagrange associé aux points x, y :

P pol d’interp. de Lagrange associé à x, y

⇐⇒ (deg(P ) 6 3, P (−2) = 4, P (0) = 0, P (1) = 0, P (2) = 4) (1)

Il n’y a plus qu’à trouver le polynôme qui satisfait toutes les propriétés de (1) (l’existence
et l’unicité du théorème du cours garantit qu’il existe et est unique). Le polynôme P1
est de degré 4, il est donc éliminé. Le polynôme P2 a un terme constant non nul : il ne
s’annule pas en 0, il est donc éliminé. Reste le polynôme P3 , on vérifie qu’il convient,
c’est donc lui.

Exercice 2. (Existence et unicité)

1. Montrez qu’il existe une infinité de polynômes de degré 2 dont le graphe passe par
les points (0, 0) et (1, 0).
Correction : Cherchons les polynômes de degré 2 p(x) = ax2 + bx + c tels que
p(0) = 0 et p(1) = 0. Ce qui est équivalent au système linéaire
(
c=0
a+b+c=0

En le résolvant, on obtient p(x) = ax(x−1) sans condition sur a, ce qui correspond

bien à une infinité de polynômes de degré 2.

1
 2. Montrez qu’il n’existe pas de polynôme de degré 2 passant par les points (0, 1),
(1, 4), (2, 15) et (3, 40).
Correction : Comme dans la question précédante, on cherche p(x) = ax2 + bx + c
tels que p(0) = 1, p(1) = 4, p(2) = 15 et p(3) = 40. Ce qui est équivalent au
système linéaire 

 c=1

 a+b+c=4

 4a + 2b + c = 15


9a + 3b + c = 40
En le résolvant, on trouve qu’il n’y a pas de solution, ce qui conclut la question.

Exercice 3. (Construction... Malin ou bourrin ?)

Remarque : C’est un bon exercice ici, maintenant que vous avez du recul d’essayer les
différentes façons de calculer un polynôme d’interpolation.
Calculer les polynômes d’interpolation de Lagrange aux points suivants :
a. x = [−1, 2, 3] et y = [4, 4, 8]
Correction : On calcule la base de Lagrange associée à x :
1 −1 1
L0 (X) = (X −2)(X −3), L1 (X) = (X +1)(X −3), L2 (X) = (X +1)(X −2)
12 3 4
et alors Pa (X) = 4L0 (X) + 4L1 (X) + 8L2 (X).
IMPORTANT : Il n’est pas demandé/nécessaire/souhaitable de développer les po-
lynômes de la base de Lagrange ni même de développer Pa , vous allez ajouter des
erreurs et le résultat final sera faux.
b. x = [−2, −1, 0, 1] et y = [0, −2, −4, 0]
Correction : Ici on voit que le polynôme a 2 racines : −2 et 1. Cela signifie
qu’il peut être factorisé par (X + 2)(X − 1), c’est à dire qu’il existe un polynôme
Q tel que Pb (X) = Q(X)(X + 2)(X − 1). Comme on sait que deg(Pb ) 6 3, alors
nécessairement Q est de degré inférieur ou égal à 1 : Q(X) = aX + b.
On cherche maintenant a et b en utilisant les autres valeurs :

Pb (−1) = −2, Pb (0) = −4

ce qui équivaut à (
−2(−a + b) = −2
−2b = −4

ce qui donne b = 2, a = 1 soit Pb (X) = (X + 2)2 (X − 1). Bien sûr, on vérifie a

posteriori que Pb convient bien.
c. x = [−1, 0, 1, 2] et y = [6, 2, 0, 0]

2
 Correction : Ici on procède de la même manière que précédemment en remar-
quant que 1 et 2 sont racines de Pc . On obtient par le même raisonnement que
précédemment
Pc (X) = −(X − 2)(X − 1).
REMARQUE : On peut évidemment calculer Pb et Pc en calculant les polynômes
de degré 3 de la base de Lagrange, mais il n’est pas nécessaire de calculer TOUS
les polynômes de la base : seuls les polynômes où P ne s’annule pas sont utiles (en
l’occurence L2 et L3 pour Pb , L1 et L2 pour Pc ).
d. x = [−1, 0, 1] et y = [1, 0, 1]
Correction : Ici un simple coup d’oeil permet de constater que X 2 convient, par
unicité, on sait donc que Pd (X) = X 2 .
e. x = [−3, −1, 2, 10] et y = [−3, −1, 2, 10]
Correction : Encore plus simple que précédemment, ici Pe (X) = X.

Exercice 4. (Utilisation de la caractérisation)

Soit P un polynôme. Montrer que son polynôme d’interpolation aux noeuds xi ∈ R,
0 6 i 6 n, est le reste de la division euclidienne de p par le polynôme πn (x) = (x −
x0 ) (x − x1 ) . . . (x − xn ).
Correction : Cet exercice vous démunit en général. Dans ce cas, revenons en à la base :
que doit-on démontrer ? On doit démontrer que le reste de la division euclidienne de P
par πn (appelons-le R, on en reparlera plus tard) est LE polynôme d’interpolation de
P aux noeuds xi , i = 1 . . . n, c’est à dire, en utilisant la caractérisation du polynôme
d’interpolation :
deg(R) 6 n, ∀i = 1 . . . n, R(xi ) = P (xi ).
Ca paraît pas mais on a beaucoup avancé en disant ça, car on sait maintenant comment
partir !
Rappelons maintenant comment est défini R :

deg(R) < deg(πn ) = n+1, et il existe un polynôme Q tel que P (X) = Q(X)πn (X)+R(X).

Ceci est la définition du reste de la division euclidienne de P par πn .

On sait donc que deg(R) 6 n, parfait, et maintenant, on évalue P (xi ) pour tout i :

P (xi ) = Q(xi )πn (xi ) + R(xi ).

Or, la définition de πn dit que pour tout i, πn (xi ) = 0. On a donc bien P (xi ) = R(xi )
pour tout i et la preuve est finie !

3
Exercice 5. (Vandermonde et interpolation de Lagrange...)
Pour (x0 , . . . , xn ) ∈ Rn+1 , on considère la matrice
1 x0 x20 xn0
 
...

 1 x1 x21 ... xn1 

 .. .. .. .. 
V (x0 , . . . , xn ) =  . . . ... .
 

 .. .. .. .. 
. . . ... .
 
 
1 xn x2n ... xnn
Y
1. Montrer que det(V (x0 , . . . , xn )) = (xj − xi ).
(i,j), 06i<j6n

Correction : C’est un exercice classique d’algèbre linéaire que vous avez peut être
déjà fait, et si ce n’est pas le cas, il est clair que vous le referez. Il y a plusieurs
manières de s’y prendre, je vous ne propose une par récurrence :
Pour i = 2 . . . n + 1, on remplace la colonne Ci par Ci − x0 Ci−1 , on obtient alors
(le déterminant est inchangé par ces opérations) :
 

 1 0 0 ... 0 


 1 x1 − x0 x21 − x0 x1 . . . xn1 − x0 xn−1
1


 .. .. .. .. .. 
V (x0 , . . . , xn ) =  . . . . . .
 

.. .. .. .. .. 
. . . . .
 
 
1 xn − x0 x2n − x0 xn . . . xnn − x0 xn−1
 
 
n

Ainsi, en développant par rapport à la première ligne, on obtient :

 
 x −x
 1 0 x1 (x1 − x0 ) . . . xn−1
1 (x1 − x0 ) 

 .. .. .. .. 

. . . . 
V (x0 , . . . , xn ) = 
 
.. .. .. .. 
. . . .
 
 
xnn−1 (xn − x0 )
 
 xn − x0 xn (xn − x0 ) . . . 

ce qui donne, par multi-linéarité :

 

 1 x1 . . . x1n−1 

 .. .. .. .. 

. . . . 
V (x0 , . . . , xn ) = (x1 − x0 ) . . . (xn − x0 )  .
 
.. .. .. ..
. . . .
 
 

 1 xn . . . n−1
xn



On conclut par récurrence.

2. Soit (y0 , . . . , yn ) ∈ Rn+1 . Montrer qu’il existe un unique polynôme P ∈ Rn tel que
P (xi ) = yi si et seulement si xi 6= xj pour tout (i, j), i 6= j.
Correction : C’est la preuve qui a été faite en amphi. Je la refais ici.
Soit donc (y0 , . . . , yn ) ∈ Rn+1 , l’existence et l’unicité d’un tel polynôme est équi-
valente à l’existence et l’unicité de coefficients a0 , . . . , an tels que (en cherchant un

4
 tel polynôme P sous la forme P (X) = a0 + a1 X + · · · + an X n et en écrivant que
pour tout i = 0 . . . n, P (xi ) = yi ) :


 a0 + a1 x0 + a2 x20 + · · · + an xn0 = y0
a0 + a1 x1 + a2 x21 + · · · + an xn1 = y1



 ...
a0 + a1 xn + a2 x2n + · · · + an xnn = yn



c’est à dire l’existence et l’unicité d’un vecteur (a0 , . . . , an ) ∈ Rn+1 tel que (en
ré-écrivant le système sous forme matricielle
V (x0 , . . . , xn )(a0 , . . . , an )T = (y0 , . . . , yn )T .
Or, d’après ce qui précède, V (x0 , . . . , xn ) est inversible si et seulement si les xi
sont deux à deux distincts. On a donc existence et unicité d’un tel polynôme si et
seulement si les xi sont deux à deux distincts.
Exercice 6. (Construction. . .)
Calculer le polynôme P de degré inférieur ou égal à 4 tel que :
1. P (−2) = 11, P (−1) = 1, P (0) = 1, P (1) = 5, P (2) = 31.
Correction : À moins d’avoir envie de se fader le calcul de l’inverse d’une matrice
de Vandermonde de taille 5 ou de calculer les 5 polynômes de la base de Lagrange
associée à ces noeuds, le mieux est sans doute ici d’utiliser la base de Newton. On
obtient en faisant le tableau des différences divisées
1
P (X) = 11−10(X+2)+5(X+2)(X+1)+(X+2)(X+1)X+ (X+2)(X+1)X(X−1).
2
2. P (−1) = 4, P 0 (−1) = −4, P (0) = 0, P (1) = 0, P 0 (1) = 0.
Correction : Un exercice un peu différent ici puisqu’il ne s’agit pas d’interpo-
lation de Lagrange : on impose aussi des valeurs aux dérivées de P aux noeuds
d’interpolation ! Quelle idée ! ! Comme souvent, deux méthodes sont envisageables
ici : la méthode "bourrin" et la méthode "malin".
La méthode bourrin consiste à chercher le polynôme sous forme indéterminée et
écrire les 5 équations vérifiées par ses coefficients. On obtiendra comme pour l’in-
terpolation de Lagrange un système linéaire de taille 5 à résoudre. Courage !
Sinon, on remarque que le polynôme que l’on cherche a le bon goût d’avoir une
racine simple : 0 et une racine double : 1 (c’est-à-dire que P ET P 0 s’annulent
en 1). On sait donc qu’on peut le factoriser par X(X − 1)2 et on le cherche donc
(puisqu’on sait qu’il est de degré inéfrieur ou égal à 4) sous la forme
P (X) = X(X − 1)2 (aX + b).
Il ne reste plus qu’à chercher a et b en utilisant les valeurs de P et de P 0 en −1.
On obtient après calcul (
a−b=1
3a + 2b = −1

5
 ce qui donne a = 45 , b = − 15 et donc finalement, P (X) = 15 X(X − 1)2 (4X − 1).
Exercice 7. (Base de Lagrange)
Soit x0 , . . . , xn (n+1) réels distincts deux à deux. Pour k ∈ {0, . . . , n}, on note
Y x − xj
Lk (x) =
x − xj
j∈{0,...,n},j6=k k

le k-ième polynôme de Lagrange.

1. Montrer que Lk est un polynôme de degré n vérifiant Lk (xi ) = δki pour tous
k, i ∈ {0, . . . , n}.
2. En déduire que la famille de polynômes {Lk }k∈{0,...,n} forme une base de Rn [X].

Correction : Cet exercice fait l’objet d’une des preuves les plus importantes du cours.
Je vous renvoie donc au cours (ou au poly).
Exercice 8. (examen 2016) (Exercice optionnel, pour aller plus loin)
Soient x0 = 0 < x1 < . . . < xn et des réels donnés yi , 0 6 i 6 n. On considère le
polynôme d’interpolation satisfaisant

P (x0 ) = y0 , P (−xi ) = P (xi ) = yi , pourtous1 6 i 6 n.

1. Montrer que le polynôme P est pair.

Correction : Cette question est un peu moins classique que le reste du TD, c’est
pourquoi cet exercice n’a pas été abordé en TD. Plusieurs d’entre vous m’en ont
demandé une correction, la voici. Je la détaille à l’extrême pour en faciliter la
compréhension. N’hésitez pas à me contacter pour toute question.
Pour simplifier les notations on va noter, pour i = 1 . . . n : x−i = −xi . On a
donc alors de l’interpolation avec 2n + 1 noeuds : x−i , xi pour i = 1 . . . n et 0. Le
polynôme que l’on cherche est donc de degré inférieur ou égal à 2n.
On rappelle par ailleurs que P n’est pas forcément de degré 2n. De plus, le fait
d’être de degré pair n’entraine pas que P soit pair. En effet, P (X) = X 2 + X + 1
n’est par exemple ni pair ni impair. Pour être pair, P doit être une somme de
polynômes pair (qui sont eux même des sommes ou produits de polynômes pairs) :
P (−X) = P (X).
On propose de commencer par se faire une idée de ce qui se passe ici en commençant
par le cas où n = 1. On a alors trois points : 0, x1 et −x1 .
On écrit les 3 polynômes de la base de Lagrange associée à 0, −x1 , x1 .
−1 −1
(X 2 − x21 )


 L0 (X) = x21
(X − x1 )(X + x1 ) = x21

1
L1 (X) = 2x2 X(X + x1 )
1
L−1 (X) = 2x12 X(X − x1 ).



1

On a alors :

P (X) = y0 L0 (X) + y1 L1 (X) + y−1 L−1 (X) = y0 L0 (X) + y1 (L1 (X) + L−1 (X))

6
puisque y−1 = y1 . On constate alors que L0 est pair et que
1 X2
L1 (X) + L−1 (X) = 2 X (X + x1 + X − x1 ) = 2
2x1 x1
est pair. P est donc finalement la somme de deux polynômes pairs. Il est donc pair.
Voyons maintenant ce qui se passe dans le cas général n > 1. On procède de la
même façon : on va calculer L0 , Lk et L−k pour chaque k = 1 . . . n.
Étape 1 : calcul de L0 :
n n n n
! !
Y X − xi Y X − x−i Y (X − xi )(X − x−i ) Y (X − xi )(X + xi )
L0 (X) = = = .
i=1
−xi i=1
−x−i i=1
xi .x−i i=1
−x2i

Où on a utilisé le fait que x−i = −xi . On a finalement :

n
Y (X 2 − x2i )
L0 (X) = qui est donc pair.
i=1
−x2i

Étape 2 : calcul de Lk et L−k , pour k ∈ {1 . . . n} fixé :

   
n n n n
! !
Y X − xi  Y X − x−i Y X − xi  Y X + xi
Lk (X) =  . = . .
x − xi
i=1,i6=k k
x − x−i
i=1 k
x − xi
i=1,i6=k k
x + xi
i=1 k

On isole alors le terme i = k dans le 2ème produit (qui correspond au terme i = −k

du polynôme) et on regroupe les autres termes dans le même produit :
 
n n
X + xk  Y (X − xi )(X + xi )  X + xk Y X 2 − x2i
Lk (X) = = .
2xk i=1,i6=k
(xk − xi )(xk + xi ) 2xk x2 − x2i
i=1,i6=k k

Il est alors clair que

n
X − xk Y X 2 − x2i
L−k (X) = .
−2xk x2 − x2i
i=1,i6=k k

Étape 3 : calcul de P :
Comme dans le cas n = 1, on écrit la décomposition de P dans la base de Lagrange :
n
X n
X
P = y0 L0 + (yk Lk + y−k L−k ) = y0 L0 + yk (Lk + L−k )
k=1 k=1

puisque y−k = yk par hypothèse. On a vu que L0 est pair, calculons maintenant

Lk + L−k pour tout k > 1 grâce aux calculs de l’étape 2.
n n
X − xk X + xk X 2 − x2i X 2 − x2i
  Y Y
Lk (X) + L−k (X) = + =
−2xk 2xk x2 − x2i
i=1,i6=k k
x2 − x2i
i=1,i6=k k

qui est donc pair. Le polynôme P est donc une somme de polynômes pairs. Il est
donc pair.

7
 2. En déduire en un minimum de calculs le polynôme d’interpolation vérifiant

P (−1) = 2, P (0) = 4, P (1) = 2.

Correction : On peut utiliser le cas que l’on a étudié (avec 3 points) pour se faire
une idée dans la question précédente. On obtient donc
−y0 2 y1
P (X) = 2 (X − x21 ) + 2 X 2 .
x1 x1
Dans le cas présent x1 = 1, y0 = 4, y1 = 2. On obtient alors

P (X) = −4(X 2 − 1) + 2X 2 = −2X 2 + 4.

On vérifie bien que P est pair et qu’il convient.

Exercice 9. (examen 1999)
1. Calculer le polynôme d’interpolation de Lagrange de la fonction f (x) = x (x2 − 1)
relativement aux points x0 = −1, x1 = 1 et x2 = 2.
Correction : Les valeurs de la fonction f aux noeuds en question sont les sui-
vantes :
f (−1) = f (1) = 0, f (2) = 6.
Il y a comme toujours deux façons de faire ici, on va présenter les deux méthodes
(elles conduisent bien sûr au même polynôme par unicité).
— On utilise les racines −1 et 1 : en effet on cherche ici un polynôme de degré
inférieur ou égal à 2, on sait qu’il admet 1 et −1 pour racine, donc P est de la
forme : P (X) = α(X − 1)(X + 1), avec α un réel à trouver. Pour le trouver, on
utilise la valeur en 2 : P (2) = 3α = 6, d’où α = 2. On a donc trouvé

P (X) = 2(X − 1)(X + 1).

— On peut sinon utiliser une manière plus systématique de calculer ce polynôme

(système, base de Lagrange ou base de Newton). En jetant un coup d’oeil à
la question suivante, on choisit la base de Newton pour éviter de devoir tout
refaire ! Le calcul des différences divisées donne le tableau suivant :
xk f [xk ] f [xk , xk+1 ] f [x0 , x1 , x2 ]

−1 0

1 0 0

2 6 6 2

On obtient ainsi dans la base de Newton :

P (X) = 0 + 0 ∗ (X + 1) + 2(X + 1)(X − 1) = 2(X − 1)(X + 1).

8
 2. Même question en rajoutant le point x3 = −2.
Correction : Toujours deux façons de faire en fonction de celle que l’on a choisie
à la question précédente.
— Si on a compris comment on construit les différences divisées, on peut utiliser
le polynôme calculé à la question précédente à l’aide des racines −1 et 1. On va
ajouter à ce polynôme un polynôme de degré 3 qui s’annule en −1, 1 et 2 (pour
garder le valeurs du polynôme calculé prcédemment en ces points : on cherche
donc β un réel tel que

Q(X) = P (X) + β(X − 1)(X + 1)(X − 2), Q(−2) = 6 − 12β

Comme on ajoute la contrainte Q(−2) = f (−2) = −6, on a nécessairement

β = 1 et donc

Q(X) = 2(X − 1)(X + 1) + (X − 1)(X + 1)(X − 2).

— On peut tout aussi bien (et c’est ,comme on l’a vu en cours, l’interêt principal
de la base de Newton) reprendre le tableau des différences divisées précédent
et rajouter une ligne avec x3 = −2 et f (x3 ) = −6. Il suffit alors de calculer
les différences divisées qui apparaissent sur la dernière ligne et de rajouter le
terme manquant à notre polynôme
xk f [xk ] f [xk , xk+1 ] f [xk , xk+1 , xk+2 ] f [x0 , x1 , x2 , f (x3 )]

−1 0

1 0 0

2 6 6 2

−2 −6 3 1 1
ce qui conduit au polynôme Q écrit dans la base de Newton :

Q(X) = 2(X − 1)(X + 1) + (X − 1)(X + 1)(X − 2).

Exercice 10. (partiel 2003)

1. Rappeler l’expression de la base de Newton de R5 [X] associée aux noeuds 1, 2, 3,
4, 5, 6.
Correction : Il s’agit du cours, elle est donnée par la famille

{1, X − 1, (X − 1)(X − 2), (X − 1)(X − 2)(X − 3), (X − 1)(X − 2)(X − 3)(X − 4),
(X − 1)(X − 2)(X − 3)(X − 4)(X − 5)}

2. Montrer qu’il s’agit bien d’une base.

9
 Correction : D’une manière générale, pour montrer qu’une famille de d vecteurs
d’un espace vectoriel de dimension d il suffit de montrer qu’elle est libre, i.e que
toute combinaison linéaire nulle de tous ses éléments est nécessairement la com-
binaison nulle. Ici, on a 6 éléments et R5 [X] est de dimension 6, il suffit donc de
montrer que cette famille est libre. Supposons trouvés des réels a0 , . . . , a5 tels que
5
X k
Y
R(X) = a0 + ak (X − i) = 0. (2)
k=1 i=1

On cherche à montrer que tous les ai sont nuls. Il est très important de comprendre
que l’équation (2) est une égalité SUR DES POLYNOMES et donc que cette égalité
est vraie pour toutes les valeurs possibles de X. On a en quelque sorte une infinité
d’équations ! ! ! On écrit cette égalité pour X = 1, ce qui annule tous les termes sauf
le premier et donne a0 = 0. De même, en évaluant en X = 2, sachant que a0 = 0,
on obtient a1 = 0. On fait de même avec X = 3, 4, 5 et on obtient a2 = a3 = a4 = 0.
Pour a5 on peut choisir d’évaluer en n’importe quelle valeur de X, par exemple
X = 0 et on obtient a5 = 0, la famille est donc libre.
3. Donner, dans cette base, l’expression du polynôme P ∈ R5 [X] tel que
P (1) = P (6) = 1, P (2) = P (3) = P (4) = P (5) = 0.
Correction : Ici on vous demande de calculer ce polynôme dans la base de Newton,
pas de question à se poser donc, on fonce sur le tableau de différences divisées...
ATTENTION : Utiliser la méthode des racines est malin bien sûr, mais HORS
SUJET puisqu’on vous demande d’exprimer le polynôme dans la base de New-
ton (sauf miracle comme dans l’exercice précédent, les expressions seront toujours
différentes)
xk f [xk ] f [xk , xk+1 ] ... ... ... f [1, 2, 3, 4, 5, 6]

1 1
2 0 -1
1
3 0 0 2

−1
4 0 0 0 6

1
5 0 0 0 0 24

1 1 1
6 1 1 2 6 24 0

On obtient donc le polynôme suivant :

1 1
P (X) = 1 − (X − 1) + (X − 1)(X − 2) − (X − 1)(X − 2)(X − 3)
2 6
1
+ (X − 1)(X − 2)(X − 3)(X − 4).
24

10
 4. Calculer P (0).
On calcule avec le polynôme précédent P (0) = 5.

Exercice 11. (Partiel 2014)

Étant donnés six réels x1 , a, b, c, d et e, on considère le tableau de différences divisées
suivant :

xk f [xk ] f [xk , xk+1 ] f [xk , . . . , xk+2 ] f [x0 , x1 , x2 , x3 ]

x0 = 0 1

x1 −1 1

x2 = −1 0 b d

2
x3 = 2 a c e 3

1. Calculer x1 , a, b, c, d et e.
Correction : En écrivant les formules des différences divisées on obtient les équa-
tions suivantes :

−2


 x1 =1




 1



 −1−x1 =b




a
3 =c





c−b
=e


4







d−e 2
=


−2 3

Ce qui finalement donne

4 19
x1 = −2, b = 1, d = 0, e= , c= , a = 19.
3 3
2. Donner, dans la base de Newton le polynôme P3 qui interpole (0, 1), (x1 , −1), (−1, 0)
et (2, a).
Correction : Le polynôme associé à ce tableau de différences divisées est donné
par :
2 2
P (X) = 1 + X + dX(X − x1 ) + X(X − x1 )(X + 1) = 1 + X + X(X + 2)(X + 1).
3 3

11
3. On considère les fonctions suivantes définies sur R :
( (
2 + 9 x2 si x > 0 0 si x > −1
f1 : x 7→ , f2 : x 7→
0 sinon −3 x2 − x3 sinon

Pour α et β deux réels, on définit la fonction f : x ∈ R 7→ α f1 (x) + β f2 (x).

Montrer que P3 est le polynôme d’interpolation de f en x0 , x1 , x2 et x3 si et
seulement si α = 21 et β = 41 .
Correction : En utilisant la caractérisation équivalente, P3 est le polynôme d’in-
terpolation de f en x0 , x1 , x2 et x3 si et seulement si

deg(P3 ) 6 3, ∀i = 0 . . . 3, P3 (xi ) = f (xi ) = αf1 (xi ) + βf2 (xi ).

Détaillons P3 (x0 ) = f (x0 ). On sait que x0 = 0, que P3 (x0 ) = 1 (par défintion) et

f (0) = αf1 (0) + βf2 (0) = 2α.

On a donc 2α = 1 soit α = 21 .
1
On obtient avec la valeur en x1 que β = 4 et on vérifie que ces valeurs collent bien
avec les égalités en x2 et x3 .
4. À partir de maintenant, on pose α = 12 et β = 14 .
Calculer dans la base de Newton le polynôme P4 qui interpole les points (0, 1),
(x1 , −1), (−1, 0), (2, a) et (1, 6).
P4 est-il le polynôme d’interpolation de f en 0, x1 , −1, 2 et 1 ?
Correction : Cela revient à rajouter x5 = 1 au problème d’interpolation de La-
grange étudié à la question 1. Il faut donc calculer la valeur en 1 de P4 : P4 (1) = 6
et ajouter une ligne au tableau de différence divisées :

xk P [xk ] f [xk , xk+1 ] f [xk , . . . , xk+2 ] f [xk , xk+1 , xk+2 , xk+3 ] f [x0 , x1 , x2 , x3 , x4 ]

x0 = 0 1

x1 = −2 −1 1

x2 = −1 0 1 0

19 4 2
x3 = 2 19 3 3 3

10 2
x4 = 1 6 13 3 3 0

Le polynôme d’interpolation P4 est alors donné par la formule du cours dans la

base de Newton :
2
P4 (X) = 1 + X + X(X + 2)(X + 1).
3

12
 Vous remarquerez peut-être qu’il s’agit finalement de P3 . En effet, si on évalue P3
en 1 on trouve bien 6. C’est un petit miracle.
Reste à savoir si P4 est le polynôme de f en 0, x1 , −1, 2 et 1. On a bien deg(P4 ) =
3 6 4. Par ailleurs, on sait déjà (par construction) que P4 a les mêmes valeurs
que f en 0, x1 , −1 et 2. Reste à vérifier donc la valeur en 1. En 1, P4 vaut 6 par
construction, et f (1) = 11/2 6= 6. Donc P4 n’est pas le polynôme d’interpolation
de f aux points 0, x1 , −1, 2 et 1.

Exercice 12. (Convergence uniforme)

On considère a = x0 < x1 < . . . < xn = b et le polynôme d’interpolation Pn tel que

∀i = 0, . . . n, Pn (xi ) = exi .

Montrer que la suite de polynômes d’interpolation Pn converge uniformément vers la

fonction exponentielle lorsque n tend vers l’infini, c’est-à-dire que

sup |Pn (x) − ex | −→ 0

x∈[a,b] n→∞

Correction : Pas besoin ici de connaître quoi que ce soit sur la convergence uniforme.
On vous demande juste de montrer que

sup |Pn (x) − ex | −→ 0.

x∈[a,b] n→∞

Pour cela, bien sûr on va utiliser la seule formule du cours qui permette d’obtenir une
estimation sur la différence entre un polynôme d’interpolation et sa fonction (ici la
fonction exponentielle). Ici, Pn est en effet LE polynôme d’interpolation de Lagrange
associé à x0 . . . , xn de la fonction exponentielle. Le théorème vu en cours assure que

kf (n+1) k∞,[a,b]
∀x ∈ [a, b], |Pn (x) − f (x)| 6 |πx0 ...,xn |(x)
(n + 1)!
Qn
où πx0 ...,xn (X) = (X − xi ).
i=0
Or, dans le cas de la fonction exponentielle, on connaît les dérivées successives : c’est
elle même. On a donc ici kf (n+1) k∞,[a,b] = eb . On a de plus (c’est toujours le cas même
si en général cette estimation est très mauvaise) supx∈[a,b] |Πx0 ...,xn |(x) 6 (b − a)n+1 .
Finalement,
(b − a)n+1
sup |Pn (x) − f (x)| 6 eb −→ 0.
x∈[a,b] (n + 1)! n→∞

13