
книги / Матричное исчисление с приложениями в теории динамических систем
..pdfАСИМПТОТИЧЕСКОЕ РАСЩЕПЛЕНИЕ УРАВНЕНИЙ МНОГОМЕРНОЙ НЕСТАЦИОНАРНОЙ ЛИНЕЙНОЙ СИСТЕМЫ ВЫБОРОМ СТРУКТУРЫ МАТРИЦЫ УПРАВЛЕНИЯ
§26Л. Постановка задачи
Вэтой главе рассматривается система автоматического управле ния, представленная системой уравнений вида
^ = A (t)x + H(t)u + h(tt х ), |
(26ЛЛ) |
где х — вектор фазового состояния {* е /?'*}, и — вектор управлений {и G /5/}, A(t), Н(1) — матрицы с размерами соответственно п х п,
п х I |
и элементами, дифференцируемыми по t на интервале |
Д = |
[?0,Г], где Г — число, превосходящее tQ, или символ «>, любое |
нужное число раз, h(t, х) — столбцовая матрица (типа л х 1), эле менты которой — нелинейные функции отношений Xj такие, что рав номерно по t на [i0,T]
Нт |
hit, х) |
= 0. |
(26Л.2) |
JC-»0 |
11*11 |
|
|
Связь между вектором управлений и вектором состояния рассмат риваем в виде
u = B(t)x, |
(26Л.З) |
где B(i) — матрица с размерами п х I и элементами, дифференци руемыми по t любое нужное число раз.
Рассмотрим следующую задачу: можно ли соответствующим вы бором матрицы K (t) линейного преобразования
x = K (t)y |
(26.1.4) |
и матрицы B(t) закона управления |
(26.1.3) добиться того, |
чтобы в новых переменных у система |
уравнений управляемого |
процесса имела |
бы в линейной |
части квазидиагональную |
м атри |
ц у , а именно |
|
|
|
^ = |
A (f)y + <p(f, У), |
А ( 0 = diag{Aa(*)}- |
(26.1.5) |
Д л я реш ения поставленной задачи, следуя не раз использованному
в предыдущ их разделах |
приему, введем в рассмотрение вспомога |
||||
тельную дифференцируемую систему |
|
|
|||
|
^ 7 = |
Л(т)х + |
Н(х)В(т, г)х, |
|
(26.1.6) |
коэф ф ициенты |
которой |
зависят |
от так называемого |
м едленного |
|
врем ени т = е t. |
Будем |
предполагать, что в векторно-матричном |
|||
уравнении (26.1.6) А ( т) |
— матрица с размерами л х |
п, |
непрерыв |
||
н ая н а [0, JL ], |
В ( х, г) — матрица с размерами I х п, |
непрерывная |
по обеим аргументам и регулярная по е в окрестности точки е = 0 .
П ри е = 1 уравнение (26.1.6) эквивалентно однородной системе |
||
(26 .1 .1)— (26.1.3) |
при |
h (t, х ) = 0, так что все формальные постро |
ения, которые |
будут |
получены применительно к уравнению |
(26.1.6) и будут выполняться тождественно по е, будут применимы |
и к исходной однородной системе, если в этих соотношениях при нять е = 1. Учитывая это, построим сначала аппарат формального расщ епления (26.1.6) путем соответствующего выбора матрицы ли нейного преобразования (26.1.4) и матрицы В ( т, е).
§26.2. Алгоритм формального расщепления уравнения (26.1.6)
Формальное расщепление уравнения (26.1.6) определяется сле дую щ ей теоремой.
Те о р е м а 26.2.1. Замена переменных
х = К у, |
(26.2.1) |
где |
|
K = l 6 m^ZQ, |
(26.2.2) |
причем — невырож денная дифференцируемая мат рица* по р яд ка п, удовлет воряю щ ая уравнению
е |
= (А + Н В + |
em + l Am + 1) K t”* - |
К<т^А^т\ |
(26.2.3) |
|
Z — невырож денное реш ение мат ричного уравнения |
|
||||
|
= A (")Z - |
ZA<m>, |
Z (/„) = |
Е п, |
(26.2.4) |
♦ Смысл верхнего индекса будет понятен из дальнейшего.
т то
А п + 1 — 2 |
Бт “ 1 (АГ, т _ а + ',1 л 1а} — |
|
|
|
v в I а = v |
|
|
|
|
|
- H B aK [m- « + v]) + |
ат |
M^>. |
(2 6 .2 . 1 0 ) |
|
|
|
|
Здесь М ^> — матрица, обратная К^т\ a D 1' 1 — квадратная матри
ц а, |
не зависящ ая от |
X IrJ, A lrl |
и |
В г при r > i. |
Так, |
||||
|
|
D [0] = |
^ |
, |
Я*1* = |
/с‘Чл1Ч |
и т .д. |
||
|
Д алее будем пользоваться обозначениями: |
|
|||||||
|
^ |
( Л |
^ |
* 2 ... К р), |
|
|
|
||
|
|
- |
|
( * J « |
|
Л ™ ) , |
(26.2.11) |
||
|
Л 1*1 = |
diag |
(Л}*1, А ^ 1, ..., A lpk]), |
|
|||||
где |
К а, |
( о = |
|
1, 2 , . . . . р; |
Л = 1, 2 , . . . , т ) |
— субматрицы м ат |
риц К и К\® — вообще говоря, прямоугольные матрицы, размеры
которых уточняю тся ниже, а Л ^ ] — субматрицы квазидиагональ-
ной матрицы |
— квадратные матрицы, размеры которых такж е |
уточняю тся ниже. |
|
§26.3. Исследование уравнения (26.2.8)
Имеем
(А + Н В 0) К = К А . |
(26.3.1) |
М атричное уравнение (26.3.1) при заданных Л и Н связы вает м еж ду собой три матрицы К ,А и В0. Аналогичным образом при задан
ных А и В0 это уравнение связывает тройку матриц К , А, Н .
В конкретных задачах, особенно в задачах синтеза, возникает потребность конструировать одну или две матрицы из этих троек так , чтобы третья матрица имела бы вполне определенный вид. Т ак возн и кает проблема: можно ли по каждой матрице из одной или другой тройки определить две другие из той же тройки и, если м ож но, то как это сделать? В соответствии с такой постановкой рас смотрим различны е возможные ситуации, которые могут возникать при реш ении практических задач, не вдаваясь пока в существо этих задач . Подробное рассмотрение всевозможных комбинаций, мысли мы х при такой постановке, потребовало бы слишком много места и вн им ания, и потому выходит за рамки настоящей книги. Мы огра-
ничимся лишь рассмотрением некоторых, на наш взгляд, наиболее интересных для приложений случаев, проведя это рассмотрение до конца, а в остальных случаях ограничимся лишь формулировкой алгебраической задачи.
Матрица BQзадана. Требуется построить матрицы К и Л. В рас сматриваемом случае матрицы А, Я, и В0 являются известными и требуется построить К и Л так, чтобы матрица Л имела заданную структуру. Мы здесь будем требовать, чтобы Л имела квазидиаго-
нальную структуру. |
|
|
|
|
|
Таким образом, |
в |
данном |
случае |
квадратная |
матрица |
А + ИВ0 = U является известной |
матрицей. Допустим, |
что квад |
|||
ратная матрица U |
порядка п |
имеет |
собственные |
значения |
|
X,, Х2,..., ^п> являющиеся |
функциями от |
/ и эти собственные |
значения могут быть разбиты на некоторое число р непересека-
ющихся |
на |
рассматриваемом |
интервале |
[t0 |
Т] |
групп |
||
Чо)(0> ^4о,(0 . •••> 4 о)(0 |
(°'= Ь 2, •••> Р; |
2 , к о = п) |
так> чт0 |
|||||
|
|
|
|
|
а= 1 |
|
|
|
|
|
|XW(0-X</)(0l г о |
О |
|
|
(26.3.2) |
||
|
(о ^ |
s, i = |
1, 2, ..., |
/ |
1,2,..., |
^ ) . |
|
|
|
|
|
В этих предположениях могут быть построены матрицы К А, удов летворяющие уравнению (26.3.1), и при этом А будет иметь квазидиагональный вид А = diag(A,, Л2, ...» Ар).
Метод построения этих матриц с подробным обоснованием изло жен в гл. 5 настоящей книги и коротко сводится к следующим по строениям. Каждой группе сг собственных значений матрицы U ста вится в соответствие многочлен
р*»
Ао(*) = П |
П |
(a=s 1.2»—.р)- |
(26.3.3) |
5*1 |
у»1 |
|
|
5^0 |
|
|
|
Ранг матрицы А0(Я), которая получается Да(Х) после замены X из на U , оказывается в точности равным числу касобственных значений матрицы £/, включенных в группу сг. Поэтому матрицу Aa(U) можно представить в виде произведения матрицы Катипа п х ка с калиней но независимыми столбцами на матрицу М0отипа кд х п с калинейнб независимыми строками:
д ,{V) = К„ма .
О
Далее вводится матрица |
|
|
М , = (М0 КаУ 1М0 . |
(26.3.5) |
|
О |
о |
|
(Здесь мы использовали то обстоятельство, что М0 Ка — невырож-
о
денная квадратная матрица порядка ка.) Построенные таким обра
зом матрицы Ка и Ма (<т= 1, 2 |
, р) связаны между собой соот |
|
ношениями |
|
|
МаК = Ь ,Е к , |
= {‘ g “ g ; |
(26.3.6) |
Ek — единичная матрица порядка ка.
о
Пусть К JL М — блочные матрицы порядка п, составленные из построенных указанным выше способом субматриц Ка , Ма , так
что
К = ( Kl K2 ... Кр), М = col (M t Мг ... Мр).
В силу (26.3.6) очевидно
КМ = М К = Еп. |
(26.3.7) |
Наконец, субматрицы А0 квазидиагональной матрицы А опреде ляются соотношениями
А0 = MaUKa . |
(26.3.8) |
Рассмотрим другой случай. Пусть теперь заданы субматрицы Ла(а = 1, 2, ..., р) матрицы Л или же собственные значения этих
матриц, а именно
Ма)>•••> Чо) |
(<* = !» 2> |
Р)- |
(26.3.9) |
о |
|
|
|
Требуется определить В0 и |
К , связанные |
с Л |
соотношением |
(26.3.1).
Алгоритм построения матриц Ка, М0, Ла(ст= 1, ..., р), приве
денный в п. 26.3.1, можно использовать в рассматриваемом случае для построения К и М , а затем и BQ. В самом деле, в случае за
дания собственных значений субматриц Л0 вопрос решается авто
матически, поскольку собственные значения этих субматриц как раз и образуют в совокупности все множество собственных значе ний матрицы А + НВ0 = U, причем, при разбиении собственных
значений матрицы U на непересекающиеся группы, группа «сг» бу дет как раз состоять из собственных значений субматрицы А а.
Р |
к * |
~ *‘/ 4 , |
|
о |
Л г - Х ^ |
|
|||
|
|
|||
-П |
П |
|
|
|
|
|
|
||
5 = 1 |
j = 1 |
|
Лр - |
Х р Е |
s=о |
|
0 |
||
|
|
|
/
П (Л, -
Я |
П ( Д 2 - х^ |
2) |
|
|
|
= П |
1 |
|
|
|
|
5 = 1 |
|
|
5 = 0 |
|
|
|
О |
П ( Л р |
\ |
|
у = 1 |
Т аким образом, |
|
|
{ Р |
Л5 |
|
п |
П (Л. - Ц ‘ >Е0 |
О |
5—1 У= 1
р
П ГКЛг-^Ч,)
5=1 У= 1 5*0
Л *(Л ) = |
(26.3.12) |
Я |
* |
П П(Л„-*$»£.)
5=1 У= 1 5*0
|
Я |
** |
о |
П |
П (Лр - X^l£t ) |
5=1 У= 1
5 * 0
Вводя обозначение
р |
*» |
Л„(Л ц) =П |
П (А-- х р Е „ ) (и = 1, 2....р), |
5 = 1 |
У = 1 |
5 * 0 |
|
Согласно (26.3.10) имеем
|
|
|
|
|
|
|
|
'м / |
|
Д„(to = К А „ (А )М = (K L К 2 ... |
|
Д„(Л») |
|
Щ |
|||||
К .) |
|
•« |
|||||||
|
|
|
|
|
|
|
|
м п |
|
|
|
|
|
о |
|
|
о/ |
р) |
|
|
|
|
|
/ |
О |
|
\ |
|
|
|
|
|
= (К 1 ... К „ ... К ) Д с ( Д . Ж |
= |
К оА а( А с) М 0. |
||||
|
|
|
|
|
* |
|
|
|
|
|
|
|
|
|
* |
|
|
|
|
|
|
|
|
|
О |
|
|
|
|
И так, |
|
ДAU) = К аА „(А а) М а. |
|
|
|
(26.3.15а) |
|||
|
|
|
|
|
|
||||
Отсюда легко получить и такие соотношения: |
|
|
|
|
|||||
|
|
|
M aA a(U ) = A a(A a)M a, |
|
|
(26.3.16а) |
|||
|
|
|
A a(U )K a = A a(A a). |
|
|
(26.3.17а) |
|||
Т аким |
образом, если заданы субматрицы |
А а |
(or= 1, 2 , . . . , р ) или |
||||||
ж е их |
собственные значения, то могут быть |
построены |
м атрицы |
||||||
М 0 , К 0 (or = 1, 2 , . . . , р ), а значит, и блочные матрицы М |
и ЯГ, ко |
||||||||
торы е |
связаны |
с |
квазидиагональной матрицей |
А |
и |
м атрицей |
|||
U = А + Н В 0 соотношением (26.3.1). |
|
|
|
|
|
||||
Теперь займемся построением В0 при заданных |
А |
и Н . И з |
|||||||
(26.3.1) |
имеем |
H B Q= К А М — А. Отсюда, |
если Н имеет псевдооб- |
||||||
ратную матрицу Н + , |
|
|
|
|
|
|
|||
|
|
|
В0 = Н +( К А М - А ). |
|
|
(26.3.18а) |
|||
§ 26.4. Определение К ^ , |
и Вк |
|
|
|
|
|
|||
Т еперь приступим к построению К [к], |
А 1*1 |
и |
В к, |
удовлетворя |
|||||
ю щ их уравнениям |
(26.3.9). Имеем |
|
|
|
|
|
|||
|
UK.W = К ^ А + K A W + H B kK + |
|
|
|
(26.4.1) |
||||
|
|
|
U = A + H B 0 ( * = 1 , 2 , . . . ) . |
|
|
||||
|
|
|
|
|
|
||||
К ак |
было показано выше, в каждом |
из рассмотренных выше |
|||||||
случаев |
могут |
быть построены |
АГ^ = К , |
A 101 |
s=A и |
В 0> которые |