Obrazy fazowe: Interpretacja układu dynamicznego bez rozwiązywania równań

Pole wektorowe to cała historia

Weź pod uwagę 2D autonomiczny układ równań różnicowych zwyczajnych: dx/dt = f(x, y) i dy/dt = g(x, y). W każdym punkcie płaszczyzny (f, g) jest wektorem prędkości, który mówi dokładnie, w jaki sposób i z jaką prędkością stan poruszałby się, gdyby tam był. Obraz wykresowy (phase portrait) rysuje siatkę tych wektorów, wraz z rodziną trajektorii, które są wszędzie styczne do nich, zamiast szkicować x(t) i y(t) w zależności od czasu. Tracisz informacje o czasie, ale zyskujesz coś bardziej wartościowego: całe jakościowe zachowanie systemu, widoczne gołym okiem, bez potrzeby znajdowania rozwiązania analitycznego.

dx/dt = a*x + b*y
dy/dt = c*x + d*y        (linear case, matrix form: dX/dt = A X)

Stałe punkty i ich klasyfikacja według wartości własnych

Stały punkt to lokalizacja, w której pole wektorowe jest dokładnie zerowe – zarówno f jak i g zanikają, więc stan umieszczony tam nigdy nie przesuwa się. W pobliżu stałego punktu można zlinearizować system za pomocą jego macierzy Jacobiego, a wartości własne tej macierzy klasyfikują zachowanie trajektorii w pobliży.

Wartości własne macierzy Jakobiego określają charakter ruchu wokół danego punktu. Wartość własna dodatnia oznacza, że trajektoria jest skierowana w kierunku, który prowadzi do większych wartości funkcji, a wartość własna ujemna oznacza odwrotność.

both eigenvalues real, negative        -> stable node (sink)
both eigenvalues real, positive        -> unstable node (source)
real eigenvalues, opposite signs       -> saddle
complex pair, negative real part       -> stable spiral
complex pair, positive real part       -> unstable spiral
purely imaginary pair                  -> center (neutral orbits)

Układy nieliniowe: cykle ograniczone i twierdzenie Poincarégo-Bendixsona

Układy nieliniowe mogą również posiadać izolowany zamknięty orbitę okresową, zwaną cyklem ograniczenia, do której bliskie trajektorie spiralują (lub oddalają się od niej) bez jej dotykania – oscylator Van der Pol jest przykładem z podręcznika. Twierdzenie Poincarégo-Bendixsona mówi, że w ograniczonym obszarze płaszczyzny, każda trajektoria ciągłego 2D systemu autonomicznego, która nie zbliża się do punktu stałego, musi zbliżyć się do okresowej orbity. Trajektorie w płaszczyźnie nie mogą przecinać się samych siebie, więc po prostu nie ma miejsca na nieskończone rozciąganie i składanie, które wymaga chaosu. Jest to powód strukturalny, dla którego fazy 2D zawsze wyglądają uporządkowanie, podczas gdy system 3D, taki jak atractor Lorenzowski, który posiada dodatkową wymiar do zaginania się w nim, może być naprawdę chaotyczny.

Czytanie stabilności i rozgałęzień z portretu fazowego

Nulklinie – krzywa, gdzie dx/dt = 0 oraz krzywa, gdzie dy/dt = 0 – przecinają się dokładnie w punktach staczooonych, a ich szkicowanie ręczne jest klasycznym sposobem na zbudowanie portretu fazowego bez komputera: pole wektorowe wskazuje jedynie pionowo na linii dx/dt = 0 oraz jedynie poziomo na linii dy/dt = 0. Wraz z ciągłą zmianą parametrów układu, punkty staczooone mogą pojawiać się, znikać, zbiegać się lub zmieniać stabilność; rozgałęzienie Hopf to moment, w którym stabilny spiralny obwód staje się niestabilny i odrzuca mały obwód limitowego wokół siebie, mechanizm odpowiedzialny za początek oscylacji w modelach drapieżnik-ofiara, jak również w układach elektronicznych.

Frequently asked questions

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Phase Portrait i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.