Méthodes composites d’intégration

Positionnement du problème

Soit f une fonction continue sur un intervalle [a, b] ⊂ R.

Objectif : Approcher la valeur de l’intégrale
Z b
I (f ) = f (x)dx
a

1 / 12
 Méthodes composites d’intégration
Positionnement du problème

Soit f une fonction continue sur un intervalle [a, b] ⊂ R.

Objectif : Approcher la valeur de l’intégrale
Z b
I (f ) = f (x)dx
a

Les méthodes composites sont obtenues lorsqu’on décompose (subdivise) l’intervalle [a, b] en n
sous-intervalles [xk , xk+1 ], 0 ≤ k ≤ n − 1, de largeur uniforme
b−a
h= = xk+1 − xk ,
n

xk = a + kh, pour tout 0 ≤ k ≤ n: nœuds d’intégration.

h : pas de subdivision.
n le nombre de sous-intervalles.
n + 1 : nombre de points d’intégration.
n
[a, b] = ∪k=0 [xk , xk+1 ].
Z x1 Z x2 Z xn n−1 Z
X xk+1
I (f ) = f (x)dx + f (x)dx + · · · f (x)dx = f (x)dx, (relation de Chasles)
1 / 12 x0 x1 xn−1 k=0 xk
 Sur chaque sous-intervalle [xk , xk+1 ] (0 ≤ k ≤ n − 1), on calcule une approximation de l’intégrale
exacte I (f ) par les méthodes simples précédentes.

R. gauche R. droite R. milieu

s s s xk + xk+1
IRg (f ) = (xk+1 − xk )f (xk ) IRd (f ) = (xk+1 − xk )f (xk+1 ) IRm (f ) = (xk+1 − xk )f ( )
2

Trapèze Simpson

(xk+1 − xk )
 
s (xk+1 − xk ) s xk + xk+1
IT (f ) = (f (xk ) + f (xk+1 )) IS (f ) = f (xk ) + 4f ( ) + f (xk+1 )
2 6 2

2 / 12
 Méthodes composites du rectangle
Rectangle à gauche

La méthode composite du rectangle à gauche consiste à approximer l’intégrale

I (f ) par la somme de l’aire de tous les rectangles dont les côtés sont (xk+1 − xk )
et f (xk ).

Soit f une fonction continue sur [a, b], alors la n−1 Z

X xk+1
formule composite du rectangle à gauche est : c
I (f ) = f (x)dx w IRg (f )
k=0 xk
n−1 n−1
c
X X
IRg (f ) = f (xk )(xk+1 − xk ) = h f (a + kh),
k=0 k=0

avec xk = a + kh, (k = 0, 1, · · · , n − 1).

3 / 12
 Méthodes composites du rectangle
Rectangle à gauche

La méthode composite du rectangle à gauche consiste à approximer l’intégrale

I (f ) par la somme de l’aire de tous les rectangles dont les côtés sont (xk+1 − xk )
et f (xk ).

Soit f une fonction continue sur [a, b], alors la n−1 Z

X xk+1
formule composite du rectangle à gauche est : c
I (f ) = f (x)dx w IRg (f )
k=0 xk
n−1 n−1
c
X X
IRg (f ) = f (xk )(xk+1 − xk ) = h f (a + kh),
k=0 k=0

avec xk = a + kh, (k = 0, 1, · · · , n − 1).

Proposition
Si f ∈ C 1 ([a, b]), alors l’erreur du rectangle à gauche est
majorée par

c c (b − a)2 0
ERg (f ) = I (f ) − IRg (f ) ≤ max f (x)
2n x∈[a,b]

3 / 12
 Méthodes composites du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du rectangle à gauche en décomposant
l’intervalle d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :

4 / 12
 Méthodes composites du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du rectangle à gauche en décomposant
l’intervalle d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1. 1 1
Z
1
I (f ) = dx = [ln(1 + x)]0 = ln(2) = 0.693
0 1+x

4 / 12
 Méthodes composites du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du rectangle à gauche en décomposant
l’intervalle d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1. 1 1
Z
1
I (f ) = dx = [ln(1 + x)]0 = ln(2) = 0.693
0 1+x
b−a
2. Ici a = 0, b = 1, n = 4 et h = n = 1
4 = 0.25

4−1
c
X
IRg = 0.25 f (xk )
k=0

= 0.25(f (x0 ) + f (x1 ) + f (x2 ) + f (x3 ))

= 0.25(f (0) + f (0.25) + f (0.5) + f (0.75)) = 0.759

4 / 12
 Méthodes composites du rectangle à gauche
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du rectangle à gauche en décomposant
l’intervalle d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1. 1 1
Z
1
I (f ) = dx = [ln(1 + x)]0 = ln(2) = 0.693
0 1+x
b−a
2. Ici a = 0, b = 1, n = 4 et h = n = 1
4 = 0.25

4−1
c
X
IRg = 0.25 f (xk )
k=0

= 0.25(f (x0 ) + f (x1 ) + f (x2 ) + f (x3 ))

= 0.25(f (0) + f (0.25) + f (0.5) + f (0.75)) = 0.759

s c
3. ERg (f ) = I (f ) − IRg = |ln(2) − 0.759| ' 0.065
4 / 12
 Méthodes composites du rectangle
Rectangle à droite

La méthode composite du rectangle à droite consiste à approximer l’intégrale

I (f ) par la somme de l’aire de tous les rectangles dont les côtés sont (xk+1 − xk )
et f (xk+1 )

Soit f une fonction continue sur [a, b], alors la n−1 Z

X xk+1
c
formule composite du rectangle à droite est : I (f ) = f (x)dx w IRd (f )
k=0 xk
n−1 n−1
c X X
IRd (f ) = f (xk+1 )(xk+1 − xk ) = h f (a + (k + 1)h),
k=0 k=0

avec xk+1 = a + (k + 1)h.

5 / 12
 Méthodes composites du rectangle
Rectangle à droite

La méthode composite du rectangle à droite consiste à approximer l’intégrale

I (f ) par la somme de l’aire de tous les rectangles dont les côtés sont (xk+1 − xk )
et f (xk+1 )

Soit f une fonction continue sur [a, b], alors la n−1 Z

X xk+1
c
formule composite du rectangle à droite est : I (f ) = f (x)dx w IRd (f )
k=0 xk
n−1 n−1
c X X
IRd (f ) = f (xk+1 )(xk+1 − xk ) = h f (a + (k + 1)h),
k=0 k=0

avec xk+1 = a + (k + 1)h.

Proposition
Si f ∈ C 1 ([a, b]), alors l’erreur du rectangle à droite est
majorée par

c c (b − a)2 0
ERd (f ) = I (f ) − IRd (f ) ≤ max f (x)
2n x∈[a,b]

5 / 12
 Méthodes composites du rectangle à droite
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x
1 Approcher la valeur de I (f ) par la méthode composite du rectangle à droite en décomposant
l’intervalle d’intégration en 4 parties.
2 Quelle est l’erreur d’intégration pour cette méthode.

Solution :

6 / 12
 Méthodes composites du rectangle à droite
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x
1 Approcher la valeur de I (f ) par la méthode composite du rectangle à droite en décomposant
l’intervalle d’intégration en 4 parties.
2 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
b−a
1. Ici a = 0, b = 1, n = 4 et h = n = 1
4 = 0.25

4−1
c
X
IRd = 0.25 f (xk+1 )
k=0

= 0.25(f (x1 ) + f (x2 ) + f (x3 ) + f (x4 ))

= 0.25(f (0.25) + f (0.5) + f (0.75) + f (1)) = 0.634

6 / 12
 Méthodes composites du rectangle à droite
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x
1 Approcher la valeur de I (f ) par la méthode composite du rectangle à droite en décomposant
l’intervalle d’intégration en 4 parties.
2 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
b−a
1. Ici a = 0, b = 1, n = 4 et h = n = 1
4 = 0.25

4−1
c
X
IRd = 0.25 f (xk+1 )
k=0

= 0.25(f (x1 ) + f (x2 ) + f (x3 ) + f (x4 ))

= 0.25(f (0.25) + f (0.5) + f (0.75) + f (1)) = 0.634

2.
c c
ERd (f ) = |I (f ) − IRd (f )| = |ln(2) − 0.634| ' 0.059

6 / 12
 Méthodes composites du rectangle
Rectangle au milieu

La méthode composite du rectangle du milieu consiste à approximer l’intégrale I (f ) par la somme de

xk +xk+1
l’aire de tous les rectangle dont les cotés sont (xk+1 − xk ) et f ( 2 )

n−1 Z
X xk+1
c
Soit f une fonction continue sur [a, b], alors la formule I (f ) = f (x)dx w IRm (f )
xk
composite du rectangle au milieu est: k=0

n−1 n−1
c
X X
IRm (f ) = h f (mk ) = h f (a + (k + 21 )h),
k=0 k=0

avec mk = (xk + xk+1 )/2 = a + (k + 21 )h.

7 / 12
 Méthodes composites du rectangle
Rectangle au milieu

La méthode composite du rectangle du milieu consiste à approximer l’intégrale I (f ) par la somme de

xk +xk+1
l’aire de tous les rectangle dont les cotés sont (xk+1 − xk ) et f ( 2 )

n−1 Z
X xk+1
c
Soit f une fonction continue sur [a, b], alors la formule I (f ) = f (x)dx w IRm (f )
xk
composite du rectangle au milieu est: k=0

n−1 n−1
c
X X
IRm (f ) = h f (mk ) = h f (a + (k + 21 )h),
k=0 k=0

avec mk = (xk + xk+1 )/2 = a + (k + 21 )h.

Proposition
Si f ∈ C 2 ([a, b]), alors l’erreur du rectangle au milieu est
majorée par

c c (b − a)3 00
ERm (f ) = |I (f ) − IRm (f )| ≤ max |f (x)|
24n2 x∈[a,b]

7 / 12
 Méthode composite du trapèze
La méthode composite de trapèze consiste à approximer l’intégrale I (f ) par la somme de l’aire de
tous les trapèzes dont les bases sont f (xk ) et f (xk+1 ) et la hauteur est (xk+1 − xk )

n−1 Z xk+1
Soit f une fonction continue sur [a, b], alors la for- X
mule composite de trapèze est: I (f ) = f (x)dx w ITc (f )
k=0 xk
n−1
c
X xk+1 − xk
IT (f )= (f (xk ) + f (xk+1 ))
k=0
2
n−1
!
h X
= f (a) + f (b) + 2 f (a + kh) ,
2 k=1

avec xk = a + kh.

8 / 12
 Méthode composite du trapèze
La méthode composite de trapèze consiste à approximer l’intégrale I (f ) par la somme de l’aire de
tous les trapèzes dont les bases sont f (xk ) et f (xk+1 ) et la hauteur est (xk+1 − xk )

n−1 Z xk+1
Soit f une fonction continue sur [a, b], alors la for- X
mule composite de trapèze est: I (f ) = f (x)dx w ITc (f )
k=0 xk
n−1
c
X xk+1 − xk
IT (f )= (f (xk ) + f (xk+1 ))
k=0
2
n−1
!
h X
= f (a) + f (b) + 2 f (a + kh) ,
2 k=1

avec xk = a + kh.

Proposition
Si f ∈ C 2 ([a, b]), alors l’erreur de trapèze est majorée par

c c (b − a)3 00
ET (f ) = |I (f ) − IT (f )| ≤ max |f (x)|
12n2 x∈[a,b]

8 / 12
 Méthode composite du trapèze
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du trapèze en décomposant l’intervalle
d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :

9 / 12
 Méthode composite du trapèze
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du trapèze en décomposant l’intervalle
d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1. 1 1
Z
1
I (f ) = dx = [ln(1 + x)]0 = ln(2) = 0.693
0 1+x

9 / 12
 Méthode composite du trapèze
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du trapèze en décomposant l’intervalle
d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1. 1 1
Z
1
I (f ) = dx = [ln(1 + x)]0 = ln(2) = 0.693
0 1+x
b−a
2. Ici a = 0, b = 1, n = 4 et h = n = 1
4 = 0.25

4−1
c 0.25 X
IT = (f (x0 ) + 2 f (xk ) + f (x4 ))
2 k=1
0.25
= (f (x0 ) + 2(f (x1 ) + f (x2 ) + f (x3 )) + f (x4 ))
2
0.25
= (f (0) + 2(f (0.25) + f (0.5) + f (0.75)) + f (1)) = 0.697
2

9 / 12
 Méthode composite du trapèze
Exercice
Soit l’intégrale suivante: Z 1 1
I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite du trapèze en décomposant l’intervalle
d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
1. 1 1
Z
1
I (f ) = dx = [ln(1 + x)]0 = ln(2) = 0.693
0 1+x
b−a
2. Ici a = 0, b = 1, n = 4 et h = n = 1
4 = 0.25

4−1
c 0.25 X
IT = (f (x0 ) + 2 f (xk ) + f (x4 ))
2 k=1
0.25
= (f (x0 ) + 2(f (x1 ) + f (x2 ) + f (x3 )) + f (x4 ))
2
0.25
= (f (0) + 2(f (0.25) + f (0.5) + f (0.75)) + f (1)) = 0.697
2

3. ETC (f ) = |I (f ) − ITc | = |ln(2) − 0.697| ' 0.003

9 / 12
 Méthode composite de Simpson

On décompose l’intervalle [a, b] en n sous-intervalles avec n = 2p, p ∈ N∗ , En introduisant les

nœuds xk = a + kh pour k = 0, 1, · · · , n = 2p (i.e. on applique la méthode de Simpson simple sur
des intervalles [x2k , x2k+2 ] de largeur 2h et le point milieu est x2x+1 .

La fonction f est tabulée, et les xk sont équidistants n−1 Z

X xk+1
Soit f une fonction continue sur [a, b], alors la formule I (f ) = f (x)dx w ISc (f )
composite de Simpson est: xk
k=0

p−1 p−1

c h X X
IS (f ) = f (x ) + f (xn ) + 2
0 f (x2k ) + 4 f (x2k+1 )
3 k=1 k=0

p−1 p−1

h X X
= f (a) + f (b) + 2 f (a + 2kh) + 4 f (a + (2k + 1)h)
3 k=1 k=0

10 / 12
 Méthode composite de Simpson

On décompose l’intervalle [a, b] en n sous-intervalles avec n = 2p, p ∈ N∗ , En introduisant les

nœuds xk = a + kh pour k = 0, 1, · · · , n = 2p (i.e. on applique la méthode de Simpson simple sur
des intervalles [x2k , x2k+2 ] de largeur 2h et le point milieu est x2x+1 .

La fonction f est tabulée, et les xk sont équidistants n−1 Z

X xk+1
Soit f une fonction continue sur [a, b], alors la formule I (f ) = f (x)dx w ISc (f )
composite de Simpson est: xk
k=0

p−1 p−1

c h X X
IS (f ) = f (x ) + f (xn ) + 2
0 f (x2k ) + 4 f (x2k+1 )
3 k=1 k=0

p−1 p−1

h X X
= f (a) + f (b) + 2 f (a + 2kh) + 4 f (a + (2k + 1)h)
3 k=1 k=0

Proposition
Si f ∈ C 4 ([a, b]), alors l’erreur de Simpson est majorée par

C c (b − a)5 (4)
ES (f ) = |I (f ) − IS (f )| ≤ max |f (x)|
2880n4 x∈[a,b]

10 / 12
 Méthode composite de Simpson
Exercice

Soit l’intégrale suivante: Z 1 1

I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite de Simpson en décomposant l’intervalle
d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :

11 / 12
 Méthode composite de Simpson
Exercice

Soit l’intégrale suivante: Z 1 1

I (f ) = dx
0 1+x

1 Calculer la valeur exacte de I (f ).

2 Approcher la valeur de I (f ) par la méthode composite de Simpson en décomposant l’intervalle
d’intégration en 4 parties.
3 Quelle est l’erreur d’intégration pour cette méthode.

Solution :
2. Ici a = 0, b = 1, n = 2p, p = 2, et h = b−a = 1 = 0.25
n 4

2−1 2−1
!
c h X X
IT = f (a) + f (b) + 2 f (a + 2kh) + 4 f (a + (2k + 1)h)
3 k=1 k=0

= 0.693

3. ESc (f ) = |I (f ) − ISc | = |ln(2) − 0.693| ' 0.0001

11 / 12
 Exercice
Travail en asynchrone

Z 1 2
Soit la fonction I (f ) = e −x dx
0
1 En décomposant l’intervalle d’intégration en 4 parties, approcher la valeur de I (f )
par les méthodes composites du
a) rectangle du milieu.
b) trapèze.
c) Simpson.
2 En utilisant les inégalités de l’estimation de l’erreur d’intégration, quel est le
nombre de sous-intervalles n à considérer garantissant une erreur inférieure à 10−5
pour chacune de ces méthodes?

12 / 12