USTHB 2020-2021 Semestre 1 Analyse numérique

Faculté de Mathématiques 2ème année Lic Maths

Série d’exercices n◦ 4 : Intégration numérique

Exercice 1 :

Calculer Arctg (3) par les méthodes d’intégration des trapèzes et de Simpson pour n = 6.
Z x
1
Indication : Arctg (x) = 1+t2
dt.
0
Puis calculer l’erreur relative dans chaque cas.

Exercice 2 :

On lance une fusée verticalement du sol et l’on mesure pendant les premières 80 secondes l’accélération

a:
t en s 0 10 20 30 40 50 60 70 80

a en m/s2 30 31.63 33.44 35.47 37.75 40.33 43.29 46.70 50.67

Calculer la vitesse V de la fusée à l’instant t = 80 s, par la méthode des trapèzes puis par Simpson.

Exercice 3 :

Déterminer le nombre de subdivisions nécessaires des intervalles d’intégration pour évaluer à 0.5×10−6

près, les intégrales suivantes par la méthode indiquée.

Z π Z 1
1
a) cos xdx Simpson. b) 1+ex
dx trapèzes.
−π 0
Exercice 4 :
2
√
Z
Calculer xdx par la formule des rectangles à gauche, à droite et avec point au milieu, en
1
décomposant l’intervalle d’intégration en dix parties. Estimer puis calculer l’erreur commise.

Exercice 5 :

Soit f une fonction continue définie sur l’intervalle [−1, 1]. Notons par p le polynôme de degré deux

qui interpole f en les points −1, 0 et 1.

Z 1
a) Exprimer p (x) dx en fonction de f (−1), f (0) et f (1).
−1

b) Vérifier que l’expression obtenue coı̈ncide avec une formule d’intégration numérique dont on

donnera le nom et la valeur du pas de discrétisation.

1/6
Solution de l’exercice 1 :
3
b−a
Z
1 3 1
On a d’après l’indication Arctg (3) = 1+x2
dx. Donc a = 0, b = 3, n = 6, h = = 6
= 2
et
0 n
1
f (x) = .
1 + x2
Alors on a le tableau suivant :
i 0 1 2 3 4 5 6
1 3 5
xi 0 2
1 2
2 2
3
1 4 1 4 1 4 1
yi = f (xi ) = 1
1 + x2i 5 2 13 5 29 10
Z 3
1
Calculons Arctg (3) = 1+x2
dx par les méthodes des trapèzes et de Simpson. On a
0
 5
 1    
h P 2 4 1 4 1 4 1
T6 (f ) = f (a) + 2 f (xi ) + f (b) = 1+2 + + + + +
2 i=1 2 5 2 13 5 29 10
 
1 1467 1 18817
= 1+ + = ' 1.2478117,
4 377 10 15080

et
" #
h P5 P 5
S6 (f ) = f (a) + 2 f (xi ) + 4 f (xi ) + f (b)
3 i=1,i pair i=1,i impair
1  
2 1 1

4 4 4

1
= 1 + 2 2 + 5 + 4 5 + 13 + 29 +
3 10
 
1 7 9392 1 9403
= 1+ + + = ' 1.2470822.
6 5 1885 10 7540

Pour l’erreur commise, on a Arctg 3 = 1.2490458

|I (f ) − T6 (f )| |Arctg (3) − 1.2478117| |1.2490458 − 1.2478117|

Er (T6 ) = = =
|I (f )| |Arctg (3)| 1.2490458
0.0012341
= ' 9.8803423 × 10−4 ' 0.099%
1.2490458

et

|I (f ) − S6 (f )| |Arctg (3) − 1.2470822| |1.2490458 − 1.2470822|

Er (S6 ) = = =
|I (f )| |Arctg (3)| 1.2490458
0.0019636
= ' 1.5720801 × 10−3 ' 0.157% .
1.2490458

Solution de l’exercice 2 :

On sait que l’accélération a est la dérivée de la vitesse V , donc

Z t
V (t) = V (0) + a (s) ds.
0

Ce qui implique que

Z 80
V (80) = 0 + a (s) ds.
0

2/6
 Calculons I = V (80) par la méthode des trapèzes. Ici, d’après le tableau des valeurs, h = 10.
7
!
h X
V (80) ' T = a (0) + 2 a (si ) + a (80)
2 i=1
   
= 5 30 + 2 31.63 + 33.44 + 35.47 + 37.75 + 40.33 + 43.29 + 46.70 + 50.67

= 3089.45 m/s.

Calculons V (80) par la méthode de Simpson

!
h X X
V (80) ' S = a (s0 ) + a (sn ) + 2 a (si ) + 4 a (si )
3 i pair i impair
10     
= 30 + 50.67 + 2 33.44 + 37.75 + 43.29 + 4 31.63 + 35.47 + 40.33 + 46.70
3
' 3087.17 m/s.

Solution de l’exercice 3 :

a) Soit
Z π
I= cos xdx.
−π

Le pas d’intégration est

b−a 2π
h= = .
n n
D’autre part l’erreur théorique sur la méthode de Simpson est donnée par

b−a 4  
|E (h)| ≤ h sup f (4) (x)
180 x∈[a,b]
 4
2π 2π  
≤ car sup f (4) (x) = sup |cos x| ≤ 1.
180 n x∈[a,b] x∈[−π,π]

Ainsi pour que |E (h)| ≤ 0.5 × 10−6 il suffit que n vérifie

 4
2π 2π
≤ 0.5 × 10−6 .
180 n

Donc
2π 24 π 4
n4 ≥ ' 1.0880700 × 108
180 0.5 × 10−6
Ainsi n vérifie
√
4
n≥ 1.0880700 × 108 ' 102.132 56.

On prendra par exemple n = 104, car pour la méthode de Simpson, le nombre de subdivisions de

l’intervalle [a, b] doit toujours être pair.

3/6
 b) Soit
Z 1
1
I= dx.
0 1 + ex
Le pas d’intégration est
b−a 1
h= = .
n n
D’autre part l’erreur théorique sur la méthode des trapèzes est donnée par

b−a 2
|E (h)| ≤ h sup |f 00 (x)|
12 x∈[a,b]
 2  
1 1 e−1
≤ e
12 n (e + 1)3

car
 x ex − 1 
 
00
sup |f (x)| = sup e
 ≤e e−1 .
3
x∈[a,b] (ex + 1)  (e + 1)3
x∈[0,1]

Ainsi pour que |E (h)| ≤ 0.5 × 10−6 il suffit que n vérifie

 2  
1 1 e−1
e 3 ≤ 0.5 × 10−6 .
12 n (e + 1)

Ce qui implique que

s   
1 1 e−1
n≥ e ' 123.05673.
12 0.5 × 10−6 (e + 1)3

On prend par exemple n = 124.

Solution de l’exercice 4 :
b−a 2−1 1 i √
On a a = 1, b = 2, n = 10, h = = 10
= 10
, xi = a + ih = 1 + 10
, 0 ≤ i ≤ 10 et f (x) = x.
n
Alors on a les tableaux suivants :

i 0 1 2 3 4 5 6 7 8 9 10

10+i 11 12 13 14 15 16 17 18 19
xi = 10
1 10 10 10 10 10 10 10 10 10
2
√ q
11
q
12
q
13
q
14
q
15
q
16
q
17
q
18
q
19
√
f (xi ) = xi 1 10 10 10 10 10 10 10 10 10
2

i 1 2 3 4 5 6 7 8 9 10

xi−1 +xi 19+2i 21 23 25 27 29 31 33 35 37 39

2
= 20 20 20 20 20 20 20 20 20 20 20
q q q q q q q q q q
xi−1 +xi
21 23 25 27 29 31 33 35 37 39
f 2 20 20 20 20 20 20 20 20 20 20

4/6
 Pour l’approximation de l’intégrale on a
10
X  q q 
1 11 19
• rectangles à gauche Rg = h f (xi−1 ) = 10
1+ 10
+ ··· + 10
' 1.1981187.
i=1
10
X
1
q
11
q
19
√ 
• rectangles à droite Rd = h f (xi ) = 10 10
+ ··· + 10
+ 2 ' 1. 2395401.
i=1
10
X q q q 
xi−1 +xi
1 21 23 39
• rectangles avec Rm =h f 2
= 10 20
+ 20
+ ··· + 20
i=1
point au milieu ' 1.2190124.

Pour l’estimation de l’erreur commise on a

b−a b−a
|Ea (Rg )| ≤ h sup |f 0 (x)| , |Ea (Rd )| ≤ h sup |f 0 (x)| ,
2 x∈[a,b] 2 x∈[a,b]

b−a 2
|Ea (Rm )| ≤ h sup |f 00 (x)| .
24 x∈[a,b]

√ 1 −1
On a f (x) = x, f 0 (x) = √ et f 00 (x) = √ donc
2 x 4x x
1 1
sup |f 0 (x)| ≤ et sup |f 00 (x)| ≤ .
x∈[1,2] 2 x∈[1,2] 4

Ce qui implique

1 1 1
|Ea (Rg )| ≤ = 0.025, |Ea (Rd )| ≤ et |Ea (Rm )| ≤ ' 1.0416667 × 10−4 .
40 40 9600

La valeur exacte de l’intégrale est

Z 2 2 √
√ 2 √

4 2−2
I (f ) = xdx = x x = ' 1.2189514.
1 3 1 3

L’erreur absolue commise est donc

|Ea (Rg )| ' |1.2189514 − 1.1981187| = 0.0208327 ≤ 0.025,

|Ea (Rd )| ' |1.2189514 − 1. 2395401| = 0.0205887 ≤ 0.025,

1
|Ea (Rm )| ' |1.2189514 − 1.2190124| = 0.0000610 ≤ .
9600

Solution de l’exercice 5 :

a) On pose x0 = −1, x1 = 0 et x2 = 1. Les polynômes auxiliaires de Lagrange associés sont

(x − x1 ) (x − x2 ) (x − 0) (x − 1) 1
`0 (x) = = = x (x − 1) ,
(x0 − x1 ) (x0 − x2 ) (−1 − 0) (−1 − 1) 2
(x − x0 ) (x − x2 ) (x − (−1)) (x − 1)
`1 (x) = = = − (x + 1) (x − 1) ,
(x1 − x0 ) (x1 − x2 ) (0 − (−1)) (0 − 1)
(x − x0 ) (x − x1 ) (x − (−1)) (x − 0) 1
`2 (x) = = = x (x + 1) .
(x2 − x0 ) (x2 − x1 ) (1 − (−1)) (1 − 0) 2
5/6
 L’expression du polynôme d’interpolation de Lagrange est

f (x) ' p (x) = f (x0 ) `0 (x) + f (x1 ) `1 (x) + f (x2 ) `2 (x)

1 1
= x (x − 1) f (−1) − (x + 1) (x − 1) f (0) + x (x + 1) f (1) .
2 2

On intègre le polynôme p sur l’intervalle [−1, 1]

Z 1 Z 1
f (x) dx ' p (x) dx
−1 −1
Z 1 Z 1 Z 1
1 1
= f (−1) x (x − 1) dx − f (0) (x + 1) (x − 1) dx + f (1) x (x + 1) dx
2 −1 −1 2 −1
1 4 1
= f (−1) + f (0) + f (1) .
3 3 3

b) L’expression obtenue est la formule d’intégration de Simpson avec h = 1 et n = 2.

6/6