Recherche

opérationnelle

Pr.
Abdessamad
Kamouss
143

Intégration
Positionnement du
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites
Analyse Numérique

Pr. Abdessamad Kamouss

API - 2
ENSAM Casablanca
2021 - 2022
 Analyse numérique
Recherche
opérationnelle

Pr.
Abdessamad
Kamouss
144

Intégration
Positionnement du
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites
Chapitre V
Intégration Numérique
 Sommaire du chapitre
Recherche
opérationnelle

Pr.
Abdessamad
Kamouss
145

Intégration
Positionnement du
problème
Méthodes d’intégration
1 Positionnement du problème
simples
Méthodes d’intégration
composites
2 Méthodes d’intégration simples
3 Méthodes d’intégration composites
 1. Positionnement du problème

Recherche
opérationnelle

Pr.
Abdessamad
Kamouss Le but de ce chapitre est d’aborder le calcul général de l’intégrale
146
d’une fonction f(x) sur un domaine fini délimité par des bornes
Intégration
Positionnement du
finies a et b (les cas des bornes infinies n’est donc pas couvert
problème
Méthodes d’intégration
ici) :
simples
Méthodes d’intégration Z b
composites

I= f (x) dx
a
L’intégration numérique (ou quadrature) consiste à intégrer (de
façon approchée) une fonction sur un intervalle borné [a , b],
c’est-à dire calculer l’aire sous la courbe représentant la fonction, à
partir d’un calcul ou d’une mesure en un nombre fini n de points.
 1. Positionnement du problème

Recherche
opérationnelle

Pr.
Abdessamad Motivations :
Kamouss
147 Dans certains cas très limités, une telle intégrale peut être
Intégration calculée analytiquement (à la main). Cependant, ce n ?est que
Positionnement du
problème très rarement possible, et le plus souvent un des cas suivants
Méthodes d’intégration
simples se présente :
Méthodes d’intégration
composites
Le calcul analytique est long, compliqué et rébarbatif
Le résultat de l’intégrale est une fonction compliquée qui
fait appel à d’autres fonctions elles-même longues à
évaluer
Cette intégrale n’a pas d’expression analytique
...
 1. Positionnement du problème

Recherche
opérationnelle

Pr.
Abdessamad
Kamouss
148
Applications :
Intégration
Positionnement du Calcul probabiliste
problème
Méthodes d’intégration
simples
Traitement du signal
Méthodes d’intégration
composites
Mécanique
Électricité
Systèmes dynamiques
...
 1. Positionnement du problème

Recherche
opérationnelle
Afin de réaliser des approximations de l’intégrale de la
Pr. fonction f sur l’intervalle [a; b], on procède à une estimation
Abdessamad
Kamouss de l’aire délimité entre la courbe de cette fonction et l’axe des
149
abscisses sur cet intervalle.
Intégration
Positionnement du
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites
 1. Positionnement du problème

Recherche
opérationnelle

Pr.
Abdessamad
Afin de réaliser des approximations de l’intégrale de la
Kamouss fonction f sur l’intervalle [a; b], on procède à une estimation
150
de l’aire délimité entre la courbe de cette fonction et l’axe des
Intégration
Positionnement du
abscisses sur cet intervalle.
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites
 2. Intégration simple

Recherche
opérationnelle

Pr.
Abdessamad
Kamouss On note Pn un polynôme d’interpolation de f de degré n.
151
Alors
Intégration
Positionnement du Z b Z b
problème
Méthodes d’intégration
simples
I= f (x) dx ≈ Pn (x) dx
Méthodes d’intégration
a a
composites
A partir de cette approche, on peut développer plusieurs
méthodes :
3 méthodes des rectangles pour n = 0.
Méthode des trapèzes pour n = 1.
Méthode de Simpson pour n = 2.
 2. Intégration simple

Recherche
opérationnelle

Pr.
Abdessamad
Définition : Méthode de degré N
Kamouss
152 C’est une méthode pour laquelle la formule approchée est exacte
Intégration
pour tout polynôme de degré au plus N et inexacte pour au moins
Positionnement du
problème
un polynôme de degré N + 1.
Méthodes d’intégration
simples
Méthodes d’intégration
composites Exemple
Trouver A, B et C pour que la formule d’intégration numérique
suivante soit exacte pour tous les polynômes de degré au plus 2.
Z 1 1
I= f (x) dx ≈ A.f (0) + B.f + C.f (1)
0 3
Quel est le degré de précision de la formule obtenue ?
 2. Intégration simple

Recherche
opérationnelle
Méthode du rectangle à gauche
Pr. On note P0 le polynôme d’interpolation de f en a.
Abdessamad
Kamouss Alors P0 (x) = f (a)
153

Intégration
Positionnement du
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites

Z b Z b Z b
I= f (x) dx ≈ P0 (x) dx = f (a) dx
a a a
Donc
I ≈ f (a)(b − a)
 2. Intégration simple

Recherche
opérationnelle
Méthode du rectangle à droite
Pr. On note P0 le polynôme d’interpolation de f en b.
Abdessamad
Kamouss Alors P0 (x) = f (b)
154

Intégration
Positionnement du
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites

Z b Z b Z b
I= f (x) dx ≈ P0 (x) dx = f (b) dx
a a a
Donc
I ≈ f (b)(b − a)
 2. Intégration simple

Recherche
opérationnelle
Méthode du rectangle au milieu
a+b
Pr.
Abdessamad
On note P0 le polynôme d’interpolation de f en .
Kamouss   2
155 a+b
Alors P0 (x) = f
Intégration 2
Positionnement du
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites

Z b Z b Z b  
a+b
I= f (x) dx ≈ P0 (x) dx = f dx
a a a 2
Donc  
a+b
I≈f (b − a)
2
 2. Intégration simple

Recherche
opérationnelle
Méthode du trapèze
Pr. On note P1 le polynôme d’interpolation de f en a et b.
Abdessamad
Kamouss
f (b) − f (a) bf (a) − af (b)
156
Alors P1 (x) = x+
b−a b−a
Intégration
Positionnement du
problème
Méthodes d’intégration
simples
Méthodes d’intégration
composites

Z b Z b Z b
f (b) − f (a) bf (a) − af (b)
I= f (x) dx ≈ P1 (x) dx = x+ dx
a a a b−a b−a
Donc
b−a

I≈ f (a) + f (b) .
2
 2. Intégration simple

Recherche
opérationnelle
Méthode de Simpson
a+b
Pr.
Abdessamad
On note P2 le polynôme d’interpolation de f en x0 = a, x1 =
Kamouss
2
157 et x2 = b.
2
Intégration Alors P2 (x) = ∑ f (xi )Li (x)
Positionnement du
problème
i=0
Méthodes d’intégration
simples
Méthodes d’intégration
composites

Z b Z b
I= f (x) dx ≈ P2 (x) dx
a a
Donc   !
b−a a+b
I≈ f (a) + 4f + f (b) .
6 2
 3. Intégration composite

Recherche
opérationnelle
Introduction globale Afin de réduire l’erreur d’approximation
Pr. commise dans les méthodes d’intégration simples et généraliser
Abdessamad
Kamouss ces méthodes, on procède à la subdivision de l’intervalle [a, b] en
158
plusieurs sous intervalles :
Intégration
Positionnement du
problème
a = x0 < x1 < x2 < . . . < xn = b.
Méthodes d’intégration
simples
Méthodes d’intégration
composites

Par la suite on utilise le faite que :

Z b n−1 Z xi+1
I= f (x) dx = ∑ f (x) dx.
a i=0 xi
 3. Intégration composite

Recherche
opérationnelle
Méthodes des rectangles :
Pr. On considère une subdivision uniforme de l’intervalle [a, b] avec
Abdessamad
b−a
Kamouss
159
xi+1 − xi = h = , pour tout i ∈ {0, 1, . . . , n − 1}.
n
Intégration
On procède par la suite à une interpolation de f par un polynôme
Positionnement du
problème
de degré 0 sur chaque sous intervalle [xi , xi+1 ].
Méthodes d’intégration
simples
Méthodes d’intégration
composites
 3. Intégration composite

Recherche
opérationnelle

Pr.
Abdessamad Méthodes des rectangles composites à gauche :
Kamouss
160
Sur chaque intervalle [xi , xi+1 ], on interpole f par le polynôme
P0i (x) = f (xi ).
Intégration
Positionnement du Donc
problème
Méthodes d’intégration Z xi+1 Z xi+1
simples
Méthodes d’intégration
composites
f (x) dx ≈ f (xi ) dx = f (xi )(xi+1 − xi ) = hf (xi ).
xi xi

Par la suite :
n−1 n−1 n−1
I≈ ∑ hf(xi ) = h ∑ f(xi ) = h ∑ f(a + ih).
i=0 i=0 i=0
 3. Intégration composite

Recherche
opérationnelle

Pr.
Abdessamad Méthodes des rectangles composites à droite :
Kamouss
161
Sur chaque intervalle [xi , xi+1 ], on interpole f par le polynôme
P0i (x) = f (xi+1 ).
Intégration
Positionnement du Donc
problème
Méthodes d’intégration Z xi+1 Z xi+1
simples
Méthodes d’intégration
composites
f (x) dx ≈ f (xi+1 ) dx = f (xi+1 )(xi+1 − xi ) = hf (xi+1 ).
xi xi

Par la suite :
n−1 n−1 n
I≈ ∑ hf(xi+1 ) = h ∑ f(xi+1 ) = h ∑ f(a + ih).
i=0 i=0 i=1
 3. Intégration composite

Recherche
opérationnelle

Pr. Méthodes des rectangles composites au milieu :

Abdessamad
Kamouss [xi , xi+1 ], on interpole f par le polynôme
Sur chaqueintervalle 
162
xi + xi+1
P0i (x) = f .
Intégration
Positionnement du
2
problème Donc
Méthodes d’intégration
simples
Z xi+1 Z xi+1    
Méthodes d’intégration
xi + xi+1 xi + xi+1
composites
f (x) dx ≈ f dx = hf .
xi xi 2 2

Par la suite :
n−1   n−1   
xi + xi+1 1
I ≈ ∑ hf = h ∑ f a+ i+ h .
i=0 2 i=0 2
 3. Intégration composite

Recherche
opérationnelle
Méthode des trapèzes composites :
Pr. Sur chaque intervalle [xi , xi+1 ], on interpole f par le polynôme
Abdessamad
Kamouss
f (xi+1 ) − f (xi ) xi+1 f (xi ) − xi f (xi+1 )
163 P1i (x) = x+ .
h h
Intégration
Donc
Z xi+1 Z xi+1
Positionnement du
problème h
Méthodes d’intégration f (x) dx ≈ P1i (x) dx = (f (xi ) + f (xi+1 )) .
simples
Méthodes d’intégration
xi xi 2
composites
 3. Intégration composite

Recherche
opérationnelle
Méthode des trapèzes composites :
Pr.
Abdessamad Par la suite :
Kamouss
164 n−1
h
Intégration
I ≈ ∑ 2 (f(xi ) + f(xi+1 ))
Positionnement du
i=0
problème
Méthodes d’intégration
simples  n−1 n 
h
Méthodes d’intégration
composites ≈ 2 ∑ f(xi ) + ∑ f(xi )
i=0 i=1

 n−1 
h
≈ 2 f(a) + 2 ∑ f(xi ) + f(b)
i=1

En final,
n−1
h 
I≈ f(a) + 2 ∑ f(a + ih) + f(b) .
2 i=1
 3. Intégration composite

Recherche
opérationnelle Méthode de Simpson composite :
Pr. Sur chaque intervalle [xi , xi+1 ], on interpole f aux points xi ,
Abdessamad
Kamouss
xi + xi+1
165
ci = et xi+1 par le polynôme
2
Intégration
Positionnement du
problème
P2i (x) = f (xi )Li (x) + f (ci )Lci (x) + f (xi+1 )Li+1 (x).
Méthodes d’intégration
simples
Méthodes d’intégration
composites
Donc
R xi+1 R xi+1
xi f (x) dx ≈ x
i
P2i (x) dx 
h
≈ 6 f (x i ) + 4f (ci ) + f (x i+1 )
 3. Intégration composite

Recherche
opérationnelle Méthode de Simpson composite :
Pr. Par la suite :
Abdessamad
Kamouss
n−1  
166 h
I ≈ ∑ 6 f(xi ) + 4f(ci ) + f(xi+1 )
Intégration i=0
Positionnement du
problème
Méthodes d’intégration
 n−1 n−1 n 
h
simples
Méthodes d’intégration
≈ 6 ∑ f(xi ) + 4 ∑ f(ci ) + ∑ f(xi )
composites i=0 i=0 i=1

 n−1 n−1 
h
≈ 6 f(a) + 2 ∑ f(xi ) + 4 ∑ f(ci ) + f(b)
i=1 i=0

En final,

n−1 n−1    
h 1 
I≈ f(a) + 2 ∑ f(a + ih) + 4 ∑ f a + i + h + f(b) .
6 i=1 i=0 2
 3. Intégration composite

Recherche
opérationnelle Erreurs d’approximation :
Pr. Méthodes des rectangles à gauche / à droite :
Abdessamad
Kamouss
167
h 0
(∃η ∈ [a, b]) e(f ) = f (η)(b − a).
Intégration
2
Positionnement du
problème
Méthode des rectangles au milieu :
Méthodes d’intégration
simples
h2 00
Méthodes d’intégration
composites (∃η ∈ [a, b]) e(f ) = f (η)(b − a).
24
Méthode des trapèzes :
h2 00
(∃η ∈ [a, b]) e(f ) = − f (η)(b − a).
12
Méthode de Simpson :
h4
(∃η ∈ [a, b]) e(f ) = − f (4) (η)(b − a).
2880