Réponses aux exercices du chapitre 6
Numéro 6. Un polynôme p(x) de degré 2 est tel que p(0) = 1, p(1) = 7 et p(2) = 17.
a) Calculer p0 (1) exactement sans calculer p(x).
b) Pouvez-vous faire de même pour calculer p00 (x), où x est arbitraire ?
c) Déterminer p(x) par la méthode de votre choix et comparer sa dérivée en x = 1 avec
la valeur trouvée en a).
1
Solution
a) Les formules à 3 points d'interpolation sont exactes pour les polynômes de degré 2. On
peut donc les utiliser pour obtenir p0 (1). Pour ce calcul, on utilise la formule suivante :
f (x + h) − f (x − h)
f 0 (x) =
2h
p(2) − p(0)
En x = 1 et avec h = 1, on a : p0 (1) = =8
2
b) On ne peut évaluer p00 (x) exactement que si x ∈ {0, 1, 2}, car ce sont les points
d'interpolation de f (À moins de calculer explicitement p(x)).
c) On a que x0 = 0, x1 = 1, x2 = 2, f (x0 ) = 1, f (x1 ) = 7 et f (x2 ) = 17.
Par la méthode de Lagrange on a que p2 (x) = f (x0 )L0 (x) + f (x1 )L1 (x) + f (x2 )L2 (x)
et donc :
(x − x1 )(x − x2 ) (x − 1)(x − 2) x2 − 3x + 2
L0 (x) = = =
(x0 − x1 )(x0 − x2 ) (0 − 1)(0 − 2) 2
2
(x − x0 )(x − x2 ) (x − 0)(x − 2) x − 2x
L1 (x) = = = = 2x − x2
(x1 − x0 )(x1 − x2 ) (1 − 0)(1 − 2) −1
(x − x0 )(x − x1 ) (x − 0)(x − 1) x2 − x
L2 (x) = = =
(x2 − x0 )(x2 − x1 ) (2 − 0)(2 − 1) 2
On a donc
x2 − 3x + 2
2
2 x −x
p2 (x) =1× + 7 × (2x − x ) + 17 ×
2 2
2 2 2
x 14x 17x 3x 28x 17x
p2 (x) = − + − + − +1
2 2 2 2 2 2
p2 (x) = 2x2 + 4x + 1
p0 (x) = 4x + 4
p0 (1) = 4(1) + 4 = 8
Or, par a) on a que toutes les méthodes à 3 points donneront cette même réponse.
2
1 √
Numéro 8. Soit f (x) = (x 1 − x2 +arcsin(x)). Pour x0 = 0,5, calculer des approximations
2
de f 0 (x0 ) par la formule de diérence arrière d'ordre 2 (voir le tableau 6.11) avec h = 0,1,
0,05, 0,023, 0,0125 et 0,00625. Dans chacun des cas, calculer l'erreur exacte en valeur absolue.
Cette erreur se comporte-t-elle en O(h2 ) ?
3
Solution
Pour x0 = 0,5, la valeur exacte de la dérivée est 0,866 025 4038.
La formule de diérence arrière d'ordre 2 est donnée par :
3f (x2 ) − 4f (x1 ) + f (x0 ) h2 f 000 (ξ2 )
f 0 (x2 ) = +
2h 3
Pour h = 0,1 :
On a x2 = 0,5, x1 = 0,4 et x0 = 0,3. Par conséquent, l'équation devient :
3f (0,5) − 4f (0,4) + f (0,3)
f 0 (0,5) ≈ ≈ 0, 870543101
0,1 × 2
Pour h = 0,05 :
On a x2 = 0,5, x1 = 0,45 et x0 = 0,4. Par conséquent, l'équation devient :
3f (0,5) − 4f (0,45) + f (0,4)
f 0 (0,5) ≈ ≈ 0, 867222895
0,05 × 2
Pour h = 0,025 :
On a x2 = 0,5, x1 = 0,475 et x0 = 0,45. Par conséquent, l'équation devient :
3f (0,5) − 4f (0,475) + f (0,45)
f 0 (0,5) ≈ ≈ 0, 864334832
0,025 × 2
Pour h = 0,0125 :
On a x2 = 0,5, x1 = 0,4875 et x0 = 0,475. Par conséquent, l'équation devient :
3f (0,5) − 4f (0,4875) + f (0,475)
f 0 (0,5) ≈ ≈ 0, 866104133
0,0125 × 2
Pour h = 0,00625 :
On a x2 = 0,5, x1 = 0,49375 et x0 = 0,4875. Par conséquent, l'équation devient :
3f (0,5) − 4f (0,49375) + f (0,4875)
f 0 (0,5) ≈ ≈ 0, 866045266
0,0125 × 2
4
Si on calcule les erreurs exactes en valeurs absolues, on obtient :
1 p 2
f (x) = x 1 − x + arcsin (x) x0 = 0,5
2
n hn en en /en+1
1 0,1 0, 004517696856111 3, 772
2 0,05 0, 001197491175961 3, 870
3 0,025 0, 000309428655761 3, 930
4 0,0125 0, 000078729449561 3, 964
5 0,00625 0, 000019861825961 −−
Le ratio en /en+1 tend vers 4 ce qui correspond à une erreur en O(h2 ).