Université Ziane Achour de Djelfa Cours Méthodes numériques

Faculté des Science et de la Technologie Licence 0 2 Elt-Aut

Département de Génie Electrique Mr. Ali Teta

Chapitre III: Intégration numérique

III.1 Définition :
Dans cette section on s'intéresse au calcul approximatif d'une intégrale d’une fonction continue sur
l’intervalle [𝑎, 𝑏] définie par l’équation (3.1) en utilisant les méthodes numériques.
𝑏
𝐼(𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 (3.1)
Dans ces méthodes numériques, la fonction f, est remplacée par une somme finie constituée de n sous-
intervalles selon :
𝑏 (𝑏−𝑎)
∫𝑎 𝑓(𝑥)𝑑𝑥 = ∑𝑛𝑖=0 𝑓(𝑥𝑖 ) (3.2)
𝑛
On note pour tout 𝑖 = 0, … , 𝑛 :
𝑥𝑖+1 = 𝑥𝑖 + ℎ (3.3)
(𝑏−𝑎)
ℎ= (3.4)
𝑛
III.2 Méthode du point milieu
Dans cette méthode on subdivise l’intervalle [𝑎, 𝑏] de la fonction f en n sous-intervalles égaux [𝑥𝑖 , 𝑥𝑖+1 ].
Le schéma explicatif de cette méthode est représenté par la figure (01).

Figure 01 : La méthode du point milieu

On peut remarquer à travers la figure (01):

𝐼0 (𝑓) ≃ ℎ. 𝑓(𝑚0 )
𝐼1 (𝑓) ≃ ℎ. 𝑓(𝑚1 ) (3.5)
𝐼2 (𝑓) ≃ ℎ. 𝑓(𝑚2 )
𝐼3 (𝑓) ≃ ℎ. 𝑓(𝑚3 )
La formule de l’intégration par la méthode du point milieu est :
𝑏 (𝑥𝑖+1 +𝑥𝑖 )
𝐼(𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≃ ℎ. ∑𝑛−1
𝑖=0 𝑓(𝑚𝑖 ) avec 𝑚𝑖 = (3.6)
2

L’erreur d’intégration par cette méthode est :

𝑏 1 (𝑏−𝑎)3
𝐸𝑃 = |∫𝑎 𝑓(𝑥)𝑑𝑥 − ℎ. ∑𝑛−1
𝑖=0 𝑓(𝑚𝑖 ) | ≤ 24 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |𝑓(𝑥)′′ | (3.7)
𝑛2

12
Université Ziane Achour de Djelfa Cours Méthodes numériques
Faculté des Science et de la Technologie Licence 0 2 Elt-Aut
Département de Génie Electrique Mr. Ali Teta

Exemple : On applique la méthode du point milieu pour le calcul approché de l’intégral (), pour sous-
intervalles égaux n=5.
𝑏 1
∫𝑎 𝑓(𝑥)𝑑𝑥 = ∫0 𝑥 2 𝑑𝑥 (3.8)
On a:
(𝑏−𝑎) (1−0)
ℎ= = = 0.2 (3.9)
𝑛 5
Par conséquent, les sous-intervalles ([𝑥𝑖 , 𝑥𝑖+1 ] 𝑝𝑜𝑢𝑟 𝑖 = 0, ,5) sont constitués de :
{[0, 0.2], [0.2, 0.4], [0.4, 0.6], [0.6, 0.8] , [0.8, 1]} (3.10)
Les points milieux (𝑚𝑖 𝑝𝑜𝑢𝑟 𝑖 = 0, ,4) de ces intervalles sont :
{0.1, 0.3, 0.5, 0.7, 0.9 } (3.11)
La formule du point milieu est :
𝑏
𝐼(𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≃ ℎ. ∑𝑛−1
𝑖=0 𝑓(𝑚𝑖 ) (3.12)
AN:
𝐼(𝑓) ≃ 0.2(𝑓(𝑚0 ) + 𝑓(𝑚1 ) + 𝑓(𝑚2 ) + 𝑓(𝑚3 ) + 𝑓(𝑚4 ))
𝐼(𝑓) ≃ 0.2(𝑓(0.1) + 𝑓(0.3) + 𝑓(0.5) + 𝑓(0.7) + 𝑓(0.7)) (3.13)
𝐼(𝑓) ≃ 0.33
L’erreur d’intégration par cette méthode (𝐸𝑃 ) est :

1 (𝑏−𝑎)3
𝐸𝑃 ≤ 24 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |𝑓(𝑥)′′ | (3.14)
𝑛2
AN:
1 (1−0)3
𝐸𝑃 ≤ 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |2| (3.15)
24 52
1 (1−0)3
𝐸𝑃 ≤ 24 . 2 ⇒ 𝐸𝑃 ≤ 0.003 (3.16)
25

III.3 Méthode du Trapèze

On considère la fonction f continue et dérivable sur l’intervalle [𝑎, 𝑏]. Dans cette méthode sur chaque
sous-intervalle[𝑥𝑖 , 𝑥𝑖+1 ], la fonction f est substituée par la droite joignant les points (𝑥𝑖 , 𝑓(𝑥𝑖 )),
et (𝑥𝑖+1 , 𝑓(𝑥𝑖+1 )).

Figure 02 : La méthode du Trapèze

13
Université Ziane Achour de Djelfa Cours Méthodes numériques
Faculté des Science et de la Technologie Licence 0 2 Elt-Aut
Département de Génie Electrique Mr. Ali Teta

A travers le schéma explicatif représenté par la figure (02), on peut remarquer :

(𝑥1 −𝑥0 ) ℎ
𝐼0 (𝑓) ≃ . [𝑓(𝑥0 ) + 𝑓(𝑥1 )] = 2 . [𝑓(𝑎) + 𝑓(𝑥1 )]
2
(𝑥2 −𝑥1 ) ℎ
𝐼1 (𝑓) ≃ . [𝑓(𝑥1 ) + 𝑓(𝑥2 )] = 2 . [𝑓(𝑥1 ) + 𝑓(𝑥2 )] (3.17)
2
(𝑥3 −𝑥2 ) ℎ
𝐼2 (𝑓) ≃ . [𝑓(𝑥2 ) + 𝑓(𝑥3 )] = 2 . [𝑓(𝑥2 ) + 𝑓(𝑥3 )]
2
(𝑥4 −𝑥3 ) ℎ
𝐼3 (𝑓) ≃ . [𝑓(𝑥3 ) + 𝑓(𝑥4 )] = 2 . [𝑓(𝑥3 ) + 𝑓(𝑏)]
2
L’intégrale par la méthode de Trapèze est approchée par la formule :
𝑏 ℎ
𝐼(𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≃ . [𝑓(𝑎) + 𝑓(𝑏)] + ℎ. ∑𝑛−1 𝑖=1 𝑓(𝑥𝑖 ) (3.18)
2
Avec l’erreur :
𝑏 1 (𝑏−𝑎)3
𝐸𝑇 = |∫𝑎 𝑓(𝑥)𝑑𝑥 − 𝐼(𝑓)| ≤ 12 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |𝑓(𝑥)′′ | (3.19)
𝑛2

Exemple : On applique la méthode des Trapèzes pour le calcul approché de l’intégral (1.47), pour sous-
intervalles égaux n=4.
𝑏 1
∫𝑎 𝑓(𝑥)𝑑𝑥 = ∫0 𝑥 2 𝑑𝑥 (3.20)
On a:

(𝑏−𝑎) (1−0)
ℎ= = = 0.25 (3.21)
𝑛 4
Par conséquent, les sous-intervalles ([𝑥𝑖 , 𝑥𝑖+1 ] 𝑝𝑜𝑢𝑟 𝑖 = 0, ,4) sont constitués de :
{[0, 0.25], [0.25, 0.5], [0.5, 0.75], [0.75, 1] } (3.22)
La formule du l’intégration par la méthode du Trapèze est :
𝑏 ℎ
𝐼(𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≃ 2 . [𝑓(𝑎) + 𝑓(𝑏)] + ℎ. ∑𝑛−1
𝑖=1 𝑓(𝑥𝑖 ) (3.23)
AN:
0.25
𝐼(𝑓) ≃ . [𝑓(0) + 𝑓(1)] + 0.25 ∑3𝑖=1 𝑓(𝑥𝑖 ) (3.24)
2
0.25
𝐼(𝑓) ≃ . [(0)2 + (1)2 ] + 0.25((0.25)2 + (0.5)2 + (0.75)2 ) ≃ 0.342 (3.25)
2
L’erreur par cette méthode est :
𝑏 1 (𝑏−𝑎)3
𝐸𝑇 = |∫𝑎 𝑓(𝑥)𝑑𝑥 − 𝐼(𝑓)| ≤ 12 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |𝑓(𝑥)′′ | (3.26)
𝑛2

AN:
1 (1−0)3
𝐸𝑇 ≤ 12 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |2| ⇒ 𝐸𝑇 ≤ 0.010 (3.27)
42

14
Université Ziane Achour de Djelfa Cours Méthodes numériques
Faculté des Science et de la Technologie Licence 0 2 Elt-Aut
Département de Génie Electrique Mr. Ali Teta

III.4 Méthode de Simpson

La méthode de Simpson est basée sur l’interpolation, de chaque sous-intervalle [𝑥𝑖 , 𝑥𝑖+1 ] par un
polynôme de degré deux. La fonction f est substituée par ce polynôme du second degré qui définit un
arc de parabole passant par les points d’ordonnées 𝑓(𝑥𝑖 ), 𝑓(𝑥𝑖+1 )𝑒𝑡 𝑓(𝑥𝑖+2 ).

Figure 03 : La méthode de Simpson

A travers le schéma représenté par la figure () :

ℎ
𝐼0 (𝑓) ≃ . [𝑓(𝑎) + 4𝑓(𝑥1 ) + 𝑓(𝑥2 )]
3
ℎ
𝐼1 (𝑓) ≃ 3 . [𝑓(𝑥2 ) + 4𝑓(𝑥3 ) + 𝑓(𝑏)] (3.28)

L’intégrale par la méthode de Simpson est approchée par la formule :

𝑏 ℎ
𝐼(𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≃ 3 . (𝑓(𝑎) + 𝑓(𝑏) + 4 ∑𝑛−1 𝑛−2
𝑖=1,3,5.. 𝑓(𝑥𝑖 ) + 2 ∑𝑗=2,4,6.. 𝑓(𝑥𝑗 )) (3.29)

Avec l’erreur :

𝑏 1 (𝑏−𝑎)5
𝐸𝑆 = |∫𝑎 𝑓(𝑥)𝑑𝑥 − 𝐼(𝑓)| ≤ 180 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |𝑓 (4) (𝑥)| (3.30)
𝑛4

Exemple : On applique la méthode de Simpson pour le calcul approché de l’intégral (1.47), pour sous-
intervalles égaux n=4.
𝑏 1
∫𝑎 𝑓(𝑥)𝑑𝑥 = ∫0 𝑥 2 𝑑𝑥 (3.31)

On a:

(𝑏−𝑎) (1−0)
ℎ= = = 0.25 (3.32)
𝑛 4

Par conséquent, les sous-intervalles ([𝑥𝑖 , 𝑥𝑖+1 ] 𝑝𝑜𝑢𝑟 𝑖 = 0, ,4) sont constitués de :
15
Université Ziane Achour de Djelfa Cours Méthodes numériques
Faculté des Science et de la Technologie Licence 0 2 Elt-Aut
Département de Génie Electrique Mr. Ali Teta

{[0, 0.25], [0.25, 0.5], [0.5, 0.75], [0.75, 1]} (3.33)

La formule du l’intégration par la méthode de Simpson est :
𝑏 ℎ
𝐼(𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 ≃ 3 . (𝑓(𝑎) + 𝑓(𝑏) + 4 ∑𝑛−1 𝑛−2
𝑖=1,3,5 𝑓(𝑥𝑖 ) + 2 ∑𝑗=2,4,6 𝑓(𝑥𝑗 )) (3.34)
ℎ
𝐼(𝑓) ≃ 3 . (𝑓(𝑥0 ) + 𝑓(𝑥4 ) + 4(𝑓(𝑥1 ) + 𝑓(𝑥3 )) + 2𝑓(𝑥2 )) (3.35)
AN:
0.25
𝐼(𝑓) ≃ . (0 + 1 + 4(0.062 + 0.562) + 2(0.25)) ≃ 0.333 (3.36)
3
L’erreur par cette méthode est :
𝑏 1 (𝑏−𝑎)5
𝐸𝑆 = |∫𝑎 𝑓(𝑥)𝑑𝑥 − 𝐼(𝑓)| ≤ 180 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |𝑓 (4) (𝑥)| (3.37)
𝑛4
AN:
1 (1−0)3
𝐸𝑆 ≤ 2880 𝑚𝑎𝑥𝑥∈[𝑎,𝑏] |0| ⇒ 𝐸𝑇 = 0 (3.38)
42

16