Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Матричное исчисление с приложениями в теории динамических систем

..pdf
Скачиваний:
7
Добавлен:
19.11.2023
Размер:
36.1 Mб
Скачать

АСИМПТОТИЧЕСКОЕ РАСЩЕПЛЕНИЕ УРАВНЕНИЙ МНОГОМЕРНОЙ НЕСТАЦИОНАРНОЙ ЛИНЕЙНОЙ СИСТЕМЫ ВЫБОРОМ СТРУКТУРЫ МАТРИЦЫ УПРАВЛЕНИЯ

§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)

Смысл верхнего индекса будет понятен из дальнейшего.

Q — некот орая невырожденная квазидиагональная мат рица по­ рядка п, приводит (26.1.6) к виду

| = ( Л М - 0 - ' | ) у + !* + ' г и . + | ^ .

(26.2.5)

Лл

Если матрицы Q, у и М = К~1 представить в виде блочных мат­

риц, а именно

 

 

 

 

(ЙГЛ

S2 — ding ( й |, Q 2,

й р ) » у У2

й г

II

 

• • •

 

 

W

Л

то легко видеть, что уравнение (26.2.5) распадается на следующие р уравнений:

dп

(Лй"> - »а-йГ)Уо - £,и + 1М аА т+ 1К у

(26.2.6)

( 0- = 1 , 2 , ..., р).

П р и м е ч а н и е 26.2.1. Введение матриц Z и Q (структуры ко­ торых мы уточним позже) с точки зрения квазидиагонализации не диктуется необходимостью. Эти матрицы в преобразовании (26.1.7) играют роль "резервных” множителей, с помощью которых преоб­ разованию (26.1.7) можно придать те или иные свойства, облегча­ ющие решение основной задачи. ■

М атрицы

А ^ п\

а также

матрицу В

представим в форме

конечных сумм:

 

 

 

 

 

т

 

 

т

KSm>=

К + 2

ъкК ш ,

Л<™> = А +

2 е*А[*].

 

* = 1

 

 

*= 1

 

 

 

т

(26.2.7)

 

 

 

 

B = £ a + 2 ^ B t . k= 1

После подстановки сумм (26,2.7) в (26.2.3) условие тождественного относительно е выполнения соотношения (26.2.3) приводит к равен­ ствам

(А + НВ0) К = К А ,

(26.2.8)

(А + НВ0)К'Н = /С|41Л + ЛГДШ +

+ H B kK + D U - »

( * - 1 , 2 , . . . , ш ) ,

(26.2.9)

т то

А п + 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 на непересекающиеся группы, группа «сг» бу­ дет как раз состоять из собственных значений субматрицы А а.

Если заданными являются сами субматрицы

A0 (or =

1 , 2 , . . . , р ) , то, во-первых, можно задачу свести к

только что

рассмотренному случаю, если определить собственные значения субматриц A 0 (cr = 1 , .... р ). Но интереснее провести построение^

и В 0 без определения собственных значений субматриц Лс, а непос­ редственно через эти субматрицы. Попробуем это. Имеем

рк>

A„(i/) = П п (U-X{pE).

 

 

 

5 = 1

j = I

 

 

 

 

S^CF

 

 

Подставляя сюда значение U =

К A M , получим

 

/

Р

Iг

\

 

 

Л о = *

П

П (Л-*'/>£) М = К А 0(А )М .

(26.3.10)

 

5 в 1

1

 

 

 

5^0

 

 

 

 

Отсюда

 

 

 

 

 

 

 

Ав(Л)

M A a(U )K .

(26.3.11)

Имеем

 

 

 

 

 

р

Д«(Л) = П П -*$*>£).

5 = 1 У = 1

s*a

Но

(А,

> 1

Jt4 ** II

0

Учитывая это, имеем

р

Д„(Л) = П п (л -

0 '

[Я*,

0 '

 

 

 

-Х$0

 

 

 

Лв

0

7

 

 

Р)

\

о

\

Ч , -

 

 

л 2 -

х<;>£4>

 

 

 

 

 

 

О

Л 0 -

№ Е к

 

р

J

%

=

s - 1 /«= 1

Р

к *

~ *‘/ 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

 

будем иметь

 

Д0 (Л ,)

0 )

До(Л) =

Д .(л „ )

о

л > „ )

Матрицы (26.3.13) обладают свойством

Д о(Л^) = 0

(ц ^ о г ).

Р *s

(26.3.14)

(26.3.15)

В самом деле, при р ф а произведение ]~J

(Л ц — № рЕк ) содер-

 

5= 1 / = 1

 

5*0

 

ж ит

множитель J~[ (Л^ — )Sj^Ek ), который

согласно теореме Га-

 

>=i

 

мильтона— Кэли равен нулю. При р = о

 

 

До(Л о) ^ 0 .

(26.3.16)

так

как

 

р*.

д0(лс) = п П (л» - XY 'E 0 * 0

j = 1

5 * 0

в силу того, что матрица Д0 № Е к при s ф а невырождена и от-

лична от нуля.

 

 

С учетом (26.3.15) и (26.3.16)

имеем

 

/ О

0

\

Д 0 (Л ) =

Д0 (Л 0)

(26.3.17)

о

о

 

где

рА«

до(ло)= П П (А. - XY E 0 * °-

J = I / = 1

Согласно (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> которые

Соседние файлы в папке книги