Mathématiques
Précédent
Suivant
Autres méthodes et estimations d'erreur

Nous allons maintenant considérer une méthode dite à deux points, aussi appelée méthode des trapèzes.

Comme précédemment, si est la fonction à intégrer sur un intervalle , on découpe cet intervalle en sous-intervalles égaux et pose

est appelé le pas (de la subdivision).

Dans chaque intervalle de la subdivision, on considère la fonction affine qui prend respectivement les valeurs et pour et . En d'autres termes, on remplace la fonction à intégrer par une fonction polynomiale par morceaux. Ici les polynômes sont de degré 1 (donc définis par leurs valeurs en 2 points).

On obtient la formule des trapèzes en remplaçant l'intégrale de par celle de , d'un point de vue plus géométrique, on remplace l'aire sous le graphe de par la somme des aires des trapèzes de base et de hauteurs (à gauche) et à droite

On peut observer cela grâce à l'application suivante :

Plus précisement, on a

et

On veut maintenant estimer l'erreur commise avec cette méthode,

On cherche l'équation de la droite telle que

On obtient alors la fonction affine par morceaux définie sur l'intervalle par

(vérifier que la fonction coïncide avec aux points d'abscisses et ).

Pour tout , il existe tel que

On obtient alors une estimation d'erreur en considérant un majorant de la dérivée seconde de sur en supposant qu'il existe (ce qui est le cas si la fonction est deux fois continûment dérivable).

On a (par le changement de variables )

On obtient alors en recollant les morceaux

On constate donc que la formule des trapèzes est d'ordre 2.

Il existe de très nombreuses autres méthodes qu'on ne détaillera pas ici .Citons néanmoins la méthode de Simpson.

Méthode : La méthode de Simpson

On prend des polynômes d'interpolation de degré 2. Si

est le polynôme d'interpolation sur , on vérifie aisément que

(écrire et faire le calcul.)

on en déduit alors l'approximation

où l'on a posé

(Petit) miracle

La formule ( ) est encore exacte pour les polynômes de degré 3 et la formule de Simpson est d'ordre 3 .

On peut donner une majoration de l'erreur , pour des fonctions 4 fois continuement dérivables :

. La valeur précise de la constante (2880) n'est pas très significative car il est souvent très délicat d'estimer (qui suppose en outre beaucoup de régularité pour ). Pour une fonction régulière, il est néanmoins évident que la méthode de Simpson, d'ordre 4, converge plus rapidement que les méthodes précédentes d'ordre 2 (trapèzes, point milieu) pour un coût sensiblement identique.

Remarque : Lien entre les deux méthodes

On aurait pu également présenter la méthode des trapèzes comme une "combinaison linéaire" des méthodes des points droit et gauche ou bien la méthode de Simpson comme une combinaison de la méthode des trapèzes et du point mileu, en cherchant les coefficients pour lesquelles on obtient la meilleure formule de quadrature (pour une présentation dans cette direction, on lira le livre de D. Monasse, Math. et Informatique, Vuibert) . Nous avons préféré une approche liée aux polynômes d'interpolation.

Enfin, il existe d'autres formules de quadrature partant d'une problématique différente : comment placer les points de quadrature pour avoir la meilleure formule possible. Plus précisement, la méthode de Gauss consiste à chercher les points de quadrature de sorte que la formule soit exacte pour les polynomes de degré inférieur ou égal à p. Cette formule étant exacte pour "beaucoup" de polynômes , il est raisonnable d'espérer qu'elle ne devrait pas être mauvaise pour des fonctions.

En pratique, si on cherche l'exactitude pour des polynômes de degré 2, on obtient la formule de quadrature suivante :

et la formule d'erreur est comparable à celle de la méthode de Simpson

(on peut renvoyer à nouveau au livre de D. Monasse).

Légende :
Apprendre
S'évaluer
S'exercer
Observer
Simuler
Réalisé avec Scenari (nouvelle fenêtre)