2 курс / Дифференциальные уравнения / Федорюк М., Обыкновенные дифф.уравнения, 1985
.pdf§ 101 УРАВНЕНИЯ, НЕ РАЗРЕШЕННЫЕ ОТНОСИТЕЛЬНО Vм 161
При а > 0 это семейство эллипсов, при а < 0 —• семейство гипербол. В частности, при а==—1 получаем ортогональные семейства равнобочных гипербол
Пример 17. Рассмотрим конфокальное семейство кривых второго порядка
1
где а>0 , Ь>0. Фокусы этих кривых расположены в точках (drYlb —a|f 0), Исключая С из данного уравнения и уравнения
2 2
находим дифференциальное уравнение семейства
Если заменить у' на — 1/у', то это уравнение не изменится. Следовательно, данное семейство софокусных кри-» вых второго порядка образует ортогональную систему.
Пусть 0 < а < & ; тогда семерство существует при С> > — Ъ. Если —6 < С < — а , то кривые являются софокусными гиперболами, ортогональное семейство — софокус-
ные эллипсы. Если —а<С<*>9 то кривые семейства— эллипсы, ортогональные кривые — гиперболы. При значениях С« —6, С*=*—а кривые второго норядка вырождаются.
Пример 18. Найдем изогональные траектории семейства прямых у=**Сх. Дифференциальные уравнения семейства есть ху' = у, уравнение траекторий есть
— ydx^tg a(xdx+y dy).
Переходя к полярным координатам г, <р, получаем
Здесь О 0 — произвольная постоянная. |
Изогональные |
траектории — логарифмические спирали |
(а^я/2), орто- |
гональные траектории — окружности г —С,
ГЛАВА 3
ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ
§ I. Теорема существования и единственности
I. Формулировка теоремы. Рассмотрим линейную систему из п уравнений
е нормальной форме. Здесь у,Ы — неизвестные, ajh{x), fj(x) — известные функции. Все эти функции предполагаются комплекснозначными. Положим
тогда система примет вид
(1)
Поставим задачу Коши
(2)
для системы (1), где у0 — заданный я-вектор. Теорема существования и единственно -
сти для линейных систем. Пусть вектор-функ* ция f (х) и матрица-функция А(х) непрерывны наот* резке I в [а, 6], точка xk e 7, Тогда
§ 1] |
ТЕОРЕМА СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ |
163 |
||
1°. Решение задачи Коши (1), |
(2; существует на всем |
|||
отрезкеL |
|
|
|
|
2\ Решение задачи Коши (1), (2) единственно: если |
||||
у{х), |
z(x)— решения системы(1) с одинаковыми дан-* |
|||
ными Коши (2), |
то у(х) = г(х) |
на всем отрезке I. |
|
|
З а м е ч а н и е . |
В отличие от |
теоремы существования |
и единственности для нелинейных систем (гл. 2, § 1), эта теорема — глобальная: решение существует на всем отрезке /.
2, Доказательство теоремы. Задача Коши (1), (2) эк-* вивалентна системе интегральных уравнений
\ |
(3) |
где обозначено |
|
Применим метод последовательных приближений: |
|
х |
|
у°(х) - g{x)t yl{x) = g{x) + J A |
|
i... |
(4) |
Так как матрица-функция АЫ и вектор-функция |
f{x) |
непрерывны на отрезке /, то все последовательные приближения у0 (х), yt (x)... — непрерывные на отрезке / вектор-функции. Докажем, что последовательные прибли-
жения уЧх) равномерно сходятся на отрезке k
^()^у(х).
Как и в доказательстве основной теоремы (гл. 2» § 1), рассмотрим ряд
...+(/+ 1 И-/(*))+... (5)
Его частичные суммы равны у0 (я), у1 {х), . . . , ун{х)£.. »v так что из равномерной сходимости ряда (5) следует равномерная сходимость последовательности [yh{x)\. Как я
164 ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ [ГЛ. 3
в доказательстве основной теоремы, достаточно доказать сходимость ряда
I.*0 (*)«с +I Vх (*) - У0 (*)Цс+ . •. + 1 / + 1 (*) - УкЛ +•
|
|
(6) |
Напомним, что |
если |
у{х) есть л-вектор, то его норма |
в пространстве |
С (/) |
определяется так: |
*)lie |
|
max (m&x\yk{x)\)f f |
|
|
) |
а норма (н X и)-матрицы А(х) определяются по формулам (гл. 2, §§ 2, 3):
Так как вектор-функция |
g(a:) и матрица-функция А Ы |
непрерывны при j e / , то |
|
где Си Сг —- постоянные. |
|
Оценим нормы в С(1) |
последовательных приближе- |
ний. Имеем Wy°(x)h^Ci. Далее, |
|
X |
II |
Iж о |
1 |
(мы воспользовались неравенством <1М(*)в11()Я и неравенством (3) из § 5, гл. 2)
Далее,
1*0 |
I |
§ 1] |
ТЕОРЕМА СУЩЕСТВОВАНИЯ И ЕДИНСТВЕННОСТИ |
165 |
(используем оценку (7))
Докажем по индукции, что выполняются оценки
(8)
При А: == 1 это неравенство доказано; совершим переход по индукции от к к к + 1. Имеем
(испольвуем оценку (8))
5 к\
Тем самым оценка (8) доказана. Так кад \х—хо\ ^ Ъ —а, то
Правая часть неравенства не зависит от х; поэтому неравенство верно для максимума по х&1 левой части (т. е. для нормы в С)
Так как числовой ряд |
^ -гг[Сг(Ь— |
сходится (его |
|
к\ |
|
сумма равна С^ |
J, то сходится ряд |
(6) и тем са« |
мым равномерная сходимость последовательности [ук (х)} доказана. В силу равномерной сходимости предельная вектор-функция у (х) непрерывна на отрезке /.
166 |
ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ |
[ГЛ. 3 |
Переходя к пределу при &-*«> в соотношении (4), что возможно в силу равномерной сходимости последовательности [ук(я)}, получаем, что у {х) удовлетворяет системе интегральных уравнений (3) при всех х&1. Тем самым существование решения задачи Коши (1), (2) доказано. Единственность вытекает из основной теоремы и теоремы о продолжении решений (гл. 2, § 6).
Рассмотрим задачу Коши для нелинейной системы из п уравнений
4JT-/('i*)f *(*о)-*°* |
(9) |
||
где х — (#i, •.., xn)i f |
— (fii *. •, fn). |
x) непре- |
|
Теорема. |
Пусть вектор-функция f{t, |
||
рывна при teI=[tOl |
t1]lx&Rn no совокупности пере* |
||
менных и |
|
|
|
a* |
t* ™\ |
|
|
|
|
|
(10) |
всех ), к. Тогда решение задачи Коши (9) существу* ет и единственно на всем отрезке L
Наметим доказательство. Для последовательных приближений
t
в силу леммы 3 (гл. 2, § 4) и оценок (9) справедливо неравенство
х
в со индукции, как и выше, доказывается неравенство
где А — некоторая постоянная. Далее см. доказательство предыдущей теоремы.
3. Линейное уравнение л-го порядка. Рассмотрим уравнение
у{п) + at{x)y^l) .+ „, + «л Ыу - fix) |
(11) |
§2] |
ФУНКЦИИ ОТ МАТРИЦ |
167 |
и поставим для него задачу Коши |
|
|
|
Уо, У'Ш =Уи . . ., »(П"!)(Жо) « l/n-i. (12) |
|
Т е о р е м а с у щ е с т в о в а н и я и е д и н с т в е н н о - |
||
сти . |
Пусть коэффициенты аДя), «.., ап(х) и правая |
|
часть fix) уравнения (11) непрерывны |
на отрезке I *=* |
|
•=*[а, |
Ь], Товда |
|
Г.Существование. Решение задачи Коши (И),
(12)существует на всем отрезке хц1.
2°. Единственность. Решение задачи Коши (11) (12) единственно на всем отрезке I.
Доказательство. Заменим уравнение (11) эквивалентной системой уравнений, полагая
Уо (*) - У(«). Уг(х) - Уо {х), ... , Уп-i (х) - Уп-2 (*)# так что jj(«)ej/U ) W. Из уравнения (11) находим
Ул-1 (^) — —«1fa)Уп-1 («) — -- . —fln(^) Уо И + / («)• Мы получили систему вида (2), где
|
1 |
... |
О |
\ |
/О |
|
о |
... |
о |
\ |
/о |
4W-I : |
: |
: |
Ь /(*)- \ |
||
|
о |
... |
1 |
I |
\ о |
Коши (12) принимают вид: у;(х0) =1/> |
|
§ 2. Функции отматриц и однородные линейные системы с постоянными коэффициентами
1. Матричная экспонента. Рассмотрим задачу Коши для линейной однородной системы из п уравнений
- у° (1)
с постоянными коэффициентами. Сведем эту задачу к системе интегральных уравнений
168 ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ [ГЛ.3
и применим метод последовательных приближений:
х
№) -У'^'ЧУ* (*) - У0 + J Ду*-1 СО л-
*«
ВЫЧИСЛЯЯ последовательные приближения, получаем
У1 (*) - У0 + (* - х0) Ау°, у* (х) = у* + (х- хо)Ау* +
В § 1 доказано, что последовательность {ук{%)} равномерно сходится к решению задачи Коши (1), (2) на любом отрезке /. Поэтому решение имеет вид
у{х) = е(*-хдАу\ |
(2) |
где введено естественное обозначение
(3)
Матрица-функция eix~*o)A называется матричной экспо* нентойи обладает рядом свойств скалярной экспоненты. Справедливы формулы
1°. е*АхА |
( + |
2°. {е**)-1 - |
е~*\ |
3°. -£
Докажем Iе. Решение системы (1) с данными Коши У (0) •=У0 имеет вид у (х) — ехАу°. По теореме единственности у {хх + х2) — z (д^), где z (ж)— решение системы (1) с данными Коши z (0) «- у (х2). Имеем
V.
так что e(*i+xi)Aye =, ехЛх*лу° и Г следует из того, что это тождество справедливо для любого вектора у0* Свойство 2° следует из 1°. Далее,
§ 2] |
|
ФУНКЦИИ |
ОТ МАТРИЦ |
|
|
169 |
|
так |
как ехЛу° |
— решение |
системы (1). |
Поскольку |
это |
||
тождество справедливо для любого вектора #°, |
то первое |
||||||
из |
равенств |
3° доказано; |
второе следует |
из |
перестано- |
||
вочности матриц А и е*А. |
|
|
|
|
|
||
|
Матричная экспонента |
&А |
бесконечно дифференци- |
||||
руема на всей оси —<» < х < °°, |
так как всякое решение |
||||||
системы (1) обладает этим свойством. |
|
|
|
||||
|
Рассмотрим ряд из матричных функций порядка |
пХп: |
|||||
|
Six) |
- ВМ +В2(х)+...+ Вк(х)+... |
|
(4) |
Этот ряд называется сходящимся (абсолютно сходящим* ся, равномерно сходящимся), если все ряды из компонент матриц Вк(х):
(Вк(х))я +Шх))я +...+ (Вк(х))я +...
сходятся (абсолютно сходятся, равномерно сходятся). Абсолютная сходимость ряда (4) в точке х0 эквивалентна сходимости ряда из норм:
а равномерная сходимость ряда (4) на интервале / эквивалентна сходимости ряда из норм в С(1):
ВВ1 (яг)Вс + Шх)\\с + ... + \\ВкШс + . . . < « .
Доказываются эти факты точно так же, как и для рядов
из |
вектор-функций (гл. |
2, § 2). Теоремы из анализа |
о |
дифференцирований и |
интегрировании функциональ- |
ных рядов справедливы и для матричных рядов. Ряд (3) для матричной экспоненты сходится абсолютно на всей
оси |
х и равномерно на любом конечном отрезке |
а^х^ |
||
«^ Ь, что следует из доказательства теоремы |
существовав |
|||
ния |
и единственности. |
|
|
|
|
2. Вычисление |
матричной экспоненты* |
Пусть |
А — |
(п |
X п)-матрица, тогда |
|
|
|
го |
Лемма . Если |
Т — невырожденная (п X п)-матрица, |
||
|
|
|
|
(6)
170 ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ [ГЛ. 3
Доказательство вытекает из тождества
Т-*А*Т = Т-1АТ'Т-1АТ...Т-*АТ-Т~*АТ Хк раз).
Эта |
лемма |
необыкновенно упрощает вычисление сте- |
пеней |
матриц. |
Ясно, что в качестве матрицы Т нужно |
взять такую, которая приводит матрицу А к простейшему виду. Пусть матрица А приводится к диагональному виду, т. е.существует матрица Т такая, что
|
(h |
0\ |
Т^АТ _ ЛН \ |
• |
|
Тогда, в силу (6), |
(}* |
° |
|
Из (6) следует, что
T-ieAT
в если матрица А приводится к диагональному виду, то
\Т-\ (8)
В частности,отсюда следует,что
(9)
8десь SpА — след матрицы А, т. е. сумма ее диагональных элементов;
Из линейной алгебры [7,17] известно, что
п
«ЭрЛ, так что S p 4 « 2^i -
Формула (9) дляопределителя матричной экспоненты (доказанная пока только для таких матриц, которые приводятся к диагональному виду) справедлива для любых квадратных матриц.
Формулы (2), (8) позволяют получить новое доказательство формулы длярешений системы (1) (гл.1, § 8). Всякое решение системы (1) имеет вид у —ехАа1 где о —