Université de Tlemcen Module Analyse numérique

Faculté de Technologie NS421 et GB422

Département de GEE, GBM

Chapitre 4: Intégration numérique

1) Introduction
Soit f une fonction définie sur un intervalle fixé [a, b] une partie de R
𝑏
Dans ce paragraphe, on s’intéresse à 𝐼 (𝑓) = ∫𝑎 𝑓(𝑥)𝑑𝑥 . Si l’intégrale définie est connue, on
utilise alors la formule de Newton-Leibnitz :𝐼 (𝑓) = 𝐹 (𝑎) − 𝐹 (𝑏) le problème se ramène à
calculer la primitive 𝐹 𝑑𝑒 𝑓 pour des valeurs données de la variable. Par ailleurs dans la pratique
la fonction 𝑓résulte de mesures physiques et donc données souvent tabulairement, et la notion
même de la primitive perd alors tout son sens. Dans ce cas on aura intérêt à utiliser les méthodes
approchées. D’où on se propose de chercher des approximations de𝐼 (𝑓) au moyen d’expression
de la forme : ∑𝑖=𝑛𝑖=0 𝐴i f(xi)=𝐴(𝑓)

a) Les abscisses xi sont données et équidistants

b) Les abscisses xi ne sont pas données à priori

2) Méthodes d’intégration
Soit f une fonction définie continue sur un intervalle [a, b]. Calculer 𝐼(𝑓) revient à calculer
l’aire limitée par les droites 𝑦 = 0, 𝑥 = 𝑎, 𝑥 = 𝑏 𝑒𝑡 𝑙𝑎 𝑐𝑜𝑢𝑟𝑏𝑒 𝑦 = 𝑓(𝑥)
Nous notons par :
𝑏
𝐼 (𝑓) = ∫ 𝑓 (𝑥)𝑑𝑥
𝑎
et 𝐴(𝑓) est la valeur approchée de 𝐼(𝑓)

Le calcul de l’intégrale est donné par les formules suivantes :

a) Formule du rectangle à gauche

le calcul de l’intégral est donné par la formule suivante :

𝑰𝑹𝑮 = 𝐼 (𝑓) ≅ 𝐴(𝑓) = (𝒃 − 𝒂)𝒇(𝒂)

b) Formule du rectangle à droite

Le calcul de l’intégrale est donné par la formule suivante :

𝐼𝑅𝐷 = 𝐼 (𝑓) ≅ 𝐴(𝑓) = (𝑏 − 𝑎)𝑓(𝑎 )

c) Formule du point du milieu

Le calcul de l’intégrale est donné par la formule suivante :

𝑎+𝑏
𝐼𝑃𝑀 = 𝐼 (𝑓) ≅ 𝐴(𝑓) = (𝑏 − 𝑎)𝑓 ( 2
)
Université de Tlemcen Module Analyse numérique
Faculté de Technologie NS421 et GB422
Département de GEE, GBM

d) Méthode des trapèzes

Le calcul de l’intégrale 𝐼 (𝑓)est donné par la formule suivante :

𝑓(𝑎)+𝑓(𝑏)
𝐼 ( 𝑓 ) ≅ 𝐴 (𝑓 ) = (𝑏 − 𝑎 ) ( )
2

Remarque
En général la précision fournie par un seul trapèze n’est pas suffisante. On divise l’intervalle [𝑎, 𝑏] en
𝑏−𝑎
segments partiels égaux [𝑎, 𝑎 + ℎ], [𝑎 + ℎ, 𝑎 + 2ℎ], … [𝑎 + (𝑛 − 1)ℎ, 𝑎 + 𝑛ℎ = 𝑏] avec ℎ = 𝑛
, et
à chacun desquels on applique la formule des trapèzes.
𝑏
𝑏 − 𝑎 𝑓(𝑎 ) + 𝑓(𝑏)
𝐼 (𝑓) = ∫ 𝑓(𝑥)𝑑𝑥 ≅ 𝐴(𝑓) = [ + 𝑓(𝑎 + ℎ ) + 𝑓 (𝑎 + 2ℎ ) + ⋯ 𝑓 (𝑎 + (𝑛 − 1)ℎ )]
𝑎 𝑛 2

e) Méthode de Simpson n=2

𝑎+𝑏
Soit f une fonction définit sur un intervalle [𝑎, 𝑏], Posons : 𝑥0 = 𝑎, 𝑥1 = 2
, 𝑥2 = 𝑏 . On interpole
la fonction f par un polynôme de Lagrange de degrés 2.

𝐿(𝑥) = 𝑓(𝑥0 )𝐿0 (𝑥) + 𝑓(𝑥1 )𝐿1 (𝑥) + 𝑓 (𝑥2 )𝐿2 (𝑥) ,
𝐿𝑖 𝑒𝑠𝑡 𝑙𝑒 𝑝𝑜𝑙𝑦𝑛ô𝑚𝑒 𝑑𝑒 𝐿𝑎𝑔𝑟𝑎𝑛𝑔𝑒 𝑎𝑢 𝑝𝑜𝑖𝑛𝑡 𝑥𝑖 𝑖 = 0,1,2 𝑑é𝑓𝑖𝑛𝑖𝑡 𝑎𝑢 𝑐ℎ𝑎𝑝𝑖𝑡𝑟𝑒 2
𝑏 𝑏
Donc 𝐼(𝑓) = ∫𝑎 𝑓 (𝑥)𝑑𝑥 = ∑20 ∫𝑎 𝑓 (𝑥𝑖 )𝐿𝑖 (𝑥)𝑑𝑥 = ∑20 𝑓(𝑥𝑖 ) 𝐼 (𝐿𝑖 ) . Après tous calculs faits l’intégrale
𝐼 (𝑓) est donnée par la formule suivante connue sous le nom de « formule de Simpson »
𝑏 𝑏−𝑎 𝑎+𝑏
𝐼 (𝑓) = ∫𝑎 𝑓 (𝑥) 𝑑𝑥 ≅ 𝐴(𝑓) = 6
[𝑓(𝑎) + 4𝑓(𝑐) + 𝑓(𝑏) ] 𝑐 =
2

REMARQUE :
Comme dans le cas de la méthode des trapèzes, il est recommandé de décomposer l’intervalle [𝑎, 𝑏]
en sous intervalles partiels égaux et appliquer la formule de Simpson à chaque intervalle.

f) Méthode de Newton (𝒏 = 𝟑)
Dans la méthode de Simpson nous avons divisé l’intervalle [𝑎, 𝑏] en deux parties égales, suite à ceci
nous avons interpolé la fonction 𝑓 par un polynôme de degré 2. La méthode de Newton consiste à
diviser l’intervalle [𝑎, 𝑏] en trois parties égales qui sont

𝑏−𝑎
[𝑎, 𝑎 + ℎ], [𝑎 + ℎ, 𝑎 + 2ℎ], [𝑎 + 2ℎ, 𝑎 + 3ℎ = 𝑏] , ℎ = dans ce cas nous interpolons f par un
3
polynôme de Lagrange de degré 3 (voir chapitre 2).

Donc :

𝑏 𝑏 𝑖=2
𝐼(𝑓) = ∫ 𝑓(𝑥) 𝑑𝑥 ≅ 𝐴(𝑓) = ∫ ∑(𝑓(𝑥𝑖 )𝐿𝑖 (𝑥))𝑑𝑥
𝑎 𝑎 𝑖=0
Université de Tlemcen Module Analyse numérique
Faculté de Technologie NS421 et GB422
Département de GEE, GBM
𝑏−𝑎 2𝑎+𝑏 𝑎+2𝑏
Ceci entraine que :𝐼 (𝑓) ≅ 𝐴(𝑓) = [𝑓 (𝑎) + 𝑓 ( )+𝑓( ) + 𝑓 (𝑏)] cette formule est la
8 3 3
formule de Newton.

g) Méthode de Newton –Cotes (𝒏 > 𝟑)

C’est exactement le même principe que la méthode de Simpson ou de Newton avec un nombre de
points plus que 4. C'est-à-dire on interpole la fonction 𝑓 avec un polynôme de Lagrange de degré
supérieur strictement à 3.

3) Méthode de Quadrature
Dans tout ce qui précède nous avons fait un calcul approché de 𝐼(𝑓) par une valeur 𝐴(𝑓) définie
sellant chaque cas

Théorème

Soit f une fonction n fois continument dérivable sur[𝑎, 𝑏] , telle que 𝑓 𝑛+1 existe sur [𝑎, 𝑏] sont
connues, 𝐴(𝑓) est l’approximation de 𝐼 (𝑓) d’ordre 𝑛 𝑑𝑒 𝐼 (𝑓) alors,

𝑏 𝑛
𝑀𝑛+1
|𝐼(𝑓) − 𝐴(𝑓)| ≤ ∫ |∏(𝑥 − 𝑥𝑖 )𝑑𝑥| 𝑀𝑛+1 = max|𝑓 𝑛+1 (𝑥)|
(𝑛 + 1)! 𝑎 [𝑎,𝑏]
𝑖=0

a) Le calcul de l’erreur dans la formule des trapèzes

(𝑏 − 𝑎)3
|𝐼(𝑓) − 𝐴(𝑓)| ≤ 𝑀2
12𝑛2

Remarque :

Pour appliquer la méthode des trapèzes au calcul d’une intégrale avec une erreur inferieur à ∈, on
choisit le nombre 𝑛 de trapèzes grâce à la majoration :

(𝑏−𝑎)3
𝑀2 12𝑛2
≤∈

b) Le calcul de l’erreur dans la formule de Simpson

(𝑏 − 𝑎)5
|𝐼(𝑓) − 𝐴(𝑓)| ≤ 𝑀4
180𝑛4

Le calcul de n se fait de la même que a)