Intégration numérique

Z b
Introduction Soit à calculer I = f (x)dx. Le calcul de l'intégrale d'une fonction n'est pas toujours facile : En eet, si
a
on a l'expression explicite de la fonction f , il est parfois dicile, voire impossible, de trouver sa primitive. Et si la fonction
f n'est dénie que dans un nombre ni de points ?
Peut-on-calculer son intégrale ?, sinon, peut-on avoir une estimation de son intégrale ? ? ? ?

Réponse :
Grâce aux Techniques d'intégration numérique, on peut, donner une approximation d'une fonction Connue juste en un
nombre ni de points ou une fonction explicitement donnée et dont il est dicile, voire, impossible de calculer sa primitive.

Idée Supposons que nous connaissons les valeurs de f en un nombre ni de points (xi , i = 0, . . . , n).

L'idée naturelle est d'utiliser le polynôme d'interpolation Pn de f aux points (xi ) (xi sont appelés points d'appui). On
écrit alors :

f (x) = Pn (x) + ε(x) où Pn (x) est le polynôme d'interpolations de f aux points xi , i = 0, . . . , n et ε(x) est l'erreur
d'interpolation. On a
Z b Z b Z b
f (x)dx = Pn (x)dx + ε(x)dx
a a a
Z n
bX Z b
= Lj (x)f (xj )dx + ε(x)dx
a j=0 a

n
!
X Z b Z b
= Lj (x)dx f (xj ) + ε(x)dx
j=0 a a

Xn Z b
= ωj f (xj ) + ε(x)dx
j=0 | a {z }
E : erreur d'intégration

n
X
≈ ωj f (xj )
j=0

Z β k−1
Exemple : On cherche à approcher f (x)dx par une série de la forme ωj f (xj ) avec [α, β] = [αi , αi+1 ]. On
X [
α j i=0
choisit ti ∈ [αi , αi+1 ] et on approxime f (x) ≈ P0 (x) = f (ti ), ∀x ∈ [αi , αi+1 ].
Z αi+1 Z αi+1 Z αi+1
f (x)dx ≈ P0 (x)dx = f (ti ) dx = (αi+1 − αi )f (ti )
αi αi αi

Z β k−1
X Z αi+1 k−1
X
f (x)dx = f (x)dx = (αi+1 − αi )f (ti )(Somme de Rieman)
α i=0 αi i=0

Les choix les plus courants sont :

La méthode des rectangles à gauche :

Z β k−1
ti = αi donc (αi+1 − αi )f (αi ).
X
f (x)dx =
α i=0

La méthode des rectangles à droite :

Z β k−1
ti = αi+1 donc (αi+1 − αi )f (αi+1 )).
X
f (x)dx =
α i=0
La méthode du point milieu :
Z β k−1
αi+1 + αi
ti = αi+1 donc .
X
f (x)dx = (αi+1 − αi )f (
α i=0
2

1
Méthode des trapèzes Sur [αi , αi+1 ], on interpole f (x) par un polynôme de degré 1 aux points αi , αi+1 .

Z αi+1 Z αi+1
f (x)dx ≈ P1 (x)
αi αi
Z αi+1 1
X
≈ Lj (x)f (xj ) (x0 = αi , x1 = αi+1 )
αi j=0
αi+1  
x − αi+1 x − αi
Z
≈ f (αi ) + f (αi+1 ) dx
αi αi − αi+1 αi+1 − αi
α α
(x − αi+1 )2 i+1 f (αi+1 ) (x − αi )2 i+1
 
f (αi )
≈ +
αi − αi+1 2 αi αi+1 − αi 2 αi
f (αi ) (αi − αi+1 )2 f (αi+1 ) (αi+1 − αi )2
≈ +
αi+1 − αi 2 αi+1 − αi 2
2
f (αi ) (αi+1 − αi )
≈ − (αi − αi+1 ) + f (αi+1 )
2 2

β k−1
X Z αi+1 k−1
(αi+1 − αi ) αi+1 − αi
Z X
f (x)dx = f (x)dx ≈ f (αi ) + f (αi+1 )
α i=0 αi i=0
2 2
Dénition :
Nous dirons qu'une méthode d'intégration numérique est d'ordre N si elle est exacte pour tout polynôme de degré ≤ N .
Par exemple, pour la méthode d'intégration des trapèzes l'ordre est égale à 1.

0.1 Exercices

Exercice I:
Z π/2
1)Déterminer par la méthode des trapèzes f (x)dx sur la base des données suivantes :
0
f (0) = 0 ; f (π/8) = 0.382683 ; f (π/4) = 0.707107 ; f (3π/8) = 0.923880 et f (π/2) = 1.
2) Ces points d'appui sont ceux donnant sinx, comparer alors les résultats obtenus avec la valeur exacte.

Exercice II:
1)Écrire le polynôme de Lagrange P (x) d'une fonction f donnée aux points : −1 ; −1/3; 1/3 et 1.
2) Par intégration du polynôme obtenu, déduire la formule d'intégration approchée suivante :
Z 1
1 3 1 3 1 1
f (x) ' f (−1) + f (− ) + f ( ) + f (1)
−1 4 4 3 4 3 4
.

Exercice III:
√
Z 1
Soit f (x) = 1 + x, et I = f (x)dx.
0
1) En utilisant les formules de quadratures, donner une approximation de I .
2) Calculer l'erreur d'intégration numérique.