Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

2 курс / Дифференциальные уравнения / Федорюк М., Обыкновенные дифф.уравнения, 1985

.pdf
Скачиваний:
844
Добавлен:
21.05.2015
Размер:
21.78 Mб
Скачать

§ 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 где о —