Solver Równań Różniczkowych Pierwszego Rzędu

O Solver Równań Różniczkowych Pierwszego Rzędu

Solver równań różniczkowych pierwszego rzędu przyjmuje równanie różniczkowe zwyczajne w postaci dy/dx = f(x, y), automatycznie klasyfikuje jego strukturę (o zmiennych rozdzielonych, liniowe, autonomiczne, dokładne lub ogólne) i generuje zarówno symboliczne rozwiązanie w postaci zamkniętej (tam, gdzie to możliwe), jak i wysokiej dokładności rozwiązanie numeryczne. Wizualizacja pola kierunkowego na żywo z nałożoną krzywą rozwiązania sprawia, że geometryczne znaczenie równania staje się natychmiast oczywiste — rozwiązania są dokładnie krzywymi stycznymi do każdej strzałki.

Co to jest równanie różniczkowe zwyczajne pierwszego rzędu?

Równanie różniczkowe zwyczajne pierwszego rzędu angażuje nieznaną funkcję y(x) i tylko jej pierwszą pochodną y'(x). Standardowa postać jawna to:

W połączeniu z warunkiem początkowym y(x₀) = y₀, definiuje to zagadnienie początkowe (IVP). Twierdzenie Picarda-Lindelöfa gwarantuje unikalne rozwiązanie w pewnym sąsiedztwie x₀, o ile f jest lipschitzowsko ciągła względem y w pobliżu (x₀, y₀). Geometrycznie IVP poszukuje unikalnej krzywej przechodzącej przez (x₀, y₀), której nachylenie w każdym punkcie odpowiada wartości f w tym punkcie — czyli krzywej stycznej do pola kierunkowego.

Sześć klas rozpoznawanych przez solver

Metoda postaci zamkniętej: Liniowe o stałych współczynnikach

Gdy prawa strona upraszcza się do dy/dx = a·y + b ze stałymi a i b, czynnik całkujący μ(x) = e^(-a·x) daje dokładne rozwiązanie. Rozwiązanie ogólne to:

Zastosowanie warunku początkowego y(x₀) = y₀ wyznacza stałą C i daje unikalne rozwiązanie szczególne. Ta jedna klasa obejmuje ogromną liczbę problemów akademickich:

• Wzrost wykładniczy — dy/dx = k·y, rozwiązanie szczególne y(t) = y₀·e^(k·t).
• Zanik wykładniczy — dy/dx = -k·y, okres połowicznego rozpadu ln 2 / k.
• Prawo stygnięcia Newtona — dy/dx = -k·(y - T_otoczenia), temperatura ciała dąży wykładniczo do temperatury otoczenia.
• Ładowanie obwodu RC — dV/dt = (1/RC)·(V_wej - V), napięcie na kondensatorze zbliża się do źródła.
• Klirens leku — farmakokinetyka pierwszego rzędu z szybkością eliminacji k.

Jak czytać pole kierunkowe

W każdym punkcie siatki (x, y) narzędzie rysuje krótki odcinek linii, którego nachylenie jest równe f(x, y). Trzy przydatne obserwacje:

• Punkty równowagi to punkty, w których f(x, y) = 0 — pole kierunkowe jest poziome. Dla równań autonomicznych są to punkty stałe y* spełniające f(y*) = 0; pobliskie trajektorie albo do nich dążą (stabilne), albo się oddalają (niestabilne).
• Izokliny to krzywe, na których f(x, y) jest równe stałej c, więc wszystkie strzałki wzdłuż krzywej mają to samo nachylenie c.
• Krzywe rozwiązań nigdy się nie przecinają (gdy f jest lipschitzowskie) — jest to wizualnie oczywiste, ponieważ dwie przecinające się krzywe wymagałyby różnych nachyleń w punkcie przecięcia.

Metoda numeryczna: Klasyczna metoda Rungego-Kutty (RK4)

Dla danych (x_n, y_n), kolejna wartość jest obliczana poprzez uśrednienie czterech oszacowań nachylenia:

RK4 posiada lokalny błąd obcięcia O(h⁵) i błąd globalny O(h⁴), co daje w przybliżeniu sześciocyfrową dokładność przy domyślnej liczbie kroków dla równań niesztywnych. Solver całkuje na zewnątrz od punktu początkowego w obu kierunkach x i zatrzymuje się czysto, jeśli wartość y przekroczy 10¹⁵ — co jest typowe dla rozwiązań dążących do nieskończoności w skończonym czasie, jak dy/dx = y².

Jak korzystać z tego kalkulatora

1. Wprowadź prawą stronę równania w polu dy/dx = .... Użyj x i y jako zmiennych, * do mnożenia, ^ lub ** dla potęg oraz standardowych funkcji, takich jak sin, cos, exp, log, sqrt. Stałe pi i e są rozpoznawane.
2. Określ warunek początkowy (x₀, y₀) — unikalna krzywa rozwiązania przejdzie przez ten punkt.
3. Wybierz zakres x, w którym ma zostać wykreślone pole kierunkowe i krzywa rozwiązania. Zakres y jest dopasowywany automatycznie na podstawie zintegrowanego rozwiązania.
4. Kliknij Rozwiąż i zwizualizuj. Najpierw uruchamiany jest klasyfikator; jeśli równanie pasuje do wzorca postaci zamkniętej (liniowe ze stałymi współczynnikami), otrzymasz odpowiedź symboliczną. Pole kierunkowe i krzywa rozwiązania są generowane zawsze.
5. Przełącz pole kierunkowe w celu skupienia się na krzywej rozwiązania lub powtórz animację rysowania krzywej, aby zobaczyć, jak postępuje proces całkowania od punktu początkowego.

Przykład: Prawo stygnięcia Newtona

Filiżanka kawy o temperaturze 80 °C stygnie w pomieszczeniu o temperaturze 20 °C. Szybkość wymiany ciepła jest proporcjonalna do różnicy temperatur:

Jest to równanie liniowe o stałych współczynnikach (a = -0.1, b = 2). Postać zamknięta to:

Po 30 minutach: T(30) = 20 + 60·e⁻³ ≈ 22.99 °C. Widok pola kierunkowego sprawia, że zachowanie graniczne jest oczywiste — każda krzywa rozwiązania, niezależnie od temperatury początkowej, dąży asymptotycznie do poziomej linii T = 20.

Typowe zastosowania

• Dynamika populacji — modele wzrostu wykładniczego, logistycznego, efektu Allee.
• Farmakokinetyka — absorpcja i eliminacja leków, obliczenia okresu półtrwania.
• Wymiana ciepła — prawo stygnięcia Newtona, modele pojemności cieplnej.
• Obwody RC i RL — liniowe stany nieustalone pierwszego rzędu.
• Rozpad promieniotwórczy — łańcuchy rozpadu pojedynczych izotopów.
• Zbiorniki mieszające — stężenie substancji rozpuszczonej przy dopływie i odpływie.
• Spadający obiekt z oporem — analiza prędkości granicznej dv/dt = g - kv.

Często zadawane pytania

Co to jest równanie różniczkowe zwyczajne pierwszego rzędu?

Równanie różniczkowe zwyczajne (ODE) pierwszego rzędu to równanie postaci dy/dx = f(x, y), które angażuje nieznaną funkcję y(x) i jej pierwszą pochodną. Rozwiązanie równania oznacza znalezienie funkcji y(x), której pochodna zgadza się z prawą stroną. Z warunkiem początkowym y(x₀) = y₀ rozwiązanie jest unikalne przy łagodnych założeniach regularności (twierdzenie Picarda-Lindelöfa).

Co to jest pole kierunkowe?

Pole kierunkowe (lub pole nachyleń) przedstawia mały odcinek linii w każdym punkcie siatki (x, y), którego nachylenie jest równe f(x, y). Krzywe rozwiązań równania różniczkowego są dokładnie tymi krzywymi, które są styczne do tych odcinków w każdym punkcie. Pole kierunkowe daje natychmiastową wizualną intuicję dotyczącą globalnego zachowania rozwiązań bez symbolicznego rozwiązywania równania.

Które klasy równań różniczkowych pierwszego rzędu rozwiązuje to narzędzie?

Narzędzie automatycznie klasyfikuje równanie do jednej z kategorii: całkowalne (zależne tylko od x, rozwiązywane przez bezpośrednie całkowanie), liniowe o stałych współczynnikach y' = a·y + b (podawana pełna postać zamknięta), autonomiczne (zależne tylko od y), o zmiennych rozdzielonych (rozpada się na g(x)·h(y)), liniowe o zmiennych współczynnikach (P(x)·y + Q(x)) lub ogólne. Dla każdej klasy generowane jest precyzyjne rozwiązanie numeryczne Rungego-Kutty oraz wizualizacja pola kierunkowego.

Jaka metoda numeryczna jest używana?

Stosowana jest klasyczna metoda Rungego-Kutty czwartego rzędu (RK4) z 300 podkrokami w każdym kierunku od punktu początkowego. RK4 posiada lokalny błąd obcięcia O(h⁵) i jest standardowym narzędziem dla niesztywnych równań różniczkowych w tej skali. Solver wykrywa rozbieżność (przepełnienie lub NaN) i czysto przerywa całkowanie, dzięki czemu wykres pozostaje poprawny.

Czym jest metoda czynnika całkującego dla liniowych równań różniczkowych?

Dla liniowego równania różniczkowego pierwszego rzędu y' + P(x)·y = Q(x), obie strony mnoży się przez czynnik całkujący μ(x) = e^∫P(x) dx. Lewa strona staje się pochodną dokładną d/dx[μ·y], więc y(x) = (1/μ(x)) · (∫μ(x)·Q(x) dx + C). Gdy P i Q są stałymi, sprowadza się to do postaci zamkniętej y = -b/a + C·e^(a·x), którą narzędzie zwraca automatycznie.

Czy to narzędzie obsługuje równania sztywne lub układy równań różniczkowych?

Ten solver jest przeznaczony dla niesztywnych skalarnych równań różniczkowych pierwszego rzędu. Problemy bardzo sztywne (gdzie rozwiązanie ma wiele skal czasowych różniących się o wiele rzędów wielkości) mogą wymagać metod niejawnych, takich jak wsteczny Euler lub Rosenbrock; układy sprzężone wymagają solvera wektorowego. W takich przypadkach należy użyć dedykowanego pakietu, takiego jak solve_ivp z SciPy lub specjalistycznych solverów do równań sztywnych.

Dalsza lektura

• Równanie różniczkowe zwyczajne — Wikipedia
• Pole kierunkowe — Wikipedia
• Metody Rungego-Kutty — Wikipedia
• Czynnik całkujący — Wikipedia
• Twierdzenie Picarda-Lindelöfa — Wikipedia