Méthode d'Euler

Un article de Wikipédia, l'encyclopédie libre.
Aller à : navigation, rechercher

En mathématiques, la méthode d'Euler, nommée ainsi en l'honneur du mathématicien Leonhard Euler, est une procédure numérique pour résoudre par approximation des équations différentielles du premier ordre avec une condition initiale. C'est la plus simple des méthodes de résolution numérique des équations différentielles.

Équation différentielle[modifier | modifier le code]

La méthode d'Euler est une méthode numérique élémentaire de résolution d'équations différentielles du premier ordre, de la forme

\forall x\in I, u'(x)=f(x,u(x))

I\ est un intervalle de \mathbb R et f\ une fonction réelle sur I\times\mathbb R.

Étant donnée une condition initiale (a,u(a))\in I\times\mathbb R, la méthode fournit pour tout point b\in I une suite (u_n(b))_{n\in\mathbb N} d'approximations de la valeur u(b)\ que prend, lorsqu'elle existe, la solution de l'équation qui correspond à cette condition initiale. Divers jeux de conditions sur f\ peuvent assurer la convergence de cette suite.

u_n(b)\ s'obtient en calculant n\ valeurs intermédiaires (y_i)_{i\in[0,n]} de la solution approchée aux points (x_i)_{i\in[0,n]} régulièrement répartis entre a\ et b\ , donnés par

x_i = a + i\frac {b-a}{n}.

En étendant cette notation à x_0 = a,\ y_0 = u(a) et x_{n} = b,\ y_{n} = u_n(b) et en utilisant l'approximation de la dérivée

u'(x_i) \simeq \frac {u(x_{i+1})-u(x_i)}{x_{i+1}-x_i}

On en déduit la relation suivante :

\frac {y_{i+1}-y_i}{x_{i+1}-x_i}=f(x_i,y_i)

Les valeurs intermédiaires sont alors données par la relation de récurrence

y_{i+1} = y_i + (x_{i+1}-x_i)f(x_i, y_i),\ i\in[0,n-1].

qui est le schéma d'Euler explicite.

En remarquant, que l'on peut aussi approcher la dérivée en x_{i+1} par la même relation

u'(x_{i+1}) \simeq \frac {u(x_{i+1})-u(x_i)}{x_{i+1}-x_i}

on en déduit la relation de récurrence

y_{i+1} = y_i + (x_{i+1}-x_i)f(x_{i+1}, y_{i+1}),\ i\in[0,n-1].

qui est le schéma d'Euler implicite. On notera que dans ce schéma, le terme y_{i+1} apparaît des deux côtés de l'équation, ce qui contraint à utiliser des méthodes de résolution numérique du type de la relation de Newton-Raphson pour déterminer y_{i+1} à chaque itération si la fonction f\ est non-linéaire.

Intégration d'une fonction[modifier | modifier le code]

L'intégration d'une fonction continue sur un segment peut être vue comme un cas particulier où la fonction f\ est continue et ne dépend que de x\ . On démontre alors, en utilisant la continuité uniforme de f\ sur [a,b]\ (théorème de Heine), que la suite (u_n(b))_{n\in\mathbb N} est de Cauchy, et donc converge par complétude de \mathbb R.

Méthode d'Euler pour une fonction[modifier | modifier le code]

Pour calculer des valeurs approchées d'une primitive G de f sur I = [x0, xn], on divise I en n intervalles et on choisit h =\frac{x_n - x_0}{n}

Pour la valeur initiale y0 on a F(x0) = G(x0) = y0 ce qui permet de placer le premier point A0 (x0 ; y0).

Pour les n valeurs x1 = x0 + h, x2 = x1 + h, …, xn = xn-1 + h, on calcule de proche en proche, en relation avec la propriété de la dérivée citée ci-dessus, les n valeurs approchées F(x1), F(x2)…,  F(xn) de G.

En effet G est dérivable en x0 et G'(x0) = f(x0) :

F(x0 + h) = F(x0) + h f(x0) donc F(x1) = y0 + hG'(x0) \approx G(x1) ; soit y1 = y0 + h f(x0) et F(x1) = y1 \approx G(x1).

On recommence avec x1 :

F(x1 + h) = F(x1) + h f(x1) donc F(x2) = y1 + hG'(x1) \approx G(x1) + hG'(x1) \approx G(x2) ; soit y2 = y1 + h f(x1) et F(x2) = y2 \approx G(x2).

Puis y3 = y2+ h f(x2) = F(x3) \approx G(x3).

Et ainsi de suite n itérations jusqu'à yn = yn-1 + h f(xn-1) = F(xn) \approx G(xn).

Exemple : f(x)= \frac{x}{2}[modifier | modifier le code]

Étant donnée la fonction f:x\mapsto\frac{x}{2} et les valeurs initiales x_0=1 et y_0 = F(x_0)= \frac{1}{4}.

Le calcul des valeurs F(x1), F(x2), F(x3)… permet d'obtenir la représentation graphique de F par les segments [A0A1], [A1A2], [A2A3]…

Integration x div 2.gif

La fonction f a pour primitive G:x\mapsto\frac{x^2}{4} avec x_0 = 1 et y_0 = G(x_0)= \frac{1}{4}.

La courbe (C) représentative de G est ici placée sur le même graphe pour visualiser le calcul des tangentes.

La fonction affine est une approximation de la primitive G.

Voir aussi[modifier | modifier le code]

Articles connexes[modifier | modifier le code]

Lien externe[modifier | modifier le code]

Méthode d'Euler. Construction point par point d'une courbe intégrale avec GéoPlan. sur le site de P. Debart