Chapitre 3: Intégration numérique

Radhia Bessi

Décembre 2020

1
Plan

1 Introduction
Formule de rectangle
Formule du point milieu
Degré de précision
2 Formules de Newton-Cotes
Formule des trapèzes
Formule de Sipmson
Fomule de Newton-Cotes à (n + 1) points.

2
Introduction

L'objectif de l'intégration numérique

Z est de calculer une
b
approximation de l'intégrale I = f (x)dx d'une
a
fonction f continue sur un intervalle [a, b].
L'expression analytique de f peut être connue dans
certains cas, comme elle peut être inconnue. On
supposera dans ce chapitre que f est connue en
(n + 1) points x0 , ..., xn de [a, b].

3
La procédure de l'intégration numérique est est donc de chercher à
Z b
remplacer f (x)dx par une somme nie
a
Z b n
X
f (x)dx = λi f (xi ) + En (f ),
a i=1

où En (f ) est le terme de l'erreur commise en remplaçant I(f ) par

n
λi f (xi ).
X

i=1
De telles formules sont appelées formules de quadratures d'intégration.
Les xi et λi sont appelés respectivement les n÷uds et les poids de la
formule. n
La formule In (f ) = λi f (xi ) = I(f ) − En (f ) est dite
X

i=1
formule fermée si x0 = a et xn = b.
formule ouverte si x0 6= a ou xn 6= b.
formule exacte si I(f ) = In (f ).

4
Formule de rectangle
Formule de rectangle à gauche
La formule de rectangle simple (ou élémentaire) à gauche consiste à
remplacer f par la constante f (a) donnant
Z b
f (x)dx = f (a)(b − a) + E(f ).
a

Cette formule est exacte pour les fonctions polynômes constants

-La formule de rectangle à gauche composite sur un intervalle [c, d] consiste
à appliquer la formule élémentaire à une subdivision
a0 = c < a1 < ... < ap = d
Z d p−1 Z
X ai+1 p−1
X
f (x)dx ' f (x)dx = f (ai )(ai+1 − ai ).
c i=0 ai i=0
5
Formule du point milieu
-La formule du point milieu élémentaire consiste à remplacer f par la
2 ) donnant
constante f ( a+b
Z b
a+b
f (x)dx = f ( )(b − a) + E(f ).
a 2

Cette formule est exacte pour les fonctions polynômes de degré ≤ 1, mais
elle n'est pas exacte pour tout polynôme de de degré exactement 2.
On dit qu'elle est de degré de précision N = 1
-La formule de rectangle à gauche composite sur un intervalle [c, d] consiste
à appliquer la formule élémentaire à une subdivision
a0 = c < a1 < ... < ap = d
Z d p−1 Z ai+1 p−1
X X ai + ai+1
f (x)dx = f (x)dx = f( ))(ai+1 − ai ) + E(f ).
c i=0 ai i=0
2
6
Degré de précision

On dit qu'une formule de quadrature est de degré de

précision N si elle est exacte pour tout polynôme de
degré inférieur ou égal à N et non exacte pour au
moins un polynôme de degré exactement N + 1.
Z b Xn
Pour montrer que f (x)dx = λi f (xi ) + En (f ),
a i=1
est de degré N , il sut de vérier En (Pi ) = 0 pour
tout i = 0, 1, ..., N et En (PN +1 ) 6= 0, pour toute base
(P0 , P1 , P2 , ..., PN , PN +1 ) de PN +1 , avec degré de
PN +1 = N + 1.
où PN désigne l'espace vectoriel des polynômes réels de
degré inférieur ou égal à N .
7
Formules de Newton-Cotes

On considère a = x0 < x1 < ...xn = b une partition équidistante de [a, b].

(xi = a + i b−a n ). Dans la formule dite de Newton-Cotes, l'idée est
d'approcher f par son polynôme d'interpolation Pn aux points d'abscisses
x0 , x1 ..., xn . SI f = Pn + Rn , alors
Z b Z b Z b
f (x)dx = Pn (x)dx + Rn (x)dx = In (f ) + En (f ).
a a a

Ainsi : Z b Z b
f (x)dx ' In (f ) = Pn (x)dx.
a a

et l'erreur d'intégration est

Z b
En (f ) = Rn (x)dx.
a

8
Exemple : Formule des trapèzes (n=1)
On suppose que f est une fonction connue au moins en x0 = a et x1 = b.
On remplace f par le polynôme de degré 1 passant par
(x0 , f (x0 ) = (a, f (a)) et (x1 , f (x1 ) = (b, f (b)). On a f = P1 + R1 (x), avec
P1 (x) = f (x0 )L0 (x) + f (x1 )L1 (x) = f (a)L0 (x) + f (b)L1 (x), pour
x−b x−a
L0 (x) = , L1 (x) = .
a−b b−a
En intégrant sur [a, b] les polynômes L0 et L1 , trouve
Z b Z b
b−a
f (x)dx = (f (a) + f (b)) + R1 (f ).
a 2 a
On prend alors comme valeur approchée de l'intégrale est
b
b−a
Z
f (x)dx ' (f (a) + f (b)). Cette formule est appelée la formule de
a 2
quadrature élémentaire des trapèzes.

9
Erreur de la formule des trapèzes

On peut évaluer l'erreur d'intégration pour f une fonction C 2 de [a, b]. En

eet, d'après la formule d'erreur d'interpolation, pour tout x ∈ [a, b], il
existe ξ(x) ∈ [a, b] tel que f (x) = P1 (x) + f ”(ξ(x))
2! (x − a)(x − b). Donc
Z b Z b Z b
f ”(ξ(x))
f (x)dx = P1 (x)dx + (x − a)(x − b)dx.
a a a 2!

10
Rappel : Formule de la valeur moyenne

Soit f et g deux fonctions continues sur [a, b]. On suppose que g garde un
signe constant sur [a, b]. Alors il existe c ∈ [a, b] tel que
Z b Z b
f g = f (c) g.
a a

11
Dans la formule des trapèzes, on
Z b
f ”(ξ(x))
E(f ) = (x − a)(x − b)dx.
a 2!

Comme (x − a)(x − b) garde un signe constant sur [a, b] et f est une

fonction C 2 sur [a, b], d'après la formule de la moyenne, il existe η ∈ [a, b]
tel que
Z b Z b
f ”(ξ(x)) f ”(η) f ”(η)
E(f ) = (x−a)(x−b)dx = (x−a)(x−b)dx = − (b
a 2! 2! a 12

Ainsi, alors si f ∈ P1 , E(f ) = 0 et E(f ) 6= 0 si f (x) = x2 . La formule des

trapèzes est de degré de précision N = 2.

12
Formule de Simpson (n=2)

On remplace f par le polynôme d'interpolation de P2 , degré 2 passant par

a+b
les points (x0 , f (x0 )), (x1 , f (x1 ) et (x2 , f (x2 )), pour x0 = a, x1 = et
2
b b
b−a
Z Z
a+b
x2 = b. On a f' P2 = [f (a) + 4f ( ) + f (b)] Si on note
a a 6 2
h= b−a
2 , on obtient
b
b−a
Z
a+b
f (x)dx ' [f (a) + 4f ( ) + f (b)].
a 6 2

C'est la formule de quadrature de Simpson.

2 ) on a
Cette formule est au moins de degré 2. Pour P (x) = (x − a+b 3
Z b
P (x)dx = 0 et b−a
6 (f (a) 2 ) + f (b)) = 0 .
+ 2f ( a+b
a
E(x 7→ P (x) = (x − a+b 3
2 ) = 0, la formule est de degré ≤ 3.
E(x 7→ P (x) = (x − a+b 4
2 ) 6= 0. Donc la formule de simpson est
exactement de degré N=3.

13
Erreur d'intégration de la formule de Simpson
Soit P3 le polynôme de Hermite vériant
a+b a+b a+b a+b
P3 (a) = f (a), P3 (b) = f (b), P3 ( ) = f( ), et P30 ( ) = f 0(
2 2 2 2
Alors,
f ( 3)(ξ(x)) a+b 2
∀x ∈ [a, b] ∃ξ(x) ∈ [a, b]/f (x)−P3 (x) = (x−a)(x−b)(x− ) .
3! 2
Ainsi
b b b
f 4 (ξ(x))
Z Z Z
a+b 2
f (x)dx − P3 (x)dx = (x − a)(x − b)(x − ) dx.
a a a 3! 2
Comme le degré de P3 ≤ 3, donc
b
b−a
Z
a+b
P3 (x)dx = (P3 (a) + 4P3 ( ) + P3 (b)). Par suite
a 6 2
f (4) (η) b
Z
a+b 2
E2 (f ) = (x − a)(x − b)(x − ) dx.
6 a 2
En posant t = x−a
hZ et h = b−a
2 , on obtient
(4)
2
f (η) 4
E2 (f ) = 3! h t(t − 1)2 (t − 2)dt.
0
14
Formule de Newton-Cotes : Cas général

En utilisant la forme de Lagrange de Pn , il vient

Z b n
(n)
X
f (x)dx = λi f (xi ) + En (f ),
a i=0

Z b
où λi(n) = Li (x)dx et pour (L0 , L1 , ..., Ln ) est la base de
(n) (n) (n) (n)

a
Lagrange associée aux points a = x0 < x1 < ..., < xn = b pour

n , i = 0, ..., n).
xi = a + i b−a
Clairement cette formule est exacte pour les polynômes Li , i = 0, ..., n,
donc elle est de degré de précision N ≥ n.
Z b n
(n)
X
f (x)dx = λi f (xi ) + En (f ).
a i=0

est la formule de de quadrature de Newton-Cotes.

15
Evaluaton de l'erreur de quadrature de Newton-Cotes

On montre que la formule de Newton-Cotes à (n + 1) points est de degré de

précision N qui est égal :
si n est impair

N =n
N =n+1 si n est pair

Puisque qu'on
1 On suppose que f ∈ C (n+1) ([a, b]). Alors l'erreur de la formule de

NewtonCotes pour n impair est

n
hn+2 (n+1)
Z
En (f ) = f (η) t(t − 1)...(t − n)dt,
(n + 1)! 0

où η ∈ [a, b] et h = b−a
n .
2 On suppose que f ∈ C (n+2) ([a, b]). Alors l'erreur de la formule de
NewtonCotes pour n pair est

n
hn+3 (n+2)
Z
En (f ) = f (η) t2 (t − 1)...(t − n)dt
(n + 2)! 0

où η ∈ [a, b] et h = b−a
n .
Formule de Newton-Cotes composite

Elle consistent à découper l'intervalle d'intégration en sous-intervalles sur

lesquels la formule de quadrature de Newton-Cotes élémentaire (ou simple)
est appliquée

17