De Raaklijn Benadering
Het Geometrische Idee
Een gewone differentiaalvergelijking dy/dx = f(x,y) wijst aan elk punt in het (x,y)-vlak een helling toe — een richtingsveld. De werkelijke oplossing y(x) is een curve die overal die toegewezen hellingen volgt.
Euler's methode zet het continue richtingsveld om in een discrete wandeling:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
Van punt (xₙ, yₙ) gaat u afstand h langs de raaklijn. Bereik een geschat volgend punt. Herhaal.
Geometrische fout: de raaklijn op (xₙ, yₙ) heeft helling f(xₙ, yₙ), maar de werkelijke curve heeft op elk punt in het interval [xₙ, xₙ + h] een ander helling. De Euler-stap gebruikt de helling op het linkereindpunt overal — 'de helling die was'. De fout per stap groeit als h².
Geaccumuleerde Fout
Over N stappen om een vast eindpunt x = a te bereiken, met h = a/N:
- Lokale afkappingsfout per stap: O(h²)
- Aantal stappen: N = a/h
- Globale fout: O(h²) × (a/h) = O(h) — nauwkeurigheid van eerste orde
Euler's methode is van eerste orde: het halveren van h halveert de globale fout.
Euler's Methode Uitvoeren
Beschouw dy/dx = y, met beginvoorwaarde y(0) = 1. Werkelijke oplossing: y(x) = eˣ, dus y(1) = e ≈ 2,71828.
Pas Euler's methode toe met h = 0,5, van x = 0 tot x = 1 (2 stappen):
Stap 1: y₁ = y₀ + h·f(x₀, y₀) = 1 + 0,5·(1) = 1,5. Nieuw punt: (0,5, 1,5).
Stap 2: y₂ = y₁ + h·f(x₁, y₁) = 1,5 + 0,5·(1,5) = 2,25. Nieuw punt: (1,0, 2,25).
Euler geeft 2,25 versus werkelijke waarde 2,71828. Fout: 0,468. Relatieve fout: ~17%.
Afleiden van Euler's Stabiliteitsgebied
Voor de testvergelijking dy/dx = λy (waarbij λ een complex getal is), geeft Euler's methode:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
De versterkingsfactor per stap: z = 1 + hλ.
Stabiliteitsvoorwaarde: de berekende oplossing blijft begrensd als en slechts als |z| ≤ 1, d.w.z. |1 + hλ| ≤ 1.
Dit is een geometrische voorwaarde in het complexe hλ-vlak: het punt hλ moet binnen de cirkel met straal 1 liggen, gecentreerd op (-1, 0).
Euler's stabiliteitsgebied: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
Voor reële, negatieve λ (een dalende ODE zoals dy/dx = -2y): hλ moet in het interval (-2, 0] op de reële as liggen. Met λ = -2 en h = 0,5: hλ = -1. Dit ligt precies op de stabiliteitsgrens — de methode is marginaal stabiel, wat de kwalitatieve mislukking in het eerdere voorbeeld verklaart.
Met h = 1 en λ = -2: hλ = -2, wat buiten het stabiliteitsgebied ligt. De oplossing oscilleert met groeiende amplitudo.
Het Stabiliteitsgebied Vinden
Runge-Kutta 4 (RK4) heeft een groter stabiliteitsgebied dan Euler, wat een van de redenen is waarom het voor de meeste problemen de voorkeur heeft.
Voor reële negatieve λ staat RK4 toe dat hλ tot ongeveer -2,785 op de reële as (versus Euler's -2 limiet).
Voor stijve vergelijkingen met eigenwaarden λ van zeer verschillende groottes — zeg λ₁ = -1 en λ₂ = -1000 — vereist stabiliteit dat hλ₂ binnen het gebied blijft. Voor RK4 op de reële as: h·(-1000) ≥ -2,785, dus h ≤ 0,002785.
Deze minuscule stapsgrootte, bepaald door de stijfe eigenwaarde λ₂, maakt de simulatie duur ondanks dat de trage component λ₁ h = 2 zou kunnen gebruiken.
Vaste Punten & Attractiebakken
Euler's methode toegepast op dy/dx = f(y) defineert een discrete afbeelding: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
Een vast punt van deze afbeelding: y zodat g(y) = y. Voor Euler op dy/dx = f(y) voldoen vaste punten aan f(y) = 0 — evenwichten van de ODE.
Stabiliteit van een vast punt: als |g'(y)| < 1, convergeren nabijgelegen iteraties naar y. Als |g'(y*)| > 1 divergeren ze.
g'(y) = 1 + h·f'(y). Op een vast punt y: |1 + h·f'(y)| < 1 voor stabiliteit.
Dit is precies de Euler-stabiliteitsvoorwaarde met λ = f'(y*) — de linearisatie van de ODE op het evenwicht.
Attractiebakken: de verzameling beginvoorwaarden die onder de Euler-afbeelding naar y* convergeren. Voor niet-lineaire systemen bepaalt de bakgrens waar de simulatie betrouwbaar de ODE-evenwicht volgt versus divergeert naar een ander attractor.
De simulatiefeedbacklus is een discreet dynamisch systeem. Het kwalitatieve gedrag — convergentie, oscillatie, divergentie — hangt af van de stapsgrootte h ten opzichte van de geometrie van het richtingsveld van de ODE.
Geometrie Verbinden met Simulatieontwerp
De geometrie van numerieke simulatie komt neer op drie vragen:
1. Waar ligt het stabiliteitsgebied? Voor Euler: de schijf |1 + hλ| ≤ 1. Groter voor RK4. Onbegrensd (geheel linker halve vlak) voor impliciete methoden.
2. Waar liggen de eigenwaarden van de ODE? De eigenwaarden λ van de Jacobische matrix van f op elk punt bepalen welk stabiliteitsgebied hλ moet bevatten.
3. Welke h houdt alle hλ in het gebied? De maximaal toegestane h = (stabiliteitsgebiedstraal) / max|λ|.
Voor stijve systemen: max|λ| is enorm, wat tiny h voor expliciete methoden afdwingt. Impliciete methoden zijn duur per stap maar staan grote h toe — ze verhandelen kostenstap voor stabiliteit.
Hamming's inzicht vertaalt: de keuze van numerieke methode codeert een gok over de geometrie van het eigenwaardespectrum van de ODE. Die gok expliciet maken is de eerste ontwerpbeslissing in elke simulatie.