next up previous contents
Next: Equations différentielles Up: Intégration Previous: Formule de quadrature de

Méthodes de Richardson et Romberg

Ces méthodes visent à accélérer la convergence sans coût supplémentaire.

Soit V(h) la valeur de l'intégrale numérique obtenue pour le pas h. On a alors

displaymath2643

L'hypothèse de Richardson est que l'on peut faire une approximation de type Taylor, ce qui nous donne

displaymath2644

ce qui est vrai pour les trapèzes. On en déduit

displaymath2645

On introduit alors

displaymath2646

On est ainsi passé d'une précision en O(h) à une précision en tex2html_wrap_inline2665 . On peut réitérer le raisonnement. Si

displaymath2647

alors

displaymath2648

displaymath2649

Plus généralement, l'algorithme de Richardson est ainsi construit tex2html_wrap2681 On peut montrer que tex2html_wrap_inline2669 .

Dans le cas des trapèzes, ces hypothèses s'appliquent très bien car on a un théorème qui indique le comportement de la formule de quadrature.

Théorème : Si tex2html_wrap_inline2671 alors

displaymath2650

tex2html_wrap_inline2673 est la valeur fournie par la formule des trapèzes et les valeurs tex2html_wrap_inline2675 des constantes.

La méthode de Romberg utilise l'algorithme de Richardson pour tex2html_wrap_inline1665 ce qui donne tex2html_wrap2683



Jean-Michel JOLION
Fri Sep 10 13:28:20 METDST 1999