Université Sidi Mohamed Ben Abdellah Année 2023-2024
Faculté Polydisciplinaire de Taza Filières : SMA (S4)
Travaux Dirigés : Analyse numérique
suite de correction de TD n◦ 4
Solution de l’exercice 3 :
1. Le polynôme d’interpolation de Lagrange de degré n sur l’ensemble des n + 1 points
{(xi , yi )}ni=0 s’écrit
n n
X Y x − xj
pn (x) = y i
i=0
xi − xj
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 )
x(x − 1) (x + 1)x
=e − (x + 1)(x − 1) + e =
2 2
= (e − 1)x2 + 1
2. Il suffit de changer les coefficients yi dans l’expression précédente :
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 1, x et x2 (i.e.
p(x) = α + βx + γx2 ) et qui interpole les trois points (−1, −1), (0, 0) et (1, −1) :
p(−1) = 1, α − β + γ = −1
p(0) = 0, ⇔ α=0
p(1) = −1, α + β + γ = −1
d’où α = 0, β = 0 et γ = −1. Le polynôme cherché est donc le polynôme p(x) = −x2 .
Solution de l’exercice 4 :
Le polynôme d’interpolation de Lagrange de degré n sur l’ensemble des n+1 points {(xi , yi )}ni=0
s’écrit
n n
X Y x − xj
pn (x) = yi
i=0
xi − xj
j=0
j̸=i
1. Ici n = 3 donc on a
(x − x1 ) (x − x2 ) (x − x3 ) (x − x0 ) (x − x2 ) (x − x3 )
P (x) = y0 + y1
(x0 − x1 ) (x0 − x2 ) (x0 − x3 ) (x1 − x0 ) (x1 − x2 ) (x1 − x3 )
(x − x0 ) (x − x1 ) (x − x3 ) (x − x0 ) (x − x1 ) (x − x2 )
+ y2 + y3 =
(x2 − x0 ) (x2 − x1 ) (x2 − x3 ) (x3 − x0 ) (x3 − x1 ) (x3 − x2 )
x(x − 1)(x − 2) (x + 1)(x − 1)(x − 2) (x + 1)x(x − 1)
= + − (x + 1)x(x − 2) + =
−6 2 2
1 1 2
= − x3 + x2 + x + 1.
6 2 3
1
2. Par construction
Q(−1) = P (−1)
Q(0) = P (0)
Q(1) = P (1)
donc le polynôme Q(x) − P (x) s’annule en −1, en 0 et en 1 , ceci signifie qu’il existe un
polynôme R(x) tel que
Q(x) − P (x) = R(x)(x + 1)x(x − 1).
Puisque P (x) a degré 3 et Q(x) a degré 2 , le polynôme Q(x) − P (x) a degré 3 , donc le
polynôme R(x) qu’on a mis en facteur a degré 0 (i.e. R(x) est une constante).
Si on n’a pas remarqué ça, on peut tout de même faire tous les calculs : dans ce cas
n = 2 donc on a
(x − x1 ) (x − x2 ) (x − x0 ) (x − x2 ) (x − x0 ) (x − x1 )
Q(x) = y0 + y1 + y2
(x0 − x1 ) (x0 − x2 ) (x1 − x0 ) (x1 − x2 ) (x2 − x0 ) (x2 − x1 )
x(x − 1)
= − (x + 1)(x − 1) + (x + 1)x
2
1 1
= x2 + x + 1
2 2
Ainsi, directement on obtient :
1 1 1 1 2 1 1 1
Q(x)−P (x) = x2 + x+1+ x3 − x2 − x−1 = x3 − x = x x2 − 1 = λx(x+1)(x−1)
2 2 6 2 3 6 6 6
1
avecλ = 6