Solver Układów Równań Różniczkowych

O Solver Układów Równań Różniczkowych

Narzędzie Solver Układów Równań Różniczkowych to kompleksowy przybornik do analizy powiązanych układów liniowych i nieliniowych. Wprowadź macierz współczynników 2×2 lub 3×3, a narzędzie przeprowadzi pełną analizę wartości i wektorów własnych, zapisze rozwiązanie ogólne i szczególne w LaTeX, sklasyfikuje punkt równowagi w zerze jako siodło, węzeł, ognisko lub centrum, oraz narysuje interaktywny portret fazowy z animowaną trajektorią. Dla nieliniowych układów płaskich możesz wpisać dowolne prawe strony \(f(x,y)\) i \(g(x,y)\), aby otrzymać portret fazowy wygenerowany metodą RK4.

Co to jest układ równań różniczkowych zwyczajnych?

Układ równań różniczkowych zwyczajnych wiąże kilka nieznanych funkcji jednej zmiennej — zazwyczaj czasu \(t\) — poprzez ich pochodne. W najbardziej zwięzłej formie:

Gdy \(\mathbf{F}(t, \mathbf{x}) = A\mathbf{x}\) dla stałej macierzy \(A\), układ jest liniowy i autonomiczny. To tutaj teoria jest najbardziej elegancka: całe długoterminowe zachowanie jest określone przez wartości własne macierzy \(A\).

Metoda wartości własnych dla układów liniowych

Dla \(\mathbf{x}' = A\mathbf{x}\) standardowa procedura to:

1. Oblicz wielomian charakterystyczny \(\det(\lambda I - A) = 0\).
2. Znajdź wartości własne \(\lambda_1, \lambda_2, \dots\).
3. Dla każdej wartości własnej znajdź wektor własny \(v\), rozwiązując \((A - \lambda I) v = 0\).
4. Złóż rozwiązanie ogólne jako kombinację liniową: \(\mathbf{x}(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} + \cdots\).
5. Wyznacz stałe \(c_i\), podstawiając warunek początkowy \(\mathbf{x}(0)\) do rozwiązania ogólnego.

Trzy przypadki dla układów 2×2

Płaszczyzna ślad-wyznacznik

Dla macierzy 2×2 o śladzie \(T = a_{11} + a_{22}\) i wyznaczniku \(D = a_{11} a_{22} - a_{12} a_{21}\), cała klasyfikacja mieści się w jednym schemacie:

Dlatego panel wyników eksponuje \(T\), \(D\) oraz \(\Delta = T^2 - 4D\) — te trzy liczby wystarczą, aby nazwać punkt równowagi.

Układy nieliniowe i portret fazowy

Większość rzeczywistych równań różniczkowych jest nieliniowa i nie posiada rozwiązania w postaci zamkniętej. Narzędzie obsługuje je, integrując równania numerycznie metodą Rungego–Kutty 4. rzędu (RK4), która posiada lokalny błąd obcięcia \(O(h^5)\) i jest standardem dla gładkich pól wektorowych.

Portret fazowy zawiera:

• Pole wektorowe próbkowane na siatce 13×13, pokazujące kierunek przepływu w każdym punkcie.
• Trajektorię od Twojego warunku początkowego, narysowaną na czerwono z animowanym pomarańczowym punktem pokazującym upływ czasu.
• Kilka linii prądu wychodzących z pierścienia punktów startowych, dających globalny obraz dynamiki.
• Dla liniowych układów 2×2, osie wektorów własnych (przerywany cyjan) — są to kierunki niezmiennicze, wzdłuż których rozwiązania przesuwają się wykładniczo.

Jak korzystać z tego solvera

1. Wybierz tryb — Liniowy 2×2, Liniowy 3×3 lub Nieliniowy 2D — za pomocą zakładek na górze formularza.
2. Wprowadź współczynniki lub równania. Kliknij dowolny Szybki przykład, aby wstępnie wypełnić układ kanoniczny (węzeł stabilny, centrum, siodło, wahadło itp.).
3. Wprowadź warunek początkowy \((x_0, y_0)\) i przedział czasu \(T\). Typowe wartości \(T\) to 6–20 dla oscylatorów i 3–6 dla szybko wygasających układów stabilnych.
4. Kliknij Rozwiąż. Pojawi się pełna strona wyników z klasyfikacją, wartościami i wektorami własnymi, rozwiązaniem w postaci zamkniętej (tryby liniowe), animowanym portretem fazowym i wykresem szeregów czasowych.
5. Powtórz trajektorię używając przycisku pod portretem fazowym, jeśli chcesz ponownie zobaczyć ruch wzdłuż krzywej rozwiązania.

Przykład — tłumiony oscylator harmoniczny

Oscylator tłumiony \(\ddot{x} + 2\zeta \omega \dot{x} + \omega^2 x = 0\) można zapisać jako układ 2D, przyjmując \(y = \dot{x}\):

Dla \(\omega = 1\) i \(\zeta = 0.2\) (słabe tłumienie), macierz to \(A = \begin{pmatrix} 0 & 1 \\ -1 & -0.4 \end{pmatrix}\). Ślad \(T = -0.4\), wyznacznik \(D = 1\), wyróżnik \(\Delta = 0.16 - 4 = -3.84 < 0\), co daje stabilne ognisko z wartościami własnymi \(-0.2 \pm 0.9798\,i\). Trajektoria spiralnie zbiega się do zera, a szeregi czasowe pokazują wykładniczo zanikające sinusoidy.

Zastosowania

• Układy mechaniczne — połączone układy masa-sprężyna, wahadła, żyroskopy.
• Obwody elektryczne — sieci RLC, filtry na wzmacniaczach operacyjnych, sterowanie w przestrzeni stanów.
• Dynamika populacyjna — model drapieżnik-ofiara Lotki–Volterry, gatunki konkurujące, epidemiologia (SIR, SIS).
• Kinetyka chemiczna — sieci reakcji, oscylatory Biełousowa–Żabotyńskiego.
• Neuronauka — model neuronu FitzHugh–Nagumo, redukcje Hodgkina–Huxleya.
• Teoria sterowania — zlinearyzowane modele obiektów, projektowanie obserwatorów, marginesy stabilności.

Wskazówki i uwagi

• Jeśli Twoja trajektoria szybko "ucieka" poza wykres, zmniejsz przedział czasu T — układ niestabilny może przekroczyć zakres widoku w kilka jednostek czasu.
• W przypadku wielokrotnej wartości własnej solver automatycznie znajduje uogólniony wektor własny \(w\), rozwiązując \((A - \lambda I)w = v\), oszczędzając Ci ręcznych obliczeń.
• W układach nieliniowych strzałki pola wektorowego ujawniają również punkty równowagi poza początkiem układu jako cyjanowe kropki — obserwuj regiony o zerowej wielkości pola.
• Dla układów 3×3 nie jest generowany portret fazowy (3D jest trudne do pokazania na płaskiej stronie), ale obliczenia rozwiązania ogólnego i werdykt stabilności nadal obowiązują.
• Warunki początkowe i przedziały czasu są niezależne od klasyfikacji: ich zmiana przesuwa tylko czerwoną trajektorię, nie zmieniając werdyktu opartego na wartościach własnych.

Często zadawane pytania

Co to jest układ równań różniczkowych zwyczajnych?

Układ równań różniczkowych zwyczajnych (ODE) to zestaw powiązanych równań opisujących pochodne kilku nieznanych funkcji jednej zmiennej niezależnej, zazwyczaj czasu. Klasyczna forma to \( \mathbf{x}'(t) = F(t, \mathbf{x}(t)) \), gdzie \( \mathbf{x} \) jest wektorem stanów, a \(F\) polem wektorowym. Układy liniowe można zapisać jako \( \mathbf{x}' = A\mathbf{x} + \mathbf{b} \), a ich dynamika zależy głównie od wartości własnych macierzy \(A\).

Jak wartości własne klasyfikują punkt równowagi układu liniowego 2×2?

Dla układu 2×2 \( \mathbf{x}' = A\mathbf{x} \) punkt zerowy klasyfikuje się na podstawie śladu \(T\) i wyznacznika \(D\) macierzy \(A\): \(D < 0\) to siodło (niestabilne); \(D > 0\) przy \(T^2 > 4D\) to węzeł (stabilny dla \(T < 0\), niestabilny dla \(T > 0\)); \(D > 0\) przy \(T^2 < 4D\) to ognisko (stabilne dla \(T < 0\), niestabilne dla \(T > 0\), centrum dla \(T = 0\)). Granica \(T^2 = 4D\) to węzeł zdegenerowany.

Jak wygląda rozwiązanie w postaci zamkniętej, gdy wartości własne są zespolone?

Jeśli \(A\) ma zespolone sprzężone wartości własne \( \alpha \pm i\beta \) z wektorem własnym \( v = p + iq \), rzeczywiste rozwiązanie ogólne to \( \mathbf{x}(t) = e^{\alpha t} \left[ c_1 (p \cos\beta t - q \sin\beta t) + c_2 (p \sin\beta t + q \cos\beta t) \right] \). Człon \(e^{\alpha t}\) steruje amplitudą, a funkcje trygonometryczne odpowiadają za obrót.

Co się dzieje, gdy macierz ma wielokrotną wartość własną?

Jeśli macierz ma powtarzającą się wartość własną \(\lambda\), ale tylko jeden liniowo niezależny wektor własny \(v\), wyznacza się uogólniony wektor własny \(w\) z równania \( (A - \lambda I) w = v \). Rozwiązanie ogólne ma formę \( \mathbf{x}(t) = c_1 v e^{\lambda t} + c_2 (tv + w) e^{\lambda t} \). Jeśli przestrzeń własna jest 2D, rozwiązanie upraszcza się do \( \mathbf{x}(t) = (c_1 v_1 + c_2 v_2) e^{\lambda t} \).

Czy to narzędzie może rozwiązywać układy nieliniowe symbolicznie?

Tryb nieliniowy rozwiązuje układ numerycznie metodą RK4 i rysuje portret fazowy. Większość takich układów nie ma rozwiązań analitycznych. Możesz zbadać lokalne zachowanie poprzez linearyzację (obliczając macierz Jacobiego w punkcie stałym) i wprowadzając ją jako macierz \(A\) w trybie liniowym.

Co to jest portret fazowy?

Portret fazowy to geometryczny obraz rozwiązań układu 2D na płaszczyźnie \(x\)–\(y\). Każde rozwiązanie tworzy krzywą (trajektorię). Zbiór trajektorii i pole wektorowe pokazują strukturę układu: czy rozwiązania dążą do punktu, oscylują, czy rozbiegają się. Portret fazowy pozwala natychmiast zrozumieć dynamikę układu.

Dalsza lektura

• Układ równań różniczkowych — Wikipedia
• Portret fazowy — Wikipedia (EN)
• Wartości własne i wektory własne — Wikipedia
• Metody Rungego–Kutty — Wikipedia
• Oscylator Van der Pola — Wikipedia

Analiza matematyczna:

• Kalkulator konwolucji
• Kalkulator pochodnych
• Kalkulator pochodnych kierunkowych
• Kalkulator podwójnych całek
• Kalkulator pochodnej niejawnej
• Kalkulator Całek
• Kalkulator odwrotnej transformaty Laplace'a
• Kalkulator transformaty Laplace'a
• Kalkulator Granic
• Kalkulator pochodnych cząstkowych
• Kalkulator pochodnych jednej zmiennej
• Kalkulator szeregu Taylora
• Kalkulator całki potrójnej
• Kalkulator promienia zbieżności
• Kalkulator krzywizny
• Kalkulator wrońskianu
• Kalkulator metody Rungego-Kutty (RK4)
• Kalkulator współczynników szeregu Fouriera
• Kalkulator Objętości Bryły Obrotowej Nowy
• Kalkulator Powierzchni Obrotowej Nowy
• Kalkulator Sumy Riemanna Nowy
• Kalkulator Reguły Trapezów Nowy
• Kalkulator Reguły Simpsona Nowy
• Kalkulator Całki Niewłaściwej Nowy
• Kalkulator Reguły L'Hospitala Nowy
• Kalkulator Szeregu Maclaurina Nowy
• Kalkulator Szeregów Potęgowych Nowy
• Kalkulator Testu Zbieżności Szeregów Nowy
• Kalkulator Sumy Szeregów Nieskończonych Nowy
• Kalkulator Średniego Tempa Zmian Nowy
• Kalkulator Chwilowego Tempa Zmian Nowy
• Kalkulator Pochodnych Powiązanych Nowy
• Kalkulator Optymalizacji (Rachunek Różniczkowy) Nowy
• Kalkulator Gradientu Wielozmiennowy Nowy
• Kalkulator Dywergencji Nowy
• Kalkulator Rotacji Nowy
• Kalkulator Całki Krzywoliniowej Nowy
• Kalkulator Całki Powierzchniowej Nowy
• Kalkulator Metody Newtona Nowy
• Solver Równań Różniczkowych Pierwszego Rzędu Nowy
• Solver Równań Różniczkowych Drugiego Rzędu Nowy
• Kreślarka Pola Kierunków i Nachyleń Nowy
• Kalkulator Metody Eulera Nowy
• Kalkulator Równania Różniczkowego Bernoulliego Nowy
• Solver Układów Równań Różniczkowych Nowy
• Kalkulator Szybkiej Transformaty Fouriera (FFT) Nowy
• Kalkulator Transformaty Z Nowy
• Kalkulator Całkowania Numerycznego Nowy