Прямой метод Ляпунова: устойчивость без решения уравнений

Александр Михайлович Ляпунов и его вклад

В 1892 году молодой петербургский математик Александр Ляпунов защитил диссертацию «Общая задача об устойчивости движения», перевернувшую теорию дифференциальных уравнений. Ляпунов предложил анализировать устойчивость, не решая уравнения в явном виде — а находя специальную вспомогательную функцию (функцию Ляпунова), поведение которой говорит об устойчивости всей системы.

Аналогия из физики: система устойчива, если её «энергия» убывает вдоль траекторий. Ляпунов обобщил это наблюдение на произвольные динамические системы, заменив настоящую энергию её математическим аналогом — функцией Ляпунова.

Определения устойчивости

Рассматривается система x' = f(x, t), f(0, t) = 0 (точка равновесия в начале координат).

Устойчивость по Ляпунову: Равновесие x* = 0 устойчиво, если для любого ε > 0 существует δ > 0 такое, что |x(t₀)| < δ влечёт |x(t)| < ε для всех t ≥ t₀.

Асимптотическая устойчивость: Дополнительно: |x(t)| → 0 при t → ∞.

Неустойчивость: Существует ε > 0 такое, что для любого δ > 0 найдётся начальное условие с |x(t₀)| < δ и |x(T)| > ε при некотором T.

Интуиция: устойчивое равновесие — «дно ямы»; малые возмущения не приводят к уходу далеко. Асимптотически устойчивое — «дно воронки»: система возвращается к равновесию. Неустойчивое — «вершина холма»: малое возмущение уводит систему прочь.

Функции Ляпунова: определение и смысл

Функция Ляпунова для системы x' = f(x) в окрестности x* = 0 — это функция V(x) такая, что:

1. V(x) > 0 при x ≠ 0, V(0) = 0 (положительно определённая),
2. V̇(x) = ∇V · f(x) = Σᵢ (∂V/∂xᵢ) fᵢ(x) ≤ 0 (производная вдоль траектории неположительна).

Производная V̇ = dV/dt вдоль траектории показывает, как V меняется со временем без явного нахождения x(t). Если V̇ ≤ 0, функция V не возрастает — траектории не могут «вылезти» за уровнь V = const, а значит, система устойчива.

Теорема Ляпунова об устойчивости: Если существует функция Ляпунова V с V̇ ≤ 0, то x* = 0 устойчиво. Если V̇ < 0 (строго), то x* асимптотически устойчиво.

Теорема о неустойчивости: Если V не является положительно определённой, но V̇ > 0 в окрестности нуля на том множестве, где V > 0, то x* неустойчиво.

Примеры функций Ляпунова

Линейная система x' = Ax с Re λᵢ < 0: Ищем V = xᵀ P x (квадратичная форма). Из условия V̇ = xᵀ(PA + AᵀP)x < 0 нужно PA + AᵀP = −Q для некоторой положительно определённой Q. Это уравнение Ляпунова — линейное матричное уравнение на P. При устойчивом A и положительно определённой Q всегда существует единственное положительно определённое P.

Гармонический осциллятор ẍ + ω²x = 0: Используем полную энергию V = ẋ²/2 + ω²x²/2. V̇ = ẋẍ + ω²xẋ = ẋ(ẍ + ω²x) = 0. Производная равна нулю — уровни энергии сохраняются. Это центр: устойчивый, но не асимптотически.

Нелинейный маятник ẍ + sin x = 0: V = ẋ²/2 + (1 − cos x). V̇ = ẋẍ + ẋ sin x = ẋ(ẍ + sin x) = 0. Снова центр — консервативная система.

Маятник с трением ẍ + bẋ + sin x = 0: V = ẋ²/2 + (1 − cos x). V̇ = ẋ(ẍ + sin x) = ẋ(−bẋ) = −bẋ² ≤ 0. При b > 0 V̇ < 0 (кроме x = 0) — асимптотически устойчиво. Трение «диссипирует» энергию, и маятник возвращается в равновесие.

Развёрнутый пример: нелинейный осциллятор

Система: ẋ = y − x³, ẏ = −x − y³.

Ищем функцию Ляпунова: Пробуем V = (x² + y²)/2.

V̇ = x·ẋ + y·ẏ = x(y − x³) + y(−x − y³) = xy − x⁴ − xy − y⁴ = −x⁴ − y⁴ ≤ 0.

При (x, y) ≠ (0, 0): V̇ = −x⁴ − y⁴ < 0. Следовательно, начало координат глобально асимптотически устойчиво — все траектории стремятся к нулю независимо от начальных условий!

Принцип LaSalle

Иногда V̇ ≤ 0, но не строго (V̇ = 0 на каком-то множестве). Принцип Барбалата–LaSalle уточняет вывод.

Формулировка: Пусть V̇ ≤ 0 и множество M = {x : V̇ = 0} не содержит полных траекторий, кроме самого x* = 0. Тогда x* = 0 асимптотически устойчиво.

Применение: Маятник с трением: V̇ = −bẋ² = 0 при ẋ = 0. Это не только начало координат, но и все точки с x = ±π, ẋ = 0. Однако при ẋ = 0, x ≠ 0 имеем ẍ = −sin x ≠ 0 — система немедленно покидает это множество. По принципу LaSalle: нижнее равновесие асимптотически устойчиво. ✓

Вопрос для размышления: Почему для линейной системы x' = Ax с устойчивой матрицей A всегда существует квадратичная функция Ляпунова? Можно ли для нелинейных систем систематически строить функции Ляпунова?