Oscylator Van der Pol: Rytm Regulujący Swoją Amplitudę

Oscylator, który sam się tłumia i odśmient

Standardowy oscylator z tłumieniem straci energię podczas każdej cykle i jego amplituda maleje do zera. Oscylator Van der Pol, wprowadzony przez holenderskiego inżyniera Balthasara van der Pola w latach 20., podczas badań nad obwodami przenośnicowymi, robi coś niezwykłego: współczynnik tłumienia zależy od amplitudy. Małe oscylacje otrzymują dodatkową energię; duże oscylacje są tłumione. Wynik to system, który osiąga i utrzymuje konkretną amplitudę, niezależnie od tego, gdzie zaczął.

x'' - mu * (1 - x^2) * x' + x = 0 x zmienna stanu (np. napięcie w obwodzie) mu parametr tłumienia/nieliniowości (mu >= 0) x' prędkość, x'' przyspieszenie Spójrz na współczynnik tłumienia, -mu*(1 - x^2). Gdy |x| < 1, wyrażenie (1 - x^2) jest dodatnie, ale ogólny znak jest ujemny-nieujemny – energia jest wprowadzana i małe oscylacje rosną. Gdy |x| > 1, wyrażenie (1 - x^2) odwraca się na ujemne, tłumienie staje się dodatnie i duże oscylacje zanikają. Przejście przy |x| = 1 to regulacja amplitudy wbudowana bezpośrednio w równanie.

x'' - mu * (1 - x^2) * x' + x = 0

x     the state variable (e.g. circuit voltage)
mu    damping/nonlinearity parameter (mu >= 0)
x'    velocity, x''  acceleration

Cykl graniczny

Dokładnie taka sama zasada odnosi się do cyklu granicznego: niezależnie od warunku początkowego -- czy to mała nawiązkę lub duży przesunięcie -- tor w przestrzeni fazowej (x rysowany przeciwko x') spiraliuje na tę samą zamkniętą ścieżkę. Ta zamknięta ścieżka jest cyklem granicznym: to izolowana okresowa ścieżka, do której są atrakcyjne (lub odpychające) torów wokół niej, a nie jedno z nieskończonego rodziny okresowych ścieżek jak u bezwadnego wahadła. System Van der Pol ma dokładnie jeden stabilny cykl graniczny, a twierdzenie Poincare-Bendixson gwarantuje, że dowolne ograniczone przepływy 2-D bez punktu stałego w swojej obszarze zatrapiającego muszą zbliżać się do okresowej ścieżki -- co dokładnie dzieje się tutaj, ponieważ jedyny punkt stały, źródło, jest nieprzewodniczący i odpycha torów na zewnątrz wokół cyklu.

Z sinusoidy do oscylacji relaksacyjnej

Parametr mu kontroluje nie tylko rozmiar, ale także kształt tego cyklu granicznego. Dla mu = 0 równanie redukuje się do prostego oscylatora harmonicznego x'' + x = 0, a technicznie „cykl graniczny” degeneruje się w jedno z zwykłych okrągłych torów oscylatora harmonicznego -- gładkiej sinusoide. Gdy mu rośnie ponad około 1, tłumienie staje się silnie nieliniowe i kształt sygnału zmienia się na wolne, łagodne wzrosty przerywane almost instantaneous, sharp discharge -- oscylacją relaksacyjną podobną do ząbka noża. Fizycznie to jest wolny proces napięcia, gdzie siła powrotowa dominuje, a następnie szybkie zwolnienie po przejściu toru przez próg, w którym tłumienie zmienia znak. Relaksacyjne oscylatory lampy neonowej i własna samonapędzająca się rytmika serca są klasycznymi przykładami tego zachowania na dwóch skalach czasowych.

mu near 0     -> near-sinusoidal limit cycle, one timescale
mu around 1-3 -> visibly distorted, asymmetric oscillation
mu >> 1       -> relaxation oscillation: slow build, fast snap, two separated timescales

Dlaczego to ważne poza jednym przewodnikiem do atmosfery

Van der Pol stworzył równanie, aby wyjaśnić stabilne, samonadzwaniające się oscylacje triodowych schematów używanych w wcześniejszych radiach -- schematów, które, zostawione samodzielnie, osiągają jedno głośne, czyste amplitudę niezależnie od tego, jak były włączone. Ale ta sama struktura -- wprowadzanie energii przy małej amplitudzie i rozpraszać przy dużych amplitudach -- pojawia się wszędzie tam, gdzie system potrzebuje utrzymywać rytm samodzielnie: model FitzHugh-Nagumo neuronu odbijającego się jest bliskim krewnym, komórki pacemakerowe są modelowane za pomocą równań typu Van der Pol, a to pozostaje standardowym przykładem w podręcznikach z dynamiki nieliniowej dokładnie dlatego, że to najprostszy system, który wykazuje cykl graniczny, którego można obliczyć ręcznie.

Numeryczne całkowanie

Zapisane jako układ pierwszego rzędu z y = x', równanie przybiera postać x' = y, y' = mu * (1 - x^2) * y - x, gotowe do użycia w RK4 lub dowolnym standardowym integratorze równań różniczkowych. Jedyną troską jest sztywność: dla dużych wartości mu szybkie zmiany fazowe "snap" zdarzają się na znacznie krótszym czasie niż wolne budowanie, więc stały integrator jawny musi mieć wystarczajną małą wielkość kroku do rozstrzygnięcia snap bez wybuchu, lub adaptacyjny/intymny integrator, jeśli mu jest bardzo wysokie. Demonstracja na tej stronie używa stałego małego kroku RK4, który jest dokładny i stabilny dla średnio sztywnej zakresu wartości mu, w którym przejście od sinusalnego do relaksacyjnego zachowania jest najbardziej widocznie interesujące.

Często zadawane pytania

Wypróbuj na żywo

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