книги из ГПНТБ / Абгарян К.А. Матричные и асимптотические методы в теории линейных систем
.pdfЗОО |
П Р И Б Л И Ж Е Н Н О Е И Н Т Е Г Р И Р О В А Н И Е |
У Р А В Н Е Н И Й |
[ГЛ |
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 |
, . .. ). |
|
* 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. Уравнениями такого типа описывают ся, в частности, малые колебания механических систем, поведение линейных объектов управления в системах авто матического управления и т. п.
Анализ и синтез процессов, описываемых системой диф ференциальных уравнений второго порядка с переменными коэффициентами, особенно если дифференциальная систе ма имеет высокий порядок, связаны с преодолением немалых трудностей. Эти затруднения в значительной мере могут быть сняты, если предварительно произвести «диагонализацию» исходной системы, т. е. соответствующей заменой пе ременных преобразовать эту систему к системе, матрицы коэффициентов которой имеют диагональную или по край ней мере квазидиагональную форму. Настоящая глава по священа изложению некоторых алгоритмов таких канони ческих преобразований.
