École Supérieure Privée d’Ingénierie et de Technologies

Analyse Numérique
Correction série d’exercices No 2 Interpolation et approximation polynomiale

Niveau : 3ème année Année universitaire : 2022-2023

Exercice 1 On considére les points (−2, 4) ; (0, 0) ;(1, 0) et (2, 4). Parmi les polynômes suivants,
lequel est le polynôme d’interpolation P aux quatre points et justifier votre réponse.
(1) P1 (X) = X 4 + 32 X 3 + 3X 2 + 38 X
(2) P2 (X) = 43 X 2 − 4
3
(3) P3 (X) = 31 X 3 + X 2 − 43 X
(4) P4 (X) = 61 X 4 + X 3 + 23 X 2 + X

Correction :
On ne demande pas ici de calculer le polynôme mais de l’identifier, on va donc utiliser la ca-
ractérisation du polynôme d’interpolation de Lagrange associé aux points.


 P(-2)=4 ;
P(0)=0 ;

P polynôme d’interpolation de Lagrange aoocié à xi ⇔ deg(P ) ≤ 3 et

 P(1)=0 ;
P(2)=4.

Il n’y a qu’à trouver le polynôme qui satisfait toute les propriètés.
Existance et unicité du polynôme :

- Le polynôme P1 est de degré 4 donc éliminé

- Le polynôme P2 a un terme constant non nul il ne s’annule pas en 0 donc éliminé
- Le polynôme P3 on vérifie qu’il convient et P4 ne vérifie pas P (1) = 0

Exercice 2 (Examen Mai 2019)

Partie I : Interpolation polynomial

(a) Justifier l’existence d’un unique polynôme P2 ∈ R2 [X] interpolant les points (−2, 16), (0, −4)
et (2, 8).
(b) Déterminer l’expression du polynôme P2 par une méthode (vue en cours) de votre choix

1
 Partie II : Approximation au sens des moindres carrées
Dans l’objectif d’étudier le chemin de freinage d’un véhicule, correspondant à la distance
parcourue en mètres (m) du début du freinage jusqu’à l’arrêt total du véhicule, en fonction
de la vitesse en Kilomètres par heure (Km/h) de ce dernier, 12 expériences indépendantes
ont été réalisées. Les résultats obtenus sont présentés dans le tableau ci-dessous. On note
par X = (xi )1≤i≤12 et Y = (yi )1≤i≤12 , où xi , et yi , désignent, respectivement, la vitesse du
véhicule et le chemin de freinage associés à l’éxpérience i.

i 1 2 3 4 5 6 7 8 9 10 11 12
xi 40 50 60 70 80 90 100 110 120 130 140 150
yi 9 11 20 27 39 45 58 78 79 93 108 12


a
(a) Déterminer les coefficients Z = de la droite f (t, Z) = a + bt, qui ajuste au mieux
b
les points (xi ; yi )1≤i≤12 au sens des moindres carrées. On donne les valeurs des sommes
suivantes :
12
X 12
X 12
X 12
X
xi = 1140; x2i = 122600; yi = 691; xi yi = 80840
i=1 i=1 i=1 i=1

(b) Rouler à une vitesse de 105Km/h, le conducteur de ce véhicule pourrait-il éviter un

obstacle survenant à une distance de 60m ? Justifir votre réponse.

Correction :
Partie I : Interpolation polynomiale
(a) Les abscisses des points (−2, 16); (0, −4) et (2, 8) sont deux à deux distincts donc il existe
un unique polynôme P2 ∈ R2 [X] passant par ces points. un unique polynôme P2 ∈ R2 [X]
passant par ces points.
(b) Méthode de Lagrange :
On considère les polynôme (Li )0≤i≤2 de Lagrange associés aux (−2, 16); (0, −4) et (2, 8).

x(x − 2) x2 − 4 x(x + 2)
L0 (x) = , L1 (x) = , L2 (x) =
8 −4 8

Alors P2 (x) = 16L0(x) − 4L1 (x) + 8L2 (x) = 4x2 − 2x − 4

Méthode de Newton :
Le polynôme de Newton est donné par :

P2 (x) = α0 w0 (x) + α1 w1 (x) + α2 w2 (x)

avec 
 w0 (x) = 1
w1 (x) = x + 2
w2 (x) = x(x + 2)


Détermination des coefficients α0 , α1 et α2 par la méthode des differences divisées : On a :

x0 = −2, y0 = 16
x1 = 0, y1 = −4

2
 x2 = 2, y2 = 8
f [x1 ,x2 ]−f [x0 ,x1 ]
alors α0 = 16 = f [x0 ], α1 = f [x0, x1] = −10, α2 = x2 −x0 =4
Ainsi α0 = 16, α1 = −10 et α2 = 4 d’où :

P2 (x) = 16w0 (x) − 10w1 (x) + 4w2 (x) = 4x2 − 2x − 4

Partie II : Approximation
  au sens des moindres carrées
a
(a) Le vecteur Z = de la droite f (t, Z) = a + bt, qui ajuste au mieux les points
b
(xi ; yi )1≤i≤12 au sens des moindres carrées
12 
X 2
F (Z, X) = f (xi , Z) − yi
i=1
 
1 x1
 1 x2 
Il est donné par la relation suivante : Z ∗ = (t AA)−1t AY avec A =   On a :
 
.. ..
 . . 
1 x12

12
 
X
 12 xi   
t
 i=1
 t 12 1140
AA =  12 12
 = AA =
 X X  1140 122600
xi x2i
 
i=1 i=1

Cherchons (t AA)−1 ?
t
On sait que (t AA)−1 = det(t1AA) com(t AA) avec det(t AA) = 12 × 122600 − (1140)2 = 171600
   
t 122600 −1140 t −1 1 122600 −1140
et com( AA) = alors ( AA) = 171600 D’autre
−1140 12 −1140 12
 
691
part t AY = Ainsi
80840
 
∗ t −1t −43, 36
Z = ( AA) AY = Donc f (t, Z) = −43, 36 + 1.06t
1, 06

(b) Une valeur estimée du chemin de freinage du véhicule à une vitesse de 105km/h est donné
par f (10.5; Z ∗ ) = −43.36 + 1.06 × 10.5 = 68.21. Le conducteur du véhicule ne pourra pas
éviter l’obstacle.

Exercice 3 (1) Construire le polynôme P d’interpolation de Lagrange aux points (−1, e); (0, 1)
et (1, e).
(2) Sans faire de calcul, donner l’expression du polynôme de Lagrange Q qui interpole les trois
points (−1, −1); (0, 0) et (1, −1).
(3) Trouver le polynôme de l’espace vectoriel Vect(1, X, X 2 ) qui interpole les trois points (−1, −1); (0, 0)
et (1, −1).

3
Correction :

1) Le polynôme P d’interpolation de Lagrange de degré n qui interpole (n+1) points {(xi , yi ); i =

0, ..., n} s’écrit :
n n
X Y x − xj
Pn (x) = yi
xi − xj
i=0 j=0,j̸=i

ici n = 2 donc on a :

(x − x1 )(x − x2 ) (x − x0 )(x − x2 ) (x − x0 )(x − x1 )

P (x) = y0 + y1 + y2
(x0 − x1 )(x0 − x2 ) (x1 − x0 )(x1 − x2 ) (x2 − x0 )(x2 − x1 )
ex(x − 1) e(x + 1)x
= − (x + 1)(x − 1) + = (e − 1)x2 + 1
2 2

2) Il suffit de changer les coefficients yi dans l’expression précedente :

x(x − 1) (x + 1)x
Q(x) = − − = −x2
2 2

3) Il s’agit de trouver un polynôme P (x) qui soit combinaison linèaire de deux polynômes as-
sigès (ie : P (x) = α + βx + γx2 ) et qui interpole les 3 points (−1, −1); (0, 0) et (1, −1).


 P(-1)=1 ;
P(0)=0 ;
P(1)=-1 ;


ce qui donne α = 0; β = 0 et γ = −1
le polynôme cherché est donc P (x) = −x2
1
Exercice 4 Soit la fonction f (x) = 1+x2
(1) Déterminer le polynôme d’interpolation de Newton aux points −2, −1, 0 et 1.
(2) Donner la valeur approchée de f au point x = 0.5
(3) Calculer l’erreur à ce point.
(4) Estimer l’erreur sur l’intervalle [−2, 1].

Correction :

1) En utilisant la méthode de Newton,

P3 (x) = β0 + β1 (x − x0 ) + β2 (x − x0 )(x − x1 ) + β3 (x − x0 )(x − x1 )(x − x2 ),

avec
β0 = y0 = f (x0 ) = 0.2
y1 −y0 f (x1 )−f (x0 )
β1 = [y0 , y1 ] = x1 −x0 = x1 −x0 = 0.3

4
 f (x2 )−f (x1 )
[y1 ,y2 ]−[y0 ,y1 ] x2 −x1
−0.3
β2 = [y0 , y1 , y2 ] = x2 −x0 = x2 −x0 = 0.25
[y2 ,y3]−[y1 ,y2 ]
[y1 ,y2 ,y3 ]−[y0 ,y1 ,y2 ] x3 −x1
−0.25
β3 = [y0 , y1 , y2 , y3 ] = x3 −x0 = x3 −x0 = −0.25
d’où

P3 (x) = β0 + β1 (x + 2) + β2 (x + 2)(x + 1) + β3 (x + 2)(x + 1)x

= 0.2 + 0.3(x + 2) + 0.1(x + 1)(x + 2) − 0.2x(x + 1)(x + 2)
= −0.2x3 − 0.5x2 + 0.2x + 1
2) P (0.5) = 0.95 ≃ f (0.5)

3) E(0.5) =| f (0.5) − P (0.5) |=| 0.8 − 0.95 |= 0.015

4) La fonction f est de Classe C 4 sur [−2, 1] et on a :

−2x
f ′ (x) = (1+x2 )2
.

6x2 −2
f (2) (x) = (1+x2 )3
.
24x(1−x2 )
f (3) (x) = (1+x2 )4
.
5x4 −10x2 +1
f (4) (x) = 24 (1+x2 )5 .
on cherche maintenant à maximiser f (4) sur [−2, 1]. On pose le changement de variable
X = x2 , ce qui revient à maximiser la fonction :
5X 2 − 10X + 1
F (X) = 24 sur [0, 4]
(1 + X)5
et on a
−15x2 + 50X − 15
F ′ (X) =
(1 + X)6
qui s’annule en 1/3 et 3, voir le TVA
1
X 0 3 3 4

F ′ (X) − 0 + 0 −
3
24 8
F (X)
− 81
8
984
3125

On a alors pour tout x ∈ [−2, 1] :

f (4) (t)
E(x) ≤ sup | || x − x0 || x − x1 || x − x2 | ≤ | x + 2 || x + 1 || x || x − 1 |
t∈[−2,1] 4!