2 9 2 П Р И Б Л И Ж Е Н Н О Е И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И И [ГЛ. XI
с невырожденной и дифференцируемой на [^0, Г] матрицей К приводит систему (1.4) к векторно-матричному уравнению
с непрерывной на U0, Т 1 матрицей U тогда и только тогда, когда
где X ( 0 — единственное решение уравнения (1.5), С — постоянная невырожденная матрица порядка п, а Z (t) — непрерывно дифференцируемая и невырожденная на [/0, Т ] матрица порядка п.
Д о к а з а т е л ь с т в о . Замена переменных (1.6 ) при водит систему (1.4) к матричному уравнению
|
|
I |
|
|
ХС ^ Ш |
= A~lH |
J G(t — t',t')T(t')X{t')C [Z{t')y(t') — |
|
|
|
|
Z(t)y(t)\dt', |
которое |
допускает |
решение |
|
|
Отсюда |
|
Z (і) у (і) = |
const. |
|
dy_ |
|
|
|
|
z~' |
(1.9) |
|
|
di |
|
|
|
|
В силу свойств матрицы Z матрица
преобразованного уравнения (1.9) непрерывна на [/0, Т]. Пусть, далее, К (і) — матрица преобразования, которое систему (1.4) приводит к уравнению (1.7). Покажем, что тог
да К ( 0 представима в форме (1 .8 ).
Матрица этого преобразования удовлетворяет уравнению
(*«- + К и - А - 'В /< )у =
t
= А~'Н { G(t — 1', t') Т (О К (t’) у (t') di'.
Имеем
5 1] |
И Н Т Е Г Р О - Д И Ф Ф Е Р Е Н Ц И А Л Ь Н А Я С И СТЕ М А |
293 |
где Y — фундаментальная матрица системы (1.7), а с — столбцовая матрица произвольных постоянных. Учитывая это, получаем
-=£- = |
A~lBK - KU + А~'НІУ~\ |
( 1. 10) |
где |
|
|
/ ( 0 - j |
G (t — t', t') T (t') К (t') Y (?) dt'. |
|
Принимая во внимание (1.5) и (1.10), а также соотноше ния
d ( X ~ ] K Y ) |
d X ~ l |
KY + |
X ~l |
dt |
Y 4 |
- X~]K |
dY |
|
dt |
dt |
“' ' |
1 |
|
' |
1 ' ‘ |
' ' |
dt |
|
= — X ~ l (A~xBX + |
A~lHI) X ~ xKY + |
X ~lA~lBKY - |
|
Hr/riu |
, v-1 Л - 1 и n/-h |
|
|
— = 0. |
|
— X KLJY |
+ |
X ~lA~lH IY~lY + X~lK |
Поэтому |
|
X ~xKY — C = const. |
|
|
|
|
|
|
|
|
Отсюда, полагая Y |
= |
Z~l, получаем |
|
|
|
|
Теорема доказана. |
|
|
К = |
XCZ. |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 .2 . О методике построения приближенного решения |
уравнений. |
Интегро-дифференциальная |
система |
(1 .4 ) со |
держится в следующем семействе систем уравнений более общего вида:
dx |
В |
|
(т) “■ |
|
А № ~дГ = |
|
|
ц = |
j |
G(t — t',i')v { t\T ')d t\ . |
(іл 1 ) |
V = |
Т (т) X |
|
|
(т = et, е > |
0 |
, р > |
0 ). |
|
Ясно, что при е = |
1 |
(1.11) совпадает с (1.4). В силу это |
го всякое решение |
х (і, е) |
системы (1 .1 1) при |
значении |
2 9 4 |
П Р И Б Л И Ж Е Н Н О Е |
И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И Й [ГЛ. XI |
параметра е, равном единице, |
будет являться решением си |
стемы (1.4). Учитывая |
это, для |
построения приближенного |
решения нестационарной интегро-дифференциальной систе мы (1.4) поступим так. Сначала для системы (1.11) постро им формальное решение в виде бесконечного ряда по степе ням е. Частичные суммы этих рядов будем трактовать как
приближенные решения системы (1 .1 1), а при е = 1 |
— как |
приближенные решения исходной системы (1.4). |
|
Такой путь построения приближенных решений системы |
(1.4) |
является эффективным и плодотворным тогда, |
когда |
А (t), |
В (t), Н (/), Т (t), а также G (t — t', t') как |
функ |
ция от второго аргумента являются медленно меняющими ся функциями.
В дальнейшем будем различать два случая:
А) воздействие регулятора на регулируемый процесс мало, так что решения уравнений замкнутой системы близ ки к решениям уравнений при и = 0, и Б) воздействие ре гулятора на регулируемый процесс нельзя считать малым.
Приближенное решение системы (1.4) будем строить на основе формального решения системы (1 .1 1) при значении р = 1 в случае А, и ц = 0 в случае Б.
§ 2. Приведение уравнений управляемого процесса к расщепленной дифференциальной системе (метод последовательных приближений)
При довольно общих предположениях решение интегродифференциальной системы (1 .1 1) можно свести к интегри рованию некоторого числа независимых друг от друга под систем линейных дифференциальных уравнений первого по рядка. Мы здесь ограничимся рассмотрением случая р = 1.
Итак, имеем
А (т) — =, В (т) X -f ЕН (т) и,
и = і G(t — t’,%')v{t'%')dt', ' |
(2Л) |
0 = 7 ’(т) X.
Пусть собственные значения матрицы U (т) = А- 1 (т) В{х) на сегменте [О, L] разделяются на р непересекающихся
§ 2] |
М Е Т О Д П О С Л Е Д О В А Т Е Л Ь Н Ы Х П Р И Б Л И Ж Е Н И Я |
295 |
групп. Предполагая, что коэффициенты уравнений в систе ме (2 .1 ) имеют на [О, L I производные пот всех порядков, решение этой системы будем искать в виде
л - |
2 Ко (т, е) уа (0 , |
|
|
(2 .2 ) |
|
|
а= 1 |
|
|
|
|
|
d y n |
= |
~ |
~ |
(о - |
1,2, |
.. . , р), |
(2.3) |
|
Да (Т, |
е) уст |
где |
оо |
|
|
|
|
со |
|
__ |
|
|
|
|
|
Ка (Т, е) = |
2 |
e*/C„[ 4 |
(т), |
А а(т, 8 |
) = |
2 е*А^] (т). |
(2.4) |
|
φχ= 0 |
|
|
|
|
φτ= 0 |
|
Всвою очередь решения уравнений (2.3) будем строить
вформе ряда
~Уо= |
2 |
(2-5) |
|
*=о |
|
Подставим значения А0 |
и уа из |
(2.4) и (2.5) в (2.3) и |
приравняем в полученном соотношении коэффициенты при одинаковых степенях е. В результате придем к следующей
системе |
уравнений: |
|
|
|
|
|
|
|
|
|
|
^ |
0] |
А |
» |
, |
|
|
|
|
|
|
|
|
|
|
d t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ft—] |
|
|
|
|
|
|
|
|
ПН |
_ л 1 0 1 / |
- 1 |
|
л ^ у п |
|
|
|
|
d,Ja ' |
' ' л' 1. У |
|
|
(k= |
1 , 2 , ...). |
--------- ^о £-'а |
і- |
|
|
|
Ус |
|
di |
|
|
|
|
|
|
|
|
|
|
|
|
|
Пусть |
Ка1 |
— фундаментальная |
матрица |
решений урав |
нения |
|
|
|
|
dy™ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
А |
» |
, |
|
|
гак что |
|
|
|
|
dt |
|
|
|
|
|
|
|
ЛО] |
|
у [0 ] |
|
|
|
|
|
|
|
|
|
|
— |
Lc |
|
|
(са — матрица-столбец |
У с |
|
* о |
|
|
|
произвольных |
постоянных). |
Тогда |
частное |
решение |
уравнения |
|
|
|
|
|
4 І Ч |
- д н у * ] |
|
fe-i |
|
|
|
|
|
|
+ |
V л |
Уа |
|
|
|
|
dt |
|
— 2Ѵа |
Уа |
і= 0 |
Jio |
|
можно |
представить |
так: |
|
|
|
|
|
|
l/ak] = ТІ°] (/) |
f |
L O J - |
1 |
( П |
У , |
|
A |
ak[ - |
‘ ] ( V ) |
y \ i \ t ' ) d t ' . |
296 П Р И Б Л И Ж Е Н Н О Е И Н Т Е Г Р И Р О В А Н И Е У Р А В Н Е Н И Й [ГЛ. X I
Обозначив |
|
|
|
|
|
|
|
|
|
|
|
|
Y™ (t) = |
К[а0] (t) [ И 01” ' (t') 5 |
' Ла*~‘] (т') У™ (t') dt' , |
(2. 6) |
|
|
|
|
|
|
У |
|
|
і= |
0 |
|
|
|
|
будем |
иметь |
|
* (I |
|
|
|
|
|
|
|
|
і/а] = |
У[а 1 |
( 0 ( У[а 1 |
’ (П № |
(т') У?1(П dt' = |
|
|
|
|
|
|
К |
|
|
|
|
|
|
|
|
|
|
|
= |
П 0 1 |
(О J ГУ]~' (О ЛУ] (т') Ft°] (t') dt’Ca = |
( 0 |
Са, |
у? |
1= |
Y™ (t) 5 |
ѴІ0] |
' (t') [ЛУ] (т') № (t’) + |
|
|
|
|
|
|
to |
|
|
|
|
|
|
|
|
|
+ ЛУ] (т') yW (t')] dt' = |
Y\?] (t) J n |
o]- ' (t') [A ^ (т') Y™ (t') |
+ |
|
|
|
|
|
|
|
|
|
|
to |
|
|
|
|
и, вообще, |
|
|
|
|
+ а У1 (T')Hö[0 1 (t')] dt’Co = кУ] ( 0 |
c0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
filo ] ^ Y [ k]ca |
|
( Ä = l , 2 , |
...)• |
(2-7) |
|
Равенства (2.6) и (2.7) определяют iß ] через Лса0]....... A ak]l . |
|
Перейдем |
к |
построению |
/(У1, |
ЛУ1 (k = 0, 1, 2, |
...). |
|
|
Подставим (2.2) и (2.3) в уравнения (2.1) и приравняем |
нулю сумму всех слаі аемых, содержащих уа. Получим |
|
[E^ |
- |
+ KaÂa - U K 0) ~уо- |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
_ |
^ |
|
|
|
|
- |
гА~'Н |
\ G(t — 1', т') Т (т') Ко (т', е) уа (t') dt' = |
0. |
|
В |
последнее |
равенство |
подставим |
разложения |
(2.4), |
(2 .5 ) и приравняем коэффициенты при одинаковых |
степе |
нях е. Будем иметь |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L ol 4 |
0] |
= 0 , |
|
|
|
а |
ч |
о] + |
|
г[*-а]..[а] I |
/I'[*-1] |
_ |
(/е= 1 , 2 , ...). |
|
а=І ь 0 |
Уо |
Т' |
‘а |
|
|
|
|
|
