Module : Analyse numérique TD: 2020-2021

Exercices corrigés

Interpolation polynômiale
Exercice 1
Déterminer le polynôme d’interpolation de Lagrange satisfaisant au tableau ci-dessous
x 0 2 3 5
f (x) −1 2 9 87

Corrigé : Rappelons que le polynôme de Lagrange basé sur les points d’appui d’abscisses x0 , x1 , . . . , xn est de degré
n et s’écrit :
n n
X Y x − xj
Pn (x) = f (xk )Lk (x) avec Lk (x) =
xk − xj
k=0 j=0,j6=k

ici les points d’appui donnés par :

x0 = 0 f (x0 ) = −1
x1 = 2 f (x1 ) = 2
x2 = 3 f (x2 ) = 9
x3 = 5 f (x3 ) = 87
déterminons donc un polynôme de Lagrange de degré 3, celui-ci s’écrit :
3
X
P3 (x) = f (xk )Lk (x)
k=0
avec
(x − x1 )(x − x2 )(x − x3 ) (x − x0 )(x − x2 )(x − x3 )
L0 (x) = L1 (x) =
(x0 − x1 )(x0 − x2 )(x0 − x3 ) (x1 − x0 )(x1 − x2 )(x1 − x3 )
(x − 2)(x − 3)(x − 5) (x − 0)(x − 3)(x − 5)
= =
(0 − 2)(0 − 3)(0 − 5) (2 − 0)(2 − 3)(2 − 5)
1 1
= − (x − 2)(x − 3)(x − 5) = x(x − 3)(x − 5)
30 6
(x − x0 )(x − x1 )(x − x3 )
L2 (x) =
(x2 − x0 )(x2 − x1 )(x2 − x3 ) (x − x0 )(x − x1 )(x − x2 )
L3 (x) =
(x − 0)(x − 2)(x − 5) (x3 − x0 )(x3 − x1 )(x3 − x2 )
= 1
(3 − 0)(3 − 2)(3 − 5) = x(x − 2)(x − 3)
1 30
= − x(x − 2)(x − 5)
6
Finalement
P3 (x) = f (x0 )L0 (x) + f (x1 )L1 (x) + f (x2 )L2 (x) + f (x3 )L3 (x)
53 3 253
= x − 7x2 + x−1
30 30
————————————————————
Exercice 2
1
Soit f (x) = . Déterminer le polynôme d’interpolation de Lagrange pour les points d’appui d’abscisses : −2, −1,
1 + x2
0, 1, 2. Ensuite discuter l’erreur d’interpolation.
1
Corrigé : Soit f (x) = . Les points d’appui sont :
1 + x2
1
x0 = −2 f (x0 ) = 5
1
x1 = −1 f (x1 ) = 2
x2 = 0 f (x2 ) = 1
1
x3 = 1 f (x3 ) = 2
1
x4 = 2 f (x4 ) = 5

1
Le polynôme de Lagrange est de degré 4. Il s’écrit
4
X
P4 (x) = f (xk )Lk (x)
k=0
avec
1 1
L0 (x) = x(x + 1)(x − 1)(x − 2) L1 (x) = − x(x + 2)(x − 1)(x − 2)
24 8
1 1
L2 (x) = (x + 2)(x + 1)(x − 1)(x − 2) L3 (x) = − x(x + 2)(x + 1)(x − 2)
4 6
1
L4 (x) = x(x + 2)(x + 1)(x − 1)
24
Finalement,
P4 (x)= f (x0 )L0 (x) + f (x1 )L1 (x) + f (x2 )L2 (x) + f (x3 )L3 (x) + f (x4 )L4 (x)
1 1
= x(x + 1)(x − 1)(x − 2) − x(x + 2)(x − 1)(x − 2)
120 12
1 1
+ (x + 2)(x + 1)(x − 1)(x − 2) − x(x + 2)(x + 1)(x − 2)
4 12
1
+ x(x + 2)(x + 1)(x − 1)
120
1 4 3 2
= x − x +1
10 5
Calculons l’erreur théorique sur cette interpolation. celle-ci est donnée ou point x par :
1
E(x) = f (x) − Pn (x) = γn+1 (x) − f (n+1) (ξx ) où ξ ∈ I = (min xi , max xi )
(n + 1)!
Elle vérifie,
n
1 Y  
|E(x)| 6 |γn+1 (x)| Mn+1 où γn+1 (x) = (x − xk ) Mn+1 = max f (n+1) (t)
 
(n + 1)! t∈I
k=0

1
Comme ici on a 5 points d’appui, cette erreur est majorée par : |E(x)| 6 |γ5 (x)| M5
Y 5!  
On a clairement γ5 (x) = (x − xk ) = x(x2 − 1)(x2 − 4). Il reste à calculer M5 = max f (5) (t). Un calcul assez long
 
t∈I
k=0
2 4
−240x(3 − 10x + 3x ) −240 
donne : f (5) (x) = de même, on trouve f (6) (x) = −21x6 + 105x3 − 63x2 + 3 .

(1 + x2 )6 (1 + x2 )7
Ainsi l’étude de f (5) donne M5 = 100. Finalement,
1  100  5
|E(x)| 6 |γ5 (x)| M5 = x(x2 − 1)(x2 − 4) = x(x2 − 1)(x2 − 4)

5! 5! 6
————————————————————

2
 Intégration numérique
Exercice 3 Z π
2
Déterminer par la méthode des trapèzes puis par celle de Simpson f (x)dx sur la base du tableau suivant :
0
π π 3π π
x 0 8 4 8 2
f (x) 0 0.382683 0.707107 0.923880 1
Ces points d’appui sont ceux donnant sin x, comparer alors les résultats obtenus avec la valeur exacte.

Corrigé :
Z π
2
I= f (x)dx
0
xn − x0 π
1. Soit T l’approximation de I par la méthode des trapèzes, le pas h donné par h = =
n 8
 3

h P
T = f (x0 ) + f (x4 ) + 2 f (xi )
2 i=1
π
= (0 + 1 + 2(0.382683 + 0.707107 + 0.92388))
16
= 0.987116

2. Soit S l’approximation de I par la méthode de Simpson. Celle-ci s’écrit,

h
S = (y0 + y4 + 4(y1 + y3 ) + 2y2 )
3
π1
= [(0 + 1 + 4(0.38... + 0.92...) + 2 × 0.707...)]
83
= 1.000135
Z π
2
Les points d’appui donnés dans cet execice correspondent à la fonction sin x. Et sin xdx = 1. On constate donc
0
que l’approximation de I donnée par la méthode de Simpson est meilleure que celle par les trapèzes,
puisque |S − I| = 0.000135 et |T − I| = 0.012884.
————————————————————
Exercice 4
On lance une fusée verticalement du sol et l’on mesure pendant les premières 80 secondes l’accélerationγ :
t (en s) 0 10 20 30 40 50 60 70 80
γ (en m/s2 ) 30 31.63 33.44 35.47 37.75 40.33 43.29 46.70 50.67
Calcule la vitesse V de la fusée à l’instant t = 80 s, par la méthode des trapèzes puis par Simpson.

Corrigé : On sait que l’acceleration γ est la dérivée de la vitesse V , donc,

Z t Z 80
V (t) = V (0) + γ(s)ds ⇒ V (80) = 0 + γ(s)ds
0 0
| {z }
I

1. Calculons I par la méthode des trapèzes. Ici, d’après le tableau des valeurs, h = 10.
 n−1

h P
I = γ(x0 ) + γ(xn ) + 2 γ(xi )
2 i=1
1
= × 10(30 + 50, 67 + 2(31, 63 + · · · + 46, 70))
2
= 3089 m/s

2. Calculons I par la méthode de Simpson

h
V (80) = (γ(x0 ) + γ(xn ) + 4(γ(1 ) + γ(x3 ) + · · · ) + 2(γ(2 ) + γ(x4 ) + · · · ))
3
10
= (30 + 50, 67 + 4(31, 63 + 35, 47 + · · · ) + 2(33, 44 + 37, 75 + · · · ))
3
= 3087 m/s

3
 ————————————————————

Exercice 5 Z π
Calculer à l’aide de la méthode des trapèzes l’intégrale I = sin x2 dx avec le nombre de points d’appui n = 5 puis
0
n = 10.
Z π
Corrigé : Soit I = sin x2 dx
0
π
1. n = 5 donc le pas d’intégration est h = . Calculons I par la méthode des trapèzes.
5
 n−1

h P
I = f (x0 ) + f (xn ) + 2 f (xi )
2 i=1
π π 2π 3π 4π
= (0 + 1 + 2(sin(π)2 + sin(0) + 2(sin( )2 + sin( )2 + sin( )2 + sin( )2 ))
10 5 5 5 5
= 0.504431

π
2. n = 10 donc le pas d’intégration est h = .
10
π π 2π 3π 4π
I = (0 + 1 + 2(sin(π)2 + sin(0) + 2(sin( )2 + sin( )2 + sin( )2 + sin( )2 ))
20 10 10 10 10
= 0.722338

alors que la valeur ‘exacte’ est approximativement 0, 772651. Avec ce pas plus petit l’approximation numérique est
meilleure.

————————————————————

Exercice 6
Trouver le nombre n de subdivisions nécessaires
Z π de l’intervalle d’intégration [−π, π], pour évaluer à 0.5 10−3 près,
grâce à la méthode de Simpson, l’intégrale cos x dx
−π

Corrigé : Soit Z π
I= cos x dx
−π

b−a 2π
Le pas d’intégration est h = = . D’autre part l’erreur théorique sur la méthode de Simpson est donnée par
n n
−(b − a) 4 (4)
E(h) = h f (ξ)
180
−2π 2π 4
= ( ) cos(ξ)
180 n
où ξ ∈ [a, b], par conséquent,  
 −2π 2π 4 
|E(h)| 6  ( ) 
180 n
 π 16π 4 
 
 6 0.5 10−3 , donc, n4 > 1 π
−3
Ainsi pour que |E(h)| 6 0.5 10 il suffit que n vérifie 

4 −3
16π 4 . Ainsi
90 n  0.5 10 90
n vérifie n > 18.6 On prendra par exemple n = 20, car pour la méthode de Simpson, le nombre de subdivisions de
l’intervalle [a, b] doit toujours être pair.

————————————————————

Exercice 7
Soit a 6 x0 < x1 < · · · < xn−1 < nn 6 b une partition fixée de l’intervalle [a, b]. Montrer qu’il existe un unique
(n + 1)-uplet (µ0 , µ1 , . . . , µn ) de nombres réels tels que
Z b n
X
P (x)dx = µi P (xi )
a i=0
Pour tout polynôme P de degré inférieur ou égal à n.

4
 Xn
Corrigé : Le polynôme P s’écrit dans la base de Lagrange P (x) = Li (x) P (xi ) (1)
n i=0
Y x − xj
avec Li (x) = , puis on intégre (1) sur [a, b], on obtient :
j=0
xi − xj
j6=i !
Z b Z bX n X n Z b n
X
P (x)dx = Li (x)P (xi ) dx = Li (x)dx P (xi ) = µi P (xi )
a a i=0 i=0 a i=0

————————————————————