Calcul Scientifique : Intégration numérique
Imane ARMAILI
CPGE ELKHANSAA - Casablanca / MPSI
2024 - 2025
1 Introduction
2 Méthodes d’intégration
3 Méthodes d’intégration composites
4 Applications
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 2 / 21
Introduction
Introduction
√
• Exemple : f (x ) = cos(πx ) x 2 + 1 :
• Peut-on calculer I(f ) = 02 f (x )dx ?
R
• Pas d’expression analytique → méthodes numérique pour appro-
cher I(f ).
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 3 / 21
Introduction
Introduction
• L’objet de ce chapitre est de décrire quelques méthodes numé-
riques classiques permettant d’ évaluer des intégrales de fonctions
dont les valeurs sont connues en un nombre fini de points.
• Soit f une fonction réelle continue surR l’intervalle [a, b]. Le calcul
explicite de l’intégrale définie I(f ) = ab f (x )dx peut être difficile,
voire impossible.
Problème : Le calcul analytique peut être long, compliqué ou impos-
sible (pas de forme analytique).
• Généralement, il est préférable de calculer numériquement la va-
leur de l’intégrale I(f ).
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 4 / 21
Méthodes d’intégration
Méthodes d’intégration
Définition
On appelle formule de quadrature ou formule d’intégration numérique
toute formule permettant de calculer une approximation de I(f).
• Soit f : [a, b] → R une fonction continue donnée sur un intervalle
[a, b] ⊂ R.
• On se propose de calculer numériquement la quantité :
Z b
I(f ) = f (x )dx
a
• Méthodes : méthode des rectangles, méthode du point milieu et
méthode des trapèzes.
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 5 / 21
Méthodes d’intégration
méthode des rectangles
• L’approximation consiste à remplacer f (x ) par une constante (po-
lynôme de degré 0) sur l’intervalle [a, b]. On peut choisir soit la
constante f (a) (point à gauche), soit la constante f (b) (point à
droite) :
Rectangles à gauche (RAG) : I(f ) = f (a)(b − a)
Rectangles à droite (RAD) : I(f ) = f (b)(b − a)
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 6 / 21
Méthodes d’intégration
méthode du point milieu
• La méthode consiste à approcher la fonction f (x ) par la valeur
qu’elle prend au point milieu de l’intervalle [a, b]
I(f ) = f ( a+b
2
)(b − a)
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 7 / 21
Méthodes d’intégration
méthode des trapèzes
• L’approximation consiste à remplacer f (x ) par le polynôme de
degré 1 passant par les points (a, f (a)) et (b, f (b)).
• L’aire sous la courbe forme un trapèze :
f (a)+f (b)
I(f ) = 2
(b − a)
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 8 / 21
Méthodes d’intégration
méthode de Simpson
• L’approximation consiste à remplacer f (x ) par le polynôme de
degré 2 passant par les points (a, f (a)) , (b, f (b)) et (m, f (m))
avec m le milieu du segment [a,b].
I(f ) = (f (a) + 4f ( a+b
2
) + f (b)) (b−a)
6
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 9 / 21
Méthodes d’intégration composites
Méthodes d’intégration composites
• Soit f : [a, b] → R une fonction continue.
• On se propose de chercher des formules approchées pour l’inté-
grale I(f ) = ab f (x )dx .
R
• Pour cela, on choisit d’abord une subdivision de l’intervalle [a, b] :
x0 = a < x1 < x2 < ... < xn = b
• La formule de Chasles donne :
Z b n−1
X Z xi+1
I(f ) = f (x )dx = f (x )dx
a i=0 xi
• On est donc ramené au problème d’évaluer l’intégrale de f sur un
petit intervalle [xi , xi+1 ].
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 10 / 21
Méthodes d’intégration composites
• Ce calcul est effectué par des formules approchées appelées for-
mules de quadratures élémentaires (qui peuvent différer a priori
sur chacun des intervalles [xi , xi+1 ]).
Méthodes de quadrature élémentaires :
Z xi+1 li
X
I(f ) = f (x )dx ≈ (xi+1 − xi ) ωi,j f (ξi,j ),
xi j=0
Pli
où ξi,j ∈ [xi , xi+1 ], 0 ≤ j ≤ li et j=0 ωi,j = 1.
• Le problème est de choisir convenablement les points ξi,j et les
coefficients ωi,j de façon à minimiser l’erreur.
• Ceci se fera en général en évaluant l’intégrale xxi i+1 f (x )dx au
R
moyen d’une interpolation de f aux points ξi,j .
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 11 / 21
Méthodes d’intégration composites
La méthode de quadrature composée associée sera :
Z b n−1
X li
X
f (x )dx ≈ (xi+1 − xi ) ωi,j f (ξi,j )
a i=0 j=0
Définition
On dit qu’une méthode de quadrature (élémentaire ou composée) est
d’ordre N si : elle est exacte pour tout polynôme de degré inférieur ou
égal à N, elle est inexacte pour au moins un polynôme de degré N + 1.
Cas le plus simple : li = 0, quel que soit i :
• On choisit alors un seul point ξi ∈ [xi , xi+1 ] et on remplace f sur
[xi , xi+1 ] par le polynôme de degré 0 : p0 (x ) = f (ξi ). On a alors :
Z xi+1
f (x )dx ≈ (xi+1 − xi )f (ξi ),
xi
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 12 / 21
Méthodes d’intégration composites
• Ce qui donne :
Z b n−1
X
f (x )dx ≈ (xi+1 − xi )f (ξi ),
a i=0
• C’est-à-dire qu’on approxime l’intégrale par une somme de Rie-
mann relative à la subdivision (xi ). Voici les choix les plus cou-
rants :
• Méthode des rectangles à gauche (ξi = xi ) :
Z b n−1
X
f (x )dx ≈ (xi+1 − xi )f (xi ),
a i=0
• Méthode des rectangles à droite (ξi = xi+1 ) :
Z b n−1
X
f (x )dx ≈ (xi+1 − xi )f (xi+1 ),
a i=0
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 13 / 21
Méthodes d’intégration composites
• Schéma explicatif : (a = α, b = β, xi = αi )
Ces méthodes sont d’ordre 0.
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 14 / 21
Méthodes d’intégration composites
xi +xi+1
• Méthode du point milieu (ξi = 2
) :
Z b n−1
X xi + xi+1
f (x )dx ≈ (xi+1 − xi )f ( ),
a i=0 2
• Schéma explicatif : (a = α, b = β, xi = αi )
L’aire du rectangle coïncide avec l’aire du trapèze indiqué en grisé. La
formule approchée est donc exacte si f est une fonction affine, par
suite la méthode est d’ordre 1.
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 15 / 21
Méthodes d’intégration composites
Cas d’une interpolation linéaire : on choisit li = 1, ∀i, ξi,0 =
xi , ξi,1 = xi+1
• On remplace f sur [xi , xi+1 ] par la fonction linéaire p1 qui interpole
f aux points xi , xi+1 :
(x − xi )f (xi+1 ) − (x − xi+1 )f (xi )
p1 (x ) =
xi+1 − xi
• On obtient les formules suivantes, correspondant à la méthode
dite des trapèzes :
Z xi+1 Z xi+1
1 1
f (x )dx ≈ p1 (x )dx = (xi+1 − xi )( f (xi ) + f (xi+1 ))
xi xi 2 2
Z b n−1
X 1 1
f (x )dx ≈ (xi+1 − xi )( f (xi ) + f (xi+1 ))
a i=0 2 2
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 16 / 21
Méthodes d’intégration composites
• Schéma explicatif : (a = α, b = β, xi = αi )
L’ordre de cette méthode est 1.
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 17 / 21
Méthodes d’intégration composites
Cas d’une interpolation lagrangienne : on choisit li = 2, ∀i, ξi,0 =
xi , ξi,1 = xi+1 et ξi,2 = xi+2
• On remplace f sur [xi , xi+1 ] par un polynôme de degré 2 p2 qui
interpole f aux points xi , xi+1 et xi+2 :
p2 (x ) = f (xi+2 )L2 (x ) + f (xi+1 )L1 (x ) + f (xi )L0 (x )
• On obtient les formules suivantes, correspondant à la méthode
dite de Simpson :
Z xi+2 Z xi+2
h
f (x )dx ≈ p2 (x )dx = (f (xi ) + 4f (xi+1 ) + f (xi+2 ))
xi xi 3
Z b n−2 n−1
h X X
f (x )dx ≈ (f (x0 ) + f (xn ) + 2f (xi ) + 4f (xi )).
a 3 i=2 i=1
i pair i impair
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 18 / 21
Méthodes d’intégration composites
• Schéma explicatif : (a = x0 , b = x6 )
L’ordre de cette méthode est 2.
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 19 / 21
Applications
Exercice
√
Considérons la fonction suivante : f (x ) = cos(πx ) x 2 + 1. Nous
cherchons à déterminer
R2
une approximation de l’intégrale I(f ) définie
comme suit : I(f ) = 0 f (x )dx .
1) Écrire une fonction integre_RAG(f, a, b, nbr_points) qui per-
met de calculer I(f ) par la formule :
Z b n−1
X
f (x )dx ≈ (xi+1 − xi )f (xi )
a i=0
.
2) Écrire une fonction integre_RAD(f, a, b, nbr_points) qui per-
met de calculer I(f ) par la formule :
Z b n−1
X
f (x )dx ≈ (xi+1 − xi )f (xi+1 )
a i=0
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 20 / 21
Applications
3) Écrire une fonction integre_milieu(f, a, b, nbr_points) qui per-
met de calculer I(f ) par la formule :
Z b n−1
X xi + xi+1
f (x )dx ≈ (xi+1 − xi )f ( )
a i=0 2
4) Écrire une fonction integre_trapeze(f, a, b, nbr_points) qui
permet de calculer I(f ) par la formule :
Z b n−1
X 1 1
f (x )dx ≈ (xi+1 − xi )( f (xi ) + f (xi+1 ))
a i=0 2 2
5) Écrire une fonction integre_simpson(f, a, b, nbr_points) qui
permet de calculer I(f ) par la formule vue precedemment.
6) Comparer les résultats obtenus avec le résultat de la fonction
quad(f, a, b) du module [Link].
Imane ARMAILI (CPGE ELKHANSAA - Casablanca
Calcul /Scientifique
MPSI ) : Intégration numérique 2024 - 2025 21 / 21
