Cours Chapitre 3 Interpolation Polynomiale

CHAPITRE 3 : INTERPOLATION POLYNÔMIALE
SLIMANE HADJ-BRAHIM(1)
(1) USTHB University, Faculty of Mathematics, RECITS Laboratory

BP 32, El-Alia, 16111 Bab-Ezzouar, Algiers, Algeria
3.2.2 Méthode de Newton
Définition 1. Soit f une fonction dont on connait les valeurs f (x0 ), f (x1 ), . . . , f (xn )
qu’elle prend aux abscisses distinctes x0 , x1 , . . . , xn . On définit les différences di-
visées de f aux points x0 , x1 , . . . , xn par les relations de récurrence :


 δ(xi ) = f (xi )
 f (xi )−f (xi+1 )
= δ(xxi )−δ(x i+1 )



 δ(xi , xi+1 ) = xi −xi+1 i −xi+1

δ(xi ,xi+1 )−δ(xi+1 ,xi+2 )
δ(xi , xi+1 , xi+2 ) = xi −xi+2

 ..



 .

 δ(x , x , . . . , x ) δ(xi ,...,xi+p−1 )−δ(xi+1 ,...,xi+p )
i i+1 i+p = xi −xi+p
.
La dernière relation du système précédent est appelée : différence divisée d’ordre p

de la fonction f aux points xi , xi+1 , . . . , xi+p .
Proposition 1. La différence divisée δ(xi , xi+1 , . . . , xi+p ) est invariante par toute
permutation des abscisses xi , xi+1 , . . . , xi+p .
Calcul des différences divisées. Pour calculer la différence divisée d’ordre n de

la fonction f aux points x0 , x1 , . . . , xn , on forme le tableau suivant en appliquant les
formules du du système précédent colonne aprés colonne :
Adresse e-mail : [email protected].

1
2 COURS DE MONSIEUR S. HADJ-BRAHIM : GÉNIE CIVIL ”D”
xi f (xi ) δ (1) δ (2) δ (3) ··· δ (n−1) δ (n)

x0 f (x0 ) = δ0
δ0,1
x1 f (x1 ) = δ1 δ0,1,2
δ1,2 δ0,1,2,3
x2 f (x2 ) = δ2 δ1,2,3
δ2,3 δ1,2,3,4
x3 f (x3 ) = δ3 δ2,3,4
δ3,4
x4 f (x4 ) = δ4
.. ..
. . δ0,1,...,n−1
δ0,1,...,n
.. ..
. . δ1,2,...,n
xn−2 f (xn−2 ) = δn−2 δn−3,n−2,n−1
δn−2,n−1 δn−3,n−2,n−1,n
xn−1 f (xn−1 ) = δn−1 δn−2,n−1,n
δn−1,n
xn f (xn ) = δn
Théorème 1. Le polynôme Pn (x) qui prend les valeurs f (xi ), i = 0, 1, 2, . . . , n, aux

points distincts xi , i = 0, 1, 2, . . . , n, peut s’écrire :
Pn (x) = δ0 + δ0,1 .(x − x0 ) + δ0,1,2 .(x − x0 ).(x − x1 ) + δ0,1,2,3 .(x − x0 ).(x − x1 ).(x − x2 )
+ · · · + δ0,1,2,...,n .(x − x0 ).(x − x1 ).(x − x2 ) · · · (x − xn−1 ).
Le polynôme Pn (x) est appelé polynôme d’interpolation de Newton pour les différences
divisés.
Théorème 2. (Relation entre différences divisées et dérivées)

Soit [a, b] un intervalle contenant les points x0 , x1 , . . . , xn , et soit f ∈ C n−1 ([a, b]).
On suppose que la n-ième dérivée de f existe sur ]a, b[. Alors il existe α ∈]a, b[ tel
que :
f (n) (α)
δ0,1,2,...,n = .
n!
3.2.3 Erreur d’interpolation
Dans la pratique, l’interpolation polynômiale sert à remplacer une fonction f qui
est, soit inconnue, soit trop compliquée, par une fonction plus simple, en l’ocurrence
CHAPITRE 3 : INTERPOLATION POLYNÔMIALE 3
un polynôme. On dit que l’on approxime f par le polynôme d’interpolation Pn (x).

Quand on utilise une approximaion, comme c’est le cas dans de nombreuses méthodes
d’analyse numériques, il est fondamental d’étudier l’erreur d’approximation. Na-
turelement, sauf cas particulier, l’expression de l’erreur ne permet pas son calcul
exact; elle peut cependant être très utile pour en calculer une majoration. C’est
ainsi que pour l’interpolation polynômiale on démontre le résulat ci-après.
Théorème 3. Soit [a, b] un intervalle contenant les points x0 , x1 , . . . , xn , et soit

f ∈ C n ([a, b]). On suppose que la (n+1)-ième dérivée de f existe sur ]a, b[. Alors
il existe α ∈]a, b[ tel que :
n
f (n+1) (α) Y
f (x) − Pn (x) = (x − xi ).
(n + 1)! i=0
La formule précédente ne permet pas le calcul de valeur exacte de l’erreur parce

que, en général, α inconnu. Elle permet, par contre, d’en déduire une majoration
(et permet meme de choisir les x0 , x1 , . . . , xn de façon optimale lorsque ceux-ci ne
sont pas imposés). D’où le résultat suivant :
Corollaire 1. Sous les hpothèses du Théorème 3, on a :

n
Mn+1 Y
|f (x) − Pn (x)| ≤ (x − xi ) .

(n + 1)! i=0

où

Mn+1 = max f (n+1) (x) .
x∈[a,b]
Exercice 1. Construire le polynôme d’interpolation de Newton associé à la fonction

tabulée suivante :
x 4 6 8 10
f (x) 1 3 8 20
Solution. I Considérons le tableau correspondant aux différences divisées d’ordre

3 de la fonction f aux points x0 , x1 , x2 , x3 :
4 COURS DE MONSIEUR S. HADJ-BRAHIM : GÉNIE CIVIL ”D”
xi f (xi ) δ (1) δ (2) δ (3)

x0 = 4 f (x0 ) = 1
δ0,1 = 1
x1 = 6 f (x1 ) = 3 δ0,1,2 = 3/8
δ1,2 = 5/2 δ0,1,2,3 = 1/12
x2 = 8 f (x2 ) = 8 δ1,2,3 = 7/8
δ2,3 = 6
x3 = 10 f (x3 ) = 20
La formule de Newton pour les différences divisées associées à f s’écrit alors :
P3 (x) = δ0 + δ0,1 .(x − x0 ) + δ0,1,2 .(x − x0 ).(x − x1 ) + δ0,1,2,3 .(x − x0 ).(x − x1 ).(x − x2 )
3 1
= 1 + 1(x − 4) + (x − 4).(x − 6) + (x − 4).(x − 6).(x − 8)
8 12
71 9 2 1 3
= −10 + x − x + x .
12 8 12
√
Exercice 2. Soit la fonction f (x) = x calculée aux points x0 = 100, x1 = 144 et
x2 = 169.
(1) Construire le polynôme de Lagrange correspodant.
(2) Calculer f (116).
(3) Estimer l’erreur correspondante. En déduire le nombre décimalee exactes.
Solution. I
1. Le polynôme de Lagrange s’écrit :
2
X
P2 (x) = f (xi )Li (x).
i=0
Déterminons les Li (x), considérons le tableau suivant :

x − 100 −44 −69
44 x − 144 −25
69 25 x − 169
On en déduit que :
(x − 144)(x − 169) 1
x2 − 313x + 24336 ,

L0 = =
(−44)(−69) 3036
(x − 100)(x − 169) 1
x2 − 269x + 16900 ,

L1 = =−
(44)(−25) 1100
(x − 100)(x − 144) 1
x2 − 244x + 14400 .

L2 = =
(69)(25) 1725
CHAPITRE 3 : INTERPOLATION POLYNÔMIALE 5
D’où
P2 (x) = −0.0000789x2 + 3.7423514x + 4.313088.

2. Posons x = 116 dans P2 (x). Alors
f (116) ' P2 (116) = 10.7804.
3. On sait que l’erreur est estimée d’après la formule

2
M3 Y
E(x) = |f (x) − P2 (x)| ≤ (x − xi ) .

3! i=0
(3)
où M3 = maxx∈[100,169] f (x). Comme, pour tout x ∈ [100, 169], on a
3 1 15 1
f (3) (x) = 2
√ > 0 et f (4) (x) = − √ < 0.
8x x 16 x3 x
Donc f (3) (x) est une fonction décroissante sur [100, 169], on déduit que
3
M = f (3) (100) = 10−5 .
8
D’où
3 −5
E(116) = |f (166) − P2 (166)| ≤ 10 |(116 − 100)(116 − 144)(116 − 169)| ' 0.1510−1 ≤ 0.510−1 .
3!8
Conséquence : f (166) = 10.8 ± 0.510−1 .

Cours Chapitre 3 Interpolation Polynomiale

Transféré par

Droits d'auteur :

Formats disponibles

Cours Chapitre 3 Interpolation Polynomiale

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Cours Chapitre 3 Interpolation Polynomiale

Transféré par

Droits d'auteur :

Formats disponibles

CHAPITRE 3 : INTERPOLATION POLYNÔMIALE

(1) USTHB University, Faculty of Mathematics, RECITS Laboratory

3.2.2 Méthode de Newton

La dernière relation du système précédent est appelée : différence divisée d’ordre p

Calcul des différences divisées. Pour calculer la différence divisée d’ordre n de

Adresse e-mail : [email protected].

xi f (xi ) δ (1) δ (2) δ (3) ··· δ (n−1) δ (n)

Théorème 1. Le polynôme Pn (x) qui prend les valeurs f (xi ), i = 0, 1, 2, . . . , n, aux

Théorème 2. (Relation entre différences divisées et dérivées)

un polynôme. On dit que l’on approxime f par le polynôme d’interpolation Pn (x).

Théorème 3. Soit [a, b] un intervalle contenant les points x0 , x1 , . . . , xn , et soit

La formule précédente ne permet pas le calcul de valeur exacte de l’erreur parce

Corollaire 1. Sous les hpothèses du Théorème 3, on a :

Exercice 1. Construire le polynôme d’interpolation de Newton associé à la fonction

Solution. I Considérons le tableau correspondant aux différences divisées d’ordre

xi f (xi ) δ (1) δ (2) δ (3)

Déterminons les Li (x), considérons le tableau suivant :

P2 (x) = −0.0000789x2 + 3.7423514x + 4.313088.

f (116) ' P2 (116) = 10.7804.

3. On sait que l’erreur est estimée d’après la formule

Conséquence : f (166) = 10.8 ± 0.510−1 .

Vous aimerez peut-être aussi