Нелинейные системы и метод линеаризации

От линейного к нелинейному

Мир в основном нелинеен. Маятник с большими колебаниями описывается уравнением θ'' + (g/L) sin θ = 0 — нелинейным. Хищник–жертва следует уравнениям Лотки–Вольтерры — нелинейным. Уравнения Навье–Стокса гидродинамики — нелинейны. Как анализировать нелинейные системы?

Главный инструмент — линеаризация вблизи точки равновесия. Идея: в малой окрестности равновесия нелинейная система «похожа» на линейную. Поведение линейной системы мы уже умеем анализировать через собственные значения.

Точки равновесия и их поиск

Точка равновесия (неподвижная точка) системы x' = f(x) — это точка x*, где f(x*) = 0, то есть производная обращается в нуль. В равновесии система «стоит» и не изменяется во времени.

Пример — маятник: θ'' = −(g/L) sin θ. В матричной форме: x₁ = θ, x₂ = θ'. Система: x₁' = x₂, x₂' = −(g/L) sin x₁. Точки равновесия: x₂ = 0 и sin x₁ = 0, то есть x₁ = nπ. Два типа: θ = 0 (висящий маятник) и θ = π (перевёрнутый маятник).

Матрица Якоби и линеаризация

Разложим f(x) по Тейлору вблизи x*: f(x* + δx) ≈ f(x*) + J(x*) δx = J(x*) δx (поскольку f(x*) = 0).

Матрица Якоби J = ∂fᵢ/∂xⱼ оценённая в точке x* — это «производная векторного поля». Линеаризованная система: δx' = J(x*) δx.

Поведение нелинейной системы вблизи x* определяется собственными значениями J(x*) — при условии, что они «не пограничные» (не на мнимой оси).

Теорема Хартмана–Гробмана

Если все собственные значения J(x*) имеют ненулевую вещественную часть (гиперболическое равновесие), то нелинейная система вблизи x* топологически эквивалентна своей линеаризации: существует непрерывное обратимое преобразование, отображающее фазовые портреты друг в друга.

Практически: при гиперболических равновесиях тип особой точки (узел, фокус, седло) определяется линеаризацией. При наличии собственных значений с нулевой вещественной частью — нет: нужен нелинейный анализ.

Пример: маятник с трением

Уравнение: θ'' + bθ' + (g/L) sin θ = 0. Матрица системы: x' = (x₂, −bx₂ − (g/L) sin x₁).

Матрица Якоби: J = [[0, 1], [−(g/L) cos x₁, −b]].

В точке x = (0, 0)* (нижнее равновесие): J = [[0, 1], [−g/L, −b]].

tr J = −b < 0 (при b > 0), det J = g/L > 0. При малом b: фокус (затухающие колебания вокруг нижнего положения). При критическом b = 2√(g/L): узел. ✓ Физически правильно!

В точке x = (π, 0)* (верхнее равновесие): J = [[0, 1], [g/L, −b]].

det J = −g/L < 0 → седло (независимо от b). Верхнее положение маятника всегда неустойчиво — как и интуиция подсказывает.

Количественный пример: затухающий маятник

Данные: g/L = 4 (рад/с)², b = 1 (с⁻¹). J(0,0) = [[0,1],[−4,−1]].

Собственные значения: λ = (−1 ± √(1−16))/2 = −0.5 ± i√(15)/2 ≈ −0.5 ± 1.94i.

Тип: устойчивый фокус. Маятник с начальными данными θ(0) = 0.5 рад, θ'(0) = 0 совершает затухающие колебания с частотой ≈ 1.94 рад/с и временем затухания 1/0.5 = 2 с.

Бифуркация типа «вилочная»

При изменении параметра точка равновесия может потерять устойчивость и «раздвоиться». Простейший пример: x' = μx − x³.

При μ < 0: единственное равновесие x* = 0, устойчивое. При μ = 0: точка бифуркации. При μ > 0: x* = 0 становится неустойчивым, появляются два новых равновесия x* = ±√μ — устойчивые.

Физический пример: Потеря устойчивости прямой балки под продольной нагрузкой (задача Эйлера о продольном изгибе). При достижении критической нагрузки μ = μcr прямое состояние становится неустойчивым, и балка выгибается в одну из двух сторон.

Вопрос для размышления: Маятник при θ → π становится «неустойчивым перевёрнутым маятником». Как можно стабилизировать его в перевёрнутом положении с помощью вибраций точки подвеса (парадокс Капицы)? Какой тип бифуркации происходит при стабилизации?

Значение линеаризации и её ограничения

Линеаризация — мощный, но ограниченный инструмент. Она полностью описывает поведение системы вблизи гиперболических особых точек (теорема Хартмана–Гробмана). Однако для центров линеаризация ненадёжна: нелинейные члены могут превратить центр в устойчивый или неустойчивый фокус. Для таких пограничных случаев необходим метод функций Ляпунова или вычисление чисел Ляпунова — инвариантов, определяющих «слабую нелинейность».

Теорема Пуанкаре–Бендиксона и предельные циклы

В системах на плоскости (ℝ²) теорема Пуанкаре–Бендиксона утверждает: любая ограниченная траектория при t → ∞ стремится либо к особой точке, либо к предельному циклу, либо к многообразию из особых точек и траекторий. Следствие: хаос невозможен в автономных двумерных системах. Критерий Бендиксона: если div f = ∂f₁/∂x₁ + ∂f₂/∂x₂ не меняет знака в односвязной области D, то в D нет замкнутых траекторий. Это удобный способ исключить предельные циклы. Нейронные осцилляторы: модель Ходжкина–Хаксли для спайков нейрона содержит предельный цикл — устойчивое периодическое решение в 4-мерном фазовом пространстве, ответственное за ритмичные разряды.