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

книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем

.pdf
Скачиваний:
97
Добавлен:
25.10.2023
Размер:
17.28 Mб
Скачать

ЗОО

П Р И Б Л И Ж Е Н Н О Е И Н Т Е Г Р И Р О В А Н И Е

У Р А В Н Е Н И Й

[ГЛ

X I

Из (3.5) следуют равенства

 

 

 

 

,

.

dK„(т, е)

~

~

 

 

 

 

 

4 ( т ) | е -----^ -------1- Ко (т, е)К (г, е)

 

 

 

 

 

 

 

 

=

В (т)Ка(т, е ) + е (х)Я

І„

( 3 . 6 )

 

 

 

(ст=

1 , 2

, . . . .

п),

 

 

 

где столбцовая матрица /

0 типа I х

1

имеет вид

 

 

 

/

 

 

 

 

 

 

 

 

/<т=

]

G(/ — 1',х')Т (х')Ка(х', е)ехр[0ст(/',е) — Qa(t, &)]dt',

оо

аѲ„ - функция, удовлетворяющая соотношению

дѳ_

 

 

 

 

 

~0Г =

(т> б)-

 

 

 

 

После замены переменных t

f

=

s

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

/о =

I G (s, т — es) Т (т — es) Ка (т — es, е) ехр [Ѳа (/ — s, е) —

 

о

 

 

 

 

 

 

 

 

 

 

 

0 а (/“, е)] e(s.

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя разложения

 

 

 

 

 

 

 

 

 

G ( S ,

т -

es) =

G(s, т) -

es -

^

т )

+

- 1

eV

Дт3

____________

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т

E S ) — Т (т) — es

 

(Т)

,

1

г2

„ 2

 

Я)

 

Яа (т — es, е) = Ка (т, е) — es

 

(t, в)

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

]

„ _ d2/(a (т, 8)

 

-

ехР

 

S ,

е) — Ѳа ( ( , е)] =

~ 2 “ е s --------------- 5^5-------------------------• • •

 

 

 

 

 

 

~

 

 

I

 

ДА._ (т, е)

 

 

 

 

ехр

sXa(т, е) -J

 

ES~ ■

 

dr

 

 

 

1

e2s3

d 2 X g (

Т , 6 )

 

 

 

 

 

 

1

,

d X a (т, е)

 

 

 

 

 

 

 

 

 

 

 

 

 

6

Дт3

 

 

=

' .

1 +

 

2

ËS

Дт

 

-- '

„2„3

d2 M

T. e)

 

+ 4

1

 

,

<^а (т- е)

 

 

 

 

dt*

+

 

T

es-

dt

 

$ 3)

 

 

 

ИНТЕГРИРОВАНИЕ

УРАВНЕНИЙ

 

301

 

 

2с8 d2K (т , в)

 

 

 

exp [— s%a(T, e)] =

----- рг- E S'

 

dz2

+

+

 

 

О

 

 

 

 

-

1+1

>

 

^

L

+

e-

dXa (T, e) ^

 

 

 

 

di

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S3

d3X0 (г, e)

+ •••

exp[— sa,0

(T,s)],

 

 

 

 

 

Т Г

dT2

 

 

 

 

 

 

 

 

 

 

интеграл Ia можно представить в виде

 

 

/о =

£оо(Я0>

т) Т (т) Яа (т, е) +

 

 

 

 

 

I R10 (ka, т)

d

[ Т (т) К а

(т, е)]

 

 

 

 

-)- Е

 

 

dT

+ Ru (^о,

T (т) /Со (^І в) -j-

где

+ 4

~

^

а>

т) т ( х

) ^т’ е)

dl°dT1 8'

i + е2 •••>

(3-7)

 

 

 

 

 

 

 

со

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ ? 0

0(^. /) =

/? (К, () =

j

G(.S, 0 e -Xsds

 

— матрица передаточных функций регулятора (с размера­ ми / X т) с параметрами, замороженными в момент време­ ни і,

Ru {К t) =

* +,* b < L в ( - і ) ' Г — і^-0- s ^ d s

'

<jjiw

v

'

J

dt>

 

( / , / = 0 , 1 , 2 , .

о

 

 

 

 

Функциональные матрицы

 

 

(Ä.0,

т) в силу второго

соотношения (3.4) в свою очередь допускают следующие

разложения

по степеням

е:

 

 

Ri 't (^о> Т-) =

R{j (^a, "О"Ь В^сг ^ (т) /?(+] /

т) -j-

 

2 /-O ^Ri+

1 / (^o. T>) -f-

^ (T)

2 / (Ха, т) -f- 8 3

(3.8)

+ e

 

 

 

 

Учитывая (3.4), (3.7) и (3.8), приравняем в (3.6) коэф­ фициенты при одинаковых степенях е:

иКо = КоК,

 

(3.9)

UKlok] = КІк]Іо +

к Л к] +

Ф = 1 , 2

, . .. ).

 

3 0 2 п р и б л и ж е н н о е И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И Й [ГЛ. XI

Здесь DÖa’ ~ 11— столбцовая матрица, которая полностью определяется величинами, фигурирующими в равенствах (3.9) с номерами 1 , 2, ..., k — 1. Так, например,

 

 

 

= S T

-

ßo, T) TKa,

D L1] =

< 1]

+

-

A~]H jfloo ß a, T) TK[']+

 

dx

 

 

 

 

 

 

 

+

Rio ßo, T)

W T KO +

d ( T K a )

+

 

 

dx

H“ Rn ßa, T) TKo H---R2 0 ßa, T) TKo

Поскольку Ko — собственный вектор матрицы U, от­ вечающий собственному значению А,0, то первое равенство (3.9) выполняется тождественно. Остальные равенства (3.9) также выполняются тождественно, если принять (см. гл. VIII, § 6 )

=

п

(3.10)

L*] _

 

D ^ + K o Q [ohl

(3.11)

Ка —

s*o

Рекуррентные соотношения (3.10), (3.11) позволяют последовательно построить члены рядов (3.4), посредством которых представляется решение (3.3) системы (3.1), что и доказывает теорему.

Эта теорема остается справедливой и в случае любого целого положительного ц, только при этом несколько

иными будут' выражения столбцовых матриц D a- 1 1 (k =

=1, 2, ...).

3.2.Приближенное решение системы. Приближенное ре­

шение системы (3.1) можно получить на основе приведен­ ного выше формального решения путем удерживания в формальных рядах (3.4) некоторого конечного числа первых членов При этом решение системы (3.1) т-го приближения будет представлено соотношениями

хт {(, е) = 2 К{ат) (т, б) у{оП),

dyT

= l [om)ß ,e )y T \

di

 

* 4]

И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И Я ( С Л У Ч А Я Б)

3 0 3

где

 

 

 

 

 

 

 

 

 

 

КТ = К0+ У і tkK\k\

 

IT =

+

2 W

1.

 

 

 

ft=i

 

 

fc=i

 

 

Полагая

в =

1,

получим

приближенные

решения

для

исходной системы (1.4). Так, например,

 

 

 

* 1

( Т е)

S

\ К а ( \ 4 " Qad) 4 "

Р оD O

У а \

 

 

 

 

0 = I

 

 

 

 

 

 

 

 

*!°L = ß a - M

aDW)ylJ\

 

 

 

 

 

 

Р о = У

кж

 

 

 

 

 

 

 

--

К

Я-

 

 

 

 

 

 

 

S = 1

5

СТ

 

 

 

 

 

 

 

*+0

 

 

 

 

 

§ 4. Приближенное интегрирование уравнений управляемого процесса (случай Б)

Для построения приближенного решения системы (1.4) здесь мы используем систему (1.11) при р = 0. Имеем

А (т)- ^ г = ß (X) X + Н (т) и,

t

(4.1)

и = J G(t — t',x')vit\x')dt\

V= Т ( х ) х .

4.1.Построение формального решения. Введем в рассмот­ рение матрицу

U (К, т) == Л- 1 (т) В (т) -f А~'Н (т) Roo (Я, т) Т (т)

 

и определяющее уравнение

 

IU ()., т) — ХЕпI = 0

(4.2)

(Еп — единичная матрица порядка /?).

Каждый корень Яа (т) уравнения (4.2) является в то же время собственным значением матрицы ІІ{а) (т) == 1) (Яст, т), так что если р;а) (т) (j = 1 , 2 , ..., п) — собственные значе­ ния матрицы и (°\ то по крайней мере одна из этих скаляр

3 0 4 П Р И Б Л И Ж Е Н Н О Е И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И Й [ГЛ. ХГ

ных функций совпадает с А0 (т). Мы ограничимся рассмот­ рением простейшего случая, когда с функцией А0 (т) при любом т £ [0 , L ] совпадает одно и то же изолированное

собственное значение матрицы U{a) (т), например |х)ст>(т). Через Ко и М0 обозначим соответственно столбцовую и

строчную матрицы, определенные равенствами

и {а} (т) /Са (Т) =

/Со (Т) РІ° (т),

(4.3)

А40 (т) і/(а1 (т)=Л4а(т)р(,0) (т),

Ма(т) Кп(т) =

1 .

 

Будем считать, что в качестве Ко и Л40 приняты те ре­ шения уравнений (4.3), которые дифференцируемы столь­

ко же раз, сколько раз дифференцируема матрица Uw . Обозначим

5 (Я,/) = A~'(t) HRl0(X, t) Т (t).

Т е о р е м а

4.1.

Пусть

А0 корень

определяющего

уравнения (4.2) и при всех т £ [0, L]

 

 

 

1) р Г ( г ) = Аа (т ),

г іа) (т )

ф pSa) (т)

(/

= 2 , 3 ............ п );

2 ) MaSaKa Ф 1

(Sa = S(Xa,t)).

 

 

 

Тогда соответствующее этому корню формальное решение

системы (4.1) можно представить в виде

 

 

 

Ха (/, е) =

Ко (Ч е) уа,

~^L = Xa(x,г)уа,

(4.4)

где Ко и Ха соответственно столбцовая матрица и

ска­

лярная функция, имеющие формальные разложения

 

Ко (т, е) = Ка (т) +

2

^Как] (т),

 

 

 

 

/te1

_

 

оо

 

 

 

 

Аа (Ч е) = Ха (г) +

2

ЕкХ ^ (т).

(4.5)

 

 

 

 

fe=i

 

 

Д о к а з а т е л ь с т в о .

Подставим

в уравнения

(4.1)

значение вектора х, определенное равенствами (4.4). Полу­ чим

dK„ (т, в)

Л ( т ) е ----- ^ -------1- Ко (Ч е) Ао (X, е) = В (т)/<0 (т,е) + НІа,

(4.6)

S 4]

И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И Я (С Л У Ч А Й Б)

3 0 5

где

по-прежнему

 

 

 

 

I

 

 

 

ІО =

j G(t

т') Т(т') 7(а (т', е) exp [Ѳа(<!', g) — Ѳо^, e)]dt' ,

a Ѳст — функция,

удовлетворяющая

соотношению

0 =•■

= К (т, е)-

 

(3.7) и (3.8), приравня­

Имея в виду соотношения (4.5),

ем в равенстве (4.6) коэффициенты при одинаковых степе­

нях

е. Получим

 

Uw Ko =

 

 

 

 

 

 

 

 

К а К ,

 

(4.7)

 

U(0)Klok] =

Ф

к + а -

SaKo) 4 *

3+ D[k- 1]

 

 

 

 

( Ä = l , 2 ,

...).

 

(4.8)

Здесь D^a~

13 — столбцовая

матрица, известная

приизвест­

ных

Ко, К .......Kla~l\

A.ff—‘-1.

Так,

например,

 

 

= dK°

■A

lH R 10(Xa, т)

d(TKa)

 

 

 

 

dx

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

+

^

 

 

1 dkn

 

’s

 

 

 

7?n (A-cj,T) -f

R z o i ^ a ,

T ) T K a j ■

В силу (4.3) и условия 1) теоремы равенство (4.7) выпол­ няется тождественно. Покажем, что при соответствующем

выборе К1а] и X,„AJ равенства (4.8) также обращаются в тождества.

Предварительно проведем некоторые дополнительные построения.

Квадратная матрица Р а = КоМа является проекцион­

ной, соответствующей

собственному значению piö) = Я0

матрицы U(a). В силу

условия 1) теоремы проекционная

матрица, соответствующая всем остальным собственным зна­

чениям матрицы

и {а), равна

Р _ а = Еп Р аРанг квад­

ратной матрицы

равен п 1 , и потому она может быть

разложена на множители К-о

и А4_а = К-a М_а)

матрицы типа соответственно п X п — I и п — 1 х п, как

и матрица Р—а, дифференцируемые по т столько же раз. сколько раз дифференцируема U(a\

3 0 6 П Р И Б Л И Ж Е Н Н О Е

И Н Т Е Г Р И Р О В А Н И Е

У Р А В Н Е Н И Й [ГЛ.

XI

Матрицы К —а и Л1_

0 друг с другом

и с матрицами

Ко,

Ма связаны соотношениями

М- аК = Е п_ 1, М__оКа = М оК-о = 0.

Далее,

если К (а) ^ (КоК-о), /И№ =

(

), а

Л(а) =

" ( 0° Л_

) ’ Где Л_а ^

M-°U{a)K-o,

то

 

 

Д(а) = /С(а)А(а)УИ(а),

Л4(а К{а) =

/С<а)М(а) == £„

(4.9)

(см. гл. V). Заметим еще, что собственными значениями мат­

рицы Л _а служат собственные значения ц/0) (/ =

2, 3, ..., п)

матрицы и іа).

 

 

 

 

 

 

 

 

 

 

Умножим теперь Іг-е равенство (4.8)

слева на М (а), за­

менив в нем U(0) выражением (4.9). Получим

 

 

 

A{a)Qok] = Qlok]Xo +

Л4(0) (Ko-SoKo)

+

M[a)Dlo ~ '\ (4.10)

Здесь

 

 

 

 

 

 

 

 

лін

 

 

 

 

 

 

 

 

 

f М о К [ок]

 

 

 

 

 

 

Q[k]

Л4іа)К ок][

 

Vao

 

 

 

 

 

=

 

 

o [ft]

 

 

 

 

 

 

 

 

 

,M _0Klok]

 

 

 

 

 

 

 

 

 

 

 

Ч—ас

 

 

 

Так

как Л(<7) — квазидиагональная

матрица,

равенство

(4.10)

распадается на следующие два:

 

 

 

 

 

М о (Ко -

SoK o ) 4 * ] +

M 0D lok- ' ] = 0 ,

 

 

 

 

I

A_0Q1*I„ = QL^

+ м_а (Ко -

SoKo) № + м_ао[а"-1].

)

 

 

 

 

 

 

 

 

 

 

 

(4.11)

В силу условия 2) теоремы первое равенство (4.11) раз­

решимо относительно A/а] и для любого с[ко получаем

 

 

 

 

 

 

W

=

-

 

 

 

 

(4.12)

 

 

 

 

 

 

і - Ч А Л т

 

 

 

 

 

Матрица Л_а не имеет собственных значений, равных Хст.

Значит,

Л_ 0 \оЕп—\ — невырожденная

матрица и

из

второго

равенства

(4.11)

можно

определить Q ^ a:

 

 

 

Q-aa =

(КЕп_

 

•Л..

M_o^o^o^oDlo~ n

(4.13)

 

 

1~

MaSaKa

 

 

 

 

 

 

 

 

 

 

 

Последняя формула представляет субматрицу

 

мат­

рицы

 

 

 

 

 

 

 

 

С*1

 

 

 

Qo*1- В качестве другой субматрицы QLoo этой матрицы

§ 4]

И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И Й (СЛУЧАЙ Б)

3 0 7

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

Зная фаЧ легко определить и искомую столбцовую мат­ рицу /Со]:

к т = K mQm = KaQw +

(4 .1 4)

Рекуррентные соотношения (4.12) — (4.14) позволяют по­ следовательно определить члены рядов (4.5), посредством которых представляется частное решение (4.4) системы (4.1). Теорема доказана.

Эта теорема, так же как и теорема 3.1, легко обобщается на случай, когда матрицы А и В — функции от т и е, до­ пускающие на [0 , L ] разложения (сходящиеся или по край­ ней мере асимптотические) по степеням е.

4.2. Приближенное решение системы. Используя постро­ енные формальные решения, решение системы (4.1) т-ѵо приближения можно представить соотношениями

х(ат)(t, е) =

Кат) (т, е) уТ \

 

- ^ - =

С ) (т,е)і/Г ,

где

 

 

т

 

т

к (пт) = Ко + 2

ekK[k),

?4т) = К + 2

k=\

 

fc=l

которые при е =

1 служат приближенным решением и для

исходной системы (1.4).

Для примера приведем простейшие приближенные ре шения системы (1.4).

При т = 0

 

 

 

 

*?' =

 

= K f i .

 

При т = 1

 

 

 

 

=

 

 

 

 

= { * а (1 + Q & ) + К —а ( № - і - А - с Г ’

^

^

О)

 

М°°Ь0]

b l )

 

dyа

К —

 

dt

\ -M as aKn

rjQ ■

 

Глава XII

НЕКОТОРЫЕ КАНОНИЧЕСКИЕ ФОРМЫ УРАВНЕНИЙ ЛИНЕЙНЫХ ПРОЦЕССОВ

Математическая модель многих процессов, происходя­ щих в реальной действительности, представляется дифферен­ циальной системой, которая в векторно-матричной записи имеет вид

 

 

М 0

- ^

+ М 0 - щ- +

М*)<7 = ф.

(0 -1 )

где

<7 — столбцовая матрица параметров процесса glt

q2, ...

...,

qn (например,

обобщенных координат механической

системы);

L0,

La,

L2 — некоторые

квадратные матрицы

порядка п

(матрицы динамических

коэффициентов

систе­

мы); ер — столбцовая матрица, элементы которой являются, вообще говоря, функциями от t и, быть может, управляю­ щих функций, которые в свою очередь определяются зна­ чениями <7 j, q2, ..., qn. Уравнениями такого типа описывают­ ся, в частности, малые колебания механических систем, поведение линейных объектов управления в системах авто­ матического управления и т. п.

Анализ и синтез процессов, описываемых системой диф­ ференциальных уравнений второго порядка с переменными коэффициентами, особенно если дифференциальная систе­ ма имеет высокий порядок, связаны с преодолением немалых трудностей. Эти затруднения в значительной мере могут быть сняты, если предварительно произвести «диагонализацию» исходной системы, т. е. соответствующей заменой пе­ ременных преобразовать эту систему к системе, матрицы коэффициентов которой имеют диагональную или по край­ ней мере квазидиагональную форму. Настоящая глава по­ священа изложению некоторых алгоритмов таких канони­ ческих преобразований.

і]

П Р Е О Б Р А З О В А Н И Е С И С Т Е М Ы У Р А В Н Е Н И Й

3 0 9

 

 

§ 1. Преобразование системы уравнений с постоянными коэффициентами к расщепленному виду

Рассмотрим простейшую систему уравнений второго по­ рядка с постоянными коэффициентами

U

_сҢ_

+ L2q = cp (t)

(det L0 Ф 0 ).

dt2

Заметим, что система более общего вида

Lo-^г - + Li -% + LiЯ = Ф if)

( 1 . 1 )

( 1.2)

может быть приведена кгформе (1 .1 ), однако при этом коэф­ фициенты преобразованной системы уже будут зависеть от t. Так, замена переменных

q = Vz,

где квадратная матрица V определена соотношением

V = е~~2 L0 1l'1

преобразует систему (1 .2

) к виду

 

 

 

 

 

 

 

Уг = У“ Ѵ

ф

1.1. Однородная

система.

Если

собственные

значения

матрицы и =

LÖ"’L2

разбиты

на р

групп v(jCT)....... vjfj (а =

р

 

п) при условии

 

 

= 1 .......р\ 2

^ 0=

 

 

а= 1

 

 

 

 

 

 

 

І ѵ Г - ѵ У '

■о

(S Ф or),

(1.3)

то могут быть построены (см. гл. V) блочные матрицы

 

 

 

'ос,

сс,

 

 

 

 

 

 

о

 

к = (KLH.2. . . Яр),

а =

о

 

ц =

 

а „

Соседние файлы в папке книги из ГПНТБ