Линейная алгебра
Матричные уравнения: Сильвестра и Ляпунова
Как проверить устойчивость системы управления самолёта без симуляции на 1000 лет? Уравнение Ляпунова даёт ответ за один шаг алгебры.
- **Авиация:** MATLAB lyap() лежит в основе анализа устойчивости управляющих систем всех коммерческих самолётов Airbus и Boeing
- **Балансировка моделей:** Hankel singular values через Ляпунов-Грамианы - метод Model Order Reduction в автомобильных ЭБУ
- **Оценка состояния (Kalman):** фильтр Калмана решает уравнение Риккати - нелинейное обобщение уравнения Ляпунова
- **Квантовые вычисления:** уравнение Линдблада (открытые квантовые системы) имеет ляпуновскую структуру над комплексными матрицами
Предварительные знания
- Разложение Шура
- Собственные значения и спектр матрицы
- Положительная определённость
Уравнение Сильвестра
Уравнение Ляпунова AX + XA^T = -Q (Q > 0) имеет единственное решение X > 0 тогда и только тогда, когда A гурвицова (все Re(lambda_i) < 0). Это фундаментальный тест устойчивости. MathWorks (MATLAB) продаёт лицензии lyap() за 2000$/год - алгоритм Бартельса-Стюарта лежит в основе анализа устойчивости всей авиационной промышленности.
Уравнение Сильвестра AX + XB = C является фундаментальным линейным матричным уравнением: оно обобщает скалярное уравнение ax + xb = c. Условие существования единственного решения - A и -B не имеют общих собственных значений - прямо следует из векторизации через кронекерово произведение.
Когда уравнение Сильвестра AX + XB = C имеет единственное решение?
Теорема Сильвестра-Розенблюма: AX + XB = C имеет единственное решение X тогда и только тогда, когда A и -B не имеют общих собственных значений. Доказательство через диагонализацию или представление через интеграл exp(At) C exp(Bt).
Уравнение Ляпунова и устойчивость
Уравнение Ляпунова - частный случай Сильвестра при B = A^T. Возникает при анализе устойчивости линейных систем: dx/dt = Ax устойчива тогда и только тогда, когда для любой Q > 0 существует P > 0, удовлетворяющее AᵀP + PA = -Q.
Связь с устойчивостью: система dx/dt = Ax устойчива (Re lambda_i < 0) тогда и только тогда, когда для любой Q > 0 уравнение Ляпунова AP + PA^T + Q = 0 имеет решение P > 0. Функция Ляпунова V(x) = x^T P x строго убывает вдоль траекторий.
Какое условие на матрицу A гарантирует устойчивость системы dx/dt = Ax?
Линейная система dx/dt = Ax устойчива (асимптотически) тогда и только тогда, когда Re(lambda_i) < 0 для всех собственных значений A. Эквивалентно: уравнение Ляпунова AᵀP + PA = -Q имеет положительно определённое решение P для любой Q > 0.
Уравнение Риккати и LQR
Уравнение Риккати AᵀP + PA - PBR^{-1}BᵀP + Q = 0 - нелинейное обобщение уравнения Ляпунова. Решение задаёт оптимальный регулятор LQR (Linear Quadratic Regulator), минимизирующий функционал J = ∫(xᵀQx + uᵀRu)dt. Применяется в управлении самолётами, ракетами, портфельной оптимизации.
Критерий Сильвестра имеет геометрический смысл: уравнение AX + XB = C задаёт линейный оператор на пространстве матриц. Невырожденность этого оператора эквивалентна тому, что спектры A и -B не пересекаются.
Чем уравнение Риккати отличается от уравнения Ляпунова?
Уравнение Ляпунова AᵀP + PA = -Q линейно по P. Уравнение Риккати AᵀP + PA - PBR^{-1}BᵀP + Q = 0 содержит квадратичный по P член и нелинейно. Решается итерационными методами: Schur decomposition, Newton-Kleinman iteration.
Связь с теорией управления и числовыми методами
Матричные уравнения Ляпунова и Риккати - краеугольный камень современной теории управления.
- Фильтр Калмана — Связанная тема
- Оптимальное управление LQR — Связанная тема
- Model Order Reduction — Связанная тема
- Спектральная теория операторов — Связанная тема
Итоги
- Уравнение Сильвестра AX+XB=C имеет единственное решение iff спектры A и -B не пересекаются
- Алгоритм Бартельса-Стюарта: разложение Шура + обратная подстановка за O(n^3)
- Уравнение Ляпунова AP+PA^T+Q=0 (Q>0) имеет P>0 iff A гурвицова (все Re(lambda)<0)
- Функция Ляпунова V(x)=x^TPx убывает вдоль траекторий dx/dt=Ax iff AP+PA^T<0
- Грамианы управляемости и наблюдаемости - решения уравнений Ляпунова, определяют балансировку системы
- Уравнение Риккати - нелинейное обобщение Ляпунова, основа фильтрации Калмана и LQR-управления
Каков физический смысл матрицы P > 0, решающей уравнение Ляпунова AP + PA^T + Q = 0?
dV/dt = d/dt(x^T P x) = x^T(AP + PA^T)x = -x^T Q x < 0. Это прямое доказательство устойчивости через функцию Ляпунова. P > 0 гарантирует, что V(x) > 0 при x != 0, а уравнение AP+PA^T = -Q < 0 гарантирует убывание V.