Le corrigé de TD3 d’analyse numérique 1- MIP, S4

1 Exercice 2
 1

1
1) Le polynôme d’interpolation d’une fonction Y qui est connue par les points : (0, 1); 2, e
2 ; (1, e).
En utilisant la méthode de Newton, on obtient
1 1 1
P2 (x) = Y (0) + Y [0, ](x − 0) + Y [0, , 1](x − 0)(x − )
2 2 2
 1
  1 
Nous trouverons P2 (x) = 2e − 4e 2 + 2 x2 + 4e 2 − e − 3 x + 1 avec x ∈ [0, 1].

(n+1)
2) On sait que e(x) = f (x)−P2 (x) = F (x) Y (n+1)!
(ζx ) 1


ici n = 2; F (x) = (x−0)× x − 2 ×(x−1)
1 eζx


et ζx ∈ [0, 1]. Donc e(x) = x × x − 2 × (x − 1) × 3! , d’où
 ζx 
1 e  e
|e(x)| 6 1 × × 1 × =⇒ |e(x)| 6
2 3! 12
3) En appliquant directement la formule de la différence centrale on aura
Y (1) − Y (0) 1
Y 0 (0, 5) ' ; ici h =
2h 2
Y 0 (0, 5) ' e − 1.

2 Exercice 3
x 0,8 0,9 1 1,1 1,2
1) On a les points suivants :
y 0,992 0,999 1 1,001 1,008
Si on appliquera par exemple la méthode de la différence progressive en deux points nous aurions,
y(1, 1) − y(1)
y 0 (1) ' = 0, 01.
1, 1 − 1
2) On a les points :
x 0 1 2 3 4
y 0 1 4 9 16
0
a) Pour approcher de y (4), On applique soit la méthode de différence régressive, soit on calcule
le polynôme d’interpolation en tous les points du tableau. Donc par la méthode de la différence
0
régressive en deux points, on a y (4) ' y(4)−y(3)
4−3 = 16 − 9 = 7.

b) Dans le cas de la méthode de la différence régressive ; c’est à dire d’utiliser le polynôme d’in-
2
terpolation justement aux points (3, 9), (4, 16), l’erreur s’écrit comme suivante : e(x) = F (x) y (ζ(x)
2 ,
2 y 2 (ζ(x)
d( )
et x ∈ [3, 4], d’où e0 (x) = F 0 (x) y (ζ(x)
2 +F (x)· 2
dx , avec F (x) = (x−3)×(x−4) = x2 −7x+12.
Donc pour x = 4 le deuxième terme sera nul, donc
y (2) (ζ4 )
e0 (4) = F 0 (4) − y[3, 4, 4] = F 0 (4). , avec ζ4 ∈ [3, 4]
2!
y (2) (ζ4 )
Puisque F 0 (x) = 2x − 7, alors F 0 (4) = 1 d’où e0 (4) = 2! .

1
 c) On a y(x) = x2 , donc y (2) (x) = 2, et y (3) (x) = 0. On sait que

sup |y (2) (t)| sup |y (3) (t)|

0 0 t∈[3,4] t∈[3,4]
|e (x)| ≤ sup |F (t)| × + sup |F (t)| ×
t∈[3,4] 2! t∈[3,4] 3!

sup |y (2) (t)|

0
0 t∈[3,4] 2
On a F (x) = 2x − 7, donc sup |F (t)| 6 1, en d’autre 2! = 2 = 1.
t∈[3,4]
sup |y (3) (t)|
t∈[3,4]
Or ∀x ∈ [3, 4] : y (3) (x) = 0, alors 3! = 0.
0 0
Enfin ∀x ∈ [3, 4] : |e (x)| 6 1, en particulier |e (4)| 6 1.

3) Calculons y 00 (2). Si on applique l’interpolation en trois points xi−1 = 1, xi = 2, xi+1 = 3, on

obtient la formule ci- dessous.
y (xi+1 ) − 2y (xi ) + y (xi−1 )
y 00 (xi ) ' , ici h = 1
h2
on obtient y 00 (2) ' y(3) − 2y(2) + y(1) = 9 − 8 + 1 = 2.

3 Exercice 4.
(p) Pn (p)
D’après l’interpolation de Lagrange, on a f (p) (x) ' Pn (x) = t=0 f (xi ) · Li (x), en par-
(p) n (p)
ticulier f (p) (α) ' Pn (α) = i=0 f (xi ) · Li (α), et on pose Ai (α) = L(p) (α) donc f (p) (α) '
P
Pn
1=0 f (xi ) · Ai (α).
Les polynômes de Lagrange Lt ne dépendent pas de f, ils dépendent que des points sur lesquels on
interpole.
(p) Pn
Alors on prend les mêmes points avec la fonction gk (x) = xk , on aura : gk (α) ' i=0 Ai (α) · xki .
(p)
Puisque si 0 6 k 6 p − 1, la dérivée d’ordre p de gk (x) = xk est gk (x) = 0, alors
n
X
Ai (α) · xki = 0
i=0

(p)
Si p 6 k 6 n, la dérivée d’ordre p de gk (x) = xk , est gk (x) = k × (k − 1) × · · · × (k − (p − 1)).xk−p ,
Pn
d’où i=0 Ai (α) · xki ' k × (k − 1) × · · · × (k − p + 1) · αk−p . Ce qui donne le résultat :
( P
n
Ai (α) · xki ' 0 si 0 6 k 6 p − 1
(∗∗) : Pni=0 k k−p
i=0 Ai (α) · xi ' k × (k − 1) × · · · × (k − p + 1) · α si p 6 k 6 n
P3
2) Ici p = 3, n = 3, et α = 2, donc f (3) (2) = i=0 Ai (2).f (xi ). Pour obtenir les Ai (2), qu’on les
notes simplement Ai , avec i = 0, 1, 2, 3, (4 inconnus) on devrais reformuler 4 équations en utilisant
( ∗∗) :
Pour k = 0 −→ A0 + A1 + A2 + A3 = 0
Pour k = 1 −→ A0 · x0 + A1 · x1 + A2 · x2 + A3 · x3 = 0
Pour k = 2 −→ A0 · x20 + A1 · x21 + A2 · x22 + A3 · x23 = 0
Pour k = 3 −→ A0 · x30 + A1 · x31 + A2 · x32 + A3 · x33 = 3!

2
On obtient le système suivant :


 A0 + A1 + A2 + A3 = 0

 0.A + A + 2A + 3A = 0
0 1 2 3


 0.A0 + A1 + 4A2 + 9A3 = 0
0.A0 + A1 + 8A2 + 27A3 = 6


Par la méthode de Gauss, on obtient :



 A0 = −1

 A =3
1
 A2 = −3


A3 = 1


Donc,
3
X
f (3) (2) = Ai (2) · f (xi ) = A0 .f (0) + A1 · f (1) + A2 · f (2) + A3 · f (3)
i=0

= −1 × 0 + 3 × 1 + (−3) × 4 + 1 × 9
= 0.

4 Exercice 5
π π 3π π
x 0 8 4 8 2
i) Soit le tableau suivant :
f(x) 0 0,382683 0,707107 0,923880 1
i) Méthode du trapèze : Choisissons la subdivision ; a = x0 = 0, x1 = π8 , x2 = π
4 , x3 = 3π
8 , et
b = x4 = π2 . D’après la formule des trapèzes avec un pas constant, on a ;
Z b "n−1 #
X 1 π
f (x)dx ' h × f (xi ) + (f (a) + f (b)) , ici h = , n = 4
a i=1
2 8

Donc, on obtient :
Z π2       π 
π π π 3π 1
f (x)dx ' × f +f +f + f (0) + f
0 8 8 4 8 2 2
 
π 1
= × 0, 382683 + 0, 707107 + 0, 923880 + ' 0, 98
8 2

ii) Méthode de Simpson : On subdivise l’intervalle 0, π2 comme suit : a = x0 = 0, x1 =

 
π π π 3π
4 , x2 = 2 = b. D’autre part η1 = 8 , η2 = 8 . Nous appliquons maintenant la formule dont le pas
π
ici est fixé h = 4 ;
Z π2 " n−1 n
#
h X X
f (x)dx ' f (a) + f (b) + 2 f (xi ) + 4 f (ηi )
0 6 i=1 i=1

Dans notre cas ici n = 2, par application numérique on obtient :

Z π2   
π π π π 3π
f (x)dx ' × f (0) + f + 2f + 4f + 4f
0 24 2 4 8 8
π
' × [0 + 1 + 2 × 0, 707107 + 4 × 0, 382683 + 4 × 0, 923880]
24

3
 R π
Donc, 2
0
f (x)dx ' 0, 99.

- Si ces points sont donnés par la fonction

h πi
sin : 0, −→ R
2
x 7−→ sin(x)

alors la valeur exacte est

Z π2 π
π
sin(x)dx = [− cos(x)]02 = − cos − (− cos(0)) = 1
0 2

5 Exercice 6
L’erreur commise pat la méthode de Simpson est donné par ;

(b − a) (4) π − (−π) 2π
e(f ) = − f (η) × h4 où f (x) = cos(x), et h = = ,
180 n n
Rπ
alors pour évaluer l’intégrale −π cos(x)dx à 0, 5 × 10−3 près, il suffit d’avoir l’inégalité |e(f )| =
   
 (b−a) (4)    (b−a) 
− 180 f (η) × h4  6 0, 5 × 10−3 · Or f (4) (η) = | cos(η)| 6 1, alors il suffit d’avoir − 180  ×

4 2π
h  6 0, 5 × 10−3 , ce qui donne : 2π × (2π)
 4
180 n 4 6 0, 5 × 10−3 =⇒ 0,5×10 180
−3 × (2π)
4
6 n4 =⇒
q
(2π)5 ×103 4 (2π)5 ×103 4 (2π)5 ×103
0,5×180 6 n =⇒ 90 6 n4 =⇒ 90 6 n =⇒ 3, 227 6 n.
Alors,il suffit de prendre le nombre de subdivision n > 4.