ТУ - лекции Шмырова / Лекции по ТУ / 11
.pdfЛекция 11. Декомпозиция линейной управляемой системы.
Рассмотрим линейную стационарную систему
x = Ax + Bu; |
(1) |
где матрица A имеет размерность (n £ n), матрица B размерности (n £ r), вектор x
имеет размерность n, вектор u размерности r. Допустимое управление вида линейной обратной связи имеет вид
u = Cx; |
(2) |
где матрица C размерности (r £n). При подстановке (2) в (1) получаем "замкнутую" систему
x = (A + BC)x: |
(3) |
В данной лекции будет рассмотрен первый этап исследования управляемой линейной системы вида (1) на полную управляемость. Этап будет заключаться в декомпозиции или разбиении системы (1) на управляемую и неуправляемую подсистемы с помощью специальной замены переменных. Также будет показано, что размерность неуправляемой подсистемы зависит от числа неуправляемых собственных чисел.
Лемма 1. О количестве неуправляемых собственных чисел. Åñëè
rang(B; AB; : : : ; An¡1B) = m; m · n;
то тогда в системе (3) будет (n ¡ m) неуправляемых собственных чисел.
Доказательство. Рассмотрим линейную оболочку, натянутую на линейно независимые вектора матриц B; AB; : : : ; An¡1B. Обозначим ее
Lin L ¢ ©B; AB; : : : ; An¡1Bª
=
Покажем, что Lin L содержит в себе вместе с любым своим элементом, элемент домноженный слева на A, ò. å.
l 2 Lin L ) Al 2 Lin L:
Это утверждение очевидно для столбцов матриц (B; AB; : : : ; An¡2B), однако не оче- видно для столбцов матрицы An¡1B. По теореме Гамильтона-Кэли, матрица является корнем своего характеристического многочлена
f(A) = An + a1An¡1 + ¢ ¢ ¢ + an = 0;
или в другой форме
An = ¡a1An¡1 ¡ ¢ ¢ ¢ ¡ an;
и умножая последнее равенство справа на B, имеем
50
AnB = ¡a1An¡1B ¡ ¢ ¢ ¢ ¡ anB;
откуда видно, что каждый столбец матрицы AnB является линейной комбинацией столбцов матриц (B; AB; : : : ; An¡1B).
Как было сформулировано в лемме 1,
|
|
rang(B; AB; : : : ; An¡1B) = m; |
m · n; |
|
||
следовательно, базисных векторов в Lin L ровно m. Обозначим их |
|
|||||
|
|
l1; l2; : : : ; lm: |
|
|
|
|
Построим матрицу |
|
|
|
|
||
|
|
~ |
~ |
~ |
|
|
|
|
T = (l1; l2; : : : ; lm; lm+1 |
; lm+2; : : : ; ln); |
|
||
ãäå ~ |
~ |
~ |
|
|
n. |
|
lm+1; lm+2; : : : ; ln - дополнение базиса Lin L до базиса R |
|
|
||||
В системе (1) сделаем замену переменных с помощью матрицы T |
||||||
|
|
x = T y: |
|
|
|
(4) |
В новых переменных система (1) будет иметь вид |
|
|
|
|||
|
|
y = T ¡1AT y + T ¡1Bu |
|
|
(5) |
|
Построим матрицу A = T ¡1AT и рассмотрим ее структуру по столбцам. Из матрич- |
||||||
ного уравнения |
e |
|
|
|
|
|
можно записать |
e |
|
|
|
A - a1 и первый |
|
|
|
T A = AT |
|
|
|
|
|
|
уравнение,связывающее первый вектор матрицы |
e e |
|||
вектор матрицы T |
|
|
|
Tea1 = Al1:
Как было отмечено выше, Al1 2 Lin L, следовательно, этот вектор принадлежит m-мерному пространству, что влечет за собой следующую структуру вектора ea
|
0 a11 |
1 |
|
||
|
|
e. |
C |
|
|
|
B a |
m1 |
|
||
|
B |
|
C |
|
|
a1 |
= B |
|
0 |
C |
; |
|
B |
|
C |
|
|
|
B e |
C |
|
||
e |
B |
|
|
C |
|
B |
|
|
C |
|
|
|
B . |
C |
|
||
|
B |
|
|
C |
|
|
B |
|
0 |
C |
|
|
B |
|
C |
|
|
|
@ |
|
|
A |
|
где последние (n ¡ m) компоненты нулевые.
51
Второй столбец матрицы A имеет аналогичную структуру |
|||||
e |
|
|
a12 |
1 |
|
|
|
0 e. |
|
||
|
|
B a |
C |
|
|
|
|
B |
m2 |
C |
|
T a2 |
= Al1; a2 |
= B |
0 |
C |
: |
|
|
B |
C |
|
|
|
|
B e |
C |
|
|
e |
e |
B |
|
C |
|
B |
|
C |
|
||
|
|
B . |
C |
|
|
|
|
B |
|
C |
|
|
|
B |
0 |
C |
|
|
|
B |
C |
|
|
|
|
@ |
|
A |
|
Можно показать что все последующие столбцы матрицы A äî m-ãî имеют такую же |
|||||||||||||||
структуру, что и первые два. Начиная с m + 1 столбца |
|
e |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
ситуация меняется и нулевых |
|||||
компонент в этих столбцах в общем случае нет. Действительно, |
|
|
|||||||||||||
|
|
|
|
T am+1 = Alm+1; am+1 = B e |
|
C: |
|
|
|||||||
|
|
|
|
|
|
|
|
|
a1 m+1 |
C |
|
|
|||
|
|
|
|
|
|
|
|
|
B e . |
|
|
|
|||
|
|
|
|
|
|
|
|
|
0 a2 m+1 |
1 |
|
|
|
|
|
|
|
|
|
e |
|
e |
e |
|
B |
|
C |
|
|
||
|
|
|
|
|
|
B an m+1 |
C |
|
|
||||||
|
|
|
|
|
|
|
|
|
B |
|
C |
|
|
||
|
|
|
|
|
|
|
|
|
B |
|
C |
|
|
||
|
|
|
|
|
|
|
|
|
@ |
|
A |
|
|
||
ò. ê. |
lm+1 принадлежит Rn-мерному пространству. e |
|
|
|
|
|
A имеет |
||||||||
e |
|
|
|
|
|
|
|
Таким образом матрица |
e |
||||||
|
|
|
|
|
|
|
|
|
|||||||
следующую структуру |
|
|
|
|
|
|
1 |
|
|
|
|
|
|||
|
|
0 a11 |
¢ ¢ ¢ |
a1m |
|
a1m+1 |
¢ ¢ ¢ |
a1n |
|
|
|
|
|
||
|
|
B e. |
|
e. |
|
e . |
|
e. |
C |
|
|
|
|
|
|
|
A = B |
|
¢ ¢ ¢ |
|
|
|
¢ ¢ ¢ |
|
C |
= |
|
|
: |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
B |
0 |
|
0 |
am+1m+1 |
|
am+1n |
C |
à |
0 |
A3 |
! |
|
|
|
|
B |
am1 |
|
amm |
|
amm+1 |
|
amn |
C |
|
A1 |
A2 |
|
|
|
|
B |
¢ ¢ ¢ |
|
¢ ¢ ¢ |
C |
|
|
|
||||||
|
e |
B |
e |
e |
|
e |
e |
C |
|
|
|
|
|
||
|
B |
|
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
B |
. |
|
. |
e . |
|
e . |
C |
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
B |
|
0 |
|
anm+1 |
|
ann |
C |
|
|
|
|
|
|
|
|
B 0 |
¢ ¢ ¢ |
|
¢ ¢ ¢ |
C |
|
|
|
|
|
||||
|
|
B |
|
|
|
|
C |
|
|
|
|
|
|||
|
|
@ |
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
e |
|
¡ |
|
e |
|
Рассмотрим теперь матрицу B = T |
|
1B. Перепишем уравнение для B â âèäå |
||||
|
e |
T B = B; |
e |
|||
соответственно, уравнение связывающиеe |
первые столбцы матриц B è B имеет вид |
|||||
соответственно, первый столбец |
e |
|
|
|
B имеет вид |
e |
|
T b1 |
= b1 2 Lin L; |
матрицы e
52
Для остальных того, что
|
|
|
|
0 b11 |
1 |
|
|
|
|
|
|
B e. |
C |
|
|
|
|
b |
|
B bm1 |
C |
; |
|
|
|
1 |
= B |
|
C |
||
|
|
|
B |
0 |
C |
|
|
|
|
|
|
B |
C |
|
|
|
|
e |
B e |
C |
|
||
|
|
B |
|
C |
|
||
|
|
|
|
B |
|
C |
|
|
|
|
|
B . |
C |
|
|
|
|
|
|
B |
0 |
C |
|
|
|
|
|
B |
C |
|
|
|
|
|
|
B |
|
C |
|
|
|
|
|
@ |
|
A |
|
r |
¡1 |
столбцов матрицы B будет такой же вид, как и для b в следствие |
|||||
|
|
|
e |
|
|
e |
B |
e |
|
T bi = bi 2 Lin L; i = 1; : : : ; r: |
Таким образом матрица |
e |
будет иметь вид |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
1 |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
0 |
|
|
¢ ¢ ¢ |
e. |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
B e. |
|
|
C |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
b11 |
|
|
b1r |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
B = B |
¢ ¢ ¢ |
|
|
|
C |
= |
|
|
|
|
|
: |
|
|
||||||
|
|
|
|
B |
|
|
|
|
|
|
|
C |
|
à |
0 |
|
! |
|
|
|||
|
|
|
|
|
|
|
|
0 |
|
|
|
|||||||||||
|
|
|
|
B 0 |
¢ ¢ ¢ |
|
C |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
B e |
e |
C |
|
|
|
B1 |
|
|
|
|||||||||
|
|
e |
|
B bm1 |
|
|
bmr |
C |
|
|
|
|
|
|
||||||||
|
|
|
B |
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
B |
|
|
|
. |
C |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
B . |
|
|
|
C |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
B |
|
|
|
0 |
C |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
B 0 |
|
|
|
C |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
B |
|
|
|
|
|
C |
|
|
|
|
|
|
|
|
|
|||
|
|
|
@ |
|
|
¢ ¢ ¢ |
|
|
|
A |
|
|
|
|
|
|
|
|
|
|||
В новых переменных y система (1) имеет вид |
! + |
à |
|
!u: |
|
|||||||||||||||||
y = |
à y2 |
! = Ã |
0 |
|
A3 |
|
!Ã y2 |
0 |
(6) |
|||||||||||||
|
|
y1 |
|
|
|
|
A1 |
|
A2 |
|
|
|
|
y1 |
|
|
|
|
|
B1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В результате замены переменных x ! y была произведена декомпозиция системы на
управляемую и неуправляемую подсистемы. Для того, чтобы задача стабилизации имела решение, требуется, чтобы собственные числа в неуправляемой подсистеме имели отрицательные вещественные части. Допустимое управление будет иметь вид
u = ¡ C1 0 ¢;
где матрица C1 имеет размерность (r £ m). Замкнутая система будет иметь вид
y = |
à y2 |
! = Ã |
|
1 0 1 |
|
1 |
A3 |
!Ã y2 |
!; |
(7) |
||||
|
|
y1 |
|
|
A |
+ B |
C |
|
A2 |
|
y1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
откуда видно, что у системы (7) набор собственных чисел (в соответствии со свойствами определителя) совпадает с набором собственных чисел диагональных блоков A1 + B1C1 è A3. Если собственные числа матрицы A3 имеют неотрицательнын соб-
ственные числа, то тогда задача стабилизации не имеет решения. В противном случае
53
ищется стабилизирующее управление для управляемой подсистемы, обеспечивающее отрицательность вещественных частей собственных чисел матрицы (A1 + B1C1). Ïî
отношению к системе (1), стабилизирующее управление будет иметь вид
u = ¡ C1 0 ¢T ¡1x:
Лемма 2. Следствие леммы 1. Управляемая подсистема
y1 = A1y1 + B1u
является полностью управляемой.
Доказательство. В формулировке леммы 1 дано, что
rang(B; AB; : : : ; An¡1B) = m; m · n;
домножим слева блочную матрицу на T ¡1. В результате этого умножения получится матрица, ранг которой будет также m , ò. ê. det T ¡1 =6 0
T ¡1(B; AB; : : : ; An¡1B) =
= (T ¡1B; T ¡1AB; : : : ; T ¡1An¡1B):
В доказательстве леммы 1 было показано, что блок T ¡1B имеет вид
T ¡1B = Ã B1 !;
0
ãäå áëîê B1 имеет размерность (m £ r). Рассмотрим, какую структуру имеют блоки
(T ¡1B; T ¡1AB; : : : ; T ¡1An¡1B), имеем
T ¡1AB = T ¡1AT T ¡1B = Ã |
A |
A2 |
!Ã |
B1 |
! = Ã |
A B |
|
!: |
1 |
|
|
1 |
1 |
||||
0 |
A3 |
0 |
0 |
|
Подобную же структуру можно проиллюстрировать и для последующих блоков, вплоть до T ¡1An¡1B
T ¡1AiB = Ã |
i |
!; i = 1; : : : ; (n ¡ 1): |
A10B1 |
Таким образом
|
(T ¡1B; T ¡1AB; : : : ; T ¡1An¡1B) = |
:! |
||||||||
= |
ÃÃ 01 |
!Ã |
10 1 |
!: : : Ã |
1 0 |
1 |
! |
|||
|
|
B |
|
A B |
|
|
An¡1B |
|
|
|
где последние (n ¡ m) строк нулевые. Ранг этой матрицы равен m, что влечет за
собой полную управляемость в соответствии с критерием Калмана.
¥
54