Département de Mathématiques X22M040
Université de Nantes TD Méthodes Numériques
Interpolation
Exercice 1 : Exemples basiques
1. Soit f1 (x) = x3 sur [−2, 2].
(a) En utilisant les polynômes de Lagrange, construire le polynôme de degré 1 qui interpole f1 aux points 0 et
1. Que dire du résultat ?
√
(b) Quel est le polynôme de degré 3 qui interpole f1 aux points −π/2, 0, π/4 et e ?
2. On considère la fonction f2 (x) = cos(πx) sur [−1, 1].
(a) Construire le polynôme de degré 1 qui interpole f2 aux points −1 et 1. Que remarque-t-on ?
(b) Construire le polynôme de degré 1 qui interpole f2 aux points 0 et 1. Cette approximation est-elle meilleure
que la précédente ?
Exercice 2 : Polynôme de Lagrange
1. Construire le polynôme de Lagrange P qui interpole les points (−1, 1), (0, 1), (1, 2) et (2, 3).
2. Soit Q le polynôme qui interpole les points (−1, 1), (0, 1) et (1, 2). Sans faire de calcul, montrer qu'il existe un
réel λ tel que
Q(x) − P (x) = λ(x + 1)x(x − 1).
Exercice 3 : Choix des points d'interpolation
1
On considère la fonction f (x) = .
1 + x2
1. Construire le polynôme p1 qui interpole f aux points −0.5, 0 et 0.5.
√ √
2 2
2. Construire le polynôme p2 qui interpole f aux points − , 0 et .
2 2
3. Calculer les intégrales de p1 et de p2 entre −1 et 1. Conclure.
Exercice 4 : Interpolation et recherche de zéro
Pour calculer le zéro d'une fonction f (x), inversible sur un intervalle [a, b], on peut utiliser l'interpolation : considérant
une discrétisation xi de [a, b], et yi = f (xi ), on interpole l'ensemble {(yi , xi )}i=0,...,n et on obtient un polynôme P tel
que P (yi ) = xi . Alors chercher le zéro x de f revient à chercher x tel que
x = P (0).
1. Utiliser cette méthode pour évaluer l'unique racine a de la fonction f (x) = exp(x) − 2 dans l'intervalle [0, 1] avec
trois points d'interpolation x0 = 0, x1 = 21 , x2 = 1.
2. Comparer le résultat obtenu avec l'approximation du zéro de f par la méthode de Newton en 3 itérations à partir
de x0 = 0.
Exercice 5 :
Soient x0 < x1 < x2 dans [a, b] et f ∈ C 4 ([a, b]).
1. Montrer qu'il existe un unique polynôme p de degré inférieur ou égal à 3 tel que
p(x0 ) = f (x0 ), p0 (x1 ) = f 0 (x1 ), p(x1 ) = f (x1 ), p(x2 ) = f (x2 ).
2. Montrer que ∀x ∈ [a, b], ∃η ∈ [a, b] tel que
1 (4)
f (x) − p(x) = f (η)(x − x0 )(x − x1 )2 (x − x2 ).
4!
1
Exercice 6 : Estimation de l'erreur
On se donne la fonction x
f (x) = sin , x ∈ I = [0, 1].
3
1. Soit pn le polynôme interpolant la fonction f aux n÷uds x0 , x1 , . . . , xn équirépartis. Estimer l'erreur d'interpolation
En (f ) sur l'intervalle I , en fonction du degré n du polynôme et étudier son comportement quand n → ∞.
2. Trouver le nombre minimal de n÷uds équirépartis pour que En (f ) ≤ 10−4 .
Exercice 7 : (?) Diérences divisées
1
On cherche à interpoler f (x) = sur [0, 1].
1 + 2x
1. Écrire les polynômes de Lagrange de f aux points 0, 0.5 et 1. En déduire le polynôme d'interpolation de f en ces
3 points.
2. Retrouver ce polynôme en utilisant la formule des diérences divisées.
3. Pour améliorer la précision, on veut maintenant utiliser 5 points d'interpolation : 0, 0.25, 0.5, 0.75 et 1. En
évitant les calculs super us, construire le polynôme d'interpolation de f en ces 5 points.
Exercice 8 : (?) Interpolation de Hermite
1. Construire (hi ), i = 1, . . . , 4, les polynômes de degré 3 tels que
• h1 (0) = 1, h1 (1) = 0, h01 (0) = 0, h01 (1) = 0.
• h2 (0) = 0, h2 (1) = 1, h02 (0) = 0, h02 (1) = 0.
• h3 (0) = 0, h3 (1) = 0, h03 (0) = 1, h03 (1) = 0.
• h4 (0) = 0, h4 (1) = 0, h04 (0) = 0, h04 (1) = 1.
2. Montrer que (hi ), i = 1, . . . , 4, est une base de R3 [X].
3. Soit f une fonction de classe C 1 . Montrer à l'aide de la base (hi ), i = 1, . . . , 4, qu'il existe un unique polynôme
p de degré 3 tel que p(0) = f (0), p(1) = f (1), p0 (0) = f 0 (0) et p0 (1) = f 0 (1).
√
4. Construire p avec f (x) = 1 + x sur [0, 1].
Exercice 9 : (?) Interpolation sur 2 intervalles
Soit f une fonction continue de I = [0, 1] à valeurs dans R. On note I1 = [0, 12 [ et I2 = [ 12 , 1].
1. On choisit d'interpoler dans I1 en utilisant les points 1
6 et 13 . Exprimer les polynômes de Lagrange correspondants.
2. Par ailleurs, on choisit d'interpoler dans I2 en utilisant les points 12 , 34 , et 1. Exprimer les polynômes de Lagrange
correspondants.
3. En utilisant ces choix, quelle sera la fonction que l'on utilisera pour interpoler f sur I ?
4. Que dire de f en 21 ? Comment aurait-on pu résoudre ce problème de manière simple ?
5. On suppose maintenant que f est de classe C 1 . Peut-on assurer cette propriété à son interpolée avec la méthode
précédente ? Est-il possible de la modi er en conséquence ?