Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
7
Добавлен:
16.04.2015
Размер:
96.42 Кб
Скачать

Лекция 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; : : : ; A1B) = m; m · n;

то тогда в системе (3) будет (n ¡ m) неуправляемых собственных чисел.

Доказательство. Рассмотрим линейную оболочку, натянутую на линейно независимые вектора матриц B; AB; : : : ; A1B. Обозначим ее

Lin L ¢ ©B; AB; : : : ; A1Bª

=

Покажем, что Lin L содержит в себе вместе с любым своим элементом, элемент домноженный слева на A, ò. å.

l 2 Lin L ) Al 2 Lin L:

Это утверждение очевидно для столбцов матриц (B; AB; : : : ; A2B), однако не оче- видно для столбцов матрицы A1B. По теореме Гамильтона-Кэли, матрица является корнем своего характеристического многочлена

f(A) = An + a1A1 + ¢ ¢ ¢ + an = 0;

или в другой форме

An = ¡a1A1 ¡ ¢ ¢ ¢ ¡ an;

и умножая последнее равенство справа на B, имеем

50

AnB = ¡a1A1B ¡ ¢ ¢ ¢ ¡ anB;

откуда видно, что каждый столбец матрицы AnB является линейной комбинацией столбцов матриц (B; AB; : : : ; A1B).

Как было сформулировано в лемме 1,

 

 

rang(B; AB; : : : ; A1B) = 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; : : : ; A1B) = m; m · n;

домножим слева блочную матрицу на T ¡1. В результате этого умножения получится матрица, ранг которой будет также m , ò. ê. det T ¡1 =6 0

T ¡1(B; AB; : : : ; A1B) =

= (T ¡1B; T ¡1AB; : : : ; T ¡1A1B):

В доказательстве леммы 1 было показано, что блок T ¡1B имеет вид

T ¡1B = Ã B1 !;

0

ãäå áëîê B1 имеет размерность (m £ r). Рассмотрим, какую структуру имеют блоки

(T ¡1B; T ¡1AB; : : : ; T ¡1A1B), имеем

T ¡1AB = T ¡1AT T ¡1B = Ã

A

A2

B1

! = Ã

A B

 

!:

1

 

 

1

1

0

A3

0

0

 

Подобную же структуру можно проиллюстрировать и для последующих блоков, вплоть до T ¡1A1B

T ¡1AiB = Ã

i

!; i = 1; : : : ; (n ¡ 1):

A10B1

Таким образом

 

(T ¡1B; T ¡1AB; : : : ; T ¡1A1B) =

:!

=

ÃÃ 01

10 1

!: : : Ã

1 0

1

!

 

 

B

 

A B

 

 

A1B

 

 

 

где последние (n ¡ m) строк нулевые. Ранг этой матрицы равен m, что влечет за

собой полную управляемость в соответствии с критерием Калмана.

¥

54

Соседние файлы в папке Лекции по ТУ