Chapitre N°3
Chapitre N°3
Chapitre N°3
Intégration Numérique
4.1 Introduction
Les méthodes d’intégration numérique, interviennent essentiellement lorsque la pri-
mitive d’une fonction f est d’expression assez compliquée ou inconnue, ou lorsque f n’est
connue qu’en certains points distincts, par exemple si elle résulte de mesures expérimen-
tales. Dans ce chapitre, on traitera seulement les intégrales du type
Z b
I (f ) = f (x) dx; (1:4)
a
Où f est une fonction d’une variable réelle dé…nie et continue sur un intervalle [a; b]
de R borné non vide.
La plupart des formules d’intégration numérique proviennent des méthodes d’inter-
polation polynomiale, où le calcul de l’intégrale d’une fonction est approché par le calcul
exact de l’intégrale d’un polynôme, interpolant cette fonction aux points xi pour i = 0; ::n.
On obtient ainsi une forme générale appelée formule de quadratique numérique.
Z b X
n
f (x) dx ' In (f ) = if (xi ) (4:2)
a i=0
49
Les xi sont alors les points d’intégration et les coe¢ cients i sont appelés les poids
de la formule de quadrature.
Dé…nition 7 Une formule de quadrature est dite exacte sur un ensemble de fonctions E
si :
8f 2 E; R (f ) = I (f ) In (f ) = 0
Dé…nition 8 Une formule de quadrature est dite de degré de précision n si elle est exacte
pour xk ; k = 0; :::; n et non exacte pour xn+1 :
X
n
P (x) = Li (x) f (xi )
i=0
On peut alors espérer que si Pn est peu di¤érent de f , il en sera de même pour leurs
intégrales. Cette méthode assure évidemment que l’erreur est nulle si f est une fonction
polynomiale de degré n et cela quelque soit le choix des xi . L’intégrale du polynôme Pn
peut s’écrire en utilisant l’équation précédente sous la forme :
Z b X
n Z b
Pn (x)dx = f (xi ) Li (x) dx
a i=0 a
X
n
= f (xi ) ILi
i=0
Rb
Où ILi = a
Li (x) dx:
50
Remarque 5 On constate que les quantités ILi ne dépendent pas de la fonction, mais
elles dépendent seulement du choix des points xi .
Parmi ces formules, on trouve celle du trapèze et celle de Simpson qu’on va présenter
en détails par la suite.
(b a)
I1T (f ) = (f (a) + f (b)): (3.4)
2
51
Figure 4.1
52
Z b Z x1 Z x2 Z xn
I (f ) = f (x) dx = f (x) dx + f (x) dx + + f (x) dx
a x0 x1 xn 1
n 1Z
X xi+1
I (f ) = f (x) dx
i=0 xi
X1
n
xi+1 xi
I (f ) ' [f (xi ) + f (xi+1 )] I (f )
i=0
2
h X
n 1
I (f ) ' [f (xi ) + f (xi+1 )]
2 i=0
" #
h X
n 1
I (f ) ' InT (f ) = f (a) + f (b) + 2 f (xi ) (4:4)
2 i=1
Géométriquement, la valeur de l’intégrale I(f ) est approchée par la somme des aires des
Erreur d’approximation
Théorème 7 Soit f une fonction de classe C 2 sur l’intervalle [a; b], il existe alors dans
l’intervalle [a; b] tel que
1 00
In (f ) I (f ) = (b a)3 f ( )
12n2
où [a; b] est divisée en n sous intervalles égaux.
Alors, on peut écrire la borne supérieure de l’erreur commise comme suit :
1
I (f ) InT (f ) (b a)3 M2 (4:5)
12n2
00
avec M2 = M ax f (x) :
x2[a;b]
53
Figure 4.2
1
I (f ) I1T (f ) (b a)3 M2
12
54
Exemple 11 Z
2
I= sin xdx
0
Corrigé
R1
I (f ) = 0 sin xdx = [ cos x]02 = 1
Pour n = 1 :
(2 0)
I1T (f ) = 2
f (0) + f 2
I1T (f ) = 4
[0 + 1] = 0:785398:
Erreur exacte :
1
Eexa (f ) = I (f ) I1T (f )
= j1 0:785398j
= 0:21460:
Erreur théorique :
1 1 3
ETth (f ) = 12 2
0 M2
00
M2 = M ax f (x) = 1
x2[0; 2 ]
3
1
Donc, Eth (f ) = 8
12
1 = 0:32298:
1
on remarque que Eexa (f ) Eth (f )
55
Pour n = 4; on a :
" #
h X
n 1
h
I2T (f ) = f (a) + f (b) + 2 f (xi ) = [f (a) + f (b) + 2(f (x1 ) + f (x2 ) + f (x1 ))]
2 i=1
2
(b a)
h= 4
= 8
et x1 = x0 + h = a + h = 8 :; x2 = x1 + h = 8
+ 8
= 4 ; x 3 = x1 + h =
3
4
+ 8
= 8
:
I4T (f ) = 16
0 + 1 + 2(sin 8
+ sin 4
+ sin 38 = 0:98711:
Erreur exacte :
4
Eexa (f ) = I (f ) I1T (f )
= j0:3333 0:34375j
= 0:01288:
Erreur théorique :
4
ETth (f ) = 1
12n2
(b a)3 M2
1 3
= 12 42 2
0 1
= 0:06341:
4 4
) Eexa (f ) Eth (f ) :
(b a) a+b
I (f ) ' I2S (f ) = f (a) + 4f + f (b) ou bien (4:6)
6 2
h a+b (b a)
I (f ) ' I2S (f ) = f (a) + 4f + f (b) avec h =
3 2 2
56
4.2.4 Formule de Simpson généralisée
Z b Z x2 Z x4 Z xn
I (f ) = f (x) dx = f (x) dx + f (x) dx + + f (x) dx
a x0 x2 xn 2
n 2Z
X xi+2
I (f ) = f (x) dx
i=0 xi
h h
I (f ) ' [f (x0 ) + f (x2 ) + 4f (x1 )] + [f (x2 ) + f (x4 ) + 4f (x3 )]
3 3
h
+ + [f (xn 2 ) + f (xn ) + 4f (xn 1 )]
2 3 3
h 4 f (x0 = a) + f (xn = b) + 4 [f (x1 ) + f (x3 ) + + f (xn 1 )]
5
I (f ) '
3 +2 [f (x2 ) + f (x4 ) + + f (xn 2 )]
" #
h Xk X
k 1
I (f ) ' InS (f ) = f (a) + f (b) + 4 f (x2i 1 ) + 2 f (x2i ) ; n = 2k (4:7)
3 i=1 i=1
57
Figure 4.3
58
Figure 4.4
59
Erreur d’approximation
Théorème 8 Soit f une fonction de classe C 4 sur l’intervalle [a; b], il existe alors dans
l’intervalle [a; b] tel que
1
I (f ) Ins (f ) = (b a)5 f 4 ( )
180n4
où n = 2m représente les subdivisions de l’intervalle [a; b] qui doivent être égales.
Alors, on peut écrire la borne supérieure de l’erreur commise comme suit :
1
jI (f ) Ins (f )j (b a)5 M4 (5:4)
180n4
4
avec M4 = M ax f (x) :
x2[a;b]
R 1:8
Trouver la valeur approchée de 1
f (x) dx en utilisant la méthode de Simpson.
Solution
(b a)
On peut appliquer la formule de Simpson car n = 4 qui est pair, alors h = n
=
1:8 1
4
= 0:2 et xi = a + ih; i = 1; : : : ; 4:
h
I6S (f ) = 3
[f (x0 ) + f (x4 ) + 4 [f (x1 ) + f (x3 )] + 2 [f (x2 )]]
60
0:2
I6S (f ) = 3
[1:543 + 3:107 + 4 [1:811 + 2:577] + 2 [2:151]]
I6S (f ) = 1:7669:
61
Série d’exercices N 4
/Exercice 1 :
1- Déterminer par la méthode du trapèze (n = 1) puis par celle de Simpson (n = 2)
l’integrale :
Z 2
I= ln xdx
1
2- Comparer les résultats obtenus avec la valeur exacte de l’intégrale donnée, ainsi
que l’erreur théorique assosiée à chaque méthode.
Exercice 2 :
Un véhicule passe par un point A et cet instant on considère t = 0, ensuite on mesure
les premières 70 secondes sa vitesse V :
t (en s) 0 10 20 30 40 50 60 70
V (en m/s) 60 61:63 63:44 65:47 67:75 70:33 73:29 76:70
62
2- Trouver n le nombre de subdivisions nécessaires de l’intervalle d’intégration [0; ],
3
pour évaluer l’intégrale I à 10 près, grâce à la méthode des trapèzes et la méthode de
Simpson.
3- Calculer la valeur approximative de l’intégrale I en utilisant la méthode de Simpson
pour le n trouvé en question 2.
63
Bibliographie
[5] H. Grar, Méthodes numériques, cours, exercices corrigés, Université Setif 1, 2017.
[9] K, Mebarki, Analyse numérique, Cours, 2eme année mathématiques, Université Béjaia.
64