Méthodes de Runge-Kutta

Un article de Wikipédia, l'encyclopédie libre.
Aller à : navigation, rechercher
image illustrant l’analyse
Cet article est une ébauche concernant l’analyse.

Vous pouvez partager vos connaissances en l’améliorant (comment ?) selon les recommandations des projets correspondants.

Les méthodes de Runge-Kutta sont des méthodes d'analyse numérique d'approximation de solutions d'équations différentielles. Elles ont été nommées ainsi en l'honneur des mathématiciens Carl Runge et Martin Wilhelm Kutta lesquels élaborèrent la méthode en 1901.

Ces méthodes reposent sur le principe de l'itération, c'est-à-dire qu'une première estimation de la solution est utilisée pour calculer une seconde estimation, plus précise, et ainsi de suite.

Principe général[modifier | modifier le code]

Considérons le problème suivant :

que l'on va chercher à résoudre en un ensemble discret Plutôt que de chercher une méthode directe, les méthodes de Runge-Kutta proposent d'introduire les points intermédiaires afin de calculer par récurrence les valeurs avec

est le pas de temps et est dans l'intervalle . Pour chaque point intermédiaire, on note la pente correspondante

Ainsi, pour une solution exacte du problème, on a

On calculera ces intégrales par une méthode de quadrature, qu'on peut choisir différentes pour deux valeurs distinctes de i :

calculées ici pour

La méthode de Runge-Kutta d'ordre q sera donc donnée par :

On résume la méthode souvent par le tableau des différents poids de quadrature, appelé tableau de Butcher :

La méthode est consistante si

Exemples[modifier | modifier le code]

La méthode de Runge-Kutta d'ordre 1 (RK1)[modifier | modifier le code]

Cette méthode est équivalente à la méthode d'Euler, une méthode simple de résolution d'équations différentielles du 1er degré.

Considérons le problème suivant :

La méthode RK1 utilise l'équation

est le pas de l'itération. Le problème s'écrit donc :

La méthode de Runge-Kutta d'ordre 2 (RK2)[modifier | modifier le code]

La méthode RK2 du point milieu est une composition de la méthode d'Euler :

h est le pas de l'itération.

Elle consiste à estimer la dérivée au milieu du pas d'intégration :

et à refaire le pas d'intégration complet à partir de cette estimation :

Ce schéma est couramment appelé schéma prédicteur-correcteur explicite.

C'est le cas particulier pour de la méthode plus générale :

On reconnait ainsi que la méthode de quadrature utilisée pour les temps intermédiaires est celle du point milieu.

C'est une méthode d'ordre 2 car l'erreur est de l'ordre de h3.

Un autre cas courant est la méthode de Heun, correspondant au cas α=1. La méthode de quadrature repose sur la méthode des trapèzes.

La méthode de Runge-Kutta classique étage quatre (RK4)[modifier | modifier le code]

C'est un cas particulier d'usage très fréquent, noté RK4.

Considérons le problème suivant :

La méthode RK4 est donnée par l'équation :

L'idée est que la valeur suivante (yn+1) est approchée par la somme de la valeur actuelle (yn) et du produit de la taille de l'intervalle (h) par la pente estimée. La pente est obtenue par une moyenne pondérée de pentes :

  • k1 est la pente au début de l'intervalle ;
  • k2 est la pente au milieu de l'intervalle, en utilisant la pente k1 pour calculer la valeur de y au point tn + h/2 par le biais de la méthode d'Euler ;
  • k3 est de nouveau la pente au milieu de l'intervalle, mais obtenue cette fois en utilisant la pente k2 pour calculer y;
  • k4 est la pente à la fin de l'intervalle, avec la valeur de y calculée en utilisant k3.

Dans la moyenne des quatre pentes, un poids plus grand est donné aux pentes au point milieu.

La méthode RK4 est une méthode d'ordre 4, ce qui signifie que l'erreur commise à chaque étape est de l'ordre de h5, alors que l'erreur totale accumulée est de l'ordre de h4.

Ces formules sont aussi valables pour des fonctions à valeurs vectorielles.

La méthode de Runge-Kutta d'ordre quatre avec dérivée seconde[modifier | modifier le code]

Dans le cas

On en déduit et

Stabilité des méthodes[modifier | modifier le code]

En tant que méthodes à un pas et explicites, il convient de s'interroger sur leur stabilité. On peut montrer le résultat suivant[1] :

Supposons que f est k-lipschitzienne en y. Soit Les méthodes de Runge-Kutta sont stables sur , avec pour constante de stabilité avec

Références[modifier | modifier le code]

  1. Jean-Pierre Demailly, Analyse numérique et équations différentielles, EDP Sciences, coll. « Grenoble Sciences », (ISBN 2-86883-891-X)

Article connexe[modifier | modifier le code]