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

Лекция 12. Управление спектром линейной системы.

Âлекции 11 была рассмотрена декомпозиция линейной стационарной системы

x = Ax + Bu;

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

u = Cx;

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

Рассмотрим управляемую систему

y = A0y + b0u;

(1)

ãäå y - m-мерный вектор, матрица A0 размерности (m £ m) и вектор b0 размерности m имеют вид

 

 

= B

0

 

0

¢ ¢ ¢ 0

¡am

A

0

 

.

 

¢.¢ ¢

 

¡ ¡

 

B .

 

.. .. . .

 

 

 

0

1

 

0

 

0

am 1

 

 

B

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

B

0

 

0

 

0

 

 

1

 

 

B

 

¢ ¢ ¢

¡

a2

 

 

B

 

 

 

 

 

 

 

 

B

0

 

0

¢ ¢ ¢

1

¡

a

 

 

 

@

 

 

 

 

 

 

1

 

 

 

 

0

1

1

 

C

;

b

0

=

0

:

C

 

 

 

B

.

C

 

C

 

 

 

 

B

C

 

C

 

 

 

 

B

 

C

 

C

 

 

 

 

B

0

C

 

C

 

 

 

 

B

C

 

C

 

 

 

 

@

 

A

 

C

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

Матрица A0 является матрицей Фробениуса и обладает рядом свойств, в частности сразу видны коэффициента характеристического полинома матрицы Фробениуса

det (¸E ¡ A) = ¸m + a1¸1 + ¢ ¢ ¢ + am

т. е. они отображены в последнем столбце. Также, интересным является тот факт, что матрица Калмана для системы (1) имеет вид

S = (b0; A0b0; : : : ; Am0 ¡1b0) = E;

откуда видно, что система (1) является полностью управляемой в соответствии с критерием Калмана. Следовательно, у этой системы нет неуправляемых собственных чисел, и задача стабилизации имеет решение, т. е. допустимое управление вида

55

u = còy;

cò = (c1; c2; : : : ; cm)ò;

(2)

0

0

 

обеспечивает для управляемой системы собственные числа с отрицательными вещественными частями. Рассмотрим лемму, посвященную задаче построения управления для системы (1) с заранее заданными собственными числами.

Лемма 1. Для любого набора комплексных чисел (l1; l2; : : : ; lm) найдется такой

вектор-строка cò

0 , что собственные числа системы

 

 

y = (A0 + b0c0ò)y;

 

 

 

 

 

 

(3)

будут совпадать с (l1; l2; : : : ; lm).

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. Построим замену переменных

 

 

 

 

 

 

 

 

 

y = Kz;

 

 

 

 

 

 

 

 

такую, что получаемая система имеет вид

 

b0u = e0

 

 

+ e0

 

 

z = K

¡1

A0Kz + K

¡1

 

 

 

(4)

 

 

 

 

 

 

A

z

 

b

u;

и обладает тем свойством, что

 

 

 

 

0

0 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A0 = A0

; b0 = B C

:

 

 

 

 

 

 

 

ò

 

 

B

0

C

 

 

 

 

 

 

e

 

e

B

 

C

 

 

 

 

 

 

 

 

 

 

 

B

1

C

 

 

 

 

 

 

 

 

 

 

 

B

C

 

 

 

 

 

Это возможно реализовать положив

 

 

 

@

 

A

 

 

 

 

 

 

am¡3

¢ ¢ ¢

 

a1

 

 

0 1

 

 

 

0 am¡2

 

 

 

 

 

 

K = B

am 1

 

am 2

 

 

a1

 

 

1

C

 

 

¡

 

. ¡ ¢.¢ ¢

 

.

 

 

.

:

 

 

B .

 

.. ..

 

 

 

C

 

 

 

B

 

 

 

 

 

 

 

 

 

 

C

 

 

 

B

a1

 

1

 

 

 

0

 

 

0

C

 

 

 

B

 

 

 

 

 

 

C

 

 

 

B

1

 

0

 

¢ ¢ ¢

 

0

 

 

0

C

 

 

 

B

 

 

 

 

 

 

C

 

 

 

B

 

 

 

 

 

 

 

 

 

 

C

 

 

 

@

 

 

 

 

¢ ¢ ¢

 

 

 

 

 

A

 

 

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

Перепишем последнее равенствоeâ0

âèäå

0

=

0

 

A = K¡1A K

 

Aò:

 

 

KAò = A0K:

 

 

 

 

0

 

 

Компонента j11

матрицы J = KAò получается в результате умножения первой стро-

 

0

 

 

 

 

ки матрицы матрицы K на первый столбец матрицы Aò

0 .

56

 

 

 

 

 

0

0

 

1

 

 

 

 

 

 

B

0

 

C

 

j = am 1 am 2

¢ ¢ ¢

a1 1

 

 

 

:

 

 

¡

 

B

 

 

 

C

 

 

 

 

 

¢B

0

 

C

 

 

 

 

 

 

B

 

C

 

 

 

¡ ¡

 

 

B

 

.

 

C

 

 

11

 

 

B

 

 

C

 

 

 

 

 

 

B

 

a

m

C

 

 

 

 

 

 

B

¡

 

C

 

матрицы A0 e

 

e

 

 

@

 

 

A

 

 

K

 

 

 

 

 

 

 

Компонента j11 матрицы J = A0K получается в результате умножения первой строки на первый столбец матрицы

j11 = ¡

0 0 ¢ ¢ ¢ 0 ¡am

e

 

и оказывается, что

e

j11 = j11:

0 a1

1

B a2

C

B

¢ ¢ ¢

C

¢B

a1

C

B

 

C

B

 

C

B

1

C

B

C

@

 

A

Компонента j12 матрицы J получается в результате умножения первой строки мат-

K. Оказывается, что

 

 

 

 

 

 

 

 

 

 

A0 на второй

e

e

рицы матрицы K на второй столбец матрицы Aò

 

, а компоненту j12

матрицы J

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

получаем путем умножения первой строки матрицы

 

 

 

 

столбец матрицы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j12 = j12:

 

 

 

 

 

 

 

 

 

Такая же ситуация имеет место и для остальных компонент матриц

ò

Je= A0K, которые имеют

 

 

am

 

0

e

 

0

 

 

 

 

0

 

 

J = KA0 è

 

 

симметричный вид

 

 

 

 

 

 

 

 

 

 

 

0 ¡0

am 2 am 3

 

¢ ¢ ¢

1

1

 

 

 

 

B

 

 

¡

 

¡

 

 

¢ ¢ ¢

0

C

 

 

 

J = J = B

0 am¡3

am¡4

 

¢ ¢ ¢

C:

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

B

 

 

.

 

 

.

 

 

..

. .

C

 

 

 

 

e B .

 

 

 

 

 

 

C

 

 

 

 

B

0

 

1

 

 

0

 

 

 

 

0

C

 

 

 

 

B

 

 

 

 

 

 

 

C

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

C

 

 

Вектор b

 

@

 

 

 

 

 

 

 

 

¢ ¢ ¢

 

A

 

 

в системе (4) имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

e0

 

 

 

 

 

 

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

b

0

= K¡1b

0

= B

.

C

:

 

 

 

 

 

 

 

 

 

 

 

B

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

C

 

 

 

 

 

 

 

 

 

e

 

 

 

 

B

0

C

 

 

 

 

 

 

 

 

 

 

 

 

 

B

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

1

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

A

 

 

 

 

 

 

 

57

Рассмотрим допустимое управление вида

u0

= c0 z; c0 = ¡

c1

c2 ¢ ¢ ¢

cm ¢

:

 

e

e

e

e

e

 

 

ò

ò

 

 

 

 

Подставим это управление в (4), имеем

0

 

 

0

 

 

 

0

 

 

1

 

 

¢ ¢ ¢

 

 

 

 

0

 

z = B

 

 

0

 

 

 

1

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

.

 

 

 

 

.

 

 

.

 

 

¢.¢ ¢

 

 

 

 

.

 

B

 

 

 

 

 

 

 

 

 

 

..

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

0

 

 

 

0

 

 

¢ ¢ ¢

 

 

0

 

 

 

 

 

1

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

¡

(a

 

¡

c )

¡

(a

 

¡

c )

¢ ¢ ¢ ¡

(a

 

¡

c )

¡

(a

 

¡

c )

@

 

 

e

 

 

e

 

 

e

 

 

e

B

 

 

 

 

 

 

1

 

 

 

2

 

 

 

1

 

 

 

m 1

 

 

 

2

 

 

 

 

1

 

 

 

m

1

C

C

C

C

CCz: (5)

C

C

A

По виду матрицы этой системы также можно определить коэффициента характеристического многочлена (они находятся в последней строке, со знаком минус)

f(¸) = ¸m + (a1 ¡ ecm)¸1 + ¢ ¢ ¢ + (am ¡ ec1):

Как было заявлено в формулировке леммы, мы хотим показать, что для любого набора комплексных чисел (l1; l2; : : : ; lm) найдется такой вектор-строка cò

0 , ÷òî ñîá-

ственные числа системы (3) будут совпадать с (l1; l2; : : : ; lm). Для этого рассмотрим многочлен

g(¸) = (¸1 ¡ l1)(¸2 ¡ l2) ¢ ¢ ¢ (¸m ¡ lm) = = ¸m + b1¸1 + ¢ ¢ ¢ + bm;

который, в случае, если (l1; l2; : : : ; lm) являются собственными числами системы (3), является характеристическим многочленом, т. е.

g(¸) = f(¸):

Тогда справедлива система

 

a1

cm = b1;

 

>

 

¡ e ¡

 

 

>

 

 

 

 

>

 

 

 

 

>

 

 

 

 

>

 

e

 

 

<

 

= b2

;

8 a2

¡ cm 1

>

> .

>

>

>

: am ¡ ec1 = bm;

откуда

8

> ec = a ¡ b ;

> 1 1 1

>

>

>

< ec2 = a2 ¡ b2;

>

> .

>

>

>

: ecm = am ¡ bm:

58

и стабилизирующее управление будет иметь вид

u = ecòz = ecòK¡1y:

0 0

откуда вектор-строка cò

0 , обеспечивающая для системы (3) любой наперед заданный набор собственных чисел имеет вид

cò = ecòK¡1:

0 0

¥

59

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