Функция Ляпунова

Функция Ляпунова — это скалярная функция, определённая в фазовом пространстве динамической системы, которая используется для исследования устойчивости положений равновесия. Метод функций Ляпунова является одним из основных инструментов качественной теории дифференциальных уравнений и теории управления. Он позволяет делать выводы об устойчивости системы без явного решения дифференциальных уравнений, а лишь на основе анализа поведения самой функции и её производной вдоль траекторий системы.

История

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

Определение

Пусть дана автономная система дифференциальных уравнений:

\[ \dot{x} = f(x), \quad x \in \mathbb{R}^n, \quad f(0) = 0, \]

где \(x=0\) — положение равновесия. Функция \(V(x)\) называется функцией Ляпунова для данной системы, если она удовлетворяет следующим условиям в некоторой окрестности \(U\) начала координат:

1. \(V(x)\) непрерывно дифференцируема.
2. \(V(0) = 0\).
3. \(V(x) > 0\) для всех \(x \in U \setminus \{0\}\) (положительная определённость).
4. Производная функции \(V\) вдоль траекторий системы \(\dot{V}(x) = \nabla V(x) \cdot f(x)\) удовлетворяет неравенству \(\dot{V}(x) \le 0\) для всех \(x \in U\) (отрицательная полуопределённость).

Если \(\dot{V}(x) < 0\) для всех \(x \in U \setminus \{0\}\), то функция называется строгой функцией Ляпунова.

Классификация

По типу устойчивости

В зависимости от свойств функции \(V\) и её производной различают несколько типов устойчивости, которые можно установить с помощью функции Ляпунова:

• Устойчивость по Ляпунову: Если существует функция Ляпунова с \(\dot{V}(x) \le 0\), то положение равновесия \(x=0\) является устойчивым.
• Асимптотическая устойчивость: Если существует строгая функция Ляпунова (\(\dot{V}(x) < 0\) для \(x \neq 0\)), то положение равновесия является асимптотически устойчивым.
• Экспоненциальная устойчивость: Если существуют положительные константы \(c_1, c_2, c_3, p\) такие, что \(c_1 \|x\|^p \le V(x) \le c_2 \|x\|^p\) и \(\dot{V}(x) \le -c_3 \|x\|^p\), то положение равновесия является экспоненциально устойчивым.
• Неустойчивость: Теоремы Четаева и Ляпунова о неустойчивости позволяют с помощью функций, не являющихся знакоопределёнными, установить неустойчивость положения равновесия.

По виду функции

Функции Ляпунова могут быть различных видов:

• Квадратичные функции: \(V(x) = x^T P x\), где \(P\) — симметричная положительно определённая матрица. Наиболее распространённый класс для линейных систем.
• Полиномиальные функции: \(V(x)\) — полином степени выше второй.
• Кусочно-линейные функции: Используются для систем с гистерезисом или разрывными правыми частями.
• Логарифмические и экспоненциальные функции: Применяются в нелинейных системах с нелинейностями специального вида.

Применение

В теории устойчивости

Основное применение функции Ляпунова — это доказательство устойчивости или неустойчивости положений равновесия динамических систем. Метод является прямым, то есть не требует решения дифференциальных уравнений. Для линейных систем с постоянными коэффициентами \(\dot{x} = A x\) функция Ляпунова может быть найдена в виде квадратичной формы \(V(x) = x^T P x\), где матрица \(P\) является решением уравнения Ляпунова:

\[ A^T P + P A = -Q, \]

где \(Q\) — произвольная симметричная положительно определённая матрица. Если все собственные значения матрицы \(A\) имеют отрицательные действительные части, то решение \(P\) существует и является положительно определённым.

В теории управления

Функции Ляпунова широко используются при синтезе законов управления:

• Стабилизация: Построение управления \(u = u(x)\) такого, что замкнутая система \(\dot{x} = f(x, u(x))\) является асимптотически устойчивой. Часто используется метод бэкстеппинга (backstepping), где функция Ляпунова строится рекурсивно.
• Адаптивное управление: Функции Ляпунова используются для настройки параметров регулятора в реальном времени.
• Робастное управление: Оценка влияния внешних возмущений на устойчивость системы. Используется понятие функции Ляпунова для систем с возмущениями (ISS-Lyapunov function).
• Оптимальное управление: В задачах аналитического конструирования оптимальных регуляторов (АКОР) функция Ляпунова часто совпадает с функцией Беллмана.

В других областях

• Нейронные сети: Функции Ляпунова используются для анализа устойчивости динамики нейронных сетей (например, сети Хопфилда).
• Экономика: Модели экономической динамики, такие как модели роста или циклов, исследуются на устойчивость с помощью функций Ляпунова.
• Биология: Анализ устойчивости экологических систем (модели «хищник-жертва») и эпидемиологических моделей.

Примеры

Пример 1: Линейный осциллятор с трением

Рассмотрим систему: \[ \ddot{x} + \dot{x} + x = 0. \]

Введём переменные состояния: \(x_1 = x\), \(x_2 = \dot{x}\). Система в форме Коши: \[ \dot{x}_1 = x_2, \quad \dot{x}_2 = -x_1 - x_2. \]

Выберем функцию Ляпунова как полную механическую энергию: \[ V(x_1, x_2) = \frac{1}{2} x_1^2 + \frac{1}{2} x_2^2. \]

Производная вдоль траекторий: \[ \dot{V} = x_1 \dot{x}_1 + x_2 \dot{x}_2 = x_1 x_2 + x_2 (-x_1 - x_2) = -x_2^2 \le 0. \]

Функция \(V\) положительно определена, а \(\dot{V}\) отрицательно полуопределена. Следовательно, положение равновесия \((0,0)\) устойчиво по Ляпунову. Поскольку \(\dot{V}\) обращается в ноль только при \(x_2=0\), а на множестве \(x_2=0\) система не может оставаться (кроме начала координат), из принципа инвариантности ЛаСалля следует асимптотическая устойчивость.

Пример 2: Нелинейная система

Рассмотрим систему: \[ \dot{x}_1 = -x_1 + 2 x_1 x_2^2, \quad \dot{x}_2 = -x_2. \]

Выберем функцию Ляпунова: \[ V(x_1, x_2) = \frac{1}{2} x_1^2 + \frac{1}{2} x_2^2. \]

Производная: \[ \dot{V} = x_1 (-x_1 + 2 x_1 x_2^2) + x_2 (-x_2) = -x_1^2 + 2 x_1^2 x_2^2 - x_2^2 = -x_1^2 (1 - 2 x_2^2) - x_2^2. \]

В области \(|x_2| < 1/\sqrt{2}\) производная \(\dot{V} < 0\) при \((x_1, x_2) \neq (0,0)\). Следовательно, положение равновесия \((0,0)\) асимптотически устойчиво в этой области. Вне этой области устойчивость не гарантируется.

Критика и ограничения

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

Кроме того, функция Ляпунова даёт лишь достаточные условия устойчивости. Если для данной системы не удаётся найти функцию Ляпунова, это не означает, что система неустойчива — возможно, просто не была найдена подходящая функция.

Интересные факты

• В 1992 году, к 100-летию метода, была опубликована монография «The Direct Method of Lyapunov: A Centennial Survey», обобщающая развитие метода за столетие.
• Метод функций Ляпунова является основой для современных методов синтеза управления, таких как управление с помощью скользящих режимов и адаптивное управление с эталонной моделью.
• В 1960-х годах метод был распространён на системы с дискретным временем (разностные уравнения) и системы с запаздыванием.

Источники

1. Ляпунов А. М. Общая задача об устойчивости движения. — М.: Гостехиздат, 1950.
2. Малкин И. Г. Теория устойчивости движения. — М.: Наука, 1966.
3. Халанай А., Векслер Д. Качественная теория импульсных систем. — М.: Мир, 1971.
4. Khalil H. K. Nonlinear Systems. — 3rd ed. — Prentice Hall, 2002.
5. Slotine J.-J. E., Li W. Applied Nonlinear Control. — Prentice Hall, 1991.