Amar D JEMA

[email protected]

13 février 2022
2
Table des matières

Table des matières 3

1 Intégration numérique 5
1.1 Méthode des trapèzes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2 Méthode de Simpson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
4
Chapitre 1

Intégration numérique

La plupart des fonctions n’admettent pas de primitives qui s’expriment à l’aide des fonctions usuelles.
Pour calculer une intégrale, on revient donc à la définition d’aire sous la courbe, aire que l’on approche,
en utilisant par exemple l’interpolation polynômiale (Polynôme de Lagrange, de Newton..).
Soit p n est le polynôme d’interpolation qui approxime une fonction f en n + 1 points xi (i = 0, .., n),
alors
Zxn xn
Z xn
Z
f (x)d x = p n (x)d x + E n (x)d x,
x0 x0 x0

où E n est l’erreur associée à l’interpolation.

1.1 Méthode des trapèzes

Commençons par la méthode la plus simple. Sur une subdivision [a, b], a = x0 < x1 < x2 < .. < xn = b
on approche la fonction par un segment. C’est-à-dire sur l’intervalle [x0 , x1 ] notre fonction f connue
seulement en deux points et est approchée par un segment (polynôme du 1er ordre), donc l’intégrale
devient l’aire du trapèze ((x0 , 0)(x1 , 0)(x1 , f (x1 ))(x0 , f (x0 ))), d’où

x1
x1 − x0
Z
f (x)d x ≃ ( f (x0 ) + f (x1 )).
2
x0

f (x1 )

f (x0 )

x0 x1 x

F IGURE 1.1 – Méthode du trapèze

6 Intégration numérique

On peut déduire facilement que

Zb xn
Z i=n−1
X xi+1 − xi
f (x)d x = f (x)d x ≃ S n = ( f (xi ) + f (xi+1 )).
i=0 2
a x0

Dans le cas où xi+1 − xi = h ∀i = 0, .., n − 1, on a

Zb Ã
i=n−1
!
h X
f (x)d x ∼ S n,h = f (x0 ) + f (xn ) + 2 f (xi ) .
2 i=1
a

La méthode de trapèze est d’ordre un et l’erreur est donnée par

¯ b ¯
¯Z ¯ 3
¯ f (x)d x − S n,h ¯ ≤ 1 (b − a) sup | f ′′ (x)|.
¯ ¯
¯ ¯ 12 n 2
x∈[a,b]
a
¯ ¯

R1
Exemple 1.1.1. calcul de la valeur approchée de x 2 d x (on en connait la valeur exacte 1/3) par ces
0
méthodes en subdivisant [0, 1] en 10 subdivisions (pas h = 1/10 et donc n = 10), donc xi = a+(b−a)/n =
i /10 pour i = 0, .., 10. On a f (xi ) = (i /10)2 donc la valeur approchée de l’intégrale est donnée

Z1 Ã
i=9
! Ã
i=9
!
2 1 X 2 1 X 2 67
x dx ∼ 1 + 2 (i /10) = 1 + 1/50 i = = 0.335.
20 i=1 20 i=1 200
0

1.2 Méthode de Simpson

Dans la méthode de Simpson, la fonction f est approchée par un polynôme
¡ x +x ¢ de second degré défi-
nissant un arc de parabole passant par les points d’ordonnées f (xi ) , f i 2 i +1 et f (xi+1 ) . Elle s’écrit

Zb xn
Z i=n−1 xi+1 − xi ³ ³x +x ´´
X i i+1
f (x)d x = f (x)d x ≃ S n = f (xi ) + f (xi+1 ) + 4f .
i=0 6 2
a x0

Si pour tout i on a xi+1 − xi = h alors

Zb
h i=n−1
f (x)d x ≃ S nh =
X ¡ ¢
f (a + i h) + f (a + (i + 1)h) + 4f (a + (i + 1/2)h) .
6 i=0
a

La méthode de Simpson est d’ordre 4 et l’erreur est donnée par

¯ b ¯
¯Z ¯ 5
¯ f (x)d x − S h ¯ ≤ 1 (b − a) sup | f (5) (x)|.
¯ ¯
n¯
¯ 180 n 4 x∈[a,b]
¯
a
¯

R1
Exemple 1.2.1. Calcul de la valeur approchée de x 2 d x, sachant que sa valeur exacte égale à 1/3 =
0
0.33333.... En subdivisant [0, 1] en n = 10 subdivisions (pas h = 1/10), donc

Z1 ¶2 ¶2 ¶2 ¶
1 i=9
µµ µ µ
2
X i i +1 i + 1/2 333
x dx ∼ + +4 = = 0.333.
60 i=1 10 10 10 1000
0
1.3 Exercices 7

1.3 Exercices
Exercice 1. Intégrer la fonction f (x) = e x dans l’intervalle [0, 1] en utilisant la méthode des trapèzes
composée avec 4 puis avec 8 intervalles. Comparer les résultats avec la valeur exacte.
1
Exercice 2. Intégrer la fonction f (x) = p dans l’intervalle [1, 9] en utilisant la méthode de Simpson
x
composée avec 6 intervalles. Comparer les résultats avec la valeur exacte.

Exercice 3. Soit
Z1
p
I= xd x.
0

1. Donner la valeur exacte de I .

2. Donner une approximation de I en utilisant la méthode de Simpson avec six sous-intervalles.
3. Le théorème du cours permettait-il de prédire l’erreur commise ?
8 Intégration numérique

R1
Corrigé de l’exo 1 : Calcul de : I = e x d x
0
1. (4 sous-intervalles)
On pose x0 = 0, x1 = 1/4, x2 = 1/2, x3 = 3/4 et x4 = 1. On a

Z1 x i +1
3 Z 3 f (xi+1 ) + f (xi ) X 3 f (x
X X i+1 ) + f (xi )
I4 = f (x)d x = f (x)d x ≃ (xi+1 − xi ) =
i=0 x i=0 2 i=0 8
0 i

à !
1 3 1
f (xi ) + f (x4 ) = (1 + e + 2(e 1/4 + e 1/2 + e 3/4 )) = 1.727221905
X
= f (x0 ) + 2
8 i=1 8

2. (8 sous-intervalles)
On pose xi = i h = i /8, i = 0, .., 8 avec le pas h = 1/8. On a

Z1 x i +1
7 Z 7 f (xi+1 ) + f (xi ) X 7 f (x
X X i+1 ) + f (xi )
I8 = f (x)d x = f (x)d x ≃ (xi+1 − xi ) =
i=0 x i=0 2 i=0 16
0 i

à !
1 7 1
(1 + e + 2(e 1/8 + e 1/4 + e 3/8 + e 1/2 + e 5/8 + e 3/4 + e 7/8 ))
X
= f (x0 ) + 2 f (xi ) + f (x8 ) =
16 i=1 16
I 8 = 1.720518592

3. I = 2.718281828, er 4 = 0.008940077, er 8 = 0.002236764 à 10−9

Corrigé de l’exo 2 :
R9 d x
Calcul de : I = p
1 x
1. (6 sous-intervalles)
On pose xi = 1 + i h = 1 + 4i /3, avec h = 4/3 i = 0, .., 6. On a xi+1 − xi = h = 4/3 On a

Z1 5 5 2³
X X xi + xi+1 ´
I6 f (x)d x = f (x)d x ≃ f (xi ) + 4f ( ) + f (xi+1 )
i=0 i=0 9 2
0

2 X5 X5 ³ xi + xi+1 ´
= (1 + 3 + 2 f (xi ) + 4 f( )) = 16.18547766
9 i=1 i=0 2

2. I = 17.33333333, er 6 = 1.14785567.
Corrigé de l’exo 3 :