Файл: Абгарян К.А. Матричные и асимптотические методы в теории линейных систем.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 167
Скачиваний: 0
154 |
М А Т Р И Ц Ы И Л И Н Е Й Н Ы Е У Р А В Н Е Н И Я |
[ГЛ . VI |
||
при условии |
||||
Y (0 = |
Е. |
(9.8) |
||
|
||||
Тогда (9.7) принимает вид |
|
|
||
|
а- (t) = Y* (/„) X(у + |
[ Y* (s) h (s) ds. |
(9.9) |
|
|
|
}„ |
|
|
Как было показано в § 8, если X — фундаментальная матрица однородной системы (9.4), то Y* (s) X (s) = С. Из условия (9.8) в данном случае получаем
C = X(t),
так что
K*(s)X(s) = X(f).
Отсюда
К* (s) = X (t) X - 1(s), Y* (У = X (О X -1(У .
Учитывая это, вместо (9.9) будем иметь
x(t) = X (О Х -' (у а(У + j X (t) X - 1(s) А (s) ds.
§ 10. Решение одного матричного уравнения
Т е о р е м а 10.1. Решение матричного уравнения
|
4 г = |
A(t)X + XB(t) |
X (У |
= С |
(10.1) |
|
представляется в виде |
|
|
|
|
||
|
|
X = Y(t)CZ(t), |
|
|
(10.2) |
|
где Y (t) |
и Z (/) — соответственно |
решения |
матричных |
|||
уравнений |
|
|
|
|
|
|
|
- § - = A(t)Y, |
Y (у |
= £, |
|
(10.3) |
|
|
~ -- = ZB (t), |
Z ( y = £ . |
|
(10.4) |
||
Д о к а з а т е л ь с т в о . |
Дифференцируя (10.2) по /, с |
|||||
учетом (10.3) и (10.4) получим |
|
|
|
|||
^ |
= — |
CZ + к с - ^ - = л к с г + |
YCZB. |
|||
$ 10] |
РЕШЕНИЕ ОДНОГО МАТРИЧНОГО УРАВНЕНИЯ |
155 |
Очевидно, к такому же виду приводится и правая часть уравнения (10.1), если вместо X подставить (10.2). Теорема доказана.
С л е д с т в и е 1. Решение матричного уравнения
= A{t)X — ХА (t) |
X (О = С |
(10.5) |
представляется в виде |
|
|
X = Y(t)CY~'(t). |
(10.6) |
|
Действительно, в данном случае уравнения (10.3) и (10.4) |
||
приобретают вид |
|
|
W - = A ( t)Y , |
Y (tQ) = Е, |
(10.7) |
- § - = - Z A ((), |
Z (t0) = E. |
(10.8) |
Как видим, (10.8) есть уравнение, сопряженное уравнению (10.7), и поэтому
Z(t) = Y~'(t).
Отсюда в силу (10.2) следует выражение (10.6).
С л е д с т в и е 2. |
Если |
А и В |
в уравнении (10.1) - |
постоянные матрицы, |
то решение |
этого уравнения пред |
|
ставляется в виде |
|
Севѵ-‘°К |
|
X = еА{і~Іо) |
|||
Действительно, в этом случае решениями уравнений (10.3) и (10.4) являются соответственно матрицы ехр [А (t — *„)] и ехр [В (t — t0)\. .
Г л а в а VII
СИСТЕМЫ ЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ПЕРВОГО ПОРЯДКА
СПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ
Вданной главе рассматриваются линейные системы, представленные векторно-матричным уравнением
— Ux + h (t)> |
(0.1) |
где U — постоянная матрица.
§ 1. Экспоненциал матрицы
Как известно, экспоненциал скалярной величины а представляется рядом ea = 1 + а + d2 + ■• • + ap— Н
По аналогии с этим вводится понятие экспоненциала квадратной матрицы А. Под экспоненциалом матрицы А понимается матричная функция
( 1. 1)
Ряд (1.1) сходится для любой квадратной матрицы, так как сходится скалярный ряд
со
составленный для нормы этой матрицы.
Из сходимости ряда (1.1) для любой квадратной матрицы следует сходимость и ряда
оо
(1.2)
S 1] |
Э К С П О Н Е Н Ц И А Л М А Т Р И Ц Ы |
157 |
где t — скалярный множитель. Этот ряд представляет собой экспоненциал произведения At, т. е. еА‘.
Сходимость ряда (1.2) в любой конечной области комп лексной плоскости параметра t равномерная в силу равно
мерной сходимости в этой области |
ряда у |
|Л||р т |
р Отме- |
||||||
|
|
|
|
|
|
|
р=0 |
pl |
|
тим некоторые свойства экспоненциала. |
|
|
|||||||
1) |
e^u+s) — eAleAs = eAseAt |
(t, |
s £ К)- |
|
|
||||
Действительно, |
|
|
|
|
|
|
|
||
|
Aktk |
^ |
A4 |
|
|
p |
kcp—k |
|
|
A tp A s |
|
|
\ 4 |
rs> |
|
|
|||
eHle |
k\ |
s - |
|
AP X |
k\ (p — &)l |
|
|||
|
|
/=0 |
|
p^O |
|
k=0 |
|
|
|
|
k=o |
|
|
= у= |
у |
Ap (tpl+ sf _ |
g.4(.'-)-s)_ |
||
P—0
Так как eA(*+s) = eA (s+f>, то из приведенной цепочки равенств следует также коммутативность матриц ем и eAs.
Полагая s = —t, будем иметь eAte~At = еА0 = Е. Значит,
2)eAt — всегда невырожденная матрица и ее обратная матрица равна е~АІ.
3)Если А В = ВА, тоеА+в — еАев = евеА. Покажем это.
Имеем
|
|
А р В |
_ |
АР |
ч” |
В* |
|
Рл ч“ |
|
АрВР |
||
|
|
|
У - — = У У |
|
|
|||||||
|
|
е*е |
|
|
|
|
||||||
|
|
|
|
рI |
ZJ |
а\ |
|
ZJ |
ZJ |
Pl?! |
||
|
|
|
|
Li |
qI |
|
|
|
|
|
||
|
|
|
|
р=0 |
р=0 |
|
|
P — |
0 |
q = |
0 |
|
Положим p + |
q = s (s = |
0, |
1,2, |
...). Тогда, учитывая, |
||||||||
что p |
s, |
имеем |
|
|
|
|
|
|
|
|
|
|
|
|
CO |
S |
„ n ^ < j _ n |
_ |
o o |
, |
S |
|
|
|
. |
e |
л Q_ v-i |
'j |
APBS p |
|
1 |
v-i |
|
|
s- ■А рВ S-P . |
|||
~ |
|
p! (s — p)I ~ |
s=0 |
si |
p=0 |
pl (s — p)l |
||||||
|
|
s=0 p=0 |
|
|
|
|
|
|
||||
С другой |
стороны, |
|
|
|
|
|
|
|
|
|||
|
|
|
|
,.4+B |
S |
i ( |
^ |
+ |
ß )s. |
|
|
|
|
|
|
|
= |
|
|
||||||
s=0
и так как А и В перестановочны, то Р—ор\ (s- р )I ■
158
И,
С И С Т Е М Ы Л И Н Е Й Н Ы Х У Р А В Н Е Н И Й |
[ГЛ. VII |
значит,
гл+’ = £ |
т £ л ^ ‘ф' в ‘- '’= е', л |
s=0 |
Р = 0 |
4) Производная экспоненциала.
Ряд, полученный формальным дифференцированием ря да (1.2), также сходится равномерно, поэтому законно по членное дифференцирование ряда (1.2). Учитывая это, полу чаем
deA ' |
= А + |
A 4 |
A 4 |
= AeAt = eAtA. |
di |
+ |
~ZT + |
§2. Решение дифференциальной системы
вформе экспоненциала
В силу свойства 4) экспоненциала матрицы матрица еи и—/о) представляет собою решение матричного уравнения
^ Г = ѴХ, X (t0) = Е
и, значит, является фундаментальной матрицей системы
В соответствии с этим общее решение неоднородной си стемы (0.1) при условии X (ta) = с можно представить в виде
I
|
X(t) = eU{‘-'°)c -f- I" eUi‘- l4e-U(s~t°)h(s) ds, |
||
или |
tо |
|
|
t |
|
|
|
|
|
|
|
|
X(t) = euv~tAc -j- ^ eUll~s)h (s) ds. |
|
|
§ 3. |
Метод Эйлера |
|
|
Решение векторно-матричного уравнения |
|||
|
-%г = их, |
X( g = с |
(3.1) |
будем искать в виде |
|
|
|
|
X = |
Кеи, |
(3.2) |
где |
К — постоянный вектор |
(столбцовая |
матрица). |
5 3] |
М Е Т О Д Э Й Л Е Р А |
159 |
|
Подстановка (3.2) в (3.1) приводит к алгебраическому |
|
равенству |
|
|
|
и К = хк. |
|
|
Отсюда видно, что (3.2) представляет собой частное ре |
|
шение уравнения (3.1), если Xесть собственное значение мат |
||
рицы Ö, а К — собственный вектор этой матрицы, отвечаю щий собственному значению X.
Рассмотрим различные случаи, которые могут иметь место.
1. Матрица U имеет л различных собственных значений A,lt Х2, ..., Хп (п — порядок матрицы U). Каждому простому собственному значению Xf отвечает единственный (с точ ностью до произвольного множителя) собственный вектор Kj, причем собственные векторы, отвечающие различным собственным значениям, линейно независимы. Общее реше ние уравнения (3.1) может быть представлено в виде
* = s |
(3-3) |
/= I |
|
где уі — произвольные постоянные.
Чтобы убедиться в этом, достаточно показать, что соот ветствующим выбором постоянных у/ можно удовлетворить любому начальному условию, ибо, очевидно, (3.3) обращает уравнение (3.1) в тождество. Прежде всего представим (3.3)
вином виде. С этой целью введем матрицы
К= (К Д 2 . .. Кп),
Тогда вместо (3.3) будем иметь
X = Кел<у.
Отсюда при t = t0 получим
x(t0) = KeAI°y |
(3.4) |
Так как К и еЛ1° — невырожденные матрицы, то уравнение