Алгебраическое уравнение Риккати - Algebraic Riccati equation

An алгебраическое уравнение Риккати это тип нелинейного уравнения, которое возникает в контексте бесконечного горизонта оптимальный контроль проблемы в непрерывное время или же дискретное время.

Типичное алгебраическое уравнение Риккати похоже на одно из следующих:

алгебраическое уравнение Риккати с непрерывным временем (CARE):

${displaystyle A^{T}P+PA-PBR^{-1}B^{T}P+Q=0,}$

или алгебраическое уравнение Риккати с дискретным временем (DARE):

${displaystyle P=A^{T}PA-(A^{T}PB)(R+B^{T}PB)^{-1}(B^{T}PA)+Q.,}$

п неизвестно п к п симметричная матрица и А, B, Q, р известны настоящий матрицы коэффициентов.

Хотя обычно это уравнение может иметь много решений, обычно указывается, что мы хотим получить единственное стабилизирующее решение, если такое решение существует.

Происхождение названия

Эти уравнения получили название Риккати из-за их связи с Дифференциальное уравнение Риккати. Действительно, CARE проверяется инвариантными во времени решениями связанного матричного дифференциального уравнения Риккати. Что касается DARE, он подтверждается инвариантными во времени решениями матричнозначного разностного уравнения Риккати (которое является аналогом дифференциального уравнения Риккати в контексте LQR с дискретным временем).

Контекст алгебраического уравнения Риккати с дискретным временем

В бесконечном горизонте оптимальный контроль В задачах, кто-то заботится о значении некоторой интересующей переменной произвольно в далеком будущем, и нужно оптимально выбрать значение контролируемой переменной прямо сейчас, зная, что он также будет вести себя оптимально в любое время в будущем. Оптимальные текущие значения управляющих переменных задачи в любое время могут быть найдены с использованием решения уравнения Риккати и текущих наблюдений за изменяющимися переменными состояния. С несколькими переменными состояния и несколькими контрольными переменными уравнение Риккати будет матрица уравнение.

Алгебраическое уравнение Риккати определяет решение постоянной времени с бесконечным горизонтом Задача линейно-квадратичного регулятора (LQR), а также бесконечного горизонта, инвариантного во времени Линейно-квадратично-гауссовская задача управления (LQG). Это две из самых фундаментальных проблем в теория управления.

Типичная спецификация линейно-квадратичной задачи управления с дискретным временем состоит в том, чтобы минимизировать

${displaystyle sum _{t=1}^{T}(y_{t}^{T}Qy_{t}+u_{t}^{T}Ru_{t})}$

с учетом уравнения состояния

${displaystyle y_{t}=Ay_{t-1}+Bu_{t},}$

куда у является п × 1 вектор переменных состояния, ты это k × 1 вектор управляющих переменных, А это п × п матрица перехода состояний, B это п × k матрица контрольных множителей, Q (п × п) является симметричным положительный полуопределенный государственный Стоимость матрица и р (k × k) - симметричная положительно определенная матрица затрат на управление.

Индукция назад во времени можно использовать для получения оптимального управляющего решения в каждый момент времени,^[1]

${displaystyle u_{t}^{*}=-(B^{T}P_{t}B+R)^{-1}(B^{T}P_{t}A)y_{t-1},}$

с симметричной положительно определенной матрицей текущих затрат п эволюционирует назад во времени из ${displaystyle P_{T}=Q}$ в соответствии с

${displaystyle P_{t-1}=Q+A^{T}P_{t}A-A^{T}P_{t}B(B^{T}P_{t}B+R)^{-1}B^{T}P_{t}A,,}$

которое известно как динамическое уравнение Риккати с дискретным временем для этой задачи. Стационарная характеристика п, актуальная для задачи с бесконечным горизонтом, в которой Т стремится к бесконечности, можно найти, повторяя динамическое уравнение несколько раз, пока оно не сходится; тогда п характеризуется удалением временных индексов из динамического уравнения.

Решение

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

Для ЗАБОТЫ управление

${displaystyle K=R^{-1}B^{T}P}$

а матрица передачи состояния замкнутого контура

${displaystyle A-BK=A-BR^{-1}B^{T}P}$

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

Для DARE управление

${displaystyle K=(R+B^{T}PB)^{-1}B^{T}PA}$

а матрица передачи состояния замкнутого контура

${displaystyle A-BK=A-B(R+B^{T}PB)^{-1}B^{T}PA}$

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

Решение алгебраического уравнения Риккати может быть получено путем матричной факторизации или повторением уравнения Риккати. Один тип итерации может быть получен в случае дискретного времени с помощью динамичный Уравнение Риккати, которое возникает в задаче с конечным горизонтом: в задачах последнего типа каждая итерация значения матрицы актуальна для оптимального выбора в каждый период, который является конечным расстоянием во времени от конечного периода времени, и если это повторяется бесконечно далеко назад во времени, он сходится к конкретной матрице, которая имеет отношение к оптимальному выбору за бесконечный промежуток времени до последнего периода - то есть, когда существует бесконечный горизонт.

Также возможно найти решение, найдя собственное разложение более крупной системы. Для CARE мы определяем Матрица гамильтониана

${displaystyle Z={egin{pmatrix}A&-BR^{-1}B^{T}-Q&-A^{T}end{pmatrix}}}$

С ${displaystyle scriptstyle Z}$ является гамильтоновым, если он не имеет собственных значений на мнимой оси, то ровно половина его собственных значений имеет отрицательную действительную часть. Если обозначить ${displaystyle scriptstyle 2n imes n}$ матрица, столбцы которой образуют базис соответствующего подпространства, в блочно-матричной записи, как

${displaystyle {egin{pmatrix}U_{1}U_{2}end{pmatrix}}}$

тогда

${displaystyle P=U_{2}U_{1}^{-1}}$

является решением уравнения Риккати; кроме того, собственные значения ${displaystyle scriptstyle A-BR^{-1}B^{T}P}$ являются собственными значениями ${displaystyle scriptstyle Z}$ с отрицательной действительной частью.

Для СМЕЛЫ, когда ${displaystyle A}$ обратима, определим симплектическая матрица

${displaystyle Z={egin{pmatrix}A+BR^{-1}B^{T}(A^{-1})^{T}Q&-BR^{-1}B^{T}(A^{-1})^{T}-(A^{-1})^{T}Q&(A^{-1})^{T}end{pmatrix}}}$

С ${displaystyle scriptstyle Z}$ является симплектическим, если у него нет собственных значений на единичной окружности, то ровно половина его собственных значений находится внутри единичной окружности. Если обозначить ${displaystyle scriptstyle 2n imes n}$ матрица, столбцы которой образуют основу соответствующего подпространства, в блочно-матричной записи, как

${displaystyle {egin{pmatrix}U_{1}U_{2}end{pmatrix}}}$

тогда

${displaystyle P=U_{2}U_{1}^{-1}}$

является решением уравнения Риккати; кроме того, собственные значения ${displaystyle scriptstyle A-B(R+B^{T}PB)^{-1}B^{T}PA}$ являются собственными значениями ${displaystyle scriptstyle Z}$ которые находятся внутри единичного круга.

Смотрите также

• Уравнение Ляпунова

Рекомендации

1. Чоу, Грегори (1975). Анализ и управление динамическими экономическими системами. Нью-Йорк: Джон Вили и сыновья. ISBN  0-471-15616-7.

• Питер Ланкастер; Лейба Родман (1995), Алгебраические уравнения Риккати, Oxford University Press, п. 504, г. ISBN  0-19-853795-6
• Алан Дж. Лауб, Метод Шура для решения алгебраических уравнений Риккати

внешняя ссылка

• Помощь решателя CARE пакета инструментов MATLAB Control.
• Помощь решателя DARE пакета инструментов MATLAB Control.
• Онлайн-программа CARE для матриц произвольного размера.
• Решатели Python CARE и DARE.
• Функция Mathematica для решения алгебраического уравнения Риккати с непрерывным временем.
• Функция Mathematica для решения алгебраического уравнения Риккати с дискретным временем.