exponenta event banner

lqr

Конструкция линейного квадратичного регулятора (LQR)

Синтаксис

[K,S,e] = lqr(SYS,Q,R,N)
[K,S,e] = LQR(A,B,Q,R,N)

Описание

[K,S,e] = lqr(SYS,Q,R,N) вычисляет матрицу оптимального усиления K.

Для непрерывной временной системы закон обратной связи состояния u = -Kx минимизирует квадратичную функцию затрат

J (u) =∫0∞ (xTQx + uTRu + 2xTNu) dt

в зависимости от динамики системы

x˙=Ax+Bu.

Помимо усиления обратной связи по состоянию K, lqr возвращает решение S связанного уравнения Риккати

ATS + SA (SB + N) R 1 (BTS + NT) + Q = 0

и собственные значения замкнутого цикла e = eig(A-B*K). K получают из S с использованием

K = R 1 (BTS + NT).

Для дискретной модели состояния-пространства u [n] = -Kx [n] минимизирует

J=∑n=0∞{xTQx+uTRu+2xTNu}

в зависимости от x [n + 1] = Ax [n] + Bu [n].

[K,S,e] = LQR(A,B,Q,R,N) - эквивалентный синтаксис для моделей непрерывного времени с x˙=Ax+Bu динамики.

Во всех случаях при пропуске матрицы N, N имеет значение 0.

Ограничения

Данные проблемы должны удовлетворять следующим требованиям:

  • Пара (A, B) является стабилизируемой.

  • R > 0 и Q−NR−1NT≥0.

  • (Q NR 1NT, A BR − 1NT) не имеет ненаблюдаемого режима на воображаемой оси (или единичной окружности за дискретное время).

Совет

lqr поддерживает дескрипторные модели с неингулярными E. Выходные данные S из lqr - решение уравнения Риккати для эквивалентной явной модели «состояние-пространство»:

dxdt = E 1Ax + E − 1Bu

См. также

| | | | |

Представлен до R2006a