Approximation d'Euler

Approximation d'Euler

Méthode d'Euler

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.

Sommaire

Équation différentielle

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 xi + 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 yi + 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 yi + 1 à chaque itération si la fonction f\ est non-linéaire.

Intégration d'une fonction

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 aisément, 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

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}

Étant donné la fonction f(x)= \frac{x}{2} et des valeurs initiales x0 = 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(x)= \frac{x}{2} a pour primitive G(x)= \frac{x^2}{4} avec x0 = 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.

Lien externe

Augmenter n pour diminuer h et obtenir de bien meilleurs résultats :

voir Méthode d'Euler - MIAM

  • Portail des mathématiques Portail des mathématiques
Ce document provient de « M%C3%A9thode d%27Euler ».

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Approximation d'Euler de Wikipédia en français (auteurs)

Игры ⚽ Поможем написать курсовую

Regardez d'autres dictionnaires:

  • Euler number — For other uses, see Euler number (topology) and Eulerian number. Also see e (mathematical constant),Euler number (physics) and Euler–Mascheroni constant. In mathematics, in the area of number theory, the Euler numbers are a sequence En of… …   Wikipedia

  • Euler's formula — This article is about Euler s formula in complex analysis. For Euler s formula in algebraic topology and polyhedral combinatorics see Euler characteristic.   Part of a series of articles on The mathematical constant e …   Wikipedia

  • Euler spiral — A double end Euler spiral. An Euler spiral is a curve whose curvature changes linearly with its curve length (the curvature of a circular curve is equal to the reciprocal of the radius). Euler spirals are also commonly referred to as spiros,… …   Wikipedia

  • Euler-Maruyama-Schema — Das Euler Maruyama Verfahren, oft auch Euler Maruyama Schema oder stochastisches Euler Schema genannt, ist das einfachste Verfahren zur numerischen Lösung von stochastischen Differentialgleichungen. Es wurde erstmals in den 1950er Jahren durch… …   Deutsch Wikipedia

  • Euler-Maclaurin-Formel — Die Euler MacLaurin Formel oder Eulersche Summenformel (nach Leonhard Euler und Colin Maclaurin) ist eine mathematische Formel, die die Berechnung eines Integrals mit der Berechnung einer Summe von Stützstellen verbindet. Inhaltsverzeichnis 1… …   Deutsch Wikipedia

  • Euler — Leonhard Euler « Euler » redirige ici. Pour les autres significations, voir Euler (homonymie). Leonhard Euler …   Wikipédia en Français

  • Euler equations (fluid dynamics) — In fluid dynamics, the Euler equations govern inviscid flow. They correspond to the Navier Stokes equations with zero viscosity and heat conduction terms. They are usually written in the conservation form shown below to emphasize that they… …   Wikipedia

  • Euler's three-body problem — In physics and astronomy, Euler s three body problem is to solve for the motion of a particle that is acted upon by the gravitational field of two other point masses that are either fixed in space or move in circular coplanar orbits about their… …   Wikipedia

  • Euler–Maclaurin formula — In mathematics, the Euler–Maclaurin formula provides a powerful connection between integrals (see calculus) and sums. It can be used to approximate integrals by finite sums, or conversely to evaluate finite sums and infinite series using… …   Wikipedia

  • Euler method — In mathematics and computational science, the Euler method, named after Leonhard Euler, is a first order numerical procedure for solving ordinary differential equations (ODEs) with a given initial value. It is the most basic kind of explicit… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”