Tangentlinjens approximation
Den geometriska idén
En ordinär differentialekvation dy/dx = f(x,y) tilldelar en lutning till varje punkt i (x,y)-planet — ett riktningsfält. Den sanna lösningen y(x) är en kurva som överallt följer dessa tilldelade lutningar.
Eulers metod konverterar det kontinuerliga riktningsfältet till en diskret promenad:
> (xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
Från punkten (xₙ, yₙ), flytta avståndet h längs tangentlinjen. Ankomst till en ungefärlig nästa punkt. Upprepa.
Geometriskt fel: tangentlinjen vid (xₙ, yₙ) har lutningen f(xₙ, yₙ), men den sanna kurvan har en olika lutning vid varje punkt längs intervallet [xₙ, xₙ + h]. Eulers steg använder lutningen vid den vänstra slutpunkten genomgående — 'den lutning som var.' Felet per steg växer som h².
Ackumulerat fel
Över N steg för att nå en fixerad slutpunkt x = a, med h = a/N:
- Lokalt trunkeringsfel per steg: O(h²)
- Antal steg: N = a/h
- Globalt fel: O(h²) × (a/h) = O(h) — första ordningens noggrannhet
Eulers metod är första ordningen: att halvera h halverar det globala felet.
Köra Eulers metod
Överväg dy/dx = y, med initialt villkor y(0) = 1. Sann lösning: y(x) = eˣ, så y(1) = e ≈ 2,71828.
Tillämpa Eulers metod med h = 0,5, från x = 0 till x = 1 (2 steg):
Steg 1: y₁ = y₀ + h·f(x₀, y₀) = 1 + 0,5·(1) = 1,5. Ny punkt: (0,5, 1,5).
Steg 2: y₂ = y₁ + h·f(x₁, y₁) = 1,5 + 0,5·(1,5) = 2,25. Ny punkt: (1,0, 2,25).
Euler ger 2,25 mot sant värde 2,71828. Fel: 0,468. Relativt fel: ~17%.
Härleda Eulers stabilitetsregion
För testekvationen dy/dx = λy (där λ är ett komplext tal), ger Eulers metod:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
Förstärkningsfaktorn per steg: z = 1 + hλ.
Stabilitetvillkor: den beräknade lösningen förblir begränsad om och endast om |z| ≤ 1, dvs |1 + hλ| ≤ 1.
Detta är ett geometriskt villkor i det komplexa hλ-planet: punkten hλ måste ligga innanför cirkeln med radie 1 centrerad vid (-1, 0).
Eulers stabilitetsregion: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
För verkliga, negativa λ (en förfallande ODE som dy/dx = -2y): hλ måste ligga i intervallet (-2, 0] på den verkliga axeln. Med λ = -2 och h = 0,5: hλ = -1. Detta är exakt på stabilitets gränsen — metoden är marginallt stabil, vilket förklarar det kvalitativa misslyckandet i det tidigare exemplet.
Med h = 1 och λ = -2: hλ = -2, vilket placerar oss utanför stabilitetsregionen. Lösningen oscillerar med växande amplitud.
Hitta stabilitets gränsen
Runge-Kutta 4 (RK4) har en större stabilitetsregion än Euler, vilket är en anledning till att den föredras för de flesta problem.
För verkligt negativ λ tillåter RK4 hλ ner till ungefär -2,785 på den verkliga axeln (kontra Eulers -2-gräns).
För styva ekvationer med egenvärden λ vid mycket olika storleksordningar — säg λ₁ = -1 och λ₂ = -1000 — kräver stabilitet att hλ₂ förblir innanför regionen. För RK4 på den verkliga axeln: h·(-1000) ≥ -2,785, så h ≤ 0,002785.
Denna lilla stegstorlek, dikterad av det styva egenvärdet λ₂, gör simuleringen dyr även om den långsamma komponenten λ₁ kunde använda h = 2.
Fixpunkter & attraktions bassänger
Eulers metod applicerad på dy/dx = f(y) definierar en diskret karta: yₙ₊₁ = g(yₙ) = yₙ + h·f(yₙ).
En fixpunkt för denna karta: y sådan att g(y) = y. För Euler på dy/dx = f(y) uppfyller fixpunkter f(y) = 0 — jämvikter för ODE:en.
Stabilitet för en fixpunkt: om |g'(y)| < 1 konvergerar närliggande iterationer till y. Om |g'(y*)| > 1 divergerar de.
g'(y) = 1 + h·f'(y). Vid en fixpunkt y: |1 + h·f'(y)| < 1 för stabilitet.
Detta är exakt Eulers stabilitetsvillkor med λ = f'(y*) — lineariseringen av ODE:en vid jämvikten.
Attraktion bassäng: mängden av initiala villkor som konvergerar till y* under Euler-kartan. För olinjära system definierar bassäng gränsen var simuleringen tillförlitligt spårar ODE-jämvikten mot divergerar till en annan attraktor.
Simuleringsslingan är ett diskret dynamiskt system. Dess kvalitativa beteende — konvergens, oscillation, divergens — beror på stegstorlek h i förhållande till ODE:ens riktnings fälts geometri.
Koppla geometri till simuleringsutveckling
Geometrin för numerisk simulering kommer ner till tre frågor:
1. Var är stabilitetsregionen? För Euler: skivan |1 + hλ| ≤ 1. Större för RK4. Obegränsad (hela vänster halv-planet) för implicita metoder.
2. Var ligger ODE:ens egenvärden? Egenvärden λ för Jacobian av f vid varje punkt bestämmer vilken stabilitetsregion som måste innehålla hλ.
3. Vilken h håller alla hλ innanför regionen? Största tillåtna h = (stabilitets region radie) / max|λ|.
För styva system: max|λ| är enorm, vilket tvingar liten h för explicita metoder. Implicita metoder är dyra per steg men tillåter stor h — de byter per-steg-kostnad för stabilitet.
Hammings insikt översätts: valet av numerisk metod kodifierar ett spel om geometrin för ODE:ens egenvärdes spektrum. Att göra det spelet explicit är det första designbeslutet i någon simulering.