Kalkulator Metody Eulera

O Kalkulator Metody Eulera

Kalkulator metody Eulera numerycznie rozwiązuje dowolny problem początkowy pierwszego rzędu postaci \( y' = f(x, y), \; y(x_0) = y_0 \) przy użyciu klasycznej (prostej) metody Eulera. Zwraca pełną tabelę iteracji, wykreśla wielobok Eulera na tle aktywnego pola kierunkowego, porównuje rozwiązanie dla trzech różnych rozmiarów kroków, abyś mógł wizualnie obserwować zbieżność metody, a jeśli podasz dokładne rozwiązanie w postaci analitycznej — generuje analizę błędu dla każdego kroku.

Co to jest metoda Eulera?

Metoda Eulera to najprostszy algorytm aproksymacji rozwiązania problemu początkowego. Zaczynając od znanego punktu \( (x_0, y_0) \) na krzywej rozwiązania, wielokrotnie przesuwa się o mały krok o rozmiarze h wzdłuż lokalnego nachylenia \( f(x, y) \):

Geometrycznie, każdy krok to krótki odcinek prostoliniowy, którego nachylenie jest równe wartości równania różniczkowego w bieżącym punkcie. Wynikowa linia łamana — wielobok Eulera — stanowi przybliżenie prawdziwego (zazwyczaj zakrzywionego) rozwiązania.

Jak dokładna jest ta metoda?

Metoda Eulera jest metodą pierwszego rzędu. Lokalny błąd obcięcia na każdym kroku wynosi \( O(h^2) \), a błąd globalny po całkowaniu w ustalonym przedziale wynosi \( O(h) \). W praktyce:

• Zmniejszenie rozmiaru kroku o połowę w przybliżeniu zmniejsza o połowę błąd globalny.
• Błąd rośnie liniowo wraz z długością przedziału całkowania.
• Błąd jest największy tam, gdzie rozwiązanie ma dużą krzywiznę.

Wbudowane porównanie rozmiarów kroków (h, h/2, h/4) pozwala bezpośrednio zaobserwować tę liniową zbieżność: włącz opcję i sprawdź, czy trzy wartości końcowe zbliżają się do wspólnej granicy, przy czym każda kolejna wartość jest w przybliżeniu o połowę bliżej granicy niż poprzednia.

Jak czytać wykres

Wizualizacja nakłada cztery rodzaje informacji na jednej płaszczyźnie współrzędnych:

• Szare pole kierunkowe — krótkie odcinki, których nachylenie odpowiada wartości \( f(x, y) \) w danym punkcie. Można o tym myśleć jak o „przepływie dyktowanym przez równanie”. Każda krzywa rozwiązania musi być styczna do pola w każdym punkcie.
• Indygo wielobok Eulera — krokowe rozwiązanie numeryczne. Każdy segment zaczyna się w poprzednim punkcie siatki i jest skierowany wzdłuż \( f(x_n, y_n) \) na dystansie h.
• Zielona przerywana krzywa dokładna — widoczna tylko wtedy, gdy podasz rozwiązanie w postaci zamkniętej. Pionowe pomarańczowe przerywane odcinki to błędy lokalne ze znakiem \( y_n - y_{\text{dokładne}}(x_n) \).
• Pomarańczowe i zielone krzywe porównawcze — ten sam problem przeliczony ponownie dla h/2 i h/4, widoczne po włączeniu porównania rozmiarów kroków.

Jak korzystać z tego kalkulatora

1. Wprowadź prawą stronę równania różniczkowego w polu oznaczonym y' =. Użyj x i y jako zmiennych. Obsługiwane operatory to + − × ÷ ^, a obsługiwane funkcje to m.in. sin, cos, tan, asin, acos, atan, exp, ln, log, log10, log2, sqrt, abs.
2. Ustaw warunki początkowe: wartość początkową x₀, początkowe y₀ w tym punkcie, rozmiar kroku h (dodatni, aby całkować w przód, ujemny, aby całkować wstecz) oraz liczbę kroków n.
3. (Opcjonalnie) Podaj dokładne rozwiązanie y(x), jeśli je znasz. Kalkulator obliczy \( |y_n - y(x_n)| \) na każdym kroku i poda błąd maksymalny oraz końcowy.
4. Przełącz opcje wizualizacji: pole kierunkowe jest domyślnie włączone; porównanie rozmiarów kroków nakłada dwie dodatkowe krzywe dla h/2 i h/4.
5. Kliknij Uruchom. Sekcja wyników pokaże podsumowanie statystyk, wykres, panel porównania zbieżności oraz pełną tabelę iteracji. Najechaniu na wiersz towarzyszy podświetlenie odpowiedniego punktu na wykresie (i na odwrót).

Przykład krok po kroku

Rozważmy \( y' = x + y, \; y(0) = 1 \) przy h = 0.1 i 10 krokach. Dokładne rozwiązanie to \( y(x) = -x - 1 + 2e^x \). Zastosowanie wzoru Eulera daje:

Błąd końcowy wynosi około 0.249. Zmniejszenie h o połowę do 0.05 obniża błąd końcowy do około 0.13, a ponowne zmniejszenie do 0.025 obniża go do około 0.067 — czysta zbieżność liniowa, dokładnie tak, jak przewiduje teoria.

Metoda Eulera vs inne metody numeryczne

Kiedy metoda Eulera zawodzi

Prosta metoda Eulera może zachowywać się nieprawidłowo w trzech sytuacjach:

• Zbyt duży rozmiar kroku — wielobok oscyluje lub rozbiega się. Rozwiązaniem jest zmniejszenie h; porównanie h, h/2, h/4 czyni to natychmiast widocznym.
• Sztywne równania różniczkowe (stiff ODEs) — równania z jednocześnie szybko i wolno zanikającymi składowymi wymuszają bardzo małe h dla zachowania stabilności. Należy wtedy przejść na metodę niejawną (wsteczny Euler) lub metody BDF.
• Osobliwości w f(x, y) — dzielenie przez zero, sqrt z liczby ujemnej lub ln z liczby niedodatniej zatrzyma proces całkowania. Kalkulator wyraźnie raportuje krok, w którym wystąpił błąd.

Typowe zastosowania

• Fizyka — druga zasada dynamiki Newtona jako układ pierwszego rzędu, rozpad promieniotwórczy \( \dot{N} = -\lambda N \), prawo stygnięcia Newtona.
• Biologia i epidemiologia — wzrost logistyczny \( \dot{y} = r\,y(1 - y/K) \), modele przedziałowe SIR.
• Ekonomia — ciągła kapitalizacja odsetek, proste modele wzrostu Solowa.
• Chemia — kinetyka reakcji pierwszego rzędu \( \dot{c} = -k c \).
• Dydaktyka — wprowadzanie koncepcji całkowania numerycznego przed przejściem do RK4 lub solverów adaptacyjnych.

Często zadawane pytania

Co to jest metoda Eulera?

Metoda Eulera to najprostsza procedura numeryczna do rozwiązywania problemu początkowego y' = f(x, y), y(x0) = y0. Na każdym kroku wyznacza rozwiązanie jako y_{n+1} = y_n + h · f(x_n, y_n), efektywnie podążając za nachyleniem w bieżącym punkcie na krótkim dystansie h. Jest to metoda pierwszego rzędu, co oznacza, że błąd globalny wynosi O(h).

Jak dokładna jest metoda Eulera?

Metoda Eulera charakteryzuje się lokalnym błędem obcięcia O(h²) i błędem globalnym O(h). Zmniejszenie rozmiaru kroku o połowę z grubsza o połowę zmniejsza błąd globalny. Dlatego porównanie zbieżności dla h, h/2 i h/4 w tym kalkulatorze jest tak pouczające: można zobaczyć, jak błąd maleje w przybliżeniu liniowo wraz z h.

Kiedy metoda Eulera zawodzi?

Metoda Eulera może stać się niestabilna w przypadku problemów sztywnych lub gdy rozmiar kroku jest zbyt duży w stosunku do lokalnej krzywizny rozwiązania. Można wtedy zaobserwować oscylacje rozwiązania numerycznego, ucieczkę do nieskończoności lub wyraźne odejście od prawdziwego rozwiązania. Zmniejszenie h zazwyczaj pomaga; dla równań sztywnych preferowane są metody niejawne, takie jak wsteczna metoda Eulera.

Jak wybrać rozmiar kroku?

Zacznij od h, które daje od 10 do 50 kroków w interesującym Cię przedziale. Jeśli wielobok Eulera wyraźnie odbiega od pola kierunkowego lub od Twojego dokładnego rozwiązania, zmniejsz h o połowę i uruchom ponownie. Użyj wbudowanego porównania h, h/2, h/4, aby sprawdzić, czy trzy krzywe zbiegają się do siebie.

Jaka jest różnica między metodą Eulera a metodą Rungego-Kutty (RK4)?

Metoda Rungego-Kutty czwartego rzędu oblicza nachylenie w czterech punktach na krok i łączy je z wagami (1, 2, 2, 1)/6, dając błąd globalny O(h⁴) — o kilka rzędów wielkości lepszy niż O(h) metody Eulera przy tej samej liczbie kroków. Euler jest nadal cenny do nauki koncepcji całkowania numerycznego oraz w bardzo prostych zastosowaniach o niskiej precyzji.

Czy mogę użyć tego dla układów równań różniczkowych?

Ten kalkulator obsługuje pojedyncze skalarne równanie różniczkowe zwyczajne pierwszego rzędu y' = f(x, y). W przypadku układów lub równań wyższych rzędów można zapisać równanie jako układ pierwszego rzędu i użyć dedykowanego solvera dla układów, lub przekształcić równanie drugiego rzędu na dwa równania pierwszego rzędu i rozwiązywać je składnik po składniku.

Czy mogę całkować wstecz w czasie?

Tak — wprowadź ujemny rozmiar kroku h. Kalkulator będzie poruszał się od x₀ w kierunku ujemnym przez n kroków. Jest to przydatne do rekonstrukcji przeszłości na podstawie znanego stanu obecnego.

Dalsza lektura

• Metoda Eulera — Wikipedia
• Metody Rungego-Kutty — Wikipedia
• Pole kierunkowe — Wikipedia
• Równanie sztywne — Wikipedia