Устойчивость по Ляпунову

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

Определение

Формально, пусть задана система обыкновенных дифференциальных уравнений:

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

где \( f(x, t) \) удовлетворяет условиям существования и единственности решения. Пусть \( x_0(t) \) — некоторое частное решение (например, положение равновесия, когда \( f(x_0, t) \equiv 0 \)). Решение \( x_0(t) \) называется устойчивым по Ляпунову, если для любого \( \varepsilon > 0 \) существует такое \( \delta > 0 \) (зависящее от \( \varepsilon \) и \( t_0 \)), что для любого решения \( x(t) \) системы, начальное условие которого удовлетворяет неравенству \( \| x(t_0) - x_0(t_0) \| < \delta \), выполняется \( \| x(t) - x_0(t) \| < \varepsilon \) для всех \( t \ge t_0 \).

Если, помимо этого, существует такое \( \delta_0 > 0 \), что при \( \| x(t_0) - x_0(t_0) \| < \delta_0 \) выполняется \( \lim_{t \to \infty} \| x(t) - x_0(t) \| = 0 \), то решение называется асимптотически устойчивым.

Если же для некоторого \( \varepsilon > 0 \) не существует такого \( \delta \), то решение называется неустойчивым.

История

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

Классификация типов устойчивости

Устойчивость по Ляпунову

Определение, данное выше, является основным. Оно описывает, что траектории, начинающиеся вблизи исследуемого решения, остаются в его окрестности.

Асимптотическая устойчивость

Более сильное свойство: траектории не только остаются в окрестности, но и стремятся к исследуемому решению с течением времени. Это наиболее распространённый тип устойчивости в прикладных задачах (например, в теории автоматического регулирования).

Экспоненциальная устойчивость

Частный случай асимптотической устойчивости, при котором скорость сходимости к решению оценивается экспоненциальной функцией: \( \| x(t) - x_0(t) \| \le C e^{-\alpha t} \| x(t_0) - x_0(t_0) \| \) для некоторых констант \( C > 0, \alpha > 0 \).

Устойчивость по Ляпунову для неавтономных систем

Для систем, явно зависящих от времени, определение усложняется: требуется, чтобы \( \delta \) не зависело от \( t_0 \) (равномерная устойчивость) или чтобы выполнялись дополнительные условия (например, равномерная асимптотическая устойчивость).

Методы исследования устойчивости

Первый метод Ляпунова (метод линеаризации)

Основан на анализе поведения системы в малой окрестности положения равновесия. Система линеаризуется (разлагается в ряд Тейлора, отбрасываются нелинейные члены), и устойчивость исходной нелинейной системы определяется по собственным значениям матрицы Якоби (матрицы линеаризации), вычисленной в точке равновесия.

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

Второй метод Ляпунова (прямой метод)

Не требует решения системы дифференциальных уравнений. Основан на построении специальной скалярной функции — функции Ляпунова \( V(x, t) \), которая является положительно определённой (в некоторой окрестности) и её производная по времени в силу системы (полная производная) является отрицательно определённой (или знакоотрицательной). Если такую функцию удаётся найти, то положение равновесия устойчиво (или асимптотически устойчиво).

Этот метод является универсальным, но для многих систем построение функции Ляпунова — нетривиальная задача.

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

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

• Критерий Гурвица (Рауса — Гурвица): проверка знаков определителей, составленных из коэффициентов характеристического полинома.
• Критерий Михайлова: геометрический критерий, основанный на поведении годографа характеристического полинома на комплексной плоскости.

Применение

Теория автоматического управления

Устойчивость по Ляпунову является ключевым требованием к любой системе управления. Система должна возвращаться к заданному состоянию (или траектории) после возмущений. Методы Ляпунова используются для анализа устойчивости нелинейных и адаптивных систем.

Механика и физика

Исследование устойчивости положений равновесия механических систем (например, маятника, спутника на орбите), устойчивости течений жидкости (гидродинамическая устойчивость), устойчивости плазмы.

Биология и экология

Модели «хищник — жертва», эпидемиологические модели (SIR) — анализ устойчивости стационарных состояний (например, эндемического равновесия).

Экономика

Исследование устойчивости рыночного равновесия, моделей экономического роста.

Примеры

Устойчивый фокус

Система \( \dot{x} = -x - y, \dot{y} = x - y \). Собственные значения матрицы линеаризации: \( -1 \pm i \). Вещественные части отрицательны — асимптотическая устойчивость. Траектории закручиваются к началу координат.

Неустойчивый узел

Система \( \dot{x} = x, \dot{y} = 2y \). Собственные значения: 1 и 2 — положительны. Траектории расходятся от начала координат.

Центр

Система \( \dot{x} = -y, \dot{y} = x \). Собственные значения: \( \pm i \). Вещественные части равны нулю. По линеаризации вывод об устойчивости сделать нельзя. В данном случае система является консервативной, и траектории представляют собой замкнутые кривые (окружности) — устойчивость по Ляпунову есть, но асимптотической устойчивости нет.

Связь с другими понятиями

• Орбитальная устойчивость: относится к устойчивости всей траектории (например, периодического решения), а не отдельной точки.
• Структурная устойчивость: свойство системы сохранять качественный характер фазового портрета при малых изменениях правых частей уравнений.
• Устойчивость по Лагранжу: ограниченность всех решений системы (без требования стремления к равновесию).

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

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

Источники

• Ляпунов А. М. «Общая задача об устойчивости движения» (1892).
• Малкин И. Г. «Теория устойчивости движения» (1952).
• Халил Х. К. «Нелинейные системы» (2002, 3-е издание).
• В. И. Зубов. «Методы А. М. Ляпунова и их применение» (1957).