MACS1 2010-2011 Analyse numérique
feuille d’exercices n˚7
Interpolation polynomiale
Exercice 1.
Considérons une fonction f dont le graphe passe par les points M0 = (0, 0), M1 = (1, 2),
M2 = (2, 16), et M3 = (3, 36).
1) Déterminer le polynôme de Lagrange p2 coïncidant avec f aux points M0 , M1 , et M2 .
2) Déterminer la forme du polynôme d’interpolation de Newton coïncidant avec f aux points M0 ,
M1 , et M2 .
3) En se servant des calculs de question précédente, déterminer le polynôme d’interpolation p3 coïn-
cidant avec f aux points M0 , M1 , M2 , et M3 .
4) En supposant que la fonction f est assez régulière, donner l’expression de l’erreur d’interpolation
E2 = f − p2 et de l’erreur E3 = f − p3 .
5) Donner une estimation de la valeur de f (1.5) par p2 et p3 .
Exercice 2 : Interpolation de Hermite.
Soit f ∈ C 1 ([a, b], R) et x1 et x2 deux points distincts. Soit p un polynôme de degré inférieur ou
égal à trois, vérifiant p(xi ) = f (xi ) et p′ (xi ) = f ′ (xi ) pour i = 1, 2.
1) Montrer qu’un tel polynôme existe et qu’il est unique.
2) Établir la majoration d’interpolation suivante : si f ∈ C 4 ([a, b], R), alors il existe ξ ∈ [a, b] tel
que :
(x − x1 )2 (x − x2 )2 (4)
f (x) − p(x) = f (ξ).
4!
3) Trouver une base (H1 , H2 , H̃1 , H̃2 ) de P3 telle que :
p(x) = f (x1 )H1 (x) + f (x2 )H2 (x) + f ′ (x1 )H̃1 (x) + f ′ (x2 )H̃2 (x), ∀p ∈ P3 ,
1
et exprimer cette base en fonction des polynômes d’interpolation de Lagrange L1 et L2 .
4) Décrire les polynômes d’interpolation de Hermite dans le cadre général.
Exercice 3 : Convergence de l’interpolation de Lagrange.
1
Soient α > 1, f une fonction définie sur [−1, 1] par :f (x) = et pn le polynôme d’interpola-
x−α
2
tion de Lagrange de f aux n + 1 points distincts xi = −1 + hi avec i ∈ [[0, n]] et h = .
n
1) Montrer que si α > 3, alors lim ||f − pn ||∞ = 0.
n→∞
2) En pratique, nous préférons utiliser des polynômes de degré peu élevé sur chaque intervalle
[xi , xi+1 ], avec i ∈ 0, . . . , n − 1. Notons fn la fonction continue telle que fn| est un polynôme
[xi ,xi+1 ]
de degré un et fn (xi ) = f (xi ) pour tout 0 ≤ i ≤ n.
a) Soit i ∈ [[0, n − 1]]. Écrire l’approximation de Lagrange P1i de degré un de f sur l’intervalle
[xi , xi+1 ].
C
b) Puis montrer que si α ∈
/ [−1, 1], alors ||f − fn ||∞ ≤ et donc que fn converge uniformément
n2
vers f lorsque n tend vers l’infini.
Exercice 4 : Propriétés des polynômes de Legendre.
Les polynômes de Legendre sont définis par :
P0 = 1,
n! dn 2
Pn (x) = (x − 1)n .
(2n)! dxn
1) Montrer que Pn est un polynôme de degré n et déterminer son coefficient devant le terme xn .
2) En utilisant des intégrations par parties, montrer que pour m ≤ n, nous avons
1 1
dn 2 dm dn+m 2
Z Z
n
(x − 1)n m (x2 − 1)m dx = (−1)n (x − 1)n (x2 − 1)m dx.
−1 dx dx −1 dxn+m
En déduire que
Z 1 0 si n 6= m
(n!)2
Pn (x)Pm (x)dx = ,
−1 (−1)n In sinon
(2n)!
2
Z 1
où In = (x2 − 1)n dx.
−1
3) Montrer que
2n
In = − In−1 .
(2n + 1)
4) Justifier l’écriture
n−1
X
xPn−1 (x) = Pn (x) + αi Pi (x),
i=0
avec αj ∈ R et 0 ≤ j ≤ n − 1.
Z 1
5) Calculer l’intégrale xPn (x)Pn−1 (x)dx.
−1
6) En écrivant Pn+1 (x) = (x − αn )Pn (x) − λn Pn−1 (x), montrer par récurrence que Pn a la même
parité que n et vérifie
P0 (x) = 1, P1 (x) = x,
n2 .
Pn+1 (x) = xPn (x) − Pn−1 (x)
(2n − 1)(2n + 1)
Exercice 5 : Splines cubiques.
Le but de cet exercice est l’étude d’un procédé d’interpolation d’une fonction, à valeurs réelles de
classe C 2 sur un intervalle [a, b], par une fonction cubique par morceaux, appelée spline cubique.
Soit (xi )0≤i≤n+1 une subdivision de l’intervalle [a, b]. Nous appelons spline cubique une fonction S
vérifiant les conditions suivantes :
• S ∈ C 2 ([a, b], R),
• S|[x ,x ] est polynôme de degré trois pour tout i = 0, . . . , n.
i i+1
Pour construire une telle approximation, nous cherchons à définir une spline S en fonction de ses
valeurs aux points xi , mais aussi en fonction de ses dérivées secondes en xi .
1) Sur un intervalle [α, β], montrer qu’il existe un unique polynôme P de degré inférieur ou égal à
trois défini par ses valeurs P (α), P (β), P ′′ (α), et P ′′ (β).
Déterminer ensuite les valeurs des dérivées premières en α et β en fonction des données.
2) En déduire qu’il existe une unique spline cubique S interpolant f au sens suivant :
S(xi ) = f (xi ), 0 ≤ i ≤ n + 1
S ′ (a) = f ′ (a), S ′ (b) = f ′ (b).
3) En considérant, pour i ∈ [[0, n + 1]], la fonction spline Si telle que
0 si j =
6 i
Si (xj ) = ,
1 si j = i
3
et Si′ (a) = Si′ (b) = 0. Puis les splines Sa et Sb telle que Sα (xi ) = 0 pour α ∈ {a, b} et Sa′ (a) = Sb′ (b) = 1
et Sb′ (a) = Sa′ (b) = 0. Montrer qu’une fonction spline S, interpolant f sur [a, b] s’écrit
n+1
X X
S(x) = fj Sj (x) + fα′ Sα (x).
j=0 α∈{a,b}