Intégration Numériques
Quelques Exemples
Patrick Rabarison
prabarison@[Link]
Exemple 1. Méthode des trapèzes
Soit l’intervalle de référence [−1, +1] et soit φ une fonction définie
R1
sur cet intervalle. On cherche à approcher I (φ) = −1 φ(t)dt par
une formule de quadrature à deux points, soit
J(φ) = w 0φ(−1) + w 1φ(1) .
Notons p le polynôme (de degré inférieur ou égal à 1) tel que
p(−1) = φ(−1) et p(1) = φ(1). Il s’écrit :
1−t t+1
p(t) = φ(−1) + φ(1)
2 2
Nous avons ainsi obtenu w 0 = w 1 = 1.
Nous aurions aussi pu obtenir les deux coefficients par la deuxième
technique. Elle consiste à écrire que l’approximation doit être
exacte pour les monômes 1 et t , soit
Z 1
p(t)dt = w0 p(−1) + w1 p(1)
−1
En prennant p(t) = 1 puis p(t) = t, on obtient bien w0 = w1 = 1.
Méthodes des trapèzes sur un intervalle quelconque
Si maintenant nous nous plaçons sur un intervalle quelconque [a,
b], nous arrivons par un changement de variable à :
Z b Z 1
b−a b−a b+a
f (t)dt = φ(u)du avec φ(u) = f ( u −+ )
a 2 −1 2 2
Nous pouvons ainsi approcher cette dernière intégrale
Z
b−a 1 b−a b−a
φ(u)du ≈ (φ(−1) + φ(1)) = (f (a) + f (b))
2 −1 2 2
uer
Estimation de l’erreru
pour la méthode des trapèzes, on a
Z b
f ′′ (ζ(t))
E (f ) = (t − a)(t − b) dt,
a 2
or sur [a, b] le produit (t − a)(t − b) ne change pas de signe d’où il
existe η ∈ [a, b] tel que
Z b
f ”(η) h3 ′′
E (f ) = (t − a)(t − b)dt = f (ζ(η))
2 a 12
Exemple de formule avec Gauss-Legendre
Établissons la formule de Gauss-Legendre à deux points : d’une
part on a
g2 (t) = 12t 2 − 4 = 4(3t 2 − 1) ,
et donc comme racines
√ √
3 3
ζ1 = − , ζ2 = − .
3 3
D’autre par, le système pour déterminer les poids , on obtient:
w1 = w2 = 1
La formule est donc :
Z 1 √ √
3 3
f (t)dt ≈ J(f ) = f (− )+f( )
−1 3 3
Exercice. Vérifiez que cette formule est exacte pour p ∈ R3 [x].