Файл: Абгарян К.А. Матричные и асимптотические методы в теории линейных систем.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 11.04.2024
Просмотров: 168
Скачиваний: 0
1 6 0 |
С И С Т Е М Ы Л И Н Е Й Н Ы Х У Р А В Н Е Н И Й |
[ГЛ . VII |
(3.4) |
разрешимо относительно у: |
|
|
у = е~ЛіаК~]х(і0) = е~АІ°К~'с. |
|
Таким образом, действительно, если все собственные зна чения матрицы U простые, то общее решение уравнения (3.1)
представляется |
выражением (3.3). |
|
2. Число |
различных собственных значений матрицы U |
|
меньше, чем |
п, |
но каждому собственному значению крат |
ности гj отвечает ровно rf линейно независимых собственных векторов K f\ В этом случае общее решение одно родной системы можно представить в виде
/=1
- 3. Число различных собственных значений матрицы U равно или меньше, чем п, и каждому собственному значению К/ кратности Г/ отвечает один или несколько (но не более, чем г() собственных векторов (общий случай). Этот случай будет рассмотрен позже при описании другого метода.
§ 4. Преобразование Лапласа
Построим решение векторно-матричного уравнения
|
-§ - = Ux + h(t) |
|
(4.1) |
при начальном условии |
|
(4.2) |
|
|
|
|
|
с помощью преобразования Лапласа. |
|
е~Р‘ и |
|
Обе части уравнения (4.1) умножим справа на |
|||
проинтегрируем по t |
в пределах от 0 до оо. Получим |
|
|
СО |
оо |
оо |
|
или |
= U L ( x ) + L ( К ) |
-f- .г (0), |
(4.3) |
p L ( X ) |
|||
где |
|
|
|
оо |
|
оо |
|
о |
о |
§ 4] П Р Е О Б Р А З О В А Н И Е Л А П Л А СА 161
— изображения Лапласа функций х и /г. Из (4.3) находим
L(x) = ( p E - U ) - ' x(0) + ( p E - U ) - 'L(h). |
(4.4) |
||
Л е м м а 4.1. Для произвольной квадратной матрицы U |
|||
имеет место равенство |
|
|
|
L - 1l(p E -U )-'] = eU!. |
|
|
|
Д о к а з а т е л ь с т в о . |
Пусть К = |
(КхК2 ■■■Кт) — |
|
матрица, преобразующая U к форме Жордана |
|
|
|
J = diag (Ji (Xj), |
J2(^2 ) 1 • •• I |
I |
|
так что
m
и = KJM = % KoJc(K)Ma.
0=1
Тогда
(pE — U)-' = (pKM — KJMГ ' = [K (pE — J) M]~l =
Ho |
|
|
|
|
|
|
= K(pE — J)~lM. |
|
|
|
|
|
|
|
|
|
|
(p£ —/) = diagKpf*, —Л)> . . . . |
(pEkm— J m)]. |
|||||||
Учитывая это, |
получим |
m |
|
|
|
|||
(pE - |
1 |
|
|
1 |
||||
U)-1= |
2 Ko (p£fc0 - |
|
Mo- |
|||||
Отсюда |
|
|
|
|
0=1 |
|
|
|
|
|
|
|
m |
|
|
|
|
L~' [(pE - |
|
1 /)'1] |
= |
[(p£*0- |
Ja) 1Л4- (4.5) |
|||
|
2 |
|||||||
|
|
|
|
|
0=1 |
|
|
|
Имеем (см. |
(5.9.9)) |
|
|
|
|
|
||
{pEk0— Ja) |
— \{р |
ha) Eka |
Hka]—1 |
— |
|
|||
|
|
|
|
|
|
|
|
,fca-l |
E>o |
+, |
H *c |
+ |
|
+ |
**a |
||
( P - K ) s |
+ |
|||||||
p — \ |
|
( P - |
V |
2 " г |
|
ІР -К ) |
||
|
|
|
|
|
|
|
|
ß K. A. Айгаря«
162 |
|
С И С Т Е М Ы Л И Н Е Й Н Ы Х |
У Р А В Н Е Н И Й |
[ГЛ. VII |
|||||
На |
основании последнего |
соотношения |
|
|
|||||
|
|
|
|
|
|
к, |
|
|
|
|
L-'l(pEka-Jar']=%Hl-lL-1 |
|
(4.6) |
||||||
|
|
|
|
|
Ѵ=1 |
(P- К ? |
|
||
Здесь принято НІа = |
Ека- Но |
|
|
|
|||||
. —1 |
|
I |
с-к°° |
|
|
. |
|
|
|
|
7* |
|
|
е?1 |
|
|
|
||
( Р - К ? |
^ |
S і7 = к Г “р ' |
|
|
|||||
|
|
|
с — І оо |
|
|
|
|
|
|
|
|
|
|
|
|
(V— 1)1 |
(і = |
Ѵ - 1), |
|
|
|
Н. і |
, |
H i Р |
|
н Г Ѵ * - ' |
|||
|
|
|
«а |
|
Аа |
+ |
«а |
|
|
|
г\п — |
|
|
|
21 |
( К - |
Di“ |
||
|
г > |
ka). Поэтому |
|
|
|||||
(так как Нр0 — 0 при |
|
|
|||||||
L“ 1[(pEka— Jo)~]] = еНка*е%а*= |
e(X°£*«+"*o> 1= eJo |
||||||||
Подставляя |
(4.7) в |
(4.5), |
|
будем |
иметь |
|
(4.7) |
||
|
|
|
|||||||
L“ 1[(рЕ — U)-1] = 2 |
|
К У *(Ха) 'Л4а = |
KeJIM = еш. |
||||||
Лемма |
доказана. |
а=>1 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|||
Используя эту лемму, вместо (4.4) записываем |
|
||||||||
|
|
L{x) = L (еш) x{0) + |
L [еш) L (h). |
|
Применяя обратное преобразование Лапласа, отсюда на ходим
x(t) = eutx(0) + L-1 [L(eUi)L(h)]. |
(4.8) |
Для скалярных функций g и /г имеет место теорема о свертке
ZT1[L (£) L (А)] = J g (t - s) h (s) ds. |
(4.9) |
о |
|
Эта формула, очевидно, остается справедливой, когда одна из функций, например h, — векторная функция. Можно показать, что формула, аналогичная (4.9), имеет место и
5 4] П Р Е О Б Р А З О В А Н И Е Л А П Л А С А 163
в |
нашем |
случае, |
когда |
одна функция (еи>) — квадратная |
||||||||
матрица, |
а |
вторая |
(h) — вектор. |
|
|
|||||||
|
Действительно, так |
как (см. (4.7) и (4.6)) |
|
|||||||||
L (eut) = ( р Е - U)-1= |
т |
Ко (pEka- |
JaГ ! Мо = |
|||||||||
£ |
||||||||||||
|
|
|
|
|
|
|
а=1 |
|
|
|
|
|
|
|
|
|
|
т |
|
|
|
Ң У ~ 1 |
|
|
|
|
|
|
|
|
СТ=аІ |
|
Ѵ = І |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U |
/ |
/V— 1 |
|
|
|
|
|
|
|
= S |
|
E ^ ' М т Г Г Т ) , |
|||
ТО |
|
|
|
|
|
|
|
0 = 1 |
|
V = 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
г —1, г , т |
|
|
|
|
|
|
|
|
|
|
|
|
L |
lL(eut)L(h)] = |
|
|
|
|
|
|
|
|
|||
|
т |
к , 2 |
|
|
|
|
|
|
,ѵ—1 |
|
||
|
= 2 |
ЯС'ЛЫ Г1 [ і ( 17Д ПГ e’- 'j L (Л) |
||||||||||
|
о=1 |
ѵ=І |
|
|
|
|
t |
|
|
|
|
|
|
|
|
К. |
|
jy—1л/f |
P |
- 5)V — |
I |
|
|||
|
|
|
|
|
Г |
й * - |
||||||
|
= £ К. 2 Л ~ 'м . |
J |
|
|
Л |
|||||||
|
0 = 1 |
V = 1 |
|
|
|
|
t |
|
|
|
|
|
|
m |
1k° |
|
|
|
|
eXa{t~s)Mah (s)ds = |
|||||
|
-E MJ тS |
|
|
(V - |
1)1 |
|||||||
|
|
|
|
* |
|
|||||||
|
|
I |
У |
|
1 |
|
|
|
|
|
|
|
|
0=1 |
о Ѵ=1 |
|
|
|
|
|
|
|
|||
|
= f £ |
„Я*о |
-s>V «- |
|
|
|
||||||
|
Кое |
at<" S>^ ' u- s)Mo/i(S)ds = |
|
|||||||||
|
О |
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
С £ |
|
KaeJ<7{ka) u~s)Mah (s) ds = I ey ('~ s) h (s) ds. |
|||||||
|
|
|
6 CT=i |
|
|
|
|
|
|
о |
|
|
|
С учетом последнего результата соотношение (4.8) при |
|||||||||||
мет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X(0 = |
етх (0) + |
{ еи {t~s)h (s) ds. |
(4.10) |
||||||
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
Остается х (0) выбрать так, чтобы удовлетворялось гра |
|||||||||||
ничное условие (4.2). Имеем |
|
|
|
|||||||||
|
|
* (д |
= |
еи (‘°}х (0) + |
J° еи (t°~s)h (s) ds. |
I *
164 |
С И С Т Е М Ы Л И Н Е Й Н Ы Х У Р А В Н Е Н И Й |
[ГЛ VII |
Отсюда
Іь
X(0) = е~и‘°х( д — J e~Ush (s) ds.
о
Подставляя значение х (0) в (4.10), получаем
x(t) = еи {t~lo)x ( g — |
е {t~s)h (s) ds + |
[ eu {‘~s)h (s) ds. |
|
о |
|
|
d |
Отсюда |
|
|
|
X (t) = eu {t~U)x (t0) + |
J eu u~s)h (s) ds. |
||
§ 5. Интегрирование путем замены переменных |
|||
Решение однородного |
уравнения |
|
|
|
^ - = Ux |
(5.1) |
|
при начальном условии |
|
|
|
|
x(t0) = |
с |
(5.2) |
можно построить и так.
Пусть К — матрица, преобразующая матрицу U к жор-
дановой матрице J, так что |
|
|
|
|
||
U = KJK~l = |
KJM |
(М = |
К~1). |
(5.3) |
||
В уравнении (5.1) произведем замену переменных |
||||||
Получим |
|
х = Ку. |
|
|
|
|
dy |
|
|
|
|
|
|
|
= |
К~ѴКу, |
|
|
||
|
dt |
|
|
|||
или, учитывая |
(5.3), |
dy |
= Jy. |
|
|
|
|
|
|
|
(5.4) |
||
|
|
dt |
|
|
|
|
Матрица J имеет квазидиагональную структуру: J = |
||||||
== diag (Jx (kj), |
..., Jp (kp)). Столбцовую матрицу у |
разобьем |
||||
на блоки так, |
чтобы число |
строк |
/-го |
блока |
равнялось |
■)