Intégration numérique
1. La formule du point milieu (point médian)
En analyse numérique, la méthode du point médian est une méthode
permettant de réaliser le calcul numérique d'une intégrale
𝒃
𝒂+𝒃
∫ 𝒇(𝒙) 𝒅𝒙 = (𝒃 − 𝒂)𝒇 ( )
𝒂 𝟐
Pour la méthode du point milieu qui est d’ordre 0, l’erreur commise est
de la forme
(𝒃 − 𝒂)𝟐 𝒎𝒂𝒙
𝑬= ∗ 𝒙∈[𝒂,𝒃] 𝒇′(𝒙)
𝟐
Pour la méthode du point milieu qui est d’ordre 1, l’erreur commise est
de la forme
(𝒃 − 𝒂)𝟑 𝒎𝒂𝒙
𝑬= ∗ 𝒙∈[𝒂,𝒃] 𝒇"(𝒙)
𝟐𝟒
L'erreur est deux fois plus petite que celle donnée par la méthode des
trapèzes. Elle est exacte pour les polynômes de degré inférieur ou égal
à 1.
2. Méthode des trapèzes
En analyse numérique, la méthode des trapèzes est une méthode pour
𝑏
le calcul numérique d'une intégrale ∫𝑎 𝑓(𝑥)𝑑𝑥 s'appuyant sur
l'interpolation linéaire par intervalles. (Voir cours)
𝒏−𝟏
𝒃 − 𝒂 𝒇(𝒂) + 𝒇(𝒃)
𝑻𝒏 = ( + ∑ 𝒇(𝒙𝒊 ))
𝒏 𝟐
𝒊=𝟏
avec 𝑥𝑖 = 𝑎 + 𝑖ℎ
La méthode des trapèzes composite à m+1 points pour calculer
l’intégrale d’une fonction 𝑓 sur l’intervalle [a,b] s’écrit
𝒃 𝒎−𝟏
𝟏 𝟏
∫ 𝒇(𝒕)𝒅𝒕 = 𝒉 ( 𝒇(𝒂) ∑ 𝒇(𝒂 + 𝒊𝒉) + 𝒇(𝒃))
𝒂 𝟐 𝟐
𝒊=𝟏
𝑏−𝑎
avec ℎ =
𝑚
𝒃
𝒂+𝒃
∫ 𝒈(𝒕) ≈ ( ) (𝒈(𝒂) + 𝒈(𝒃))
𝒂 𝟐
Pour une fonction à valeurs réelles, deux fois continûment
différentiable sur le segment [a , b], l'erreur commise est de la forme
(𝒃 − 𝒂)𝟑 𝒎𝒂𝒙
𝑬=− ∗ 𝒙∈[𝒂,𝒃] 𝒇"(𝒙)
𝟏𝟐
Conformément aux expressions de l’erreur, la méthode des trapèzes
est souvent moins performante que celle du point milieu.
Exercice 3
21
On considère l’intégrale 𝐼 = ∫1 𝑑𝑥
𝑥
1. Calculer la valeur exacte de I.
2. Évaluer numériquement cette intégrale par la méthode des
trapèzes avec n égale à 3 sous-intervalles.
1. La valeur exacte est 𝐼 = [𝑙𝑛𝑥]𝑥=1
𝑥=2
= 𝑙𝑛2
2. La méthode des trapèzes composite à n points pour calculer
l’intégrale d’une fonction 𝑓 sur l’intervalle [a , b] s’écrit
3. Méthode de Gauss
La méthode de Gauss est une méthode de calcul numérique d'intégrales
à poids. Elle constitue une application directe de la théorie des
polynômes orthogonaux. Toutefois, elle a un intérêt théorique plus que
pratique.
4. Méthode des rectangles
(Voir cours)
5. Méthode de Simpson
La méthode de Simpson est une méthode de calcul approché
d'intégrale. Elle consiste en l'approximation suivante :
𝒂
𝒃−𝒂 𝒂+𝒃
∫ 𝒇(𝒕)𝒅𝒕 ≈ (𝒇(𝒂) + 𝟒𝒇 ( ) + 𝒇(𝒃))
𝒂 𝟔 𝟐
Cette formule est exacte pour tous les polynômes de degré inférieur
ou égal à 3 : on dit que la méthode de Simpson est d'ordre 3.
En général, pour appliquer cette méthode d'intégration, on découpe
(𝑏−𝑎)
l'intervalle [a,b] en n intervalles de longueur , et on applique la
𝑛
formule précédente sur chacun des sous-intervalles. En posant
𝑏−𝑎
ℎ=
𝑛
l’on a :
𝒃
𝟏 𝒂+𝒃
∫ 𝒈(𝒕) ≈ ( ) (𝒈(𝒂) + 𝟒𝒈 ( ) + 𝒈(𝒃))
𝒂 𝟔 𝟐
Si f est 4 fois continument différentiable sur [a, b], l'erreur
d'approximation commise est de la forme :
(𝒃 − 𝒂)𝟓 𝒎𝒂𝒙 𝟒
𝑬=− ∗ 𝒙∈[𝒂,𝒃] 𝒇 (𝒙)
𝟗𝟎
Cette expression du terme d'erreur signifie que la méthode de
Simpson est exacte (c'est-à-dire que le terme d'erreur s'annule) pour
tout polynôme de degré inférieur ou égal à 3.
Exercice 1
Quelle est l’équation de la parabole passant par les points (1,2), (3,7), (4,1).
Réponse
La parabole passant par les points (1,2), (3,7), (4,1) est donnée
par :
(𝑥 − 3)(𝑥 − 4) (𝑥 − 1)(𝑥 − 4) (𝑥 − 1)(𝑥 − 3)
𝑝2 (𝑥 ) = 2 +7 +1
(1 − 3)(1 − 4) (3 − 1)(3 − 4) (4 − 1)(4 − 3)
(𝑥 − 3)(𝑥 − 4) (𝑥 − 1)(𝑥 − 4) (𝑥 − 1)(𝑥 − 3)
𝑝2 (𝑥 ) = −7 +1
3 2 3
Puis développer, réduire et ordonner afin d’avoir la forme
𝑝2 (𝑥 ) = 𝛼𝑥 2 + 𝛽𝑥 + 𝑐
Exercice 3
Trouver le polynôme de l’espace vectoriel 𝑉𝑒𝑐{1 + 𝑥 2 , 𝑥 4 } qui interpole les points
(0,1) et (1,3).
Réponse
Trouvons un polynôme 𝑝(𝑥) qui soit combinaison linéaire des deux polynômes
assignés c'est-à-dire
𝑝(𝑥) = 𝛼(1 + 𝑥²) + 𝛽(𝑥 4 )
et qui interpole les deux points (0, 1) et (1, 3) :
