ArithmétiqueMéthode · Glossaire
méthode de Romberg
La méthode de Romberg est une technique d'intégration numérique qui combine des approximations successives par la règle des trapèzes, sur des maillages uniformes emboîtés, puis l'extrapolation de Richardson pour annuler les principaux termes d'erreur. Pour une fonction suffisamment régulière, elle fournit rapidement une estimation très précise de l'intégrale tout en réutilisant les calculs précédents.
Sommaire
Ce que vous allez apprendre
- Construire le tableau de Romberg à partir de raffinements trapézoïdaux emboîtés.
- Appliquer la formule d'extrapolation de Richardson avec le facteur correspondant à chaque colonne.
- Refaire un exemple exact pour l'intégrale de x puissance quatre entre 0 et 1.
- Reconnaître les situations où le manque de régularité ou les arrondis rendent la méthode peu fiable.
En clair
Imaginez une courbe dont vous voulez mesurer l'aire. La méthode des trapèzes fournit une première estimation, puis recommence en coupant l'intervalle en deux, en quatre, en huit morceaux. Chaque découpage réutilise les points déjà calculés.
Romberg ne se contente pas de garder l'estimation la plus fine. La méthode compare les résultats successifs et compense l'erreur qui décroît de façon prévisible. Une nouvelle comparaison peut ensuite retirer le terme d'erreur suivant. On obtient ainsi un tableau triangulaire dont la diagonale donne, dans les cas réguliers, des approximations de plus en plus précises.
Définition
La méthode de Romberg approche l'intégrale d'une fonction sur un intervalle fini. Pour un entier k positif ou nul, on note Tk la règle composite des trapèzes construite avec 2k sous-intervalles égaux. Les découpages sont emboîtés : passer de k à k + 1 conserve les anciennes évaluations et ajoute les milieux.
On place Tk dans la première colonne d'un tableau, sous le nom Rk,0. Pour chaque entier j compris entre 1 et k, la colonne suivante combine deux valeurs de la colonne précédente : . Cette extrapolation de Richardson annule successivement les termes en puissances paires du pas.
Lorsque la fonction est assez régulière pour que l'erreur des trapèzes possède ce développement asymptotique, la diagonale Rk,k converge souvent très vite. Cette accélération n'est toutefois pas une garantie universelle : une singularité, une discontinuité ou les erreurs d'arrondi peuvent détruire le comportement attendu.
Le principe
Soit f une fonction à intégrer sur l'intervalle [a, b]. Calculez d'abord les approximations trapézoïdales emboîtées T0, T1, …, où le niveau k utilise 2k sous-intervalles égaux, puis posez Rk,0 = Tk.
Pour les indices j allant de 1 à k, appliquez :
Le calcul s'arrête lorsque deux valeurs diagonales successives sont suffisamment proches pour la tolérance choisie, à condition que cette stabilisation ne soit pas seulement due aux arrondis numériques.
Quand l'utiliser
Le cadre naturel est une intégrale propre sur un intervalle fini, avec une fonction évaluable aux extrémités et aux points du découpage. Les niveaux doivent utiliser des subdivisions uniformes emboîtées, chacune divisant par deux le pas précédent.
L'accélération attendue demande surtout une fonction assez régulière pour que l'erreur de la règle des trapèzes se développe en puissances paires du pas. Plus on veut remplir de colonnes, plus il faut de régularité. Le résultat obtenu est une approximation, accompagnée en pratique d'un critère d'arrêt et d'un contrôle d'erreur.
Si l'intégrande présente une singularité à une borne, une discontinuité ou une oscillation que le maillage uniforme résout mal, les colonnes peuvent cesser de s'améliorer. Il faut alors scinder l'intervalle, changer de variable ou préférer une quadrature adaptative conçue pour ce comportement.
Un exemple, pas à pas
Approchons l'aire sous la courbe de la fonction f définie par f(x) = x4, entre 0 et 1. Le tableau de Romberg permet ici de voir chaque correction et d'obtenir une valeur que l'on peut contrôler exactement.
Données.
Intervalle : [0, 1].
Niveaux trapézoïdaux : 1, 2 puis 4 sous-intervalles égaux.
Valeur de contrôle : une primitive de x4 est x5/5.
Intervalle : [0, 1].
Niveaux trapézoïdaux : 1, 2 puis 4 sous-intervalles égaux.
Valeur de contrôle : une primitive de x4 est x5/5.
Étape 1. Avec un seul sous-intervalle, la règle des trapèzes donne T0 = (f(0) + f(1))/2 = 1/2.
Étape 2. En partageant ensuite l'intervalle en deux puis quatre parties, on obtient T1 = 9/32 et T2 = 113/512. La première extrapolation donne R1,1 = (4T1 − T0)/3 = 5/24, puis R2,1 = (4T2 − T1)/3 = 77/384.
Étape 3. La seconde extrapolation combine ces deux corrections :
Résultat et contrôle. Romberg donne exactement 1/5, soit 0,2, pour cet exemple. Le calcul direct confirme que l'intégrale vaut [x5/5] entre 0 et 1, donc 1/5. Cette exactitude après trois niveaux vient du polynôme choisi ; elle n'est pas promise pour toute fonction.
En pratique
Pour une fonction lisse sur un intervalle court, Romberg est utile lorsque chaque évaluation coûte cher : le raffinement dyadique réutilise tous les points anciens. Une quadrature adaptative devient préférable si l'erreur se concentre dans une petite zone identifiable.
Dans un programme, on construit le tableau ligne par ligne et on compare deux termes diagonaux successifs. On s'arrête lorsque leur écart passe sous une tolérance absolue ou relative choisie, tout en imposant un nombre maximal de niveaux.
Pour un calcul ponctuel avec peu de points, la méthode de Simpson peut suffire. Romberg devient intéressant quand plusieurs raffinements cohérents sont possibles et que les colonnes continuent réellement à réduire l'écart entre approximations.
À ne pas confondre
Méthode des trapèzes. Elle fournit la première colonne du tableau ; Romberg ajoute les extrapolations. Si le calcul ne combine jamais deux niveaux trapézoïdaux, il ne s'agit pas encore de Romberg.
Extrapolation de Richardson. C'est le procédé général qui élimine un terme d'erreur connu. Romberg en est une organisation particulière, fondée sur les trapèzes emboîtés et les facteurs 4j. Une extrapolation appliquée à un autre schéma n'est donc pas nécessairement une méthode de Romberg.
Méthode de Simpson. Sur trois points également espacés, la première extrapolation de Romberg coïncide avec la règle de Simpson. Les niveaux suivants les distinguent : Romberg poursuit le tableau, tandis que Simpson peut être utilisée seule ou sous forme composite.
Limites et pièges
Régularité insuffisante. Un angle, une discontinuité ou une dérivée non bornée peut empêcher l'erreur de suivre les puissances paires attendues. Le symptôme est une colonne qui stagne ou se dégrade ; il faut scinder l'intervalle ou employer une règle adaptée.
Fausse convergence. Deux diagonales presque égales ne prouvent rien si les arrondis dominent. Le seuil d'arrêt doit combiner une tolérance absolue et une tolérance relative ; si l'écart cesse de diminuer sur plusieurs niveaux, poursuivre les extrapolations n'ajoute plus de chiffres fiables.
Coût du raffinement. Le niveau k emploie 2k sous-intervalles et le niveau suivant en emploie 2k+1 : chaque nouvelle ligne double donc le maillage. Si une zone étroite cause l'erreur, un raffinement uniforme gaspille des évaluations ; une méthode adaptative les concentre au bon endroit.
Intégrales impropres ou très oscillantes. Appliquer directement le tableau peut produire des valeurs instables sans avertissement. Il faut d'abord traiter la singularité, tronquer avec une erreur contrôlée ou choisir une quadrature spécialisée. La convergence géométrique décrite dans le cas régulier n'est pas universelle.
Pour aller plus loin
La méthode des trapèzes détaille la première approximation et montre comment le raffinement du maillage construit la colonne de départ de Romberg.
La méthode de Simpson permet de comparer une règle quadratique directe avec la première extrapolation du tableau de Romberg.
Explorez les mathématiques autrement
Retrouvez nos magazines, podcasts et jeux pour explorer les mathématiques autrement.
Découvrir les offres