INTEGRATION NUMERIQUE

1
1- On lance une fusée verticalement du sol et l’on mesure
l’accélération pendant les premières 80s on obtient le tableau
suivant :
T(s) 0 10 20 30 40 50 60 70 80
ɤ (m/s2) 30 31,63 33,44 35,47 37,74 40,33 43,29 46,70 50,67

Calculer la vitesse de la fusée à l’instant t= 70s

2- calculer l’intégrale:
2 3/2
2  cos( 1  x ) 0 . 5 x

0 1  0 , 5 . sin( x )
e dx

2
Dans certaines situations, le calcul de l’intégrale
d’une fonction est difficile voir même impossible,
ceci est étant dû, soit à l’expression analytique
complexe de la fonction, ou du fait que
l’intégrande n’est pas donné sous forme analytique,
mais numériquement en un nombre de valeurs
discrètes.
D’où la nécessité d’approcher cette intégrale
numériquement.

3
Les méthodes numériques d'intégration d'une
fonction sont nombreuses et les techniques très
diverses.
Dans ce chapitre, on s’intéressera aux méthodes des
rectangles, des trapèzes et de Simpson.

4
 Méthode des rectangles
La méthode consiste à diviser l’intervalle [a,b] en n
parties et à approximer la surface de chaque tranche
par un rectangle.
y

y y y y … y
0 1 2 3 n-1

x
x x x x x
0 1 2 3 n-1

5
 La méthode des rectangles présente trois techniques de
résolution:

Point Point Point droit

gauche milieu

La seule différence réside sur la considération de x et de f(x)

pour chaque tranche
6
Somme à gauche:
Cette méthode consiste à approximer l’aire sous une
courbe à l’aide des sommes d’aires de rectangles, dont
les hauteurs sont calculées en utilisant les extrémités
gauches des sous-intervalles.
y
y
f (xi)
f (xi-1) y0 y1 y2 y3 … y
n-1

x0 x1 x2 x3 xn-1 x

x
xi-1 xi
7
Dans ce cas, pour chaque sous-intervalle de la forme [xi‑1, xi],
on construit un rectangle dont la hauteur correspond à f(xi-1).
L’intégrale est alors la somme des aires de tous les rectangles
ainsi construits. On suppose que les xi sont équidistants.
On écrit ainsi :
b n 1
A   f ( x ) dx  
i 1
Ai
a

 hf ( x 0 )  hf ( x 1 )   hf ( x n  1 )
n 1 n 1
 h  f ( x i )  h  f ( x 0  ih )
i0 i0

b  a
Avec h  x i  x i 1 
n 8
Méthode de droite :
Elle consiste à approximer l’aire sous une courbe à
l’aide des sommes d’aires de rectangles dont les
hauteurs sont calculées en considérant les extrémités
droites des sous-intervalles.
y y
f (xi)
f (xi-1)
y1 y2 y3 … … yn

x
x x0 x1 x2 x3 xn
xi-1 xi
9
Dans ce cas, pour chaque sous-intervalle de la forme [xi‑1, xi],
on construit un rectangle dont la hauteur correspond à f(xi).
L’intégrale est alors la somme des aires de tous les rectangles
ainsi construits. On écrit ainsi :
b n 1
A   f ( x)dx  
i 1
Ai
a

 h f ( x1 )  h f ( x 2 )   h f ( x n )
n n
 h  f ( xi )  h  f ( x 0  ih )
i 1 i 1

b  a
A vec h  x i  x i 1 
n 10
 Méthode du point mileu:
La méthode « du point milieu » consiste à approximer
l’aire sous une courbe à l’aide des sommes d’aires de
rectangles dont les hauteurs sont calculées en utilisant
les milieux mi des sous-intervalles.

y
y
f (mi)

f (mn)
f (m1)

f (m2)

f (m3)
x m1 m2 m3 mn x
xi-1 mi xi

11
Dans ce cas, pour chaque sous-intervalle de la forme [xi, xi+1],
on construit un rectangle dont la hauteur correspond à f(mi)
avec mi est le milieu des points xi et xi+1. L’intégrale est alors
la somme des aires de tous les rectangles ainsi construits. On
écrit ainsi :
b n 1
I   f ( x ) dx  
i 1
Ai
a

 hf ( m 1 )  hf ( m 1 )   hf ( m n 1 )
n 1 n 1
h
 h f (m i )  h  f ( x 0  ih  )
i 1 i0 2
b  a
Avec h  x i  x i 1 
n 12
 4
Calculons l’intégrale : 1 2

0
( x  1 ) dx
8

1 2
3

y  x 1
8 2

0 x4 1

4
1 2 0 1 2 3 4

A   ( x  1) dx
0
8 4
1 3 20 6,66667
A x x A
24 0
3
13
 3

Calcul par la méthode du

point gauche
2

1 2 1
y  x 1 0 x4
8
0 1 2 3 4

3 3
A  h 
i  0
f ( x i )  h 
i  0
f ( x 0  ih )

b  a 4  0
h   1
n 4
1 1 9
A  1    1   1  5 . 75
8 2 8

Très sous estimée !!!

14
 3
La méthode de point
droit:
2

1 2 0 x4
y  x 1 1

8
0 1 2 3 4
4 4
A  h 
i1
f (xi)  h 
i1
f (x 0  ih )

h  1
1 4 9 16
A   1   1   1   1  7 , 75
8 8 8 8
Trop surestimée !!!
15
 3

La méthode du point
milieu: 1

1 2 0 x4
y  x 1
8 0 1 2 3 4

3 3
h
A  h 
i1
f (m i)  h 
i1
f (x 0  ih 
2
)

h  1
1 1 1 9 1 25 1 49
A  .  1  .  1  .  1  .  1  6 , 625
8 4 8 4 8 4 8 4

16
 La Méthode des trapèzes
La méthode des trapèzes consiste à remplacer la
courbe f(x) par une ligne brisée et à calculer l’aire
de chaque trapèze. L’intégrale est alors approximée
en additionnant l’aire de chaque trapèze adjacent
permettant l’approximation de l’aire sous la courbe
d’une fonction f(x).

y0 y1 y2 y3 … … yn

x
x0 x1 x2 x3 xn
17
 La Méthode des trapèzes
Si Ai est l’aire du ieme trapèze, la hauteur de
chacun des trapèzes est donnée par h et les
bases sont données par yi-1, et yi, on aura alors :
b n
A  f (x)  
i 1
Ai
a

h n
h n 1

2
(f
i 1
i 1  fi )  ( f0  2 fi  fn )
2 i 1

18
 3

Reprenons l’exemple 2

précédent:
1

1 2
y  x 1 0 x4
8 0 1 2 3 4

1  9  1  9 3  1  3 17  1  17 
T  1              3 
2 8 28 2 22 8  2 8 
1  9 9 3 3 17 17 
T  1        3 
2 8 8 2 2 8 8 
1  27  27
T     6.75
2 2  4 19
 
0
x 2 sin( 2 x )dx

n=2
I = -1.4239 e-15
Exact = -4. 9348

20
 
0
x 2 sin( 2 x )dx

n=4
I = -3.8758
Exact = -4. 9348

21
 
0
x 2 sin( 2 x )dx

n=8
I = -4.6785
Exact = -4. 9348

22
 
0
x 2 sin( 2 x )dx

n = 16
I = -4.8712
Exact = -4. 9348

23
 Remarque:
n
h n 1
A   Ai  ( f 0  2 f i  f n )
i 1 2 i 1

h
 ( f 0  f1  f 2   f n 1 )  ( f1  f 2   f n ) 
2
1
 h ( f 0  f1  f 2   f n 1 )  h ( f1  f 2   f n ) 
2
1 n 1 n
 ( h  f i h  f i )
2 i 0 i 1

24
 Elle consiste à remplacer, entre trois points, la
fonction par l’arc de parabole passant par fi-1, fi
et fi+1. et on montre que :

h
Ai  ( f i 1  4 f i  f i 1 )
3

25
L’aire totale est obtenue en sommant toutes
les aires élémentaires Ai. On aura alors:

b n 1
A   f ( x)dx  i 1
Ai
a

h
  ( ( f i 1  4 f i  f i 1 ) )
3
26
27
 Ø Le calcul de l’intégrale est donné par:
x2 x4 x n 1
A  x0
f ( x ) dx  
x2
f ( x ) dx    
xn 3
f ( x ) dx
Ø Soit:

f ( x0 )  4 f ( x1 )  f ( x2 ) f ( x 2 )  4 f ( x3 )  f ( x 4 )
Ah h
3 3
f ( xn  2 )  4 f ( xn 1 )  f ( xn 1 )
  h
3

28
 Ø On obtient finalement

h  n 2 n3 
A   f ( x0 )  4  f ( xi )  2  f ( x j )  f ( xn1 )
3 i 1,i impaire j 2,i paire 
ba
Avec x0  a, xn  b, h 
n

29
 L’algorithme de la méthode de Simpson
• h=(b-a)/n
• S4=0
• Pour i=2 à n-1, pas=2 faire
S4==S4+f(a+(i-1)*h)

• FinPour
• S2=0
• Pour i=3 à n-2, pas=2 faire
S2==S2+f(a+(i-1)*h)

• FinPour
• Aire=(h/3.)*(f(a)+4*S4+2*S2+f(b))
• Fin 30
 Estimation des erreurs

• La méthode des des rectangle point milieu :

- En supposant que la fonction f est assez

régulière, l’erreur est donnée par l’expression:
b
ba n xi  xi 1 (b  a)3
|  f ( x)dx  
n i 1
f(
2
) |
24n 2
max(| f "( x) |
a

pour x  [ a, b]

31
 Estimation des erreurs

• La méthode des trapèzes:

- En supposant que la fonction f est assez

régulière, l’erreur est donnée par l’expression:
b
ba n (b  a)3
|  f (x)dx  
2n i1
( f (xi )  f (xi1)) |
12n 2
max(| f "(x) |
a

pour x [a, b]

32
 Estimation des erreurs

• La méthode des Simpson:

- En supposant que la fonction f est assez

régulière, l’erreur est donnée par l’expression:
b
ba n xi  xi 1
|  f ( x)dx  
6n i 1
f ( xi )  f (
2
)  f ( xi 1 ) |
a

(b  a)4 (4)
 4
max(| f ( x) |
120*24* n
pour x  [a, b] 33
 La formule de Newton-COTES
Pour une fonction f connue en un nombre de points
x i (i=0,1…..,n), on peut construire un polynôme
d’interpolation P n (x) passant par ces points et
l’intégrer. On suppose que ces points sont équidistants,
et on écrit:
b b n x i1

 f ( x ) dx   P n ( x ) dx  
i 0
 P n ( x ) dx
a a x i

b  a
où x0  a, h 
n
34
 La formule de Newton-COTES
Si on considère le polynôme Pn d’interpolation de Lagrange,
on a:
n
Pn ( x )  
i0
f ( xi ) Li ( x )

Le calcul d’intégrale est alors:

b b n
A   f ( x ) dx  P n ( x ) dx  
i0
 i f ( xi )
a a
b
Avec i  L
a
i ( x ) dx

35
 La formule de Newton-COTES

Si on considère le polynôme Pn d’interpolation de Lagrange,

on a: n
Pn ( x )  
i0
f ( xi ) Li ( x )
Le calcul d’intégrale est alors:
b b n
A   f ( x ) dx  P n ( x ) dx  
i0
 i f ( xi )
a a
b
Avec i  L
a
i ( x ) dx

C’est la formule de Newton-COTES

36
• Pour n=1, le polynôme de Lagrange s’écrit :
x  x1 x  x0
P1 ( x )  f ( x0 )  f ( x1 )
x 0  x1 x1  x 0
xa dx
On a : a  x 0 , b  x1, t  , dt  ; h  ba
ba h
x  a  t  0
   P1 ( t )  (1  t ) f ( a )  ( t ) f ( b )
x  b  t  1 
b b 1
a
f ( x ) dx  
a
P1 ( x ) dx  h  P1 ( t ) dt
0
1 1
 f ( a ) h  (1  t ) dt  f ( b ) h  tdt
0 0

2 1 2 1
t t h
 f ( a ) h (t  )  f (b ) h   f ( a )  f (b ) 
2 0
2 0
2 37
Pour n=2, le polynôme de Lagrange s’écrit :
( x  x 1 )( x  x 2 ) ( x  x 0 )( x  x 2 )
P2 ( x )  f ( x0 )  f ( x1 )
( x 0  x 1 )( x 0  x 2 ) ( x 1  x 0 )( x 1  x 2 )
( x  x 0 )( x  x 1 )
 f ( x2 )
( x 2  x 0 )( x 2  x 1 )
a  b
Soit x 0  a, x
2  b, x 1 
2
ba x  x1 dx
h  ,t  , dt 
2 h h
 x  x0  t  1

 x  x1  t  0
x  x  t  1
 2

t ( t  1) t ( t  1)
P2 ( t )  f ( x 0 )  (1  t 2 ) f ( x 1 )  f (x2 )
2 2 38
 t ( t  1) t ( t  1)
P2 ( t )  f ( x 0 )  (1  t 2 ) f ( x 1 )  f (x2 )
2 2
b 1 h 1

a
f ( x)dx  h  P2 (t )dt  f ( x0 )  t (t  1)dt
1 2 1
1 h 1
 f ( x1 )h  ( 1  t )dt  f ( x2 )  t (t  1)dt
2
0 2 1
3 2 1 3 1
h t t t
 f ( x0 ) (  )  f ( x1 )h(t  )
2 3 2 1 3 1
3 2 1
h t t
 f ( x2 ) (  )
2 3 2 1

h
f ( x )dx   f ( x0 )  4 f ( x 1 )  f ( x 2 )
b
a 3 39