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

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

Solver Równań Różniczkowych Drugiego Rzędu obsługuje liniowe równania różniczkowe zwyczajne postaci a·y″ + b·y′ + c·y = g(x) o stałych współczynnikach rzeczywistych. Narzędzie automatycznie wyprowadza równanie charakterystyczne, klasyfikuje reżim tłumienia (silne, krytyczne, słabe, brak tłumienia lub niestabilne) i generuje zarówno symboliczne rozwiązanie w formie zamkniętej, jak i wysoce dokładne rozwiązanie numeryczne. Interaktywny wynik łączy dwukrzywy wykres czasowy y(x) i y′(x) z trajektorią na płaszczyźnie fazowej (y, y′) — widokiem, który pozwala natychmiast rozpoznać reżim: spirala zbieżna dla tłumienia słabego, węzeł dla silnego, zamknięta pętla dla braku tłumienia i spirala rozbieżna dla drgań niestabilnych.

Co to jest liniowe równanie różniczkowe zwyczajne drugiego rzędu o stałych współczynnikach?

Liniowe równanie różniczkowe zwyczajne drugiego rzędu o stałych współczynnikach rzeczywistych to równanie postaci

gdzie a ≠ 0, b, c są stałymi rzeczywistymi, a g(x) jest członem wymuszającym. Dwa warunki początkowe y(x₀) = y₀ oraz y′(x₀) = y′₀ zmieniają je w zagadnienie początkowe (problem Cauchy’ego) o unikalnym rozwiązaniu w otoczeniu x₀ — co wynika z twierdzenia Picarda-Lindelöfa zastosowanego do równoważnego układu pierwszego rzędu.

Jeśli g(x) = 0, równanie jest jednorodne. W przeciwnym razie jest niejednorodne, a jego pełne rozwiązanie składa się z sumy:

gdzie y_h jest rozwiązaniem ogólnym skojarzonego równania jednorodnego (zawierającym dwie dowolne stałe), a y_p jest dowolnym rozwiązaniem szczególnym pełnego równania. Zastosowanie dwóch warunków początkowych pozwala wyznaczyć wartości obu wolnych stałych.

Równanie charakterystyczne

Przyjęcie założenia y = e^(r·x) w równaniu jednorodnym prowadzi do równania charakterystycznego (pomocniczego):

Jest to równanie kwadratowe, którego wyróżnik Δ = b² − 4ac determinuje całe zachowanie jakościowe układu:

Trzy przypadki pierwiastków i reżimy tłumienia

Metoda przewidywań (dla przypadku niejednorodnego)

Gdy g(x) przyjmuje jedną z poniższych prostych form, metoda przewidywań dostarcza rozwiązanie szczególne poprzez założenie postaci rozwiązania identycznej z wymuszeniem (z nieznanymi współczynnikami) i obliczenie ich:

• Stała g(x) = k. Przewidujemy: y_p = K. Jeśli c = 0, mnożymy przez x; jeśli dodatkowo b = 0, ponownie mnożymy przez x.
• Wielomian stopnia n. Przewidujemy: ogólny wielomian stopnia n. Mnożymy przez x lub x², jeśli wyraz stały lub liniowy rezonuje z pierwiastkami.
• Funkcja wykładnicza g(x) = A·e^(k·x). Przewidujemy: y_p = K·e^(k·x). Jeśli k pokrywa się z pierwiastkiem charakterystycznym, mnożymy przez x (pierwiastek pojedynczy) lub x² (pierwiastek podwójny) — jest to zjawisko rezonansu.
• Funkcje sinusoidalne g(x) = A·cos(ω·x) + B·sin(ω·x). Przewidujemy: y_p = K₁·cos(ω·x) + K₂·sin(ω·x). Mnożymy przez x, jeśli iω jest pierwiastkiem (rezonans częstotliwościowy).
• Iloczyny i sumy powyższych form wynikają z liniowości i reguły iloczynu.

Interpretacja płaszczyzny fazowej

Równoważny układ pierwszego rzędu to u = y, v = y′, gdzie u′ = v oraz v′ = (g(x) − b·v − c·u)/a. Wykreślenie v względem u parametrycznie względem x tworzy trajektorię na płaszczyźnie fazowej. Dla układów jednorodnych autonomicznych (brak x w g), orbity są jednoznacznie określone przez punkt startowy (y₀, y′₀) i ujawniają reżim na pierwszy rzut oka:

• Tłumienie słabe: trajektoria spiralnie zbiega do początku układu.
• Tłumienie silne: trajektoria zbliża się do początku układu wzdłuż linii niezmienniczej.
• Tłumienie krytyczne: węzeł zdegenerowany, trajektoria styczna do pojedynczego wektora własnego.
• Brak tłumienia: zamknięta elipsa wokół początku układu — wieczne oscylacje.
• Niestabilne: trajektoria spiralnie lub gwałtownie ucieka do nieskończoności.

Przykład: Tłumiony oscylator harmoniczny z wymuszeniem

Rozważmy równanie y″ + 2·y′ + 5·y = 10 z warunkami y(0) = 0, y′(0) = 0 — układ z wymuszeniem i słabym tłumieniem.

1. Równanie charakterystyczne: r² + 2r + 5 = 0 → Δ = 4 − 20 = −16 → r = −1 ± 2i.
2. Rozwiązanie jednorodne: y_h = e^(−x)·(C₁·cos 2x + C₂·sin 2x).
3. Rozwiązanie szczególne dla stałego wymuszenia g = 10: zakładamy y_p = K, skąd 5K = 10, czyli y_p = 2.
4. Zastosowanie WP: y(0) = 0 → C₁ + 2 = 0 → C₁ = −2. y′(0) = 0 → −C₁ + 2C₂ = 0 → C₂ = −1.
5. Wynik końcowy: y(x) = 2 − e^(−x)·(2·cos 2x + sin 2x) — układ oscyluje z zanikającą amplitudą, dążąc do y → 2.

Jak korzystać z tego kalkulatora

1. Wprowadź współczynniki a, b, c w górnym rzędzie. Współczynnik a musi być niezerowy (inaczej równanie staje się pierwszego rzędu).
2. Wpisz człon wymuszający g(x) lub pozostaw 0 dla problemu jednorodnego. Rozwiązania w formie zamkniętej są generowane dla stałych, wielomianów do stopnia 2 i funkcji wykładniczych A·e^(k·x) (w tym dla rezonansu).
3. Podaj warunki początkowe (x₀, y₀, y′₀). Należy określić zarówno y, jak i y′ w punkcie x₀, ponieważ równanie jest drugiego rzędu.
4. Wybierz zakres x dla wykresów. Solver wykonuje całkowanie od x₀ w obu kierunkach przy użyciu metody RK4.
5. Kliknij Oblicz i wizualizuj. Otrzymasz równanie charakterystyczne z pierwiastkami na płaszczyźnie zespolonej, klasyfikację reżimu, rozwiązania symboliczne oraz interaktywne wykresy i trajektorię fazową.

Typowe zastosowania

• Układy mechaniczne masa-sprężyna-tłumik: m·x″ + c·x′ + k·x = F(t). Tłumienie silne, krytyczne i słabe odpowiadają różnym współczynnikom tłumienia ζ = c/(2·√(m·k)).
• Obwody elektryczne RLC: szeregowe obwody RLC podlegają równaniu L·Q″ + R·Q′ + Q/C = V(t) — identyczna struktura, inne symbole.
• Wahadło (małe kąty): θ″ + (g/L)·θ = 0 daje prosty ruch harmoniczny; dodanie oporu powietrza daje drgania tłumione.
• Reakcja budynków na trzęsienia ziemi: modelowanie struktur o jednym stopniu swobody z przyspieszeniem podłoża jako wymuszeniem.
• Systemy serwo z regulatorem PID: dynamika błędu w pętli zamkniętej często sprowadza się do równania II rzędu, gdzie tłumienie decyduje o przeregulowaniu.
• Modele populacyjne z inercją: wzrost gospodarczy z opóźnieniem akumulacji kapitału lub modele ekologiczne z opóźnioną reakcją.

Metoda numeryczna — klasyczna metoda Rungego-Kutty (RK4)

Narzędzie sprowadza równanie a·y″ + b·y′ + c·y = g(x) do układu pierwszego rzędu:

z warunkami u(x₀) = y₀, v(x₀) = y′₀. Czteroetapowa metoda Rungego-Kutty jest stosowana do wektora stanu (u, v). RK4 posiada lokalny błąd obcięcia O(h⁵) i błąd globalny O(h⁴); domyślne 400 kroków w każdym kierunku zapewnia dokładność rzędu sześciu cyfr dla problemów, które nie są sztywne.

Często zadawane pytania

Co to jest liniowe równanie różniczkowe zwyczajne drugiego rzędu o stałych współczynnikach?

Liniowe równanie różniczkowe zwyczajne drugiego rzędu o stałych współczynnikach ma postać a·y″ + b·y′ + c·y = g(x), gdzie a, b, c są stałymi rzeczywistymi, a g(x) jest członem wymuszającym (niejednorodnym). Przy dwóch warunkach początkowych y(x₀) = y₀ i y′(x₀) = y′₀ rozwiązanie jest unikalne. Przypadek jednorodny g(x) = 0 zawsze pozwala na rozwiązanie w formie zamkniętej poprzez równanie charakterystyczne a·r² + b·r + c = 0; przypadek niejednorodny rozwiązuje się jako y(x) = y_h(x) + y_p(x).

Co to jest równanie charakterystyczne?

Dla równania a·y″ + b·y′ + c·y = 0, podstawienie założenia y = e^(r·x) daje a·r² + b·r + c = 0 — równanie charakterystyczne lub pomocnicze. Jego pierwiastki określają postać rozwiązania jednorodnego: dwa różne pierwiastki rzeczywiste dają y_h = C₁·e^(r₁·x) + C₂·e^(r₂·x); pierwiastek podwójny r daje y_h = (C₁ + C₂·x)·e^(r·x); sprzężone pierwiastki zespolone α ± β·i dają y_h = e^(α·x)·(C₁·cos(β·x) + C₂·sin(β·x)).

Co oznaczają pojęcia: słabe, krytyczne i silne tłumienie?

Terminologia ta pochodzi z modelu masa-sprężyna-tłumik m·x″ + c·x′ + k·x = 0. Tłumienie silne (nadkrytyczne, wyróżnik > 0, dwa pierwiastki rzeczywiste) oznacza, że układ powraca do równowagi powoli bez oscylacji. Tłumienie krytyczne (wyróżnik = 0, podwójny pierwiastek) to najszybszy powrót bez przeregulowania. Tłumienie słabe (wyróżnik < 0, pierwiastki zespolone) daje oscylacje zanikające. Brak tłumienia (b = 0, c/a > 0) daje stałe oscylacje sinusoidalne.

Na czym polega metoda przewidywań?

Dla prostych wymuszeń g(x) — stałych, wielomianów, funkcji wykładniczych, sinusów, cosinusów i ich iloczynów — rozwiązanie szczególne y_p zakłada się w tej samej postaci co g z nieznanymi współczynnikami, które wyznacza się przez podstawienie do równania i porównanie wyrazów. Założenie należy pomnożyć przez x (lub x² dla pierwiastków podwójnych), gdy g(x) rezonuje z pierwiastkiem charakterystycznym.

Co to jest płaszczyzna fazowa?

Dla równania drugiego rzędu sprowadzonego do układu 2D (y, y'), płaszczyzna fazowa przedstawia wykres y′ w funkcji y wraz ze zmianą x. Krzywe na płaszczyźnie fazowej ujawniają reżim na pierwszy rzut oka: spirale zbieżne dla tłumienia słabego, węzły dla silnego, zamknięte elipsy dla ruchu harmonicznego i spirale rozbieżne dla oscylacji niestabilnych. Jest to geometryczny odpowiednik wykresu pierwiastków równania charakterystycznego.

Jakiej metody numerycznej używa to narzędzie?

Klasyczna metoda Rungego-Kutty czwartego rzędu (RK4) jest stosowana do równoważnego układu pierwszego rzędu u = y, v = y′, gdzie u′ = v oraz v′ = (g(x) − b·v − c·u)/a. RK4 ma lokalny błąd obcięcia O(h⁵), a domyślne 400 kroków w każdą stronę zapewnia dokładność rzędu sześciu cyfr dla równań o umiarkowanej sztywności w wybranym oknie.

Dodatkowe materiały

• Równanie różniczkowe liniowe — Wikipedia
• Równanie charakterystyczne — Wikipedia
• Metoda przewidywań — Wikipedia
• Oscylator harmoniczny — Wikipedia
• Płaszczyzna fazowa (Portret fazowy) — Wikipedia
• Metody Rungego-Kutty — Wikipedia