Aproksymacja linii tangensjalnej
Podstawowa idea geometryczna
Różniczkowa równanie zwyczajne dy/dx = f(x,y) przypisuje nachylenie do każdego punktu w płaszczyźnie (x,y) — pole kierunków. Prawdziwa rozwiązanie y(x) to krzywa, która wszędzie śledzi te przypisane nachylenia.
Metoda Eulera przekształca ciągłe pole kierunków w ciągły spacer:
>(xₙ, yₙ) → (xₙ + h, yₙ + h·f(xₙ, yₙ))
Od punktu (xₙ, yₙ), przesuwaj się odległość h wzdłuż linii tangensjalnej. Dokonaj dokładnego następnego punktu. Powtórz.
Błąd geometryczny: linia tangensjalna w punkcie (xₙ, yₙ) ma nachylenie f(xₙ, yₙ), ale prawdziwa krzywa ma inne nachylenie na każdym punkcie w przedziale [xₙ, xₙ + h]. Krok metody Eulera używa nachylenia na lewym końcu przez cały czas — 'nachylenie, które było'. Błąd na krok wzrasta w miarę jak h².
Zbiorczy błąd
W ciągu N kroków do osiowania punktu stałego x = a, z h = a/N:
- Lokalny błąd obcięcia na krok: O(h²)
- Liczba kroków: N = a/h
- Błąd globalny: O(h²) × (a/h) = O(h) — dokładność pierwszego rzędu
Metoda Eulera ma dokładność pierwszego rzędu: podzielając h, podzielony jest błąd globalny.
Uruchamianie metody Eulera
Rozważ równanie dy/dx = y, z warunkiem początkową y(0) = 1. Prawdziwe rozwiązanie: y(x) = eˣ, więc y(1) = e ≈ 2.71828.
Zastosuj metodę Eulera z h = 0.5, od x = 0 do x = 1 (2 kroki):
Krok 1: y₁ = y₀ + h·f(x₀, y₀) = 1 + 0.5·(1) = 1.5. Nowy punkt: (0.5, 1.5).
Krok 2: y₂ = y₁ + h·f(x₁, y₁) = 1.5 + 0.5·(1.5) = 2.25. Nowy punkt: (1.0, 2.25).
Euler daje 2.25 wobec prawdziwego wartości 2.71828. Błąd: 0.468. Błąd względny: ~17%.
Otrzymywanie obszaru stabilności Eulera
Dla równania testowego dy/dx = λy (gdzie λ to liczba zespolona), metoda Eulera daje:
> yₙ₊₁ = yₙ + h·λ·yₙ = yₙ·(1 + hλ)
Współczynnik wzrostu na krok: z = 1 + hλ.
Warunek stabilności: obliczona rozwiązanie pozostaje ograniczone tylko wtedy, gdy |z| ≤ 1, to znaczy |1 + hλ| ≤ 1.
To jest warunek geometryczny w zespolonej hλ-planie: punkt hλ musi leżeć wewnątrz koła o promieniu 1, ośrodek w (-1, 0).
Obszar stabilności Eulera: { hλ ∈ ℂ : |1 + hλ| ≤ 1 }
Dla rzeczywistego, ujemnego λ (np. ODE zanikającej jak dy/dx = -2y): hλ musi leżeć w przedziale (-2, 0] na osi rzeczywistej. Z λ = -2 i h = 0,5: hλ = -1. To jest dokładnie na granicy stabilności - metoda jest marginalnie stabilna, co wyjaśnia niepowodzenie jakościowe w wcześniejszym przykładzie.
Z h = 1 i λ = -2: hλ = -2, co umieszcza nas poza obszarem stabilności. Rozwiązanie oscyluje z rosnącą amplitudą.
Znalezienie granicy stabilności
Metoda Runge-Kutta 4 (RK4) ma większy obszar stabilności niż metoda Eulera, co jest jednym z powodów, dla których preferowana jest w większości problemów.
Dla ujemnych rzeczywistych λ, RK4 pozwala na hλ aż do około -2,785 na osi rzeczywistej (w porównaniu z limitem Eulera -2).
Dla szybko zaniku równań, gdy wartości własne λ mają różne wielkości - na przykład λ₁ = -1 i λ₂ = -1000 - stabilność wymaga, aby hλ₂ znajdować się wewnątrz obszaru. Dla RK4 na osi rzeczywistej: h·(-1000) ≥ -2,785, więc h ≤ 0,002785.
Ten mały rozmiar kroku, wymuszony twardą wartością własną λ₂, sprawia, że symulacja jest droga nawet jeśli wolny składnik λ₁ mógłby użyć h = 2.
Fixed Points & Bazy Atrakcji
Metoda Eulera zastosowana do dy/dx = f(y) definiuje dyskretne mapowanie: y₊n₋1 = g(y₊n) = y₊n + h…f(y₊n).
Punkt ustawowy tego mapowania: y tak, że g(y) = y. Dla Eulera na dy/dx = f(y), punkty ustawowe spełniają warunek f(y) = 0 - równowagi równania różniczkowego.
Stabilność punktu ustawowego: jeśli |g'(y)| < 1, bliskie iteracje konwergują do y. Jeśli |g'(y*)| > 1, odchodzą.
g'(y) = 1 + h…f'(y). W punkcie ustawowym y: |1 + h…f'(y)| < 1 dla stabilności.
To jest dokładnie warunek stabilności metody Eulera z λ = f'(y*) - liniowaizacja równania różniczkowego w punkcie równowagi.
Baza atrakcji: zestaw warunków początkowych, które konwergują do y* pod dyskretnym mapowaniu Eulera. Dla nieliniowych systemów, brzeg bazy definiuje, gdzie symulacja będzie wiarygodnie śledzić równowagę równania różniczkowego vs odchodzić do innego atraktora.
Pętla symulacji tworzy dyskretny system dynamiczny. Jego zachowanie jakościowe - konwergencja, oscylacja, odchylenie - zależy od wielkości kroków h w stosunku do geometrii pola kierunkowego równania różniczkowego.
Połączenie geometrii z projektem symulacji
Geometria symulacji numerycznej sprowadza się do trzech pytań:
1. Gdzie znajduje się region stabilności? Dla Eulera: dysk |1 + hλ| ≤ 1. Większy dla RK4. Nieograniczony (cała lewa połowa płaszczyzny) dla metod implisity.
2. Gdzie znajdują się własne wartości ODE? Własne wartości λ macierzy Jacobiego funkcji f w każdym punkcie określają, który region stabilności musi zawierać hλ.
3. Jakie h utrzymuje wszystkie hλ wewnątrz regionu? Maksymalne dopuszczalne h = (promień regionu stabilności) / max|λ|.
Dla układów sztywnych: max|λ| jest olbrzymie, co zmusza do bardzo małego h dla metod eksperymentalnych. Metody implisity są drogie w kroku, ale pozwalają na duże h — wymieniają koszt kroku na stabilność.
Włamanie Hamminga przekłada się: wybór metody numerycznej koduje zakład o geometrii spektrum własnych wartości ODE. Ujawnienie tego zakładu to pierwsze decyzja projektowa w każdej symulacji.