Chapitre 4: Interpolation et approximation polynômiale de

fonctions

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025
Introduction
Principe : Chercher des fonctions simples (polynômes, polynômes par morceaux,
polynômes trigonométriques) passant par des points données :
(x0 , y0 ), (x1 , y1 ), · · · , (xn , yn ).
Exemple :
Soit une expérience où on enregistre la distance parcourue par un objet en fonction du
temps, les résultats sont donnés dans le tableau suivant :

t(s) 0 1 2 3 4
x(m) 0 5 15 0 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 2 / 22
Introduction
Principe : Chercher des fonctions simples (polynômes, polynômes par morceaux,
polynômes trigonométriques) passant par des points données :
(x0 , y0 ), (x1 , y1 ), · · · , (xn , yn ).
Exemple :
Soit une expérience où on enregistre la distance parcourue par un objet en fonction du
temps, les résultats sont donnés dans le tableau suivant :

t(s) 0 1 2 3 4
x(m) 0 5 15 0 3

But : avoir une forme analytique de x en fonction de t, x(t); afin de calculer la position
de l’objet au temps t = 2.5s.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 2 / 22
 Introduction
Principe : Chercher des fonctions simples (polynômes, polynômes par morceaux,
polynômes trigonométriques) passant par des points données :
(x0 , y0 ), (x1 , y1 ), · · · , (xn , yn ).
Exemple :
Soit une expérience où on enregistre la distance parcourue par un objet en fonction du
temps, les résultats sont donnés dans le tableau suivant :

t(s) 0 1 2 3 4
x(m) 0 5 15 0 3

But : avoir une forme analytique de x en fonction de t, x(t); afin de calculer la position
de l’objet au temps t = 2.5s.

=⇒ Dans ce qui suit, nous allons chercher à approximer la fonction x(t) à l’aide d’un
polynôme de degré n, s’exprimant sous la forme :

x(t) = a0 + a1 t + a2 t 2 + · · · + an t n

où ai (i = 0, 1, · · · , n) sont des coefficients à déterminer.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 2 / 22
 Introduction
Principe : Chercher des fonctions simples (polynômes, polynômes par morceaux,
polynômes trigonométriques) passant par des points données :
(x0 , y0 ), (x1 , y1 ), · · · , (xn , yn ).
Exemple :
Soit une expérience où on enregistre la distance parcourue par un objet en fonction du
temps, les résultats sont donnés dans le tableau suivant :

t(s) 0 1 2 3 4
x(m) 0 5 15 0 3

But : avoir une forme analytique de x en fonction de t, x(t); afin de calculer la position
de l’objet au temps t = 2.5s.

=⇒ Dans ce qui suit, nous allons chercher à approximer la fonction x(t) à l’aide d’un
polynôme de degré n, s’exprimant sous la forme :

x(t) = a0 + a1 t + a2 t 2 + · · · + an t n

où ai (i = 0, 1, · · · , n) sont des coefficients à déterminer.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 2 / 22
Introduction

f (x) = x · sin(πx)
Pn (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 3 / 22
Introduction

f (x) = x · sin(πx)
Pn (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 3 / 22
Introduction

f (x) = x · sin(πx)
Pn (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 3 / 22
Introduction
f (x) = x · sin(πx)
Pn (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 3 / 22
Introduction
f (x) = x · sin(πx)
Pn (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 3 / 22
 Introduction

Définition
Soit f une fonction réelle définie sur un intervalle [a, b] contenant n + 1 points
distincts x0 , x1 , · · · , xn .
Soit Pn un polynôme de degré inférieur ou égal à n.
On dit que Pn est un interpolant de f ou interpole f en x0 , x1 , · · · , xn si :

Pn (xi ) = f (xi ) pour 0 ≤ i ≤ n

Théorème
Un polynôme de degré n dont la forme générale est :

Pn (x) = a0 + a1 x + a2 x 2 + · · · + an x n (an ̸= 0)

possède très exactement n racines qui peuvent être réelles ou complexes conjuguées.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 4 / 22
 Introduction

Définition
Soit f une fonction réelle définie sur un intervalle [a, b] contenant n + 1 points
distincts x0 , x1 , · · · , xn .
Soit Pn un polynôme de degré inférieur ou égal à n.
On dit que Pn est un interpolant de f ou interpole f en x0 , x1 , · · · , xn si :

Pn (xi ) = f (xi ) pour 0 ≤ i ≤ n

Théorème
Un polynôme de degré n dont la forme générale est :

Pn (x) = a0 + a1 x + a2 x 2 + · · · + an x n (an ̸= 0)

possède très exactement n racines qui peuvent être réelles ou complexes conjuguées.


=⇒ Il existe un unique polynôme Pn qui interpole les nœuds (xi , f (xi )) i=0:n

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 4 / 22
 Introduction

Définition
Soit f une fonction réelle définie sur un intervalle [a, b] contenant n + 1 points
distincts x0 , x1 , · · · , xn .
Soit Pn un polynôme de degré inférieur ou égal à n.
On dit que Pn est un interpolant de f ou interpole f en x0 , x1 , · · · , xn si :

Pn (xi ) = f (xi ) pour 0 ≤ i ≤ n

Théorème
Un polynôme de degré n dont la forme générale est :

Pn (x) = a0 + a1 x + a2 x 2 + · · · + an x n (an ̸= 0)

possède très exactement n racines qui peuvent être réelles ou complexes conjuguées.


=⇒ Il existe un unique polynôme Pn qui interpole les nœuds (xi , f (xi )) i=0:n

Question : Comment construire le polynôme Pn ?

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 4 / 22
 Introduction

Définition
Soit f une fonction réelle définie sur un intervalle [a, b] contenant n + 1 points
distincts x0 , x1 , · · · , xn .
Soit Pn un polynôme de degré inférieur ou égal à n.
On dit que Pn est un interpolant de f ou interpole f en x0 , x1 , · · · , xn si :

Pn (xi ) = f (xi ) pour 0 ≤ i ≤ n

Théorème
Un polynôme de degré n dont la forme générale est :

Pn (x) = a0 + a1 x + a2 x 2 + · · · + an x n (an ̸= 0)

possède très exactement n racines qui peuvent être réelles ou complexes conjuguées.


=⇒ Il existe un unique polynôme Pn qui interpole les nœuds (xi , f (xi )) i=0:n

Question : Comment construire le polynôme Pn ?

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 4 / 22
Construction du polynôme Pn

On peut écrire Pn dans plusieurs bases de Rn :

1 base canonique : (1, x, x 2 , x 3 , · · · , x n )


2 base de Lagrange : L0 (x), L1 (x), L2 (x), · · · , Ln (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 5 / 22
Construction du polynôme Pn

On peut écrire Pn dans plusieurs bases de Rn :

1 base canonique : (1, x, x 2 , x 3 , · · · , x n )


2 base de Lagrange : L0 (x), L1 (x), L2 (x), · · · , Ln (x)


3 base de Newton : 1, N1 (x), N2 (x), · · · , Nn (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 5 / 22
Construction du polynôme Pn

On peut écrire Pn dans plusieurs bases de Rn :

1 base canonique : (1, x, x 2 , x 3 , · · · , x n )


2 base de Lagrange : L0 (x), L1 (x), L2 (x), · · · , Ln (x)


3 base de Newton : 1, N1 (x), N2 (x), · · · , Nn (x)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 5 / 22
Construction du polynôme Pn : base canonique

On cherche un polynôme de degré n tel que Pn (xi ) = f (xi ), i = 0, 1, · · · , n.

Cette méthode consiste à déterminer (a0 , a1 , a2 , · · · , an ) tel que :
X n
Pn (x) = ak x k
k=0
Écrivons explicitement Pn (xi ) = f (xi ) :

a0 + a1 xi + a2 xi2 + · · · + an−1 xin−1 + an xin = f (xi )

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 6 / 22
Construction du polynôme Pn : base canonique

On cherche un polynôme de degré n tel que Pn (xi ) = f (xi ), i = 0, 1, · · · , n.

Cette méthode consiste à déterminer (a0 , a1 , a2 , · · · , an ) tel que :
X n
Pn (x) = ak x k
k=0
Écrivons explicitement Pn (xi ) = f (xi ) :

a0 + a1 xi + a2 xi2 + · · · + an−1 xin−1 + an xin = f (xi )

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 6 / 22
Construction du polynôme Pn : base canonique

On cherche un polynôme de degré n tel que Pn (xi ) = f (xi ), i = 0, 1, · · · , n.

Cette méthode consiste à déterminer (a0 , a1 , a2 , · · · , an ) tel que :
X n
Pn (x) = ak x k
k=0
Écrivons explicitement Pn (xi ) = f (xi ) :

a0 + a1 xi + a2 xi2 + · · · + an−1 xin−1 + an xin = f (xi )




 a0 + a1 x0 + a2 x02 + · · · + an−1 x0n−1 + an x0n = f (x0 )
a0 + a1 x1 + a2 x12 + · · · + an−1 x1n−1 + an x1n = f (x1 )


=⇒ .. ..



 . .
a + a x + a x 2 + · · · + a x n−1 + a x n = f (x )

0 1 n 2 n n−1 n n n n

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 6 / 22
Construction du polynôme Pn : base canonique

On cherche un polynôme de degré n tel que Pn (xi ) = f (xi ), i = 0, 1, · · · , n.

Cette méthode consiste à déterminer (a0 , a1 , a2 , · · · , an ) tel que :
X n
Pn (x) = ak x k
k=0
Écrivons explicitement Pn (xi ) = f (xi ) :

a0 + a1 xi + a2 xi2 + · · · + an−1 xin−1 + an xin = f (xi )




 a0 + a1 x0 + a2 x02 + · · · + an−1 x0n−1 + an x0n = f (x0 )
a0 + a1 x1 + a2 x12 + · · · + an−1 x1n−1 + an x1n = f (x1 )


=⇒ .. ..



 . .
a + a x + a x 2 + · · · + a x n−1 + a x n = f (x )

0 1 n 2 n n−1 n n n n

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 6 / 22
 Construction du polynôme Pn : base canonique

On cherche un polynôme de degré n tel que Pn (xi ) = f (xi ), i = 0, 1, · · · , n.

Cette méthode consiste à déterminer (a0 , a1 , a2 , · · · , an ) tel que :
X n
Pn (x) = ak x k
k=0
Écrivons explicitement Pn (xi ) = f (xi ) :

a0 + a1 xi + a2 xi2 + · · · + an−1 xin−1 + an xin = f (xi )




 a0 + a1 x0 + a2 x02 + · · · + an−1 x0n−1 + an x0n = f (x0 )
a0 + a1 x1 + a2 x12 + · · · + an−1 x1n−1 + an x1n = f (x1 )


=⇒ .. ..



 . .
a + a x + a x 2 + · · · + a x n−1 + a x n = f (x )

0 1 n 2 n n−1 n n n n

Sous forme matricielle :

        −1  
1 x0 x02 · · · x0n a0 f (x0 ) a0 1 x0 ··· x0n f (x0 )
1
 x1 x12 · · · x1n  a  f (x )
  1 1 
a  1
1
   x1 ··· x1n  f (x1 )
 
  ..  =  ..  =⇒  ..  =  .. .. ..   .. 
. .. .. .. ..         ..  
.
. . . . .  .   .   .  . . . .  . 

n
1 xn xn2 · · · xnn an f (xn ) an 1 xn ··· xn f (xn )
| {z }
Matrice de Vandermonde
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 6 / 22
 Construction du polynôme Pn : base canonique

On cherche un polynôme de degré n tel que Pn (xi ) = f (xi ), i = 0, 1, · · · , n.

Cette méthode consiste à déterminer (a0 , a1 , a2 , · · · , an ) tel que :
X n
Pn (x) = ak x k
k=0
Écrivons explicitement Pn (xi ) = f (xi ) :

a0 + a1 xi + a2 xi2 + · · · + an−1 xin−1 + an xin = f (xi )




 a0 + a1 x0 + a2 x02 + · · · + an−1 x0n−1 + an x0n = f (x0 )
a0 + a1 x1 + a2 x12 + · · · + an−1 x1n−1 + an x1n = f (x1 )


=⇒ .. ..



 . .
a + a x + a x 2 + · · · + a x n−1 + a x n = f (x )

0 1 n 2 n n−1 n n n n

Sous forme matricielle :

        −1  
1 x0 x02 · · · x0n a0 f (x0 ) a0 1 x0 ··· x0n f (x0 )
1
 x1 x12 · · · x1n  a  f (x )
  1 1 
a  1
1
   x1 ··· x1n  f (x1 )
 
  ..  =  ..  =⇒  ..  =  .. .. ..   .. 
. .. .. .. ..         ..  
.
. . . . .  .   .   .  . . . .  . 

n
1 xn xn2 · · · xnn an f (xn ) an 1 xn ··· xn f (xn )
| {z }
Matrice de Vandermonde
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 6 / 22
Construction du polynôme Pn : base canonique

Exemple :
On souhaite trouver le polynôme passant par les points (1, 2), (2, 3), (4, 1).
La formule de la base canonique à 3 points est :

P2 (x) = a0 + a1 x + a2 x 2

avec,

             
1 x0 x02 a0 f (x0 ) 1 1 1 a0 2 a0 −1
   3
1 x1 x12  a1  = f (x1 ) =⇒ 1 2 4  a1  = 3 =⇒ a1  =  3 
         
1 x2 x22 a2 f (x2 ) 1 4 16 a2 1 a2 − 23

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 7 / 22
Construction du polynôme Pn : base canonique

Exemple :
On souhaite trouver le polynôme passant par les points (1, 2), (2, 3), (4, 1).
La formule de la base canonique à 3 points est :

P2 (x) = a0 + a1 x + a2 x 2

avec,

             
1 x0 x02 a0 f (x0 ) 1 1 1 a0 2 a0 −1
   3
1 x1 x12  a1  = f (x1 ) =⇒ 1 2 4  a1  = 3 =⇒ a1  =  3 
         
1 x2 x22 a2 f (x2 ) 1 4 16 a2 1 a2 − 23

Le polynôme recherché est donc :

P2 (x) = 0.9x 2 − 2.9x + 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 7 / 22
Construction du polynôme Pn : base canonique

Exemple :
On souhaite trouver le polynôme passant par les points (1, 2), (2, 3), (4, 1).
La formule de la base canonique à 3 points est :

P2 (x) = a0 + a1 x + a2 x 2

avec,

             
1 x0 x02 a0 f (x0 ) 1 1 1 a0 2 a0 −1
   3
1 x1 x12  a1  = f (x1 ) =⇒ 1 2 4  a1  = 3 =⇒ a1  =  3 
         
1 x2 x22 a2 f (x2 ) 1 4 16 a2 1 a2 − 23

Le polynôme recherché est donc :

P2 (x) = 0.9x 2 − 2.9x + 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 7 / 22
Construction du polynôme Pn : base canonique
Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points, sous la forme P3 (x) = a0 + a1 x + a2 x 2 + a3 x 3 .

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 8 / 22
Construction du polynôme Pn : base canonique
Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points, sous la forme P3 (x) = a0 + a1 x + a2 x 2 + a3 x 3 .

(x0 , f (x0 )) = (0, 1)



(x , f (x )) = (1, 2)
1 1


 (x2 , f (x2 )) = (2, 9)

(x3 , f (x3 )) = (3, 28)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 8 / 22
Construction du polynôme Pn : base canonique
Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points, sous la forme P3 (x) = a0 + a1 x + a2 x 2 + a3 x 3 .

(x0 , f (x0 )) = (0, 1)



(x , f (x )) = (1, 2)
1 1


 (x2 , f (x2 )) = (2, 9)

(x3 , f (x3 )) = (3, 28)

         
1 x0 x02 x03 a0 f (x0 ) 1 0 0 0 a0 1
1 x1 x12 x13  a1  f (x1 ) 1 1 1 1  a 1   2 
         
  =   =⇒    =  
1 x2 x22 x23  a2  f (x2 ) 1 2 4 8  a 2   9 

1 x3 x32 x33
a3 f (x3 ) 1 3 9 27 a3 28

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 8 / 22
 Construction du polynôme Pn : base canonique
Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points, sous la forme P3 (x) = a0 + a1 x + a2 x 2 + a3 x 3 .

(x0 , f (x0 )) = (0, 1)



(x , f (x )) = (1, 2)
1 1


 (x2 , f (x2 )) = (2, 9)

(x3 , f (x3 )) = (3, 28)

         
1 x0 x02 x03 a0 f (x0 ) 1 0 0 0 a0 1
1 x1 x12 x13  a1  f (x1 ) 1 1 1 1  a 1   2 
         
  =   =⇒    =  
1 x2 x22 x23  a2  f (x2 ) 1 2 4 8  a 2   9 

1 x3 x32 x33
a3 f (x3 ) 1 3 9 27 a3 28

La solution obtenue par décomposition LU est [a0 , a1 , a2 , a3 ]T = [1, 0, 0, 1]T ,

et le polynôme recherché est donc :
P3 (x) = 1 + x 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 8 / 22
 Construction du polynôme Pn : base canonique
Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points, sous la forme P3 (x) = a0 + a1 x + a2 x 2 + a3 x 3 .

(x0 , f (x0 )) = (0, 1)



(x , f (x )) = (1, 2)
1 1


 (x2 , f (x2 )) = (2, 9)

(x3 , f (x3 )) = (3, 28)

         
1 x0 x02 x03 a0 f (x0 ) 1 0 0 0 a0 1
1 x1 x12 x13  a1  f (x1 ) 1 1 1 1  a 1   2 
         
  =   =⇒    =  
1 x2 x22 x23  a2  f (x2 ) 1 2 4 8  a 2   9 

1 x3 x32 x33
a3 f (x3 ) 1 3 9 27 a3 28

La solution obtenue par décomposition LU est [a0 , a1 , a2 , a3 ]T = [1, 0, 0, 1]T ,

et le polynôme recherché est donc :
P3 (x) = 1 + x 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 8 / 22
Construction du polynôme Pn : base canonique
Exercice 2 :
Trouver le polynôme passant par les points (−1, 8), (0, 3), (1, 6).

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 9 / 22
Construction du polynôme Pn : base canonique
Exercice 2 :
Trouver le polynôme passant par les points (−1, 8), (0, 3), (1, 6).
Étant donné 3 points, on va chercher le polynôme P2 de degré 2 qui interpole ces points,
sous la forme P3 (2) = a0 + a1 x + a2 x 2 .

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 9 / 22
Construction du polynôme Pn : base canonique
Exercice 2 :
Trouver le polynôme passant par les points (−1, 8), (0, 3), (1, 6).
Étant donné 3 points, on va chercher le polynôme P2 de degré 2 qui interpole ces points,
sous la forme P3 (2) = a0 + a1 x + a2 x 2 .

(x0 , f (x0 )) = (−1, 8)


(x1 , f (x1 )) = (0, 3)

(x2 , f (x2 )) = (1, 6)


Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 9 / 22
Construction du polynôme Pn : base canonique
Exercice 2 :
Trouver le polynôme passant par les points (−1, 8), (0, 3), (1, 6).
Étant donné 3 points, on va chercher le polynôme P2 de degré 2 qui interpole ces points,
sous la forme P3 (2) = a0 + a1 x + a2 x 2 .

(x0 , f (x0 )) = (−1, 8)


(x1 , f (x1 )) = (0, 3)

(x2 , f (x2 )) = (1, 6)


         
1 x0 x02 a0 f (x0 ) 1 −1 1 a0 8
1 x1 x12  a1  = f (x1 ) =⇒ 0 0 0 a1  = 3
          
1 x2 x22 a2 f (x2 ) 1 1 1 a2 6

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 9 / 22
 Construction du polynôme Pn : base canonique
Exercice 2 :
Trouver le polynôme passant par les points (−1, 8), (0, 3), (1, 6).
Étant donné 3 points, on va chercher le polynôme P2 de degré 2 qui interpole ces points,
sous la forme P3 (2) = a0 + a1 x + a2 x 2 .

(x0 , f (x0 )) = (−1, 8)


(x1 , f (x1 )) = (0, 3)

(x2 , f (x2 )) = (1, 6)


         
1 x0 x02 a0 f (x0 ) 1 −1 1 a0 8
1 x1 x12  a1  = f (x1 ) =⇒ 0 0 0 a1  = 3
          
1 x2 x22 a2 f (x2 ) 1 1 1 a2 6

La solution obtenue est [a0 , a1 , a2 ]T = [3, −1, 4]T ,

et le polynôme recherché est donc :

P2 (x) = 4x 2 − x + 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 9 / 22
 Construction du polynôme Pn : base canonique
Exercice 2 :
Trouver le polynôme passant par les points (−1, 8), (0, 3), (1, 6).
Étant donné 3 points, on va chercher le polynôme P2 de degré 2 qui interpole ces points,
sous la forme P3 (2) = a0 + a1 x + a2 x 2 .

(x0 , f (x0 )) = (−1, 8)


(x1 , f (x1 )) = (0, 3)

(x2 , f (x2 )) = (1, 6)


         
1 x0 x02 a0 f (x0 ) 1 −1 1 a0 8
1 x1 x12  a1  = f (x1 ) =⇒ 0 0 0 a1  = 3
          
1 x2 x22 a2 f (x2 ) 1 1 1 a2 6

La solution obtenue est [a0 , a1 , a2 ]T = [3, −1, 4]T ,

et le polynôme recherché est donc :

P2 (x) = 4x 2 − x + 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 9 / 22
Construction du polynôme Pn : base de Lagrange

La base de Lagrange est (L0 (x), L1 (x), · · · , Ln (x)).

Le polynôme interpolateur Pn s’écrit sous la forme :

n
X
Pn (x) = ai Li (x)
i=0

avec,
n  
x − x0 x − x1 x − xi−1 x − xi+1 x − xn Y x − xj
Li (x) = · ··· · ··· =
xi − x0 xi − x1 xi − xi−1 xi − xi+1 xi − xn j=0
xi − xj
j̸=i

=⇒ Le terme (x − xi ) n’apparaı̂t pas explicitement dans l’expression de Li (x).

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 10 / 22
Construction du polynôme Pn : base de Lagrange

La base de Lagrange est (L0 (x), L1 (x), · · · , Ln (x)).

Le polynôme interpolateur Pn s’écrit sous la forme :

n
X
Pn (x) = ai Li (x)
i=0

avec,
n  
x − x0 x − x1 x − xi−1 x − xi+1 x − xn Y x − xj
Li (x) = · ··· · ··· =
xi − x0 xi − x1 xi − xi−1 xi − xi+1 xi − xn j=0
xi − xj
j̸=i

=⇒ Le terme (x − xi ) n’apparaı̂t pas explicitement dans l’expression de Li (x).

Chaque fonction Li (x) satisfait les conditions suivantes :

(
Li (xi ) = 1 ∀i
Li (xj ) = 0 ∀j ̸= i

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 10 / 22
Construction du polynôme Pn : base de Lagrange

La base de Lagrange est (L0 (x), L1 (x), · · · , Ln (x)).

Le polynôme interpolateur Pn s’écrit sous la forme :

n
X
Pn (x) = ai Li (x)
i=0

avec,
n  
x − x0 x − x1 x − xi−1 x − xi+1 x − xn Y x − xj
Li (x) = · ··· · ··· =
xi − x0 xi − x1 xi − xi−1 xi − xi+1 xi − xn j=0
xi − xj
j̸=i

=⇒ Le terme (x − xi ) n’apparaı̂t pas explicitement dans l’expression de Li (x).

Chaque fonction Li (x) satisfait les conditions suivantes :

(
Li (xi ) = 1 ∀i
Li (xj ) = 0 ∀j ̸= i

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 10 / 22
Construction du polynôme Pn : base de Lagrange

Écrivons explicitement Pn (xj ) = f (xj ) :

a0 L0 (xj ) + a1 L1 (xj ) + a2 L2 (xj ) + · · · + an Ln (xj ) = f (xj )



 · · · + an Ln (x0 )
a0 L0 (x0 ) + a1 L1 (x0 ) + a2 L2 (x0 ) + |{z} = f (x0 )

 | {z } | {z } | {z } | {z }

 =1 =0 =0 =0 =0


a0 L0 (x1 ) + a1 L1 (x1 ) + a2 L2 (x1 ) + · · · + an Ln (x1 )
 = f (x1 )
| {z } | {z } | {z } |{z} | {z }
 =0 =1 =0 =0 =0



 .. ..


 . .

a0 L0 (xn ) + a1 L1 (xn ) + a2 L2 (xn ) + · · · + an Ln (xn )

= f (xn )
| {z } | {z } | {z } |{z} | {z }
=0 =0 =0 =0 =1

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 11 / 22
Construction du polynôme Pn : base de Lagrange

Écrivons explicitement Pn (xj ) = f (xj ) :

a0 L0 (xj ) + a1 L1 (xj ) + a2 L2 (xj ) + · · · + an Ln (xj ) = f (xj )

 

 · · · + an Ln (x0 )
a0 L0 (x0 ) + a1 L1 (x0 ) + a2 L2 (x0 ) + |{z} = f (x0 ) 

 a0 = f (x0 )

 | {z } | {z } | {z } | {z } 


 =1 =0 =0 =0 =0



 

a0 L0 (x1 ) + a1 L1 (x1 ) + a2 L2 (x1 ) + · · · + an Ln (x1 )
 = f (x1 ) a1 = f (x1 )

| {z } | {z } | {z } |{z}
=0
| {z } =⇒

 =0 =1 =0 =0 



 .. .. 

 ..


 . . 

 .
 
a0 L0 (xn ) + a1 L1 (xn ) + a2 L2 (xn ) + · · · + an Ln (xn )

= f (xn ) a n

= f (xn )
| {z } | {z } | {z } |{z} | {z }
=0 =0 =0 =0 =1

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 11 / 22
Construction du polynôme Pn : base de Lagrange

Écrivons explicitement Pn (xj ) = f (xj ) :

a0 L0 (xj ) + a1 L1 (xj ) + a2 L2 (xj ) + · · · + an Ln (xj ) = f (xj )

 

 · · · + an Ln (x0 )
a0 L0 (x0 ) + a1 L1 (x0 ) + a2 L2 (x0 ) + |{z} = f (x0 ) 

 a0 = f (x0 )

 | {z } | {z } | {z } | {z } 


 =1 =0 =0 =0 =0



 

a0 L0 (x1 ) + a1 L1 (x1 ) + a2 L2 (x1 ) + · · · + an Ln (x1 )
 = f (x1 ) a1 = f (x1 )

| {z } | {z } | {z } |{z}
=0
| {z } =⇒

 =0 =1 =0 =0 



 .. .. 

 ..


 . . 

 .
 
a0 L0 (xn ) + a1 L1 (xn ) + a2 L2 (xn ) + · · · + an Ln (xn )

= f (xn ) a n

= f (xn )
| {z } | {z } | {z } |{z} | {z }
=0 =0 =0 =0 =1

n
X
=⇒ Pn (x) = f (xi )Li (x)
i=0

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 11 / 22
Construction du polynôme Pn : base de Lagrange

Écrivons explicitement Pn (xj ) = f (xj ) :

a0 L0 (xj ) + a1 L1 (xj ) + a2 L2 (xj ) + · · · + an Ln (xj ) = f (xj )

 

 · · · + an Ln (x0 )
a0 L0 (x0 ) + a1 L1 (x0 ) + a2 L2 (x0 ) + |{z} = f (x0 ) 

 a0 = f (x0 )

 | {z } | {z } | {z } | {z } 


 =1 =0 =0 =0 =0



 

a0 L0 (x1 ) + a1 L1 (x1 ) + a2 L2 (x1 ) + · · · + an Ln (x1 )
 = f (x1 ) a1 = f (x1 )

| {z } | {z } | {z } |{z}
=0
| {z } =⇒

 =0 =1 =0 =0 



 .. .. 

 ..


 . . 

 .
 
a0 L0 (xn ) + a1 L1 (xn ) + a2 L2 (xn ) + · · · + an Ln (xn )

= f (xn ) a n

= f (xn )
| {z } | {z } | {z } |{z} | {z }
=0 =0 =0 =0 =1

n
X
=⇒ Pn (x) = f (xi )Li (x)
i=0

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 11 / 22
Construction du polynôme Pn : base de Lagrange

Exemple :
On souhaite trouver le polynôme passant par les points (1, 2), (2, 3), (4, 1).
La formule de Lagrange à 3 points est :

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

avec,

(x − x1 )(x − x2 ) (x − 2)(x − 4) (x − 2)(x − 4)

L0 (x) = = =
(x0 − x1 )(x0 − x2 ) (1 − 2)(1 − 4) 3
(x − x0 )(x − x2 ) (x − 1)(x − 4) (x − 1)(x − 4)
L1 (x) = = =
(x1 − x0 )(x1 − x2 ) (2 − 1)(2 − 4) −2
(x − x0 )(x − x1 ) (x − 1)(x − 2) (x − 1)(x − 2)
L2 (x) = = =
(x2 − x0 )(x2 − x1 ) (4 − 1)(4 − 2) 6

(x − 2)(x − 4) (x − 1)(x − 4) (x − 1)(x − 2)

=⇒ P3 (x) = 2 · −3· +1·
3 2 6

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 12 / 22
Construction du polynôme Pn : base de Lagrange
Exercice :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 13 / 22
Construction du polynôme Pn : base de Lagrange
Exercice :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points.  
x − xj
On sait que : Li (x) = nj=0
Q
j̸=i xi − xj

L’interpolation de Lagrange donne dans ce cas :

(x − 1)(x − 2)(x − 3) (x − 0)(x − 2)(x − 3)
P3 (x) = 1· +2·
(0 − 1)(0 − 2)(0 − 3) (1 − 0)(1 − 2)(1 − 3)

(x − 0)(x − 1)(x − 3) (x − 0)(x − 1)(x − 2)

+9 · + 28 ·
(2 − 0)(2 − 1)(2 − 3) (3 − 0)(3 − 1)(3 − 2)

(x − 1)(x − 2)(x − 3)
=⇒P3 (x) = − + x(x − 2)(x − 3)
6
x(x − 1)(x − 3) x(x − 1)(x − 2)
+9 · + 14 ·
2 3
=⇒P3 (x) = 1 + x 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 13 / 22
Construction du polynôme Pn : base de Lagrange
Exercice :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points.  
x − xj
On sait que : Li (x) = nj=0
Q
j̸=i xi − xj

L’interpolation de Lagrange donne dans ce cas :

(x − 1)(x − 2)(x − 3) (x − 0)(x − 2)(x − 3)
P3 (x) = 1· +2·
(0 − 1)(0 − 2)(0 − 3) (1 − 0)(1 − 2)(1 − 3)

(x − 0)(x − 1)(x − 3) (x − 0)(x − 1)(x − 2)

+9 · + 28 ·
(2 − 0)(2 − 1)(2 − 3) (3 − 0)(3 − 1)(3 − 2)

(x − 1)(x − 2)(x − 3)
=⇒P3 (x) = − + x(x − 2)(x − 3)
6
x(x − 1)(x − 3) x(x − 1)(x − 2)
+9 · + 14 ·
2 3
=⇒P3 (x) = 1 + x 3

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 13 / 22
Construction du polynôme Pn : base de Newton

La base de Newton est (1, N1 (x), N2 (x), · · · , Nn (x)).

Cette méthode consiste à approcher une fonction f (x) à l’aide d’un polynôme
d’interpolation construit avec les différences divisées de Newton.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 14 / 22
Construction du polynôme Pn : base de Newton

La base de Newton est (1, N1 (x), N2 (x), · · · , Nn (x)).

Cette méthode consiste à approcher une fonction f (x) à l’aide d’un polynôme
d’interpolation construit avec les différences divisées de Newton.
Le polynôme interpolateur de Newton est donné par :
n
X
Pn (x) = ai Ni (x)
i=0

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 14 / 22
Construction du polynôme Pn : base de Newton

La base de Newton est (1, N1 (x), N2 (x), · · · , Nn (x)).

Cette méthode consiste à approcher une fonction f (x) à l’aide d’un polynôme
d’interpolation construit avec les différences divisées de Newton.
Le polynôme interpolateur de Newton est donné par :
n
X
Pn (x) = ai Ni (x)
i=0

avec, Ni (x) est un polynôme de degré i, s’annulant en x0 , . . . , xi−1 :

i−1
Y
Ni (x) = (x − x0 )(x − x1 )(x − x2 ) · · · (x − xi−1 ) = (x − xj )
j=0

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 14 / 22
Construction du polynôme Pn : base de Newton

La base de Newton est (1, N1 (x), N2 (x), · · · , Nn (x)).

Cette méthode consiste à approcher une fonction f (x) à l’aide d’un polynôme
d’interpolation construit avec les différences divisées de Newton.
Le polynôme interpolateur de Newton est donné par :
n
X
Pn (x) = ai Ni (x)
i=0

avec, Ni (x) est un polynôme de degré i, s’annulant en x0 , . . . , xi−1 :

i−1
Y
Ni (x) = (x − x0 )(x − x1 )(x − x2 ) · · · (x − xi−1 ) = (x − xj )
j=0

Pn (x) = a0 N0 (x) + a1 N1 (x) + a2 N2 (x) + · · · + an Nn (x)

| {z }
=1

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 14 / 22
Construction du polynôme Pn : base de Newton

La base de Newton est (1, N1 (x), N2 (x), · · · , Nn (x)).

Cette méthode consiste à approcher une fonction f (x) à l’aide d’un polynôme
d’interpolation construit avec les différences divisées de Newton.
Le polynôme interpolateur de Newton est donné par :
n
X
Pn (x) = ai Ni (x)
i=0

avec, Ni (x) est un polynôme de degré i, s’annulant en x0 , . . . , xi−1 :

i−1
Y
Ni (x) = (x − x0 )(x − x1 )(x − x2 ) · · · (x − xi−1 ) = (x − xj )
j=0

Pn (x) = a0 N0 (x) + a1 N1 (x) + a2 N2 (x) + · · · + an Nn (x)

| {z }
=1

Pn (x) = a0 + a1 (x − x0 ) + a2 (x − x0 )(x − x1 ) + · · · + an (x − x0 )(x − x1 ) · · · (x − xn−1 )

| {z } | {z } | {z }
N1 (x) N2 (x) Nn (x)
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 14 / 22
Construction du polynôme Pn : base de Newton

La base de Newton est (1, N1 (x), N2 (x), · · · , Nn (x)).

Cette méthode consiste à approcher une fonction f (x) à l’aide d’un polynôme
d’interpolation construit avec les différences divisées de Newton.
Le polynôme interpolateur de Newton est donné par :
n
X
Pn (x) = ai Ni (x)
i=0

avec, Ni (x) est un polynôme de degré i, s’annulant en x0 , . . . , xi−1 :

i−1
Y
Ni (x) = (x − x0 )(x − x1 )(x − x2 ) · · · (x − xi−1 ) = (x − xj )
j=0

Pn (x) = a0 N0 (x) + a1 N1 (x) + a2 N2 (x) + · · · + an Nn (x)

| {z }
=1

Pn (x) = a0 + a1 (x − x0 ) + a2 (x − x0 )(x − x1 ) + · · · + an (x − x0 )(x − x1 ) · · · (x − xn−1 )

| {z } | {z } | {z }
N1 (x) N2 (x) Nn (x)
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 14 / 22
 Construction du polynôme Pn : base de Newton
Pn (xi ) = a0 +a1 (xi −x0 )+a2 (xi −x0 )(xi −x1 )+· · ·+an (xi −x0 )(xi −x1 ) · · · (xi −xn−1 ) = f (xi )



 P(x0 ) = a0 = f (x0 )

P(x ) = a + a (x − x ) = f (x1 )



 1 0 1 1 0

P(x2 ) = a0 + a1 (x2 − x0 ) + a2 (x2 − x0 )(x2 − x1 ) = f (x2 )

 .. ..
. .





P(xn ) = a0 + a1 (x − x0 ) + · · · + an (xn − x0 )(xn − x1 ) · · · (xn − xn−1 ) = f (xn )


Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 15 / 22
 Construction du polynôme Pn : base de Newton
Pn (xi ) = a0 +a1 (xi −x0 )+a2 (xi −x0 )(xi −x1 )+· · ·+an (xi −x0 )(xi −x1 ) · · · (xi −xn−1 ) = f (xi )



 P(x0 ) = a0 = f (x0 )

P(x ) = a + a (x − x ) = f (x1 )



 1 0 1 1 0

P(x2 ) = a0 + a1 (x2 − x0 ) + a2 (x2 − x0 )(x2 − x1 ) = f (x2 )

 .. ..
. .





P(xn ) = a0 + a1 (x − x0 ) + · · · + an (xn − x0 )(xn − x1 ) · · · (xn − xn−1 ) = f (xn )




 a0 = f (x0 ) = f [x0 ]
a1 = f (x1 ) − f (x0 ) = f [x0 , x1 ]



x1 − x0



f [x1 , x2 ] − f [x0 , x1 ]


a2 = = f [x0 , x1 , x2 ] =⇒ ai : les différences divisées
x2 − x0
.





 ..

an = f [x1 , x2 , · · · , xn ] − f [x0 , x1 , · · · , xn−1 ] = f [x0 , x1 , · · · , xn ]




xn − x0

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 15 / 22
 Construction du polynôme Pn : base de Newton
Pn (xi ) = a0 +a1 (xi −x0 )+a2 (xi −x0 )(xi −x1 )+· · ·+an (xi −x0 )(xi −x1 ) · · · (xi −xn−1 ) = f (xi )



 P(x0 ) = a0 = f (x0 )

P(x ) = a + a (x − x ) = f (x1 )



 1 0 1 1 0

P(x2 ) = a0 + a1 (x2 − x0 ) + a2 (x2 − x0 )(x2 − x1 ) = f (x2 )

 .. ..
. .





P(xn ) = a0 + a1 (x − x0 ) + · · · + an (xn − x0 )(xn − x1 ) · · · (xn − xn−1 ) = f (xn )




 a0 = f (x0 ) = f [x0 ]
a1 = f (x1 ) − f (x0 ) = f [x0 , x1 ]



x1 − x0



f [x1 , x2 ] − f [x0 , x1 ]


a2 = = f [x0 , x1 , x2 ] =⇒ ai : les différences divisées
x2 − x0
.





 ..

an = f [x1 , x2 , · · · , xn ] − f [x0 , x1 , · · · , xn−1 ] = f [x0 , x1 , · · · , xn ]




xn − x0

=⇒Pn (x) = f [x0 ] + f [x0 , x1 ](x − x0 ) + f [x0 , x1 , x2 ](x − x0 )(x − x1 ) + · · ·

| {z } | {z }
N1 N2
+f [x0 , · · · , xn ](x − x0 )(x − x1 ) · · · (x − xn−1 ) = f (x)
| {z }
Nn
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 15 / 22
 Construction du polynôme Pn : base de Newton
Pn (xi ) = a0 +a1 (xi −x0 )+a2 (xi −x0 )(xi −x1 )+· · ·+an (xi −x0 )(xi −x1 ) · · · (xi −xn−1 ) = f (xi )



 P(x0 ) = a0 = f (x0 )

P(x ) = a + a (x − x ) = f (x1 )



 1 0 1 1 0

P(x2 ) = a0 + a1 (x2 − x0 ) + a2 (x2 − x0 )(x2 − x1 ) = f (x2 )

 .. ..
. .





P(xn ) = a0 + a1 (x − x0 ) + · · · + an (xn − x0 )(xn − x1 ) · · · (xn − xn−1 ) = f (xn )




 a0 = f (x0 ) = f [x0 ]
a1 = f (x1 ) − f (x0 ) = f [x0 , x1 ]



x1 − x0



f [x1 , x2 ] − f [x0 , x1 ]


a2 = = f [x0 , x1 , x2 ] =⇒ ai : les différences divisées
x2 − x0
.





 ..

an = f [x1 , x2 , · · · , xn ] − f [x0 , x1 , · · · , xn−1 ] = f [x0 , x1 , · · · , xn ]




xn − x0

=⇒Pn (x) = f [x0 ] + f [x0 , x1 ](x − x0 ) + f [x0 , x1 , x2 ](x − x0 )(x − x1 ) + · · ·

| {z } | {z }
N1 N2
+f [x0 , · · · , xn ](x − x0 )(x − x1 ) · · · (x − xn−1 ) = f (x)
| {z }
Nn
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 15 / 22
 Construction du polynôme Pn : base de Newton
Définition
On définit les premières différences divisées de la fonction f (x) par :

f (xi+1 ) − f (xi )
f [xi , xi+1 ] =
xi+1 − xi
a1 = f [x0 , x1 ]
Les deuxièmes différences divisées de f (x) sont définies à partir des premières
différences divisées par :
f [xi+1 , xi+2 ] − f [xi , xi+1 ]
f [xi , xi+1 , xi+2 ] =
xi+2 − xi
a2 = f [x0 , x1 , x2 ]

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 16 / 22
 Construction du polynôme Pn : base de Newton
Définition
On définit les premières différences divisées de la fonction f (x) par :

f (xi+1 ) − f (xi )
f [xi , xi+1 ] =
xi+1 − xi
a1 = f [x0 , x1 ]
Les deuxièmes différences divisées de f (x) sont définies à partir des premières
différences divisées par :
f [xi+1 , xi+2 ] − f [xi , xi+1 ]
f [xi , xi+1 , xi+2 ] =
xi+2 − xi
a2 = f [x0 , x1 , x2 ]
Les n-ièmes différences divisées de f (x) sont données par :

f [x1 , x2 , . . . , xn ] − f [x0 , x1 , . . . , xn−1 ]

f [x0 , x1 , . . . , xn ] =
xn − x0
an = f [x0 , x1 , · · · , xn ]

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 16 / 22
 Construction du polynôme Pn : base de Newton
Définition
On définit les premières différences divisées de la fonction f (x) par :

f (xi+1 ) − f (xi )
f [xi , xi+1 ] =
xi+1 − xi
a1 = f [x0 , x1 ]
Les deuxièmes différences divisées de f (x) sont définies à partir des premières
différences divisées par :
f [xi+1 , xi+2 ] − f [xi , xi+1 ]
f [xi , xi+1 , xi+2 ] =
xi+2 − xi
a2 = f [x0 , x1 , x2 ]
Les n-ièmes différences divisées de f (x) sont données par :

f [x1 , x2 , . . . , xn ] − f [x0 , x1 , . . . , xn−1 ]

f [x0 , x1 , . . . , xn ] =
xn − x0
an = f [x0 , x1 , · · · , xn ]

Remarque :
Les 0ème différences divisées sont simplement les valeurs f (xi ).
Chaque niveau utilise les résultats du niveau précédent.
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 16 / 22
 Construction du polynôme Pn : base de Newton
Définition
On définit les premières différences divisées de la fonction f (x) par :

f (xi+1 ) − f (xi )
f [xi , xi+1 ] =
xi+1 − xi
a1 = f [x0 , x1 ]
Les deuxièmes différences divisées de f (x) sont définies à partir des premières
différences divisées par :
f [xi+1 , xi+2 ] − f [xi , xi+1 ]
f [xi , xi+1 , xi+2 ] =
xi+2 − xi
a2 = f [x0 , x1 , x2 ]
Les n-ièmes différences divisées de f (x) sont données par :

f [x1 , x2 , . . . , xn ] − f [x0 , x1 , . . . , xn−1 ]

f [x0 , x1 , . . . , xn ] =
xn − x0
an = f [x0 , x1 , · · · , xn ]

Remarque :
Les 0ème différences divisées sont simplement les valeurs f (xi ).
Chaque niveau utilise les résultats du niveau précédent.
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 16 / 22
Construction du polynôme Pn : base de Newton

Calcul pratique des coefficients ai :

On construit un tableau triangulaire basé sur les nœuds (xi ) et les valeurs associées
f (xi ).

xi f [xi ] f [xi , xi+1 ] f [xi , xi+1 , xi+2 ] f [xi , xi+1 , xi+2 , xi+3 ]
x0 f [x0 ]
f [x0 , x1 ]
x1 f [x1 ] f [x0 , x1 , x2 ]
f [x1 , x2 ] f [x0 x1 , x2 , x3 ]
x2 f [x2 ] f [x1 , x2 , x3 ]
f [x2 , x3 ]
x3 f [x3 ]

Remarque :
Le polynôme de Newton utilise uniquement les termes situés sur la première ligne de
chaque colonne.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 17 / 22
Construction du polynôme Pn : base de Newton

Théorème
L’unique polynôme de degré n passant
par les (n + 1) points d’interpolation
(xi , f (xi )) pour i = 0, 1, 2, · · · , n peut s’écrire sous la forme récursive :

Pn (x) = Pn−1 (x) + an (x − x0 )(x − x1 ) · · · (x − xn−1 ) (1)

Les coefficients de ce polynôme sont les différences divisées :

an = f [x0 , x1 , x2 , · · · , xn ] pour 0 ≤ i ≤ n (2)

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 18 / 22
Construction du polynôme Pn : base de Newton

Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points.

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 19 / 22
Construction du polynôme Pn : base de Newton

Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points.

xi f [xi ] f [xi , xi+1 ] f [xi , xi+1 , xi+2 ] f [xi , xi+1 , xi+2 , xi+3 ]
0 1
1
1 2 3
7 1
2 9 6
19
3 28

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 19 / 22
Construction du polynôme Pn : base de Newton

Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points.

xi f [xi ] f [xi , xi+1 ] f [xi , xi+1 , xi+2 ] f [xi , xi+1 , xi+2 , xi+3 ]
0 1
1
1 2 3
7 1
2 9 6
19
3 28

=⇒ P3 (x) = 1 + 1(x − 0) + 3(x − 0)(x − 1) + 1(x − 0)(x − 1)(x − 2)

=⇒ P3 (x) = x 3 + 1

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 19 / 22
Construction du polynôme Pn : base de Newton

Exercice 1 :
On souhaite trouver le polynôme passant par les points (0, 1), (1, 2), (2, 9), (3, 28).
c-à-d, étant donné 4 points, on va chercher le polynôme P3 de degré 3 qui interpole ces
points.

xi f [xi ] f [xi , xi+1 ] f [xi , xi+1 , xi+2 ] f [xi , xi+1 , xi+2 , xi+3 ]
0 1
1
1 2 3
7 1
2 9 6
19
3 28

=⇒ P3 (x) = 1 + 1(x − 0) + 3(x − 0)(x − 1) + 1(x − 0)(x − 1)(x − 2)

=⇒ P3 (x) = x 3 + 1

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 19 / 22
Construction du polynôme Pn : base de Newton
Exercice 2 :
Si on souhaite ajouter un point d’interpolation, et calculer le polynôme de degré 4
passant par les points (0, 1), (1, 2), (2, 9), (3, 28), (5, 54).

xi f [xi ] f [xi , xi+1 ] f [xi , xi+1 , xi+2 ] f [xi , xi+1 , xi+2 , xi+3 ] f [xi , xi+1 , xi+2 , xi+3 , xi+4 ]
0 1
1
1 2 3
7 1 3
2 9 6 −
19 −2 5
3 28 −2
13
5 54

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 20 / 22
Construction du polynôme Pn : base de Newton
Exercice 2 :
Si on souhaite ajouter un point d’interpolation, et calculer le polynôme de degré 4
passant par les points (0, 1), (1, 2), (2, 9), (3, 28), (5, 54).

xi f [xi ] f [xi , xi+1 ] f [xi , xi+1 , xi+2 ] f [xi , xi+1 , xi+2 , xi+3 ] f [xi , xi+1 , xi+2 , xi+3 , xi+4 ]
0 1
1
1 2 3
7 1 3
2 9 6 −
19 −2 5
3 28 −2
13
5 54

On sait que :
Pn (x) = Pn−1 (x) + an (x − x0 )(x − x1 ) · · · (x − xn−1 )
Alors,
P4 (x) = P3 (x) + f [x0 , x1 , x2 , x3 , x4 ](x − 0)(x − 1)(x − 2)(x − 3)

3
=⇒ P4 (x) = P3 (x) − (x − 0)(x − 1)(x − 2)(x − 3)
5
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 20 / 22
Construction du polynôme Pn : base de Newton
Exercice 2 :
Si on souhaite ajouter un point d’interpolation, et calculer le polynôme de degré 4
passant par les points (0, 1), (1, 2), (2, 9), (3, 28), (5, 54).

xi f [xi ] f [xi , xi+1 ] f [xi , xi+1 , xi+2 ] f [xi , xi+1 , xi+2 , xi+3 ] f [xi , xi+1 , xi+2 , xi+3 , xi+4 ]
0 1
1
1 2 3
7 1 3
2 9 6 −
19 −2 5
3 28 −2
13
5 54

On sait que :
Pn (x) = Pn−1 (x) + an (x − x0 )(x − x1 ) · · · (x − xn−1 )
Alors,
P4 (x) = P3 (x) + f [x0 , x1 , x2 , x3 , x4 ](x − 0)(x − 1)(x − 2)(x − 3)

3
=⇒ P4 (x) = P3 (x) − (x − 0)(x − 1)(x − 2)(x − 3)
5
Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 20 / 22
 Erreur d’interpolation

Théorème
Soit f ∈ C n+1 ([a, b]) une fonction (n + 1)-fois continûment dérivable sur un intervalle [a, b],
et soient x0 , x1 , . . . , xn des points distincts de cet intervalle.
Alors, pour tout x ∈ [a, b], il existe un ξ ∈ [a, b] tel que :

f (n+1) (ξ)
En (x) = f (x) − Pn (x) = (x − x0 )(x − x1 ) · · · (x − xn )
(n + 1)!

où Pn (x) est le polynôme interpolateur de Newton (ou de Lagrange) de degré n.

On définit :
Mn+1 = max f (n+1) (x)
a≤x≤b

Alors l’erreur est bornée par :

Mn+1
En (x) ≤ (x − x0 )(x − x1 ) · · · (x − xn )
(n + 1)!

Remarque :
En (xi ) = 0 pour i = 0, 1, 2, · · · , n

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 21 / 22
 Erreur d’interpolation
Exemple :

Soit f (x) = cos(x) définie sur l’intervalle [0, π2 ]. On se donne les points d’interpolation
suivants x0 = 0 et x1 = π2 . Le point évalué est x = π4 .

Formule de l’erreur d’interpolation :

f (2) (x)
E1 (x) = (x − x0 ) (x − x1 )
2!
Application numérique :
  
M2 π π π
E1 ≤ −0 −
2! 4 4 2

Évaluation de M2 :

M2 = maxπ f (2) (x) = maxπ − cos(x) ⇒ M2 = 1

0≤x≤ 2 0≤x≤ 2

π2
=⇒ E1 ≤ ≈ 0.308
32

Pr. Ikrame DAQAQ Interpolation et approximation polynômiale de fonctions November 25, 2025 22 / 22