METHODES D’INTEGRATION NUMERIQUE
IV.1 INTRODUCTION
On veut calculer l’intégrale d’une fonction 𝑓 :
𝑰(𝒇∫) =
𝒃
𝒇(𝒙)𝒅𝒙
𝒂
– On ne connaît la valeur de 𝑓 qu’en certains points 𝑥0, 𝑥1, … , 𝑥𝑛 , et il n’est pas possible
d’avoir d’autres valeurs que celles-ci (c’est le cas quand la fonction 𝑓 est tabulée).
– Il est possible de calculer 𝑓(𝑥) pour un 𝑥 quelconque, mais la primitive de 𝑓 n’est pas
connue (par exemple𝑓(𝑥) = 𝑒sin ), ou bien l’expression analytique de 𝑓 est trop
(𝑥)
compliquée pour être explicitée (𝑓(𝑥) est par exemple le résultat d’un code de calcul trop
complexe).
Dans ce chapitre, on va présenter certaines méthodes numériques permettant de calculer d’une
façon approchée 𝐼(𝑓).En plus ,nous allons créer des programmes sur le logiciel bien connu
MATLAB pour utiliser ces méthodes numériques et calculer des intégrations.
On appelle 𝐽(𝑓) la valeur approchée de 𝐼(𝑓).
Pour certaines fonctions, il est impossible d'exprimer la fonction primitive en utilisant les
fonctions usuelles (algébriques, trigonométriques, exponentielles,..). C'est le cas par exemple
de la fonction 𝑒𝑠𝑖𝑛𝑥 . Cette fonction est pourtant intégrable, son graphe est donné par la figure
17.
Figure 1 : Exemple d’une fonction intégrable
Il faut signaler que ∫ 𝑓(𝑥)𝑑𝑥 représente l’aire de la surface délimitée par les droites 𝑥 =
𝑏
𝑎, 𝑥 = 𝑏, 𝑦 = 0 et la courbe 𝑦 = 𝑓(𝑥) (Figure 18).
Figure 2 : Surface délimitée
IV.2 Méthode des rectangles
La méthode des rectangles consiste à :
- diviser l'intervalle [a, b] en n segments égaux. On obtient ainsi (𝑛 + 1) points équidistants :
On pose :
𝑥𝑖 = 𝑎 + 𝑖h, 𝑖 = 0,1, …, 𝑛 avec h =(b-a)/n
- approximer la surface de chaque "tranche" par un rectangle (Figure 3).
Figure 3 : Méthode des rectangles
𝑥𝑖
∫ 𝑓(𝑥)𝑑𝑥 = (𝑥𝑖 − 𝑥𝑖−1)𝑓(𝛼𝑖) = h 𝑓(𝛼𝑖), 𝛼𝑖 ∈ [𝑥𝑖−1, 𝑥𝑖 ]
𝑥𝑖−1
La fonction est donc remplacée par une constante (polynôme de degré 0) sur chaque sous-
intervalle. On peut prendre 𝛼𝑖 = 𝑥𝑖 (point à droite) ou 𝛼𝑖 = 𝑥𝑖−1 (point à gauche), mais la
= (𝑥𝑖−1 + 𝑥 ).
1
𝛼
meilleur valeur de 𝑖 est celle du point milieu, c’est-à-dire :
2
𝛼𝑖 𝑖
En additionnant la somme des surfaces de tous les rectangles, on obtient :
𝑏 𝑥 𝑛
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓(𝑥)𝑑𝑥 = ∑ ∫ 𝑓(𝛼𝑖)𝑑𝑥
n xi
𝑎 𝑥0 𝑥𝑖−
𝑖=
𝑱(𝒇) = 𝒉 ∑𝒏𝒊= 𝒇(𝜶𝒊
1 1
𝟏
)
Dans la figure précédente, on a choisi de construire le rectangle à partir du point milieu de
chaque sous-intervalle.
IV.2.1 Exemple
- Calculer ∫5 e− x2 𝑑𝑥 en prenant n=5 .
0
Solution :
b - 5−0
𝑛=5⇒h = =1
a 5
=
𝑛
𝑥𝑖−1 + 𝑥𝑖
𝑥𝑖 𝑓(𝛼𝑖)
𝛼𝑖 =
2
𝑥0 = 𝑎 = 0
𝛼 = =0.5
𝑥0 +𝑥1
1
0.7788
2
𝑥1 = 𝑥0 + h = 1
𝛼 = =1.5
𝑥1 +𝑥2
0.1053
2 2
𝑥2 = 𝑥1 + h = 2 𝛼 =
𝑥2 +𝑥3
=2.5
3 2 0.0019
𝑥3 = 𝑥2 + h = 3
𝛼 = =3.5
𝑥3 +𝑥4
4
0.0000047
2
𝑥4 = 𝑥3 + h = 4 𝛼 = =4.5
𝑥4 +𝑥5
5 2
0.0000016
𝑥5 = 𝑥4 + h = 5 = 𝑏
5 𝑛
∫ e− x2 𝑑𝑥 = h ∑ 𝑓(𝛼𝑖) = 1 (0.7788+0.1053+0.0019+0.0000047+0.0000016) =
0.8860063
0
𝑖=1
IV.2.2 Programme matlab de la méthode des rectangles
Ce code effectue une approximation numérique de l’intégrale d’une fonction à l’aide de
méthode des rectangles, également connue sous le nom de méthode de la somme de
Riemann. Voici une explication détaillée du code :
f=@(x) x^(1/4) ; -Définit une fonction anonyme ‘f(x)’qui correspond à la fonction à
1
intégrer, dans ce cas x 4 .
a=0 ; b=5 ; - Définit les bornes d’intégration de a à b .
n=5 ; -Nombre de sous –intervalles.
h=(b-a)/n ; largeur de chaque sous-intervalle .
s=0 ; initialise une variable ‘s’ à zéro pour stocker la somme des valeurs de la
fonction du largeur d’un rectangle (h) dans la méthode des rectangles .
for i = 1 :n ; -Démarre une boucle qui itère à travers chaque sous – intervalle.
S= s+ f(a +i*h) ; -Ajoute à la somme (s) la valeur du fonction évaluée au point médian
de chaque sous-intervalle.
end ; -Terminer la boucle .
I= h*s ; -Calculer l’approximation numérique de l’intégrale .
IV.3 Méthode des Trapèzes
La méthode des Trapèzes consiste à :
- diviser l'intervalle [a, b] en n segments égaux. On obtient ainsi (𝑛 + 1) points équidistants.
x= 𝑎 + 𝑖h, 𝑖 = 0,1, … , 𝑛 avec h=(b−a)/n
On pose :
- approximer la surface de chaque "tranche" par un trapèze construit à partir des valeurs de la
fonction aux bornes de chaque sous-intervalle (Figure 4).
Figure 4: Méthode des trapèzes
𝑓(𝑥𝑖−1)
𝑥𝑖
∫ 𝑓(𝑥)𝑑𝑥
=h +
𝑥𝑖−1
𝑓(𝑥𝑖) 2
La fonction 𝑓 est donc remplacée par une droite (polynôme de degré 1) sur chaque sous-
intervalle. En additionnant la somme des surfaces de tous les trapèzes, on obtient :
𝑏 𝑛
𝑥𝑛 𝑥𝑖
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓(𝑥)𝑑𝑥 = ∑ ∫𝑓(𝑥)𝑑𝑥
𝑎 𝑥0 𝑖=1 𝑥𝑖−1
𝑏
h h h
∫� 𝑓(𝑥)𝑑𝑥 = (𝑓(𝑥0) + 𝑓(𝑥1)) + (𝑓(𝑥1) + 𝑓(𝑥2)) + ⋯ (𝑓(𝑥𝑛−1) +
𝑓�(𝑥𝑛 )) 2 2 2
𝑱(𝒇) = [𝒇(𝒂) + 𝒇(𝒃) + 𝟐 ∑𝒏−𝟏 𝒇(𝒙 )]
𝒉
𝟐 𝒊=𝟏 𝒊
IV.3.1 Exemple
- Calculer ∫5 e− x2 𝑑𝑥 en prenant n=5 .
0
b-a 5−0
h = =1
= 𝑛 5
𝑥𝑖 𝑓(𝑥𝑖)
𝑥0 = 𝑎 = 0 1
𝑥1 = 𝑥0 + h = 1 0.3678
𝑥2 = 𝑥1 + h = 2 0.0183
𝑥3 = 𝑥2 + h = 3 0.00012
𝑥4 = 𝑥3 + h = 4 0.000000112
𝑥5 = 𝑥4 + h = 5 = 0.00000000001
𝑏
D’après la formule des trapèzes, on a :
5
∫ 𝑓(𝑥)𝑑𝑥 (
h/2) [𝑓(0) + 𝑓(5) + 2(𝑓(1) + 𝑓(2) + 𝑓(3) +
= 𝑓(4))] = 0.8
0
IV.3.2 Programme matlab de la méthode des trapèzes
Ce code effectue une approximation numérique de l’intégrale d’une fonction à l’aide de
méthode des trapèzes. Voici une explication détaillée du code :
f=@(x) x^(1/4) ; -Définit une fonction anonyme ‘f(x)’qui correspond à la fonction à
1
intégrer, dans ce cas x 4 .
a=0 ; b=1 ; - Définit les bornes d’intégration de a à b .
n=100; -Nombre de sous –intervalles.
h=(b-a)/n ; largeur de chaque sous-intervalle .
s=f(a)+f(b); initialise une somme avec les valeurs de la fonction aux extrémités de
l’intervalle.
for i = 1 :n-1 ; -Démarre une boucle .
S= s+ 2*f(a +i*h) ; -Ajoute à la somme (s) la valeur du fonction évaluée sur chaque
point médian de chaque sous-intervalle.
end ; -Terminer la boucle .
I= (h/2)*s ; -Calculer l’approximation numérique de l’intégrale en utilisant la
méthode des trapèzes .
I= (h/2)*s ; -Calculer l’approximation numérique de l’intégrale en utilisant la
méthode des trapèzes .
IV.4 Méthode de Simpson
La méthode de Simpson consiste à :
- diviser l'intervalle [a, b] en n segments égaux avec 𝒏 un nombre pair(𝑛 = 2𝑚).
On
obtient ainsi (2𝑚 + 1) points équidistants
= 𝑎 + 𝑖h , 𝑖 = 0,1, … , 𝑛 avec h =
𝑥𝑖
(b-a)/n
- approximer la fonction sur chaque "tranche" par une parabole construite à partir de trois
points consécutifs (Figure 5).
Figure 5: Méthode de Simpson
b x2 x4 x2 m
∫ f ( x ) dx=∫ f ( x ) dx +∫ f ( x ) dx + …+ ∫ f ( x ) dx
a x0 x2 x 2 m−2
Ce qui donne après calcul :
𝑓(𝑥0 4𝑓(𝑥1 𝑓(𝑥2) h
𝑥2
) )
∫ 𝑓(𝑥)𝑑𝑥 h h
+ + h = (𝑓(𝑥0) + 4𝑓(𝑥1) + 𝑓(𝑥2))
= 3 3 3 3
𝑥0
m m−1
h
I ( f )= ( f ( a )+ f ( b ) ) + 4 ∑ f ( x 2i−1 ) +2 ∑ f ( x 2 i )
3 i=1 i=1
IV.4.1 Exemple (4.10)
5
𝑒𝑠𝑖𝑛𝑥 𝑑𝑥 en partageant l’intervalle [0, 5] en 4 segments .
∫ 0
- Calculer
Solution : Pour la méthode de Simpson, le nombre de segments doit être pair, c’est la raison
pour laquelle on a pris n=4 et non pas n=5 comme pour les cas précédents.
h= = = 1.25, 𝑛 = 2𝑚 = 4, ce qui donne 𝑚 = 2
𝑏−𝑎 5−0
𝑛 4
𝑥𝑖 𝑓(𝑥𝑖)
𝑥0 = 𝑎 = 0 1
𝑥1 = 𝑥0 + h =
1.25
2.5831
𝑥2 = 𝑥1 + h = 2.5
1.8193
𝑥3 = 𝑥2 + h =
1.5646
3.75
0.3833
𝑥4 = 𝑥3 + h = 5 =
𝑏
D’après la formule de Simpson, on a :
h
2 1
5
∫ 𝑓(𝑥)𝑑𝑥 = (𝑓(𝑎) + 𝑓(𝑏) + 4 ∑ 𝑓(𝑥2𝑖−1) + 2 ∑ 𝑓(𝑥2𝑖))
0 3 𝑖=1 𝑖=1
5
1.25
∫ 𝑓(𝑥)𝑑𝑥 (𝑓(𝑥0) + 𝑓(𝑥4) + 4(𝑓(𝑥1) + 𝑓(𝑥3)) + 2𝑓(𝑥2)) = 7.3387
= 3
0
IV.4.2 Programme MATLAB de la méthode de Simpson
Ce code effectue une approximation numérique de l’intégrale d’une fonction à l’aide de
méthode de Simpson . Voici une explication détaillée du code :
f=@(x) x^(1/4) ; -Définit une fonction anonyme ‘f(x)’qui correspond à la fonction à
1
intégrer, dans ce cas x 4 .
a=0 ; b=1 ; - Définit les bornes d’intégration de a à b .
n=100; -Nombre de sous –intervalles.
h=(b-a)/n ; largeur de chaque sous-intervalle .
s=f(a)+f(b); initialise une somme avec les valeurs de la fonction aux extrémités de
l’intervalle.
for i = 1 :2 :n-1 ; -Démarre une boucle .
S= s+ 4*f(a +i*h) ; -Ajoute à la somme (s) la valeur quatre fois les valeurs du fonction
évaluées sur les points médians des sous-intervalles.
end ; -Terminer la boucle .
for i= 1 :2 :n-2 ;-Démarre une boucle .
s = s +2 *f(a +i * h) ; -Ajoute a la somme actuelle (s) deux fois les valeurs du fonction
évaluées sur les points médians des sous-intervalles restants.
end ; -Terminer la boucle ;
I= (h/3)*s ; -Calculer l’approximation numérique de l’intégrale en utilisant la
méthode des trapèzes .