Chapitre 3 Intégration Numérique:

Méthode composites

x
x0 x1 x2 x3 xn

Unité Pédagogique de Mathématiques

1
 Principe des méthodes composites

𝑏
Soit 𝑓: 𝑎, 𝑏 → ℝ une fonction continue et soit: 𝐼 𝑓 = ‫ 𝑡𝑑 𝑡 𝑓 𝑎׬‬.

Objectif: approcher la valeur de 𝐼 𝑓 . 2

 Principe des méthodes composites
Solution:
1. Subdiviser l’intervalle [𝑎, 𝑏] en 𝑛 sous-intervalles 𝑥𝑖 , 𝑥𝑖+1 , 0 ≤ 𝑖 ≤ 𝑛 − 1 de largeur
𝑏−𝑎
ℎ= . 𝑥𝑖 = 𝑎 + 𝑖ℎ (0 ≤ 𝑖 ≤ 𝑛)
𝑛

ℎ: pas de subdivision
𝑛: nombre de sous-intervalles
𝑛 + 1 : nombre de points d’intégration

Ainsi on obtient via une multiple application de la relation de Chasles:

𝒃 𝒙𝟏 𝒙𝟐 𝒙𝒏+𝟏
𝑰 𝒇 = ‫𝒇 𝒂׬‬ 𝒕 𝒅𝒕= ‫𝒇 𝒙׬‬ 𝒕 𝒅𝒕+ ‫𝒇 𝒙׬‬ 𝒕 𝒅𝒕+…+ ‫𝒙׬‬ 𝒇 𝒕 𝒅𝒕
𝟎 𝟏 𝒏
3
 Principe des méthodes composites
𝑥𝑖+1
2. Approcher ‫𝑓 𝑥׬‬ 𝑡 𝑑𝑡, (0 ≤ 𝑖 ≤ 𝑛 − 1) par l’une des méthodes f(x)
𝑖 f (xi+1)
simples appliquées sur 𝑥𝑖 , 𝑥𝑖+1 : f (xi)

➢ Rectangle à droite simple: 𝑰𝒔𝑹𝒅 (𝒇) = 𝒇 𝒙𝒊+𝟏 𝒙𝒊+𝟏 − 𝒙𝒊

x
xi xi+1

f(x)
f (xi+1)
f (xi)

➢ Rectangle à gauche simple: 𝑰𝒔𝑹𝒈 (𝒇) = 𝒇 𝒙𝒊 𝒙𝒊+𝟏 − 𝒙𝒊

x
xi xi+1
4
 Principe des méthodes composites
f(x)
f (mi)
𝒙𝒊+𝟏 +𝒙𝒊
➢ Rectangle du milieu simple: 𝑰𝒔𝑹𝒎 (𝒇) =𝒇 𝒙𝒊+𝟏 − 𝒙𝒊
𝟐

xi mi xi+1 x
f(x)
f

𝒇(𝒙𝒊+𝟏 )+𝒇(𝒙𝒊 ) 𝒙𝒊+𝟏 + 𝒙𝒊

➢ Trapèze simple : 𝑰𝒔𝑻 𝒇 = 𝒙𝒊+𝟏 − 𝒙𝒊 𝒎𝒊 =
𝟐
𝟐

xi xi+1 x
f(x)

f
𝟏 𝒙𝒊+𝟏 −𝒙𝒊 𝒙𝒊+𝟏 +𝒙𝒊
➢ Simpson simple: 𝑰𝒔𝑺 𝒇 = (𝒇 𝒙𝒊 + 𝟒𝒇 + 𝒇 𝒙𝒊+𝟏 )
𝟑 𝟐 𝟐

xi mi xi+1 x
5
Les abréviations 𝑅𝑑, 𝑅𝑔, 𝑅𝑚, 𝑇, 𝑆 , c et 𝑠 désignent respectivement, Rectangle à droite, Rectangle à gauche, Rectangle au milieu, Trapèze, Simpson composite et simple.
Principe des méthodes composites

Dans tout ce qui suit nous considérons une fonction 𝒇 continue sur 𝒂, 𝒃
et 𝒙𝒊 = 𝒂 + 𝒊𝒉 , 𝟎 ≤ 𝒊 ≤ 𝒏 les points d’intégration issus de la subdivision
𝒃−𝒂
de 𝒂, 𝒃 en 𝒏 sous-intervalles avec un pas uniforme 𝒉 =
𝒏

6
 Méthodes des rectangles composites
❖ Rectangle à gauche
𝑖=𝑛−1 𝑥𝑖+1 𝑖=𝑛−1 𝑖=𝑛−1

𝐼 𝑓 = ෍ න 𝑓 𝑥 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖 𝑥𝑖+1 − 𝑥𝑖 = ℎ ෍ 𝑓 𝑥𝑖
𝑖=0 𝑥𝑖 𝑖=0 𝑖=0

𝒊=𝒏−𝟏

𝑰 𝒇 ≈ 𝑰𝒄𝑹𝒈 (𝒇) = 𝒉 ෍ 𝒇 𝒙𝒊 𝑓 𝑥
𝒊=𝟎

𝒄 désigne composite

𝑓 𝑥0 … … … … 𝑓 𝑥𝑛−1

x
x x x x x
0 1 2 3 n-1

7
 Méthodes des rectangles composites
❖ Rectangle à droite
𝑖=𝑛−1 𝑥𝑖+1 𝑖=𝑛−1 𝑖=𝑛−1

𝐼 𝑓 = ෍ න 𝑓 𝑥 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑖+1 𝑥𝑖+1 − 𝑥𝑖 = ℎ ෍ 𝑓 𝑥𝑖+1

𝑖=0 𝑥𝑖 𝑖=0 𝑖=0

𝒊=𝒏−𝟏

𝑰 𝒇 ≈ 𝑰𝒄𝑹𝒅 (𝒇) = 𝒉 ෍ 𝒇 𝒙𝒊+𝟏 𝑓 𝑥

𝒊=𝟎

𝒄 désigne composite
𝑓 𝑥1 … … … … 𝑓 𝑥𝑛

x
x x x x x
0 1 2 3 n

8
 Méthodes des rectangles composites
❖ Rectangle au milieu
𝑖=𝑛−1 𝑥𝑖+1 𝑖=𝑛−1 𝑖=𝑛−1
𝒙𝒊+𝟏 + 𝒙𝒊
𝐼 𝑓 = ෍ න 𝑓 𝑥 𝑑𝑥 ≈ ෍ 𝑓 𝑥𝑚𝑖 𝑥𝑖+1 − 𝑥𝑖 = ℎ ෍ 𝑓 𝑥𝑚𝑖 , 𝒎𝒊 =
𝟐
𝑖=0 𝑥𝑖 𝑖=0 𝑖=0

𝒊=𝒏−𝟏

𝑰 𝒇 ≈ 𝑰𝒄𝒎 (𝒇) = 𝒉 ෍ 𝒇 𝒙𝒎𝒊 f((x)

𝒊=𝟎

f (mn-1)
f (m0)

f (m2)
m m m m x
0 1 2 n-1

9
 Méthodes des rectangles composites
1 1
Exercice: Soit 𝐼(𝑓) = ‫׬‬0 𝑓(𝑥)𝑑𝑥 où 𝑓(𝑥) =
1+𝑥

1. Calculer la valeur exacte de 𝐼(𝑓)

2. Approcher la valeur de 𝐼(𝑓) par les méthodes d’intégration simples suivantes : rectangle à gauche et
rectangle à droite.
3. Approcher la valeur de 𝐼(𝑓) par les méthodes d’intégration composites suivantes : rectangle à
gauche et rectangles à droite, pour 𝑛 = 4 .

Solution:
1 1 1
1. 𝐼(𝑓) = ‫׬‬0 𝑑𝑥 = log(𝑥 + 1) 0 = log(2) ≈ 0,693
1+𝑥
2. 𝑰𝒔𝑹𝒈 (𝒇) = 𝟏 − 𝟎 𝒇(𝟎) = 𝟏,
𝟏
𝑰𝒔𝑹𝒅 (𝒇) = 𝟏 − 𝟎 𝒇(𝟏) = = 𝟎. 𝟓
𝟐

10
 Méthodes des rectangles composites
1−0 1−0 1
3. 𝑛 = 4 (4 sous intervalles et 5 points d’intégration) alors: ℎ = 𝑛
=
4
= = 0.25
4

𝑰𝒄𝑹𝒈 𝒇 = 𝟎. 𝟐𝟓 σ𝒊=𝟒−𝟏
𝒊=𝟎 𝒇 𝒙𝒊

= 𝟎. 𝟐𝟓 𝒇 𝒙𝟎 + 𝒇 𝒙𝟏 + 𝒇 𝒙𝟐 + 𝒇 𝒙𝟑

= 𝟎. 𝟐𝟓 𝒇 𝟎 + 𝒇 𝟎. 𝟐𝟓 + 𝒇 𝟎. 𝟓 + 𝒇 𝟎. 𝟕𝟓

= 𝟎. 𝟕𝟓𝟗

11
 Méthodes des rectangles composites

𝑰𝒄𝑹𝒅 𝒇 = 𝟎. 𝟐𝟓 σ𝒊=𝟒−𝟏
𝒊=𝟎 𝒇 𝒙𝒊+𝟏

= 𝟎. 𝟐𝟓 𝒇 𝒙𝟏 + 𝒇 𝒙𝟐 + 𝒇 𝒙𝟑 + 𝒇 𝒙𝟒

= 𝟎. 𝟐𝟓 𝒇 𝟎. 𝟐𝟓 + 𝒇 𝟎. 𝟓 + 𝒇 𝟎. 𝟕𝟓 + 𝒇 𝟏

= 𝟎. 𝟔𝟑𝟒

12
 Méthodes des trapèzes composites
❖ Trapèzes composites
𝑖=𝑛−1 𝑥𝑖+1 𝑖=𝑛−1 𝑖=𝑛−1 𝑖=𝑛−1
𝑓(𝑥𝑖+1 ) + 𝑓(𝑥𝑖 ) ℎ
𝐼 𝑓 = ෍ න 𝑓 𝑥 𝑑𝑥 ≈ ෍ 𝑥𝑖+1 − 𝑥𝑖 = ( ෍ 𝑓 𝑥𝑖 + ෍ 𝑓 𝑥𝑖+1 )
2 2
𝑖=0 𝑥𝑖 𝑖=0 𝑖=0 𝑖=0
𝑖=𝑛−1
ℎ
= (𝑓 𝑥0 + 2 ෍ 𝑓 𝑥𝑖 + 𝑓 𝑥𝑛 )
2
𝑖=1
f(x)
𝒊=𝒏−𝟏
𝒄 𝒉
𝑰 𝒇 ≈ 𝑰𝑻 (𝒇) = 𝒇 𝒙𝟎 + 𝟐 ෍ 𝒇 𝒙𝒊 + 𝒇 𝒙𝒏
𝟐
𝒊=𝟏

𝑰𝒄𝑹𝒈 𝒇 + 𝑰𝒄𝑹𝒅 (𝒇)

𝑰 𝒇 ≈ 𝑰𝒄𝑻 (𝒇) =
𝟐
x
x0 x1 x2 x3 xn

13
 Méthodes des trapèzes composites
1 1
Exercice: Soit 𝐼(𝑓) = ‫׬‬0 𝑓(𝑥)𝑑𝑥 où 𝑓(𝑥) =
1+𝑥

1. Calculer la valeur exacte de 𝐼(𝑓)

2. Approcher la valeur de 𝐼(𝑓) par la méthode d’intégration simple de Trapèze,
3. Approcher la valeur de 𝐼(𝑓) par les méthodes d’intégration composites des trapèzes pour 𝑛 = 4 .

Solution:
1 1 1
1. 𝐼(𝑓) = ‫׬‬0 𝑑𝑥 = log(𝑥 + 1) 0 = log(2) ≈ 0,693
1+𝑥
2. 𝑰𝒔𝑻 𝒇 = 𝟎. 𝟕𝟓

14
 Méthodes des trapèzes composites
1−0 1−0 1
3. 𝑛 = 4 (4 sous intervalles et 5 points d’intégration) alors: ℎ = 𝑛
=
4
= = 0.25
4

𝟎.𝟐𝟓
𝑰𝒄𝑻 𝒇 = 𝒇 𝒙𝟎 + 𝟐 σ𝒊=𝟒−𝟏
𝒊=𝟏 𝒇 𝒙𝒊 + 𝒇 𝒙𝟒
𝟐

𝟎.𝟐𝟓
= 𝒇 𝒙𝟎 + 𝟐(𝒇 𝒙𝟏 + 𝒇 𝒙𝟐 + 𝒇 𝒙𝟑 ) + 𝒇 𝒙𝟒
𝟐

= 𝟎. 𝟔𝟗𝟕

15
 Méthodes de Simpson composites
❖ Simpson composites

Si on ne connait pas l’expression de 𝒇 (fonction tabulée): On suppose que 𝒏 est pair

𝒏 𝒏
𝒊=𝟐−𝟏 𝒊=𝟐−𝟏
𝒉
𝑰 𝒇 ≈ 𝑰𝒄𝑺 (𝒇) = 𝒇 𝒙𝟎 + 𝟒 ෍ 𝒇 𝒙𝟐𝒊+𝟏 + 𝟐 ෍ 𝒇 𝒙𝟐𝒊 + 𝒇 𝒙𝒏
𝟑
𝒊=𝟎 𝒊=𝟏

Si on connait l’expression de 𝒇: cette formule est valable pour 𝒏 quelconque

𝒊=𝒏−𝟏
𝒉 𝒙𝒊+𝟏 + 𝒙𝒊
𝑰 𝒇 ≈ 𝑰𝒄𝑺 (𝒇) = ෍ 𝒇 𝒙𝒊 + 𝟒𝒇 + 𝒇 𝒙𝒊+𝟏
𝟔 𝟐
𝒊=𝟎

16
 Méthodes de Simpson composites
1 1
Exercice: Soit 𝐼(𝑓) = ‫׬‬0 𝑓(𝑥)𝑑𝑥 où 𝑓(𝑥) =
1+𝑥

1. Calculer la valeur exacte de 𝐼(𝑓)

2. Approcher la valeur de 𝐼(𝑓) par la méthode d’intégration simple de Simpson.
3. Approcher la valeur de 𝐼(𝑓) par les méthodes d’intégration composites de Simpson pour 𝑛 = 4 .

Solution:
1 1 1
1. 𝐼(𝑓) = ‫׬‬0 𝑑𝑥 = log(𝑥 + 1) 0 = log(2) ≈ 0,693
1+𝑥
2. 𝑰𝒔𝑺 𝒇 = 𝟎. 𝟔𝟗𝟒

17
 Méthodes de Simpson composites
1−0 1−0 1
3. 𝑛 = 4 (4 sous intervalles et 5 points d’intégration) alors: ℎ = 𝑛
=
4
= = 0.25
4
𝒉 𝒊=𝒑−𝟏 𝒊=𝒑−𝟏
On a : 𝑰 𝒇 ≈ 𝑰𝒄𝑺 (𝒇) = 𝒇 𝒙𝟎 + 𝟒 σ𝒊=𝟎 𝒇 𝒙𝟐𝒊+𝟏 + 𝟐 σ𝒊=𝟏 𝒇 𝒙𝟐𝒊 + 𝒇 𝒙𝒏
𝟑
Avec 𝒏 = 𝟐𝒑, 𝒑𝝐ℕ∗

Comme 𝑛 = 4 alors 𝑝 = 2

𝟎.𝟐𝟓
𝑰𝒄𝑺 𝒇 = 𝒇 𝒙𝟎 + 𝟒 σ𝒊=𝟐−𝟏
𝒊=𝟎 𝒇 𝒙𝟐𝒊+𝟏 + 𝟐 σ𝒊=𝟐−𝟏
𝒊=𝟏 𝒇 𝒙𝟐𝒊 + 𝒇 𝒙𝟒
𝟑

𝟎.𝟐𝟓
= 𝒇 𝒙𝟎 + 𝟒(𝒇 𝒙𝟏 + 𝒇 𝒙𝟑 ) + +𝟐𝒇 𝒙𝟐 + 𝒇 𝒙𝟒
𝟑

= 𝟎. 𝟔𝟗𝟑
18
 Erreurs d’intégration
❖ Méthodes du rectangles à gauche et à droite:
Si 𝑓 est une fonction de classe ∁1 sur 𝑎, 𝑏 alors les erreurs d’approximation de 𝑰(𝒇) par 𝑰𝒄𝑹𝒈 (𝒇)
et 𝑰𝒄𝑹𝒅 (𝒇 ) peuvent être majorées comme suit :

(𝒃 − 𝒂)𝟐 𝒄 désigne composite

𝑰 𝒇 − 𝑰𝒄𝑹𝒈 (𝒇) ≤ 𝒎𝒂𝒙 𝒇(𝟏) (𝒙)
𝟐𝒏 𝒙𝝐 𝒂,𝒃

(𝒃 − 𝒂)𝟐
𝑰 𝒇 − 𝑰𝒄𝑹𝒅 (𝒇) ≤ 𝒎𝒂𝒙 𝒇(𝟏) (𝒙)
𝟐𝒏 𝒙𝝐 𝒂,𝒃

❖ Méthodes du rectangles au milieu:

Si 𝑓 est une fonction de classe ∁2 sur 𝑎, 𝑏 alors l’erreur d’approximation de 𝑰(𝒇) par 𝑰𝒄𝒎 (𝒇) peut
être majorée comme suit :

(𝒃 − 𝒂)𝟑
𝑰 𝒇 − 𝑰𝒄𝒎 (𝒇) ≤ 𝒎𝒂𝒙 𝒇(𝟐) (𝒙)
𝟐𝟒𝒏𝟐 𝒙𝝐 𝒂,𝒃
19
 Erreurs d’intégration

❖ Méthodes du trapèze:
Si 𝑓 est une fonction de classe ∁2 sur 𝑎, 𝑏 alors l’erreur d’approximation de 𝑰(𝒇) par 𝑰𝒄𝑻 (𝒇) peut
être majorée comme suit :

(𝒃 − 𝒂)𝟑
𝑰 𝒇 − 𝑰𝒄𝑻 (𝒇) ≤ 𝒎𝒂𝒙 𝒇(𝟐) (𝒙)
𝟏𝟐𝒏𝟐 𝒙𝝐 𝒂,𝒃

❖ Méthodes de Simpson:
Si 𝑓 est une fonction de classe ∁4 sur 𝑎, 𝑏 alors l’erreur d’approximation de 𝑰(𝒇) par 𝑰𝒄𝑺 (𝒇) peut
être majorée comme suit :
(𝒃 − 𝒂)𝟓
𝑰 𝒇 − 𝑰𝒄𝑺 (𝒇) ≤ 𝒎𝒂𝒙 𝒇(𝟒) (𝒙)
𝟐𝟖𝟖𝟎𝒏𝟒 𝒙𝝐 𝒂,𝒃
20