ТУ - лекции Шмырова / Лекции по ТУ / 12
.pdfЛекция 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¸m¡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 am¡1 |
1 |
|
B am¡2 |
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 |
|
|
|
|
|
|
m¡1 |
|
|
|
2 |
|
|
|
1 |
|
|||||
|
|
m 1 |
|
|
|
2 |
|
|
|
|
m¡1 |
|
|
|
m |
1
C
C
C
C
CCz: (5)
C
C
A
По виду матрицы этой системы также можно определить коэффициента характеристического многочлена (они находятся в последней строке, со знаком минус)
f(¸) = ¸m + (a1 ¡ ecm)¸m¡1 + ¢ ¢ ¢ + (am ¡ ec1):
Как было заявлено в формулировке леммы, мы хотим показать, что для любого набора комплексных чисел (l1; l2; : : : ; lm) найдется такой вектор-строка cò
0 , ÷òî ñîá-
ственные числа системы (3) будут совпадать с (l1; l2; : : : ; lm). Для этого рассмотрим многочлен
g(¸) = (¸1 ¡ l1)(¸2 ¡ l2) ¢ ¢ ¢ (¸m ¡ lm) = = ¸m + b1¸m¡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