L'Approximation par Ligne Tangente
L'Idée Géométrique
Une équation différentielle ordinaire dy/dx = f(x,y) assigne une pente à chaque point du plan (x,y) — un champ de directions. La vraie solution y(x) est une courbe qui partout suit les pentes assignées.
La méthode d'Euler convertit le champ de directions continu en une marche discrète :
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
À partir du point (xₙ, yₙ), on se déplace de distance h le long de la ligne tangente. On arrive à un prochain point approximatif. On répète.
Erreur géométrique : la ligne tangente au point (xₙ, yₙ) a la pente f(xₙ, yₙ), mais la vraie courbe a une pente différente à chaque point de l'intervalle [xₙ, xₙ + h]. L'étape d'Euler utilise la pente au point final gauche partout — « la pente qui était ». L'erreur par étape augmente comme h².
Erreur Accumulée
Sur N étapes pour atteindre un point final fixe x = a, avec h = a/N :
- Erreur de troncature locale par étape : O(h²)
- Nombre d'étapes : N = a/h
- Erreur globale : O(h²) × (a/h) = O(h) — précision du premier ordre
La méthode d'Euler est du premier ordre : diviser h par deux divise l'erreur globale par deux.
Exécution de la Méthode d'Euler
Considérez dy/dx = y, avec condition initiale y(0) = 1. Vraie solution : y(x) = eˣ, donc y(1) = e ≈ 2,71828.
Appliquez la méthode d'Euler avec h = 0,5, de x = 0 à x = 1 (2 étapes) :
Étape 1 : y₁ = y₀ + h·f(x₀, y₀) = 1 + 0,5·(1) = 1,5. Nouveau point : (0,5, 1,5).
Étape 2 : y₂ = y₁ + h·f(x₁, y₁) = 1,5 + 0,5·(1,5) = 2,25. Nouveau point : (1,0, 2,25).
Euler donne 2,25 vs vraie valeur 2,71828. Erreur : 0,468. Erreur relative : ~17%.
Dérivation de la Région de Stabilité d'Euler
Pour l'équation de test dy/dx = λy (où λ est un nombre complexe), la méthode d'Euler donne :
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
Le facteur d'amplification par étape : z = 1 + hλ.
Condition de stabilité : la solution calculée reste bornée si et seulement si |z| ≤ 1, c'est-à-dire |1 + hλ| ≤ 1.
C'est une condition géométrique dans le plan complexe hλ : le point hλ doit se trouver à l'intérieur du cercle de rayon 1 centré en (-1, 0).
Région de stabilité d'Euler : { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
Pour λ réel et négatif (une EDO décroissante comme dy/dx = -2y) : hλ doit se trouver dans l'intervalle (-2, 0] sur l'axe réel. Avec λ = -2 et h = 0,5 : hλ = -1. C'est exactement sur la limite de stabilité — la méthode est marginalement stable, ce qui explique l'échec qualitatif dans l'exemple antérieur.
Avec h = 1 et λ = -2 : hλ = -2, nous mettant hors de la région de stabilité. La solution oscille avec une amplitude croissante.
Trouver la Limite de Stabilité
Runge-Kutta 4 (RK4) a une région de stabilité plus grande que celle d'Euler, ce qui est une raison pour laquelle il est préféré pour la plupart des problèmes.
Pour λ réel négatif, RK4 permet à hλ d'aller jusqu'à approximativement -2,785 sur l'axe réel (vs la limite -2 d'Euler).
Pour les équations raides avec valeurs propres λ à magnitudes très différentes — disons λ₁ = -1 et λ₂ = -1000 — la stabilité nécessite que hλ₂ reste à l'intérieur de la région. Pour RK4 sur l'axe réel : h·(-1000) ≥ -2,785, donc h ≤ 0,002785.
Cette minuscule taille de pas, dictée par la valeur propre raide λ₂, rend la simulation coûteuse même si la composante lente λ₁ pourrait utiliser h = 2.
Points Fixes & Bassins d'Attraction
La méthode d'Euler appliquée à dy/dx = f(y) définit une application discrète : yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
Un point fixe de cette application : y tel que g(y) = y. Pour Euler sur dy/dx = f(y), les points fixes satisfont f(y) = 0 — les équilibres de l'EDO.
Stabilité d'un point fixe : si |g'(y)| < 1, les itérés proches convergent vers y. Si |g'(y*)| > 1, ils divergent.
g'(y) = 1 + h·f'(y). À un point fixe y : |1 + h·f'(y)| < 1 pour la stabilité.
C'est exactement la condition de stabilité d'Euler avec λ = f'(y*) — la linéarisation de l'EDO à l'équilibre.
Bassin d'attraction : l'ensemble des conditions initiales qui convergent vers y* sous l'application d'Euler. Pour les systèmes non linéaires, la limite du bassin définit où la simulation suivra de manière fiable l'équilibre de l'EDO vs divergera vers un autre attracteur.
La boucle de simulation est un système dynamique discret. Son comportement qualitatif — convergence, oscillation, divergence — dépend de la taille de pas h par rapport à la géométrie du champ de directions de l'EDO.
Connecter la Géométrie à la Conception de Simulation
La géométrie de la simulation numérique se résume à trois questions :
1. Où se trouve la région de stabilité ? Pour Euler : le disque |1 + hλ| ≤ 1. Plus grand pour RK4. Non borné (demi-plan gauche entier) pour les méthodes implicites.
2. Où se trouvent les valeurs propres de l'EDO ? Les valeurs propres λ du jacobien de f à chaque point déterminent quelle région de stabilité doit contenir hλ.
3. Quel h garde tous les hλ à l'intérieur de la région ? Le h maximal autorisé = (rayon de la région de stabilité) / max|λ|.
Pour les systèmes raides : max|λ| est énorme, forçant un h minuscule pour les méthodes explicites. Les méthodes implicites sont coûteuses par étape mais permettent un h grand — elles échangent le coût par étape pour la stabilité.
L'intuition de Hamming se traduit : le choix de la méthode numérique encode un pari sur la géométrie du spectre de valeurs propres de l'EDO. Rendre ce pari explicite est la première décision de conception dans toute simulation.