Научная статья на тему 'Алгоритм пошаговой стабилизации линейного объекта управления'

Алгоритм пошаговой стабилизации линейного объекта управления Текст научной статьи по специальности «Математика»

CC BY
63
14
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по математике, автор научной работы — Леонов Г. А., Шумафов М. М.

Приведено элементарное доказательство теоремы о стабилизируемости вполне управляемой линейной системы.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Elementary proof of the stabilization theorem of the completely controllable system is given.

Текст научной работы на тему «Алгоритм пошаговой стабилизации линейного объекта управления»

УДК 62-50

алгоритм пошаговой стабилизации линейного объекта управления

© 2005 г Г.А. Леонов, М.М. Шумафов

Elementary proof of the stabilization theorem of the completely controllable system is given.

Введение

Задача стабилизации линейным объектом управления с помощью линейной обратной связи является классической и рассматривалась многими авторами [1-7]. Однако все существующие доказательства теоремы о стабилизации для векторного и вещественного случая (когда число управляющих воздействий больше единицы) достаточно громоздки. В настоящей работе предложено элементарное доказательство теоремы о стабилизации, основанное на процедуре «пошаговой» стабилизации. Эта процедура состоит в последовательном сдвиге всех комплексно-сопряженных собственных значений данной матрицы А в левую полуплоскость комплексной плоскости с последующей «поточечной» заменой вещественной части спектра матрицы А на отрицательные вещественные числа.

Постановка задачи

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

X = Ах + Ьи , X £ К" , и £ ят, (1)

где А и Ь - произвольные вещественные (" х ") - и (" х т) -матрицы соответственно; х = х(/) - вектор состояния; и = и(/) - вектор управления; т > 1.

Напомним определения понятий полной управляемости и стабилизируемости системы (1).

Определение 1. Система (1) называется полностью управляемой (пара (А,Ь) называется полностью

управляемой), если для любых векторов х0 £ К", х1 £ К" существует такое управление и(/), являющееся кусочно-непрерывной функцией, заданной на [/0,tl], что для решения х(/) системы (1) с этим управлением и начальным условием х^0) = х0 выполнено равенство х(^) = х1.

Определение 2. Система (1) называется стабилизируемой (пара (А, Ь) называется стабилизируемой), если существует такая обратная связь

и = s * х, (2)

где 5 - постоянная вещественная (" х т) - матрица,

что замкнутая система (1), (2), т. е. система

*

х = (А + Ь5 )х асимптотически устойчива.

*

Здесь знак означает транспонирование.

Задача стабилизации системы (1) сводится к следующей: даны вещественные (" х ") - и (" х т) -

матрицы А и Ь соответственно; найти вещественную

*

(" х т) - матрицу 5 такую, чтобы А + Ь5 была гурви-

цевой, т.е. чтобы вещественные части ее собственных чисел были отрицательными.

Все рассматриваемые в настоящей статье векторы и матрицы являются вещественными. Поэтому слово «вещественный» будем опускать часто.

Вспомогательные утверждения

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

Лемма 1 [4, 8]. Система (1) полностью управляема тогда и только тогда, когда ранг (" х пт) -

матрицы 11 Ь, АЬ, к, А" 1Ь II равен п.

Лемма 2 [4]. Если пара (А, Ь) полностью управляема, то для любой (п х т) -матрицы 8 пара (А + Ь5*,Ь) также полностью управляема.

Лемма 3. Пусть система (1) полностью управляема и

х = ду, у £ я" - (3)

невырожденное линейное преобразование векторной переменной х (0> - неособая матрица).

Тогда полностью управляема также и система

у = Ау + Ьи, у £ я", и £ ят, (4)

где А = д- Ад, Ь = д~1Ь, получающаяся из (1) заменой переменных (3).

Лемму 3 можно переформулировать так: если пара (А, Ь) полностью управляема, то полностью управляема также пара (А, Ь ) ,где матрицы А и Ь определяются из (4).

В дальнейшем нам понадобятся следующие две простые леммы.

Лемма 4. Пусть полностью управляемая система (1) с помощью преобразования (3) приведена к виду

(у1 =4^1 + V, у1 £ я

[у2 = Су + Ь2и, у2 £ я"-1, и £ ят , где 4 £ К ; С - матрица порядка (" -1) х " ;

Ь1 = (Ь11, к, Ь1т); Ь2 - матрица порядка (" -1) х т ;

* *

и = (иь...,ит ) ; у = (у1,у2) .

Тогда существует такое преобразование вектора управления и = Ту (у£ят), что в преобразованной

матрице-строке Ь{ = Ь1Т первый элемент Ь' ^ 0 (здесь Ь1 = (Ь^,..., Ь1т)).

Лемма 5. Пусть полностью управляемая система (1) с помощью преобразования (3) приведена к виду

у1 =41 у1 + 42 у 2 + ь1и, у1 £ К у2 = 4зу1 +44у2 + Ь2и, у2 £ я (6)

у3 = Бу + Ь3и, у3 £ я"-2, и £ ят,

(5)

где А] е Я (] = 1,4); В - матрица порядка (п-2) х п;

Ъ = (Ьп,...,Ът) (' = 1,2); Ъ3 - матрица порядка

* *

(п -2)хт ; и = (иь...,ит) ; у = (уъу2,Уз) .

Тогда существует такое преобразование вектора управления и = Ту (у е Ят), что в преобразованных

матрицах-строках: ¿1 = Ь1Т, ¿2 = Ъ2Т по крайней мере один из первых элементов ¿1'1, ¿21 отличен от нуля: Ъ^ * 0 или ¿21 * 0 (здесь Ъ' = (Ьд,...,Ъ'т), ' = 1,2).

Доказательства лемм 4 и 5 следуют из того, что в силу полной управляемости системы (1) и, следовательно, в силу леммы 3 систем (5) и (6) должно быть

Ъ1 * 0 (в системе (5)) и Ъ1 * 0 или Ъ2 * 0 (в системе

(6)) (иначе в первом уравнении системы (5), а в системе (6) в первом и втором уравнениях не входило бы управление и и тем самым системы (5) и (6) были бы не полностью управляемыми).

Пусть к - номер первой отличной от нуля координаты вектор-строки Ъ1 (в (5)) и Ъ1 или Ъ2 (в(6)). Тогда, сделав переобозначения переменных и1 = ук, ик = у1, и] = у! (] * 1, к), получим утверждения лемм 4 и 5.

Замечание. В силу лемм 4 и 5 без умаления общности всегда можно считать, что в системе (5) Ъ11 * 0 , а в системе (6) или Ъ11 * 0 или Ъ21 * 0.

Следующие три леммы позволяют по заданной (2 х 2) -матрице с комплексно-сопряженными собственными значениями построить другую матрицу того же порядка, имеющую комплексно-сопряженные собственные значения с отрицательными вещественными частями.

Лемма 6. Пусть л

M =

Р =

( а,ю е R, ю Ф 0) с матрицей подобия

(pi е R , i = 1,4), причем

(7)

N (p) =

1 lb

2 - bL+.

fr2 + b22

2 I b b2 (pb + РэЬ?) + (Р2Ь + P4b>)

2 '

(P12 + P22) •bL- (P32 + P4 ) •^

Г ( P) =

A^ (bi2 + b22)

2 2 ( Plb1 + P3b2) + ( P 2 b1 + P 4 b2)

(8)

a = pip4 - p2p3'

если bi * 0

b2 * 0;

n (p) = К ( P) = N (P) = ■ К (P) = -

pi p3 + p2 p4 2 , 2 p3 + p4

A

2 2 P32 + P4

если b = 0

Pi P3 + P2 P4

2 2 Pi2 + Pi

2 2 Pi2 + P2

если bi * 0

b2 * 0:

b2 = 0.

(9)

(i0)

Далее, существует значение векторного параметра

p = p такое, что

N(p0)<0, ^(p0) >0.

(s s Л

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(11) где числа

Доказательство. Пусть S =

Vs2 s 4 )

s^ е R , i = 1,4, подлежат определению. Ищем матрицу S, удовлетворяющую матричному уравнению: Л + BS = P-1 MP . (12)

По условию леммы 6 по крайней мере один из элементов bi (i = 1,4) матрицы B отличен от нуля. Поэтому или a) rank B = 2, или b) rank B = 1. В случае a) решение уравнения (12) очевидно:

S = B _1(M - Л), P := E =

i 0 0 i

Остается рассмотреть случай b): rank B = 1.

Представив (12) в виде

PBS = MP - РЛ, (13)

далее записав (13) как две системы линейных уравнений относительно si, s2, s3, s4 соответственно и учитывая, что rank PB = rank B = 1, получим, что условие разрешимости уравнения (12) относительно S имеет вид:

а-в], в> ъз

а) {¿2 ¿4/

(2 х 2) -матрицы (а, в, Ъ' е Я , ' = 1,4), причем а > 0 , в> 0, + |Ъ2| + |Ъ3| + Ъ4 * 0 .

Тогда существует матрица 8 порядка 2 х 2 такая, что матрица Л + подобна некоторой матрице вида

К® а,

Рз КР2 Р4 , а = а + вЫ(р), ю = рК(р), где р = (рх, р2, Рз , р4),

Pibi + P3b2 P2bi + P 4b2

Pi (P-a) - P2® - P3ß

Pi® + P2(a-a) - P4ß

Pibi + P3b2 Piß + P3(a-a) - P4® P2bi + P 4b2 P2ß + P3®+ P4(a-a)

bi® (Pi2 + P22)+ b2® (PiP3 + P2P4) = = a[biß + b2 (-a)],

b2® (P32 + P42 )+ bi®(PiP3 + P2P4 ) = = a[ß-bi (-a)],

= 0,

= 0,

или

(i4)

(i5)

(a =

Pi P4 - P2 P3

) -

Могут представиться три случая: 1) * 0, ¿2 * 0, 2)Ъ = 0, ¿2 * 0, 3) Ъ1 * 0, Ъ2 = 0 (в случае Ъ1 = Ъ2 = 0 равенства (14) и (15) очевидны).

Из (14) и (15) находим, что в случае 1) имеют место равенства (7), (8), в случае 2) - равенства (7), (9) и в случае 3) - равенства (7), (10). В равенствах (7) положим: 1) в случае 1) р1 = р4 = 0

и Р2,Рз: |Р2 <|РЗ|, Р2РЗ <(16) если Ъ^2 > 0, и р1 = р4 = 0

и р2,рз : \р2\чрз^ р2рз <0 (17)

A

х

2

i

Р0 =1

если Ътр2 < 0;

2) в случае 2)

Р 4 = 0 и Pi, P2, Р3 : Pi Р3 < 0 Р 2 Р3 < 0; (18)

3) в случае 3):

Р2 = 0 и Рь Р3> Р4 : Р1Р3 < 0 Р1Р 4 < (19)

В результате получим неравенства (11), где

/ 0 0 0 0 1 ,

: p , Р2 , Р3 , Р4 / - любое значение параметра р, удовлетворяющее неравенствам (16) - (19) в зависимости от случаев 1) - 3). Лемма 6 доказана.

(а -р\ Г Ъ1 ^

Лемма 7. Пусть Л = 1 I, B = 1 ,

U а J V Ъ2 J

где а, в, Ъг е Л , / = 1,2, а> 0, в> 0, ^ + |Ъ2| Ф 0.

Тогда существует матрица-строка S (порядка 1 х 2) такая, что матрица Л + BS подобна некоторой матри-

це вида I I, где а < 0 и о > 0.

V® а J

Доказательство проводится аналогично доказательству леммы 6. Действительно, по условию леммы 7, В Ф 0, поэтому rank PB = rank B = i (здесь P -

та же самая неособая матрица-параметр, что и в лемме 6). Пусть S = (, s2).

Тогда условие разрешимости уравнения (17) относительно S принимает вид:

Р1Ъ1 + Р3Ъ2 = Р1(а-а)-Р 2о- Р3в = Р2h + Р4Ъ2 Pi® + Р2(а-а)- Р4в = PiP + Р3 (а-а)- Р4о

Р2в + Р3о+ Р4 (а-а)' где р j (j = 1,4 ) - параметры. Возможны три случая: i) Ъ1 Ф 0, Ъ2 Ф 0, 2) = 0, Ъ2 Ф 0, 3)Ъ1 ф 0, Ъ2 = 0. В 1) и 2) положим Р2 = 0. Тогда вышеуказанное условие разрешимости примет вид:

(а - а)2 =

— + ^ +

Р4 Pi Pi Р4

2 л / Р3 ß®-(

2

(а - а)2 =

1

Pi

2 Л

+ Рз +-Рз Рз

ß®- (ß2 +®2),

Последние равенства выполнены при а < 0 и а > 0, если ^3 - достаточно большое положительное

число, а р12 =в р3 -1. Лемма 7 доказана.

а

Лемма 8. Пусть матрицы Л и В те же самые, что и в лемме 6.

Тогда существует матрица £ такая, что матрица Л + ВБ имеет комплексно-сопряженные собственные значения с отрицательной вещественной частью.

Доказательство. По лемме 6 существует матрица такая, что Л1 = Л + подобна некоторой матрице М1 с собственными значениями = а1 ± а (они же и собственные числа матрицы Л1, так как Л1 подобна М1), для которых выполняются соотношения (7), (11), где Н и К определяются равенствами (8) или (9), или (10) в зависимости от значений Ь1 и Ь2 . Таким образом:

а1 =а + вЦр0), а =вк(р0) (20)

причем

N(p0 )< 0, к(р0)> 0. (21)

Из (20), (21) видно, что а1 <а и а>1 > 0 (так как Р > 0 ). Если а1 < 0, то лемма 8 доказана.

Пусть а1 > 0. Тогда, применяя лемму 6 к матрицам Л и В найдем, как и выше, матрицу Б2 такую, что Л 2 = Л1 + ВБ 2 имеет собственные числа /и2, /и2 =а2 ± 1а2 , причем а2 = а1 +а1 к(р0)

в

а2 < а1 и а2 > 0.

Очевидно, что Л 2 = Л + В(Б1 + Б 2) и

а 2 =а+(в + а )N (р 0), а 2 = рК(р 0).

Продолжая этот процесс последовательного применения леммы 6 к матрицам Лк ((= 1,2,.) и В (здесь Л0 .•= Л), будем получать Бк+1 (к = 0,1,2,. ) такие, что Лк+1 = Л + ВБк+1, (к = 0,1,2,. ) имеет соб-

®2 =®1к(р0). Отсюда в силу (21) следует, что

м. + р3 1® -р~.

р4Ь2 р4 ) Ь2 В случае 2) последнее равенство имеет вид:

р3

а -а =-а.

р 4

Последние равенства выполняются при а < 0 и а > 0, если: а) Ь1Ь2 < 0, р3 = 0 , р4 = 1, а р1 - достаточно большое положительное число; Ь) Ь]Ь2 > 0, или Ьх = 0 , Ь2 Ф 0, рг = р4 = 1, а р3 -достаточно большое по модулю отрицательное число.

В случае 3) положим р2 = -1, р4 = 0 . Тогда условие разрешимости уравнения (13) принимает вид:

ственные числа ^k+l, /лк+1 = ak+l ± iok+1, где

ak+i =ак +ок N(p 0 )=... = а + (в+®1+к +ок )(p 0 ),

®k+1 = ®k К

(Р1 +1)®= Рз ß.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

:(р0 )= к = р (к(р0 )). (22)

Здесь а0 := а, а0 := р.

Поскольку в > 0, ау > 0, у = 1,.,к , то в силу (21), (22) существует такое значение к = ", что собственные числа /ип+1, /ип+х = ап+1 ± ¡а^ матрицы

Л "+! = Л" + ВБ"+1 = к = Л + В(£ + к + Б") будут удовлетворять неравенствам

а"+1 < 0, а"+1 > 0 (23)

Полагая Б = + „. + Б", получаем, что для собственных чисел ап+2 ± ¡а^ матрицы Л + ВБ имеют место неравенства (23). Лемма 8 доказана.

Замечание. Утверждение леммы 8 справедливо также в случае, когда В - матрица-столбец (см. лемму 7).

Теорема о стабилизации

Докажем теперь основное утверждение. Теорема (о стабилизации). Пусть А и Ь - произвольные вещественные (п х п) - и (п х т) - матрицы и пара (А, Ь) полностью управляема.

Тогда существует вещественная (п х т) - матрица

5 такая, что матрица А + Ъs * гурвицева.

Доказательство. Пусть матрица А имеет к вещественных у\,---,7к и некоторое число 1 комплексно-

сопряженных пар А1,Л1,...,А1 ,А1 ; А],А}- = а ±'в ],

в] * 0 (] = 1,1), к + 21 = п, собственных чисел (с учетом их кратности). Не умаляя общности, можно считать, что в] > 0.

Приведя матрицу А с помощью преобразования подобия А = Q0-AQ0, где Q0 - неособая матрица (det Qo * 0), к вещественной «нижней» жордановой нормальной форме [9,10], а затем перенумеровав подряд все (2 х 2) - матрицы и числа, стоящие вдоль

главной диагонали, запишем А в виде:

( С1 О 2.......О 2 о 2

Р 2 С 2....... О 2 О 2

A =

O 2 O 2.....C, O 2

O 2 O 2........ F 2 C,

O

г1 0......0 0

/ г2 .... 0 0

где

C- =

O

Ч -Ч

0 0.

. /1 г

ч 0 0

F =

1 0 0 1

1 * *

или

(24)

F = O2 =|0 0|; /1 = 1, или /1 = 0.

Пусть

~ = Q0:b =

b11 b12

Vb21 b22 у

(да > 2)

(25)

(26)

Ai =

(27)

-* I sii O si =

(да > 2),

(28)

КО О

где (2 х 2) - матрица 51'1 подлежит определению (в случае т = 1 - матрица-строка размерности 2).

Из (26) - (28) получаем, что А + Ъ имеет вид (27), если выполнено равенство

С + ¿4: = Л. (29)

В силу леммы 8 (замечания к лемме 8, если т = 1) существует 51'1 такая, что имеет место (29). Очевидно, что А1 = А + Ч* = Q 01(А + Ъ ?1* Q 01 ^ 0.

Так как пара (А, Ъ) полностью управляема, то, по

лемме 2, полностью управляема также и пара (А1, Ъ), * * ~ * 1 где А1 = А + Ьз1, 8! = 51 Q0 .

Так как А1 подобна А1 , то в силу (27)

а(Б) = а (Г) и {А2,АА2;кАг,А1 ;у1,...,ук

Здесь а обозначает спектр соответствующей матрицы.

2) Далее в (25) перестановкой 1 и 3-й, 2 и 4-й строк и 1 и 3-го, 2 и 4-го столбцов переставим местами матрицы Г1 и С2, расположенные вдоль главной диагонали. Эта операция, как легко убедиться, равносильна умножению А1 слева и справа на некоторую ортогональную матрицу Т1 , получающуюся из единичной перестановкой 1 и 3-й, 2 и 4-й строк. Матрица Т1А1Т1 имеет те же собственные числа, что и А1, но она, вообще говоря, уже не имеет нижнего треугольного вида (24), (25). Сделаем над ней преобразование подобия с матрицей подобия Р1, приводящее ее к вещественной «нижней» жордановой нормальной форме

А = Q1-1А^, где Q1 = Т^.

Пусть b = Q-1 b =

(~11 ~12 ^

b21 b22

(да > 2 ) (в случае

( в случае т = 1, Ъ11 и Ъ12 - одностолбцовые матрицы размерностей 2 и п-2 соответственно).

Так как пара (А, Ъ) полностью управляема, то, согласно лемме 5, Ъ11 * 0 .

1) Построим (п х т) - матрицу ~1 такую, чтобы А1 = А + Ъ имела вид Ч О

Ч Ч

где Г1 - гурвицева (2 х 2) - матрица; Ох - некоторая (п - 2) х 2 - матрица; - (п - 2) х (п - 2) - матрица,

получающаяся из А вычеркиванием первых двух ее строк и столбцов.

Будем искать ~1 в виде

т = 1 Ъ11 и Ъ21 - одностолбцовые матрицы размерностей 2 и п-2 соответственно).

Поскольку пара (А1,Ъ) полностью управляема, то,

по лемме 3, полностью управляема и пара ( А1 ,Ъ ) .

Поэтому в силу леммы 5 Ъ11 * 0 .

Как и в п. 1), построим (п х т) - матрицу ~2 та-

чтобы А2 = А1 + Ъ ~2*

кую,

A2 =

имела

вид

(Г 2 O ^

где Г2 - гурвицева ( 2х 2)- матрица;

02 - некоторая (п - 2) х 2 - матрица; О 2 -(п - 2) х (п - 2) - матрица, спектр которой равен

а(0.2) = а (Г) и Аз, Аз;... а А ;п,.,Гк }

тт «

Для этого, как и выше, найдем матрицу 52 вида

(28) такую, что имеет место равенство С2 + Ъ11«"1 = Г 2, где - ненулевая (2 х 2)- подматрица матрицы ~2*.

а

Учитывая, что А1 = Адо, Ь = её^ получаем:

А = А + ь?;=е1-1бА41 + ¿^гег1 ^ = = е 1-1 е о1 (а 1 + ьз г о е 1, где з2=^е^ео1-

Очевидно, что матрица

А2 = А1 + Ьзг = А + Ь^" + имеет спектр

У ( a2 ) = У (Г) и .У (Г 2) и {л

3 л 3/к л ¡,л ¡; г1,...,г k

г k}.

венными числами а j ± im j

шагов, очевидно, мы получим А(1 вида (24), в которой

все С}- заменены на Г, а у1,. и в (24). При этом будем иметь

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

,Yk те же самые, что

A? = Qi\ ...Qo(A + bs'c)Q0 ...QM

sc = Z s,

(30)

j=1

где ео >■■■ ,е-1 - неособые матрицы подобия,

¿у = 3]-е]--1 ■•• ео а = 1,1) , - (т х п)- матрица,

имеющая вид

sj =

(sj O

O O

J)

(m > 2,j = l,l) .

Г} = С, + bVsj

b() = Qj-1 .Q01 b =

( b()

h

11 b12

v b21

(j) ,(j)

22

(j = 1, l) - матрица

порядка (п х т) ; Ь^ - (2 х 2) - подматрица (при т=1

ьЦ и Ь21 - матрицы- столбцы размерностей 2 и (п-2) соответственно).

Заметим, что в силу леммы 5 ЬЦ Ф о. Поскольку

в силу (3о) А1 = А + Ьз" подобна А(,) , то ее спектр имеет вид у(А,) = {у 1 ± 1щ1,.,у , ±¡щ,}{,...,гк}, где

< о, ю} > о (а = 1,.,/).

Перейдем теперь от Аг к ее вещественной «нижней» жордановой нормальной форме А,. Сделаем в А1 ряд перестановок строк и столбцов (что равносильно умножению А1 слева и справа на некоторую ортогональную матрицу Т, получающуюся из единичной перестановкой некоторого числа соответствующих строк), запишем полученную матрицу ТА1Т в виде:

TAlT =

Продолжая этот процесс, последовательно заменим двумерные матрицы С(} = 1,1) (с собственными числами а а ± 1ра) на некоторые гурвицевы Г} (с собст-

<у а < о, со а > о). Через I

( у1 0...0 0 /1 Y2 .00

0 0 к Yk-10

0 0. /1 Yk

O

O

Г1 O2 .O2 O2

F2 Г2.

•O2 O2

О2 О2 -Г, -1 О2 «2 «2 - F2 Г^ где /1 и определяются из (25), а

(32)

Г j =

(aj -mjЛ m j а j

v

а. < 0, m. > 0 (/=1,l).

Пусть т1,---,тк - произвольные вещественные числа. Проведём в (32), аналогично как и выше, процедуру последовательной замены собственных чисел Y1,■■■,Y2 соответственно на числа т, • • •, тк . В результате получим матрицы

А1+1 = А1 + Ьз*+1, А1+к = А1+к-1 + Ьэ*+к , где

* п-1 * п-1 п-sl+1 = sl+1Ql , •••, sl+k = sl+kQl+k-1 ■■■Ql

(33)

(31)

Здесь - (2 х 2)- матрица (при т = 1 в (31) -матрица-строка размерности 2) такая, что гурвицева (в силу леммы 8);

(а) ^

такие, что спектр ст1 (А,+д) =

= {о"1 ± ¡Ю1,...,а1 ± ¡ю,} и {Т1,...,гч } и {уч+1,---,Ук }, (Ч = 1,..., к). (34)

тч-Уч.

Здесь sl+q II Olq+P O),

b (l+q) = Ql+1q- ...Qrlb = (h (l+q) u 11 b (l+q) v" 21

Pl+q

12 b (l+q) u22 У

b

(l+q)

(35)

где е1, ■, ег+к-1 - некоторые неособые матрицы

матрица-строка размерности

(b(l+q) - число; b,(l+q)

12

т -1; Ь21+Ч) - матрица-столбец размерности п -1).

Заметим, что в (35) Ь^+Ч) Ф о в силу леммы 4. Из (34) при ч=к получаем:

°(А1+к) = {^1 ± ¡Ю1,...,а1 ±ю,} и {Т1,..., Тк }. (36) Поскольку А, = А + Ьз*, то А,+к = А + Ь(зс + ^ )*, (37)

* * где SY = Z sl+q .

q=1

Положим

s = sc + SY .

(38)

(39)

Тогда с учетом (37) соотношение (36) можно переписать так:

а(А + Ья *) = {о"1 ± ¿Ю1,...,а1 ± ¡ю,} и {Т1,..., Тк },

<о, ю} >о. (4о)

Если числа тч (4 = 1, к) выбраны отрицательными,

то соотношение (4о) означает, что матрица А + Ьз * гурвицева, где стабилизирующая матрица з определяется по формулам (33), (35), (38), (39) и (3о), (31).

Тем самым, теорема о стабилизации полностью доказана.

Замечание 1. Нами доказано несколько более общее утверждение: для произвольной вещественной матрицы А, спектр которой состоит из к вещественных и I комплексно-сопряжённых пар, и произвольного набора {тц }к=1 вещественных чисел существует вещественная матрица 5 такая, что имеет место соотношение (40).

Замечание 2. В ходе доказательства теоремы о стабилизации дан конструктивный метод построения стабилизирующей матрицы 5.

Литература

1. Попов В.М. Гиперустойчивость автоматических систем. М., 1970.

2. Зубов В.И. Теория оптимального управления судном и другими подвижными объектами. Л., 1966.

3. Аксенов Г. С. К задаче стабилизации линейного объекта управления // Вестн. Ленингр. ун-та. Сер. i. 1997. № 7 С. 5-8.

4. Гелиг А.Х., Леонов Г.А., Якубович В.А. Устойчивость нелинейных систем с неединственным состоянием равновесия. М., i978.

5. Смирнов Е.Я. Стабилизация программных движений. СПб., i997.

6. Wonham W.M. // Trans. Aut. Contr. i967. AC-i2. № 6. P. 660-665.

7. Wonham W. M. Linear Multivariable Control: a Geometric Approach. New York, Heidelberg, Berlin, i979.

8. Леонов Г.А. Математические проблемы в теории управления. Мотивация к анализу. СПб., i999.

9. Гантмахер Ф.Р. Теория матриц. М., i967.

10. Демидович Б.П. Лекции по математической теории устойчивости. М., i967.

Санкт-Петербургский государственный университет, Адыгейский государственный университет

10 июня 2004 г.

i9

i Надоели баннеры? Вы всегда можете отключить рекламу.