CHAPITRE III : INTEGRATION NUMERIQUE
Le but de ce chapitre est d’aborder le calcul de l’intégrale d’une fonction 𝑓(𝑥) sur un
domaine fini délimité par des bornes finies 𝑎 et 𝑏.
II.1. Intégration Numérique :
Dans la plupart des cas, les fonctions analytiques, du fait de leurs complexités, ne sont
pas intégrables analytiquement. Dans d'autres cas, on a des fonctions qui sont évaluées
numériquement en différents points de l’intervalle où ces dernières sont données, et l’intégrale
de ces types de fonctions ne peut être obtenue que par des approches numériques. Dans ce
chapitre, on s’intéresse aux méthodes utilisées fréquemment ; à savoir la méthode des
trapèzes, la méthode de Simpson.
L’idée de ces deux méthodes est d’approximer la fonction 𝑓 à intégrer par un
polynôme 𝑃𝑛 (𝑥),.en utilisant la formule vue en interpolation de newton ::
0𝛼 𝛼(𝛼 − 1) 2 𝛼(𝛼 − 1) … … (𝛼 − 𝑛 + 1) 𝑛
𝑃𝑛 (𝑥) = 𝑦0 + ∆𝑦0 + 2
∆ 𝑦0 + ⋯ … … + ∆ 𝑦0
1! ℎ 2! ℎ 𝑛! ℎ𝑛
𝑥−𝑥0
Avec 𝛼 = .
ℎ
III.2. Méthode des trapèzes :
La formule des trapèzes est obtenue en remplaçant 𝑓(𝑥) par un polynôme de degré 1
(figure 1.) il vient :
𝑥1 𝑥0 +ℎ 1
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓(𝑥)𝑑𝑥 = ∫(𝑓0 + 𝛼∆𝑓0 )ℎ𝑑𝛼
𝑥0 𝑥0 0
𝑥−𝑥0
Avec : 𝛼 = ⟹ 𝑑𝑥 = ℎ𝑑𝛼
ℎ
Figure 1.
Analyse numérique I U.K.M.O Page 1
CHAPITRE III : INTEGRATION NUMERIQUE
l’intégrale est alors approchée par :
𝑥1
1
𝛼2 1 1 ℎ
∫ 𝑓(𝑥)𝑑𝑥 = ℎ [𝑓0 𝛼 + ∆𝑓0 ] = ℎ(𝑓0 + ∆𝑓0 ) = ℎ (𝑓0 + (𝑓1 − 𝑓0 )) = (𝑓0 + 𝑓1 )
2 0
2 2 2
𝑥0
Et l’erreur commise dans l’intervalle [𝑥0 , 𝑥1 ] sera :
1 1
𝛼(𝛼 − 1) 2 1 𝛼3 𝛼2
𝑒=∫ ∆ 𝑓0 ℎ𝑑𝛼 = ℎ∆2 𝑓0 [ − ]
2 2 3 2 0
0
∆2 𝑓0
Et on a : 𝑓 ′′ ≅ ⟹ ∆2 𝑓0 = ℎ𝑓 ′′ (𝑥)
ℎ
1
On obtient : 𝑒 = − 12 ℎ3 𝑓 ′′ (𝑥)
En répétant ce procédé pour les 𝑛 sous intervalle, on a finalement la formule des trapèzes :
𝑏 𝑛 𝑥 1 2 𝑥 𝑥𝑛 𝑥 ℎ
∫𝑎 𝑓(𝑥)𝑑𝑥 = ∫𝑥 𝑓(𝑥)𝑑𝑥 = ∫𝑥 𝑓(𝑥)𝑑𝑥 + ∫𝑥 𝑓(𝑥)𝑑𝑥 + ⋯ … … + ∫𝑥 𝑓(𝑥)𝑑𝑥 = 2 (𝑓0 +
0 0 1 𝑛−1
ℎ ℎ ℎ
𝑓1 ) + 2 (𝑓1 + 𝑓2 ) + ⋯ … . . + 2 (𝑓𝑛−1 + 𝑓𝑛 ) = 2 (𝑓0 + 2𝑓1 + 2𝑓2 + ⋯ … … … + 2𝑓𝑛−1 + 𝑓𝑛 )
III.2.1. L’erreur de la méthode :
On a trouver l’erreur pour un seul pas est :
1
𝑒 = − 12 ℎ3 𝑓 ′′ (𝑥) avec 𝑥0 ≤ 𝑥 ≤ 𝑥1
Donc l’erreur commise sur l’intervalle [𝑎, 𝑏] devient :
𝑛 3 ′′
𝑒=− ℎ 𝑓 (𝑥)
12
1 3 ′′ 1 𝑀𝑎𝑥
𝑒 = |− ℎ 𝑓 (𝑥)| = 𝑛ℎ3 |𝑓 ′′ (𝑥)|
12 12 𝑎≤𝑥≤𝑏
Exemple :
Calculer par la méthode des trapèzes la valeur approchée de l’intégrale I et l’erreur de la
1
méthode. On donne : 𝐼 = ∫0.5 𝑥𝑒 𝑥 𝑑𝑥 et ℎ = 𝑂. 1
Solution :
𝑥 0.5 0.6 0.7 0.8 0.9 1
𝑓(𝑥) 0.8244 1.0933 1.4096 1.7804 2.2136 2.7183
La valeur approchée de 𝐼 :
Analyse numérique I U.K.M.O Page 2
CHAPITRE III : INTEGRATION NUMERIQUE
1
ℎ
𝐼 = ∫ 𝑥𝑒 𝑥 𝑑𝑥 = (𝑦(0) + 2(𝑦(1) + 𝑦(2) + 𝑦(3) + 𝑦(4)) + 𝑦(5)) = 0.8286
0.5 2
L’erreur de la méthode :
𝑓(𝑥) = 𝑥𝑒 𝑥 ; 𝑓 ′ (𝑥) = 𝑥𝑒 𝑥 + 𝑒 𝑥 ; 𝑓 ′′ (𝑥) = (2 + 𝑥)𝑒 𝑥
𝑚𝑎𝑥|𝑓 ′′ (𝑥)| = |𝑓 ′′ (1)| = 8.1548
𝑛ℎ3
𝑒= 8.1548 = 0.003398
12
III.2.2. Programme matlab :
%Intégration numérique
%Méthode des trapèzes
clear all; clc
a=1;
b=2;
h=0.2;
n=(b-a)/h;
x=a:h:b
f=x.*exp(-(x).^3)
Itr=0
for i=1:n
Itr=(f(i)+f(i+1))*h/2+Itr;
end
Itr
III.3. Méthode de Simpson :
La formule de Simpson est obtenue en remplaçant 𝑓(𝑥) par un polynôme de degré 2
(la parabole, figure 2.) il vient :
𝑥0 +2ℎ 2 2
𝛼(𝛼 − 1) 2 𝛼2 𝛼3 𝛼2
∫ 𝑓(𝑥)𝑑𝑥 = ∫ (𝑓0 + ∆𝑓0 + ∆ 𝑓0 ) ℎ𝑑𝛼 = ℎ [𝑓0 𝛼 + ∆𝑓0 + ( − )∆2 𝑓0 ]
2 2 6 4 0
𝑥0 0
1 ℎ
= ℎ (2𝑓0 + 2∆𝑓0 + ∆2 𝑓0 ) = (𝑓0 + 4𝑓1 + 𝑓2 )
3 3
Analyse numérique I U.K.M.O Page 3
CHAPITRE III : INTEGRATION NUMERIQUE
Figur 2.
En répétant ce procédé pour les 𝑛 sous intervalle, on a finalement la formule de Simpson :
𝑏 𝑥𝑛
∫ 𝑓(𝑥)𝑑𝑥 = ∫ 𝑓(𝑥)𝑑𝑥
𝑎 𝑥0
𝑥2 𝑥4 𝑥𝑛
ℎ
= ∫ 𝑓(𝑥)𝑑𝑥 + ∫ 𝑓(𝑥)𝑑𝑥 + ⋯ … … + ∫ 𝑓(𝑥)𝑑𝑥 = (𝑓0 + 4𝑓1 + 𝑓2 )
3
𝑥0 𝑥2 𝑥𝑛−2
ℎ ℎ ℎ
+ (𝑓2 + 4𝑓3 + 𝑓4 ) + (𝑓4 + 4𝑓5 + 𝑓6 ) + ⋯ … + (𝑓𝑛−2 + 4𝑓𝑛−1 + 𝑓𝑛 )
3 3 3
ℎ
= (𝑓0 + 4𝑓1 + 2𝑓2 + 4𝑓3 + 2𝑓4 + 4𝑓5 + 2𝑓6 + ⋯ . . +2𝑓𝑛−2 + 4𝑓𝑛−1 + 𝑓𝑛 )
3
ℎ
= (𝑓0 + 4(𝑓1 + 𝑓3 + 𝑓5 + ⋯ . +𝑓𝑛−1 ) + 2(𝑓2 + 𝑓4 + 𝑓6 + ⋯ . +𝑓𝑛−2 ) + 𝑓𝑛
3
III.3.1. L’erreur de la méthode :
2 2 2
𝛼 ′ 𝛼 − 1)(𝛼 − 2) 3 ℎ∆3 𝑓0 ℎ∆3 𝑓0 𝛼 4
𝑒 = ℎ∫ ∆ 𝑓0 𝑑𝛼 = ∫(𝛼 3 − 3𝛼 2 + 2𝛼)𝑑𝛼 = [ − 𝛼3 + 𝛼2]
3! 6 6 4 0
0 0
ℎ∆3 𝑓0 16
= ( − 8 + 4) = 0
6 4
Donc l’erreur commise est à l’ordre suivant :
2
𝛼 ′ 𝛼 − 1)(𝛼 − 2)(𝛼 − 3) 4
𝑒 = ℎ∫ ∆ 𝑓0 𝑑𝛼
4!
0
2 2
ℎ∆4 𝑓0 ℎ∆4 𝑓0 𝛼 5 6𝛼 4 11𝛼 3
= ∫(𝛼 4 − 6𝛼 3 + 11𝛼 2 − 6𝛼)𝑑𝛼 = [ − + − 3𝛼 2 ]
24 24 5 4 3 0
0
4
ℎ∆ 𝑓0 4 ℎ∆4 𝑓0
= (− ) = −
24 15 90
Analyse numérique I U.K.M.O Page 4
CHAPITRE III : INTEGRATION NUMERIQUE
∆4 𝑓0
Et comme :𝑓 (4) (𝑥) ≅ on obtient :
ℎ4
ℎ5 𝑚𝑎𝑥
𝑒= |𝑓 (4) (𝑥)|
90 0 𝑥 ≤ 𝑥2
𝑥 ≤
et l’erreur commise sur l’intervalle [𝑎, 𝑏] devient :
𝑛 ℎ5 𝑚𝑎𝑥 (4) ℎ5 𝑚𝑎𝑥
𝑒 = ( ) 𝑥 ≤ 𝑥 ≤ 𝑥 |𝑓 (𝑥)| = |𝑓 (4) (𝑥)|
2 90 0 2 180 0 𝑥 ≤ 𝑥2
𝑥 ≤
Exemple :
Calculer par la méthode de Simpson la valeur approchée de l’intégrale I et l’erreur de la
1.1
méthode. 𝐼 = ∫0.5 𝑥𝑒 𝑥 𝑑𝑥 avec ℎ = 0.1
Solution :
𝑥 0.5 0.6 0.7 0.8 0.9 1 1.1
𝑓(𝑥) 0.8244 1.0933 1.4096 1.7804 2.2136 2.7183 3.3046
La valeur approchée de 𝐼 :
1.1 ℎ
𝐼 = ∫0.5 𝑥𝑒 𝑥 𝑑𝑥 = 3 (𝑦(0) + 4(𝑦(1) + 𝑦(3) + 𝑦(5)) + 2(𝑦(2) + 𝑦(4)) + 𝑓(6)) =6.7486
L’erreur de la méthode :
𝑓(𝑥) = 𝑥𝑒 𝑥 ; 𝑓 ′ (𝑥) = 𝑥𝑒 𝑥 + 𝑒 𝑥 ; 𝑓 ′′ (𝑥) = (2 + 𝑥)𝑒 𝑥 ; 𝑓 ′′′ (𝑥) = (3 + 𝑥)𝑒 𝑥 ;
𝑓 ′′′ ′(𝑥) = (4 + 𝑥)𝑒 𝑥
𝑚𝑎𝑥|𝑓 5 (𝑥)| = |𝑓 4 (1.1)| = 15.3212
𝑛ℎ5
𝑒= 15.3212 = 0.000078
180
III.3.2. Programme matlab :
%Intégration numérique
%Méthode de Simpson
clear all; clc
a=1;
b=2.6;
h=0.2;
n=(b-a)/h;
x=a:h:b
f=x.*exp(-(x).^3)
Itr=0
for i=1:2:n-1
Analyse numérique I U.K.M.O Page 5
CHAPITRE III : INTEGRATION NUMERIQUE
Itr=(f(i)+4*f(i+1)+f(i+2))*h/3+Itr;
end
Itr
Analyse numérique I U.K.M.O Page 6
