Осцилятор Ван дер Поля: Ритм, який сам визначає свою амплітуду

Дисперсний осцилятор Ван дер Поля

Звичайний дмпнутий осцилятор втрачає енергію з кожної фази та його амплітуда зменшується до нуля. Осцилятор Ван дер Поля, представлений голландським інженером Балтасаром ван дер Полем у 1920-х роках під час вивчення вакуумних ламп, робить щось дивне: його коефіцієнт гасіння залежить від амплітуди. Невеликі коливання отримують енергію; великі коливання пригнічуються. В результаті система стабілізується на певному конкретному рівні та залишається там без змін, незалежно від того, звідки вона почала.

x'' - mu * (1 - x^2) * x' + x = 0 x змінна стану (наприклад, напруга в ланцюгу) mu коефіцієнт гасіння/нелінійність (mu >= 0) x' швидкість, x'' прискорення Розгляньте коефіцієнт гасіння, -mu*(1 - x^2). Коли |x| 1, (1 - x^2) перевертається на негативне значення, пригнічення стає позитивним, і великі коливання зменшуються. Перетин при |x| = 1 є вбудованим саморегулюючим контролем амплітуди, що будується безпосередньо в рівнянні.

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

Цикл обмеження

Який б початковий стан ви не обрали – крихітний імпульс або велике зміщення – траєкторія в фазовому просторі (x, проти x') спірається на один і той самий замкнений цикл. Цей замкнений цикл є циклом обмеження: ізольований періодичний орбіта, до якої притягуються або відштовхуються сусідні траєкторії, а не просто один член безперервної сім’ї орбіт, як у необмеженого набору вкладених періодичних орбіт без тертя має маятник на кожній енергії. Система Ван дер Поля має лише один стабільний цикл обмеження, і теорема фон Берджіса-Бенджіоні гарантує, що будь-який замкнений 2D потік без фіксованої точки в його області утримання повинен наблизитися до періодичної орбіти – що саме відбувається тут, оскільки єдина фіксована точка, початок координат, нестабільна і відштовхує траєкторії зовні циклу.

Від синусоїди до релаксації осциляцій

Параметр му контролює не лише розмір, але й форму цього циклу обмежень. При му = 0 рівняння спрощується до простого гармонічного осцилювача x'' + x = 0, і технічно «цикл обмежень» деградує в один з звичайних колісних орбіт гармонічного осцилювача — гладка синусоїда. Коли му перевищує приблизно 1, гасіння стає сильно нелінійним, і форма хвилі спотворюється у повільний, ніжний підйом, перерваний майже миттєвим, різким розрядом — як пілообразний релаксаційний осцилятор. Фізично це повільна фаза зарядки, де домінуюча сила відновлення, слідує за швидкою фазою вивільнення, коли траєкторія перетинає поріг, де змінюється знак гасіння.

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

Чому це має значення за межами одного вакуумного трубки

Ван дер Поля побудував рівняння, щоб пояснити стабільні, самопідтримувані коливання триводельних схем, які використовувалися в ранніх радіопристроях – схеми, які, залишаючись без уваги, встановлюють один гучний і чистий амплітудний рівень незалежно від того, як їх перемикають. Але та ж сама структура — енергетична ін’єкція на невеликій амплітуді, розсіювання на великій — проявляється будь-яку систему, яка потребує автономного підтримки ритму: модель Фітцгера-Нагумо нейрона є близьким родичем, кардіостимулятори моделюються за допомогою рівнянь типу Ван дер Поля, і це залишається стандартним прикладом нелінійних динамік у підручниках саме тому, що це найпростіша система, яка демонструє цикл обмеження, який можна обчислити вручну.

Інтеграція в чисельний вигляд

Переписавши рівняння як першорядну систему, де y = x', отримане рівняння становить x' = y, y' = mu*(1 - x^2)*y - x, готове для RK4 або будь-якого стандартного чисельного розв’язника РLД. Особливість полягає у відцентності: при великих значеннях mu швидка фаза «розриву» змінюється набагато швидше, ніж повільна збірка, тому для вирішення цієї проблеми потрібно використовувати невеликий крок в явному чисельному розв’язнику або адаптивний/неявний розв’язувач при дуже високих значеннях mu. Демонстрація на цій сторінці використовує фіксований невеликий крок RK4, який є точним і стабільним для помірних значень mu, де відбувається перехід від синусоїдальної до релаксації поведінки.

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Van der Pol Oscillator і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.