Критерии устойчивости: от Рауса до Ляпунова

Зачем нужны алгебраические критерии

Собственные значения матрицы A определяют устойчивость. Но для матрицы высокого порядка вычисление корней характеристического многочлена — трудоёмкая задача. В XVIII–XIX веках математики искали алгебраические критерии, позволяющие проверить устойчивость без явного нахождения корней.

Историческая задача: Проектировщику регулятора паровой машины нужно знать, будет ли двигатель устойчив при данных параметрах. Решать полиномиальное уравнение шестого порядка он не может. Критерии Рауса (1877) и Гурвица (1895) дают ответ в виде нескольких неравенств для коэффициентов многочлена.

Критерий Рауса–Гурвица

Для многочлена P(λ) = λⁿ + a_{n-1} λ^{n-1} + ... + a₁ λ + a₀ строим матрицу Гурвица:

H_n = [[a_{n-1}, a_{n-3}, a_{n-5}, ...], [1, a_{n-2}, a_{n-4}, ...], [0, a_{n-1}, a_{n-3}, ...], ...]

Теорема Гурвица: Все корни P(λ) имеют отрицательную вещественную часть (система устойчива) тогда и только тогда, когда все главные миноры Δ₁, Δ₂, ..., Δₙ матрицы H_n положительны.

Для n = 2: P = λ² + a₁λ + a₀. Условия: a₁ > 0 и a₀ > 0.

Для n = 3: P = λ³ + a₂λ² + a₁λ + a₀. Условия: a₂ > 0, a₀ > 0, и a₂a₁ > a₀.

Пример: P(λ) = λ³ + 3λ² + 2λ + 1. a₂ = 3 > 0 ✓, a₀ = 1 > 0 ✓, a₂a₁ = 3·2 = 6 > 1 = a₀ ✓. Все условия выполнены → все корни устойчивы.

Таблица Рауса: алгоритмическая форма

Метод Рауса — более удобный алгоритм проверки (особенно при ручном счёте). Строим таблицу:

Строка 1: λⁿ, a_{n-2}, a_{n-4}, ... Строка 2: a_{n-1}, a_{n-3}, a_{n-5}, ... Строки 3–(n+1): каждый элемент = (b × элемент выше-слева − элемент выше × элемент слева) / b.

Система устойчива тогда и только тогда, когда все элементы первого столбца таблицы имеют одинаковый знак. Число смен знака в первом столбце = числу корней в правой полуплоскости.

Теорема Ляпунова о линеаризации

Прямой метод Ляпунова применяется к нелинейным системам. Но для гиперболических равновесий существует элегантная связь с линеаризацией.

Теорема (Ляпунов, 1892): Рассмотрим x' = Ax + g(x), где |g(x)| = o(|x|) (нелинейные члены исчезают быстрее линейных). Тогда:

• Если все собственные значения A имеют отрицательную вещественную часть → x* = 0 асимптотически устойчиво для нелинейной системы.
• Если хотя бы одно собственное значение имеет положительную вещественную часть → x* = 0 неустойчиво для нелинейной системы.
• Если максимальная вещественная часть собственных значений = 0 → линеаризация ничего не говорит; нужен нелинейный анализ.

Физический смысл: «Маленькая нелинейность» не влияет на тип устойчивости гиперболических равновесий. Именно поэтому линейные модели работают вблизи стабильных рабочих точек.

Критический случай: пограничная устойчивость

Теорема Ляпунова о линеаризации «молчит» при пограничных собственных значениях. Рассмотрим два примера:

Пример 1: ẋ = −x³. Линеаризация: f'(0) = 0. Исходная система: V = x²/2, V̇ = xẋ = −x⁴ < 0. Асимптотически устойчиво — несмотря на нейтральную линеаризацию.

Пример 2: ẋ = x³. Линеаризация: f'(0) = 0. V = −x²/2, V̇ = −x⁴ < 0 для V, значит x нарастает. Неустойчиво — хотя линеаризация та же.

Вывод: при нулевых вещественных частях собственных значений высшие нелинейные члены принципиально важны.

Устойчивость систем с переменными коэффициентами: парадокс Перрона

Для систем с постоянными A достаточно Re λᵢ < 0 для устойчивости. Для систем x' = A(t)x с переменными коэффициентами это не так!

Пример Перрона: A(t) = [[−1 + 1.5 cos²t, 1 − 1.5 sin t cos t], [−1 − 1.5 sin t cos t, −1 + 1.5 sin²t]]. Мгновенные собственные значения: Re λ = −0.25 < 0 всегда. Но решение системы: x(t) = e^{0.5t} (cos t, −sin t)ᵀ (растёт экспоненциально!).

Показатели Ляпунова для систем с переменными коэффициентами: σᵢ = lim_{T→∞} (1/T) ln|xᵢ(T)|. Система устойчива тогда и только тогда, когда все σᵢ < 0. Вычисление показателей Ляпунова для нелинейных систем — важная задача теории хаоса.

Вопрос для размышления: Как проектировщик системы управления может использовать критерий Рауса–Гурвица для выбора параметров ПИД-регулятора? Почему этот критерий предпочтительнее прямого вычисления собственных значений на начальных этапах проектирования?

Критерий Найквиста и запас устойчивости

Критерий Найквиста (частотный) позволяет оценить устойчивость замкнутой системы по годографу разомкнутой — без вычисления корней характеристического уравнения. Ключевое понятие: запас устойчивости по фазе (phase margin) и по усилению (gain margin). Если годограф Найквиста не обходит точку (−1, 0) — система устойчива. Запас по фазе PM ≥ 45° — практический стандарт в проектировании авиационных и робототехнических систем. Связь с частотным анализом: пиковое значение замкнутой передаточной функции |H(jω)|_max связано с запасом устойчивости: чем меньше PM, тем выше пик, тем более «колебательным» является переходный процесс.