Kalkulator Równania Bernoulliego

O Kalkulator Równania Bernoulliego

Kalkulator Równania Bernoulliego rozwiązuje jedno z najsłynniejszych nieliniowych równań różniczkowych pierwszego rzędu — równanie Bernoulliego y' + P(x)y = Q(x)yⁿ — i zamienia klasyczne akademickie wyprowadzenie w interaktywny przewodnik krok po kroku. Narzędzie linearyzuje równanie poprzez podstawienie v = y¹⁻ⁿ, wyznacza czynnik całkujący μ(x) i nakłada wynikową krzywą analityczną na numeryczne rozwiązanie RK4 oraz pole kierunków, dzięki czemu możesz zobaczyć każdy szczegół jednocześnie.

Co to jest równanie różniczkowe Bernoulliego?

Wprowadzone przez Jacoba Bernoulliego w 1695 roku, równanie Bernoulliego jest równaniem różniczkowym zwyczajnym (ODE) pierwszego rzędu postaci

Gdy n = 0, równanie jest liniowe; gdy n = 1, jest ono o zmiennych rozdzielonych. Dla każdej innej wartości rzeczywistej n równanie jest nieliniowe, ale klasyczne podstawienie v = y¹⁻ⁿ przekształca je w liniowe ODE względem v, które można rozwiązać standardową metodą czynnika całkującego.

Metoda Bernoulliego w sześciu krokach

Zaczynając od y' + P(x)y = Q(x)yⁿ:

1. Dzielenie przez yⁿ: \( y^{-n}y' + P(x)\,y^{1-n} = Q(x) \).
2. Podstawienie v = y¹⁻ⁿ: zauważ, że \( v' = (1-n)y^{-n}y' \), więc \( y^{-n}y' = v'/(1-n) \).
3. Linearyzacja: \( v' + (1-n)P(x)\,v = (1-n)Q(x) \) — liniowe ODE pierwszego rzędu względem v.
4. Czynnik całkujący: \( \mu(x) = \exp\!\left(\int (1-n)P(x)\,dx\right) \), zatem \( (\mu v)' = \mu(1-n)Q(x) \).
5. Rozwiązanie dla v(x): \( v(x) = \frac{1}{\mu(x)}\left[\mu(x_0)v_0 + \int_{x_0}^{x}\mu(t)(1-n)Q(t)\,dt\right] \).
6. Podstawienie zwrotne: \( y(x) = v(x)^{1/(1-n)} \).

Gdy zaangażowane całki są elementarne, otrzymujemy czystą postać analityczną; w przeciwnym razie kalkulator oblicza je numerycznie przy użyciu reguły Simpsona, aby wykreślić krzywą rozwiązania.

Przypadki szczególne obsługiwane automatycznie

Przykład — n = 2, typ logistyczny

Rozważmy y' + y/x = x·y² z warunkiem początkowym y(1) = 1. Tutaj P(x) = 1/x, Q(x) = x oraz n = 2, więc 1 − n = −1.

1. Podstawiamy v = y⁻¹ = 1/y. Wtedy v' = −y⁻²y' i równanie przyjmuje postać v' − (1/x)v = −x.
2. Czynnik całkujący: μ(x) = exp(∫−1/x dx) = 1/x.
3. (μ·v)' = μ·(−x) = −1. Całkując: (1/x)·v = −x + C, czyli v = −x² + Cx.
4. Stosujemy WP: dla x = 1, v = 1/1 = 1, więc 1 = −1 + C ⇒ C = 2. Stąd v(x) = −x² + 2x.
5. Podstawienie zwrotne: y(x) = 1/v(x) = 1/(2x − x²) = 1/(x(2 − x)).

Rozwiązanie analityczne y = 1/(x(2−x)) posiada pionowe asymptoty w x = 0 i x = 2 — to dokładnie ten rodzaj zjawiska, który pole kierunków czyni oczywistym na pierwszy rzut oka.

Jak korzystać z tego kalkulatora

1. Wypełnij kreator równania. Wpisz P(x) i Q(x) w niebieskie pola, a wykładnik n w małe pole w indeksie górnym. Układ odzwierciedla postać standardową y' + P(x)y = Q(x)yⁿ.
2. Ustaw warunek początkowy (x₀, y₀) oraz zakres wykresu [x min, x max]. Zakres powinien zawierać x₀.
3. Kliknij Rozwiąż. Kalkulator wykryje, czy mamy do czynienia z przypadkiem szczególnym (n = 0 lub n = 1) i wyświetli odpowiednie wyprowadzenie. W przeciwnym razie przeprowadzi pełne podstawienie Bernoulliego z równaniami w formacie MathJax.
4. Przeanalizuj wykres. Pomarańczowa krzywa to numeryczne rozwiązanie RK4. Niebieska przerywana krzywa to postać analityczna wyznaczona przez czynnik całkujący. Pole strzałek pokazuje nachylenie y' w każdym punkcie.
5. Skopiuj CSV punktów próbnych, jeśli chcesz zaimportować trajektorię do innego programu.

Wskazówki, pułapki i przypadki graniczne

• Niecałkowite n wymaga y > 0. Solver flaguje kombinacje takie jak n = 1/2 przy y₀ ≤ 0, gdzie yⁿ byłoby wartością zespoloną.
• y₀ = 0 jest często punktem osobliwym. Każde równanie Bernoulliego z Q ≠ 0 i n > 0 posiada trywialne rozwiązanie y ≡ 0, które zazwyczaj nie jest gałęzią, której szukasz.
• Unikaj osobliwości P(x) w pobliżu x₀. Wyrażenia takie jak 1/x wymagają x₀ ≠ 0; solver sprawdza to przed obliczeniami.
• Bardzo duże wykładniki (|n| > 20) są odrzucane, aby zapobiec przepełnieniu (overflow). W praktyce równania Bernoulliego z tak dużym n prawie nigdy nie występują w rzeczywistych problemach.
• Asymptoty pionowe. Jeśli RK4 staje się rozbieżne, spróbuj zawęzić zakres x do tej strony x₀, po której rozwiązanie pozostaje skończone.

Gdzie pojawiają się równania Bernoulliego

• Dynamika populacji — równanie logistyczne y' = ry(1 − y/K) to w istocie równanie Bernoulliego (n = 2 po przekształceniu).
• Kinetyka chemiczna — reakcje autokatalityczne często podlegają zależności y' ∝ y − y².
• Obwody elektryczne — niektóre obwody RL z nieliniowymi rezystorami dają postać Bernoulliego.
• Mechanika płynów — równania warstwy przyściennej po redukcji podobieństwa.
• Modele epidemiologiczne — frakcja podatna w modelu SIR może zostać zredukowana do postaci Bernoulliego.
• Wzrost gospodarczy — model Solowa-Swana ze stałą stopą oszczędności jest równaniem Bernoulliego z n = α.

Często zadawane pytania

Co to jest równanie różniczkowe Bernoulliego?

Równanie Bernoulliego to równanie różniczkowe zwyczajne pierwszego rzędu postaci y' + P(x)y = Q(x)yⁿ, gdzie P i Q są funkcjami ciągłymi, a n jest dowolną liczbą rzeczywistą. Jest to klasyczny przykład nieliniowego równania, które można przekształcić w liniowe poprzez podstawienie v = y¹⁻ⁿ.

Jak działa podstawienie v = y¹⁻ⁿ?

Pomnóż oryginalne równanie przez y⁻ⁿ, tak aby każdy składnik z y stał się y¹⁻ⁿ lub y⁻ⁿy'. Przyjęcie v = y¹⁻ⁿ daje v' = (1−n)y⁻ⁿy'. Podstawienie to przekształca równanie Bernoulliego w v' + (1−n)P(x)v = (1−n)Q(x), które jest liniowe względem v i rozwiązalne za pomocą czynnika całkującego.

Co się dzieje, gdy n = 0 lub n = 1?

Gdy n = 0, równanie jest już równaniem liniowym pierwszego rzędu, więc podstawienie nie jest wymagane. Gdy n = 1, wzór Bernoulliego wymagałby dzielenia przez 1 − n = 0, więc przypadek ten obsługujemy oddzielnie: równanie sprowadza się do y' = (Q(x) − P(x))·y, które jest równaniem o zmiennych rozdzielonych z rozwiązaniem analitycznym y = y₀·exp(∫(Q−P) dx).

Czy równania Bernoulliego zawsze można rozwiązać w postaci analitycznej?

W zasadzie tak, ale wynikowe całki zawierające czynnik całkujący mogą nie mieć elementarnych funkcji pierwotnych. W takim przypadku kalkulator oblicza je numerycznie za pomocą reguły Simpsona i rysuje krzywą rozwiązania. Sama metoda zawsze sprowadza równanie Bernoulliego do kwadratur.

Dlaczego ujemne y i niewymierne n powodują problemy?

Jeśli n nie jest liczbą całkowitą, yⁿ jest definiowane jako exp(n·ln y) i jest rzeczywiste tylko dla y > 0. Podanie ujemnego y dałoby liczbę zespoloną. Solver sygnalizuje tę sytuację i prosi o y₀ > 0 lub całkowity wykładnik, aby rozwiązanie pozostało rzeczywiste.

Co pokazuje pole kierunków?

Pole kierunków to siatka małych odcinków stycznych, których nachylenie odpowiada wartości y' w danym punkcie (x, y). Każda krzywa rozwiązania musi podążać za tymi stycznymi, dzięki czemu pole kierunków pozwala zobaczyć jakościowy kształt wszystkich rozwiązań jednocześnie, a warunek początkowy wyodrębnia konkretną krzywą.

Dalsza lektura

• Równanie różniczkowe Bernoulliego — Wikipedia (EN)
• Czynnik całkujący — Wikipedia (EN)
• Funkcja logistyczna — Wikipedia (EN)
• Pole kierunków — Wikipedia (EN)