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

Обыкновенные дифференциальные уравнения в задачах. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
688 Кб
Скачать

ГЛАВА 3.

Нелинейные дифференциальные уравнения

В этой главе мы приведем некоторые сведения из теории нелинейных дифференциальных уравнений. Мы сформулируем классические теоремы о локальной разрешимости задачи Коши для нелинейных дифференциальных уравнений, обсудим важное понятие непродолжаемого решения, приведем теоремы о непрерывной и дифференцируемой зависимости непродолжаемых решений задачи Коши, сформулируем теорему Ляпунова о существовании решения в целом . В последнем параграфе мы рассмотрим некоторые классы дифференциальных уравнений, решения которых можно получить в квадратурах с помощью элементарных методов интегрирования, а также классы дифференциальных уравнений, допускающих понижение порядка.

Ÿ 3.1. Системы нелинейных дифференциальных уравнений

Рассмотрим систему дифференциальных уравнений следующе-

ãî âèäà:

 

 

 

dy

 

 

 

 

 

 

= f(t; y);

(3:1:1)

 

 

 

 

dt

 

 

 

 

 

 

где вектор-функция f(t; y) :

G Rn+1 ! Rn задана. В разверну-

той записи система (3.1.1) имеет вид

 

8

 

dy

 

 

 

1

= f1(t; y1; : : : ; yn);

 

 

dt

 

>

: : : : : : : : : : : : : : : : :

 

>

 

 

 

 

 

 

>

 

 

 

 

 

 

>

 

 

 

 

 

 

>

 

 

 

 

 

 

<

 

 

 

 

 

 

>

 

dy

 

 

 

dtn = fn(t; y1; : : : ; yn):

 

>

 

 

 

 

 

 

>

 

 

 

 

 

 

>

>

:

51

Мы всегда будем рассматривать случай, когда вектор-функция f(t; y) является непрерывной в области G.

Определение. Решением системы (3.1.1) на интервале (a; b) называется вектор-функция

y(t) = (y1(t); : : : ; yn(t))T 2 C1(a; b)

такая, что

(t; y(t)) 2 G ïðè t 2 (a; b); dtd y(t) f(t; y(t)) íà (a; b):

Аналогичным образом дается определение решения системы (3.1.1) на отрезке [a; b] (или на промежутках [a; b), (a; b]) ñ òîé

лишь разницей, что в граничных точках рассматриваются односторонние производные.

Определение. Множество всех точек

ft; y1(t); : : : ; yn(t)g;

t 2 (a; b);

называется интегральной кривой системы (3.1.1).

Основной задачей для системы (3.1.1), которую мы будем рассматривать, является задача Коши.

Определение. Пусть (t0; y0) 2 G фиксированная точка. Задача о нахождении решения y(t) системы (3.1.1), принимающего ïðè t = t0 значение y0, называется задачей Коши:

8

dy

= f(t; y);

 

 

 

dt

(3:1:2)

: j

t=t0 = y0:

 

< y

 

Условие y(t0) = y0 называется начальным условием.

52

Очевидно, нахождение решения задачи Коши (3.1.2) геометрически эквивалентно построению интегральной кривой системы

(3.1.1), проходящей через заданную точку (t0; y0).

В дальнейшем, чтобы не загромождать изложение, мы будем интересоваться решениями задачи Коши (3.1.2), определенными

на интервале (a; b) или на отрезке [a; b], для которых точка t0 ÿâ-

ляется внутренней. Следующая теорема гарантирует локальную разрешимость задачи Коши (3.1.2).

Теорема (Дж. Пеано). Пусть замкнутый цилиндр

Ö = f(t; y) 2 Rn+1 : jt

 

t0j T; y

y0 rg

(3:1:3)

содержится в области G. Обозначим

 

 

 

 

 

max

f(t; y)

 

;

T

 

= min

T;

r

:

 

 

F

 

F = (t;y)2Ö

 

 

 

0

 

{

}

 

Тогда на отрезке fjt t0j T0g существует решение задачи Коши (3.1.2).

Доказательство теоремы Пеано можно прочитать, например, в учебниках [2, 9, 10].

Отметим, что в теореме Пеано речь идет только о существовании решения задачи Коши. Простейшие примеры показывают, что для того, чтобы гарантировать единственность решения, нуж-

ны дополнительные требования на правую часть f(t; y) системы

дифференциальных уравнений.

В качестве одного из таких примеров рассмотрим следующую задачу Коши для скалярного уравнения

8 dy

< dt = jyj;

: yjt=t0 = y0:

Для наглядности неединственности решения полезно нарисовать интегральные кривые в плоскости переменных (t; y), проходящие

в окрестности оси Ot.

53

В следующей теореме гарантируется однозначная разрешимость задачи Коши (3.1.2).

Теорема (Э. Пикар). Пусть цилиндр Ц содержится в области G и в этом цилиндре вектор-функция f(t; y) удовлетворя-

ет условию Липшица по y, т. е. существует константа L > 0 такая, что для любых двух точек (t; y1), (t; y2) 2 Ц выполнено неравенство

f(t; y1) f(t; y2) L y1 y2 :

Тогда на отрезке fjt t0j T0g, указанном в теореме Пеано, определено единственное решение задачи Коши (3.1.2).

Следствие. Пусть вектор-функция f(t; y) имеет непрерыв-

@

ные в области G производные @yj fi(t; y), i; j = 1; : : : n. Тогда на отрезке fjt t0j T0g определено единственное решение задачи Коши (3.1.2).

Доказательство теоремы Пикара можно прочитать, например, в [2, 9, 10]. Однако его нетрудно провести самостоятельно, используя схему, основанную на доказательстве следующей леммы.

Лемма 3.1.1. Решение задачи Коши (3.1.2), определенное на

некотором промежутке a; b , ïðè ýòîì t0 2

(a; b), является ре-

шением системы интегральных уравнений

 

y(t) = y0 + t

f(s; y(s)) ds:

(3:1:4)

t0

 

 

И наоборот, непрерывное на a; b решение системы интегральных уравнений (3.1.4) является решением задачи Коши (3.1.2).

Из леммы 3.1.1 вытекает, что для доказательства теоремы Пикара достаточно доказать однозначную разрешимость на отрезке

54

fjt t0j T0g системы интегральных уравнений (3.1.4) в классе

непрерывных вектор-функций.

Для доказательства разрешимости системы (3.1.4) построим последовательность вектор-функций fyk(t)g по следующему пра-

âèëó:

y0(t) y0;

 

 

yk(t) y0 + t f(s; yk 1(s)) ds;

k 1:

(3:1:5)

t0

 

 

Лемма 3.1.2. На отрезке fjt t0j T0g определены и непрерывны все вектор-функции yk(t), k = 1; 2; : : : .

Лемма 3.1.3. На отрезке fjt t0j T0g имеют место оценки

 

yk(t) yk 1

(t)

 

F Lk 1

jt

t0jk

;

k = 1; 2; : : : :

 

k!

 

 

 

 

 

 

Лемма 3.1.4. Последовательность fyk(t)g равномерно сходится на отрезке fjt t0j T0g, при этом предельная функция y(t) является непрерывным решением системы интегральных уравнений (3.1.4).

Лемма 3.1.5. Непрерывное на отрезке fjt t0j T0g решение y(t) системы (3.1.4) определено однозначно.

Замечание 1. Отрезок fjt t0j T0g называется отрезком Пеано Пикара.

Замечание 2. Последовательность приближенных решений fyk(t)g, определенных формулами (3.1.5), называется последова-

тельностью Пикара. Из лемм 3.1.3, 3.1.4 можно получить оценку погрешности для k-го приближения Пикара.

55

В следующих двух теоремах речь идет о свойствах гладкости решения задачи Коши (3.1.2) в зависимости от гладкости правой

части f(t; y).

Теорема (О. Коши). Пусть вектор-функция f(t; y) анали-

тична в окрестности точки (t0; y0). Тогда задача Коши (3.1.2) имеет единственное аналитическое решение, определенное в некоторой окрестности точки t0.

Доказательство теоремы Коши можно прочитать, например, в [9].

Теорема Коши позволяет искать решение задачи (3.1.2) в виде степенного ряда

y(t) = c0 + c1(t t0) + c2(t t0)2 + : : : ;

где в силу начального условия c0 = y0.

Отметим, что исторически доказательство теоремы Коши было первым доказательством однозначной разрешимости задачи (3.1.2), когда правая часть уравнения имеет произвольный вид.

Конечно, условие аналитичности f(t; y) является достаточно жестким и на практике выполняется очень редко.

Теорема. Пусть f(t; y) 2 Cm(G), тогда в некоторой окрест- ности точки t0 решение задачи Коши (3.1.2) имеет непрерывные производные до порядка m + 1 включительно.

Рассмотрим теперь скалярное дифференциальное уравнение произвольного порядка

dnu

= g (t; u;

du

; : : : ;

dn 1u

dtn

 

dt

dtn 1

)

;(3:1:6)

где заданная функция g : G Rn+1 ! R непрерывна.

56

Определение. Решением уравнения (3.1.6) на интервале (a; b)

называется функция

u(t) 2 Cn(a; b)

такая, что

 

 

 

du(t)

 

 

 

 

dn 1u(t)

 

 

 

 

 

(t; u(t);

 

 

 

 

; : : : ;

 

 

 

 

 

) 2 G;

 

 

t 2 (a; b);

 

dt

 

 

 

dtn 1

 

 

nu(t)

 

 

 

 

 

 

 

 

 

du(t)

 

 

dn 1u(t)

 

 

d

g

(t; u(t);

 

 

 

; : : : ;

 

 

 

)

íà (a; b):

 

dtn

 

dt

 

dtn 1

 

Определение. Пусть (t0; u0; u1; : : : ; un

1) 2 G фиксирован-

ная точка. Следующая задача

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dnu

 

 

 

 

 

 

 

 

 

du

 

 

 

dn

1u

 

 

 

8

 

dtn

 

= g (t; u; dt

; : : : ; dtn 1 ) ;

 

 

 

> u t=t0 = u0;

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

du

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

= u

;

 

 

 

 

 

 

 

 

(3:1:7)

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

t=t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<

 

 

 

 

: : : : : : : : : : : : : : : : : : : :

 

 

 

 

: : :

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

d

n

 

1

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

n

1

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

= u

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

t=t0

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

называется задачей Коши для уравнения (3.1.6).

Из следующей леммы вытекает, что решение задачи Коши (3.1.7) можно эквивалентным образом свести к решению задачи Коши вида (3.1.2).

Лемма 3.1.6. Пусть u(t) решение дифференциального уравнения (3.1.6) на интервале (a; b), тогда вектор-функция

y(t) = (u(t);

du(t)

; : : : ;

dn 1u(t)

)

T

dt

dtn 1

 

57

является решением системы дифференциальных уравнений

dy

= f(t; y); y = (y1; : : : ; yn)T ;

(3:1:8)

dt

 

 

ãäå

f(t; y) = (y2; : : : ; yn; g(t; y1; : : : ; yn))T ;

на интервале (a; b). И наоборот, пусть y(t) = (y1(t); : : : ; yn(t))T решение системы дифференциальных уравнений (3.1.8) на интервале (a; b), тогда вектор-функция y(t) имеет вид

y(t) = (y1(t);

dy1(t)

 

dn 1y1(t)

)

T

 

; : : : ;

 

;

dt

dtn 1

ïðè ýòîì y1(t) удовлетворяет уравнению (3.1.6) на интервале

(a; b).

Упражнение. Сформулируйте теоремы об однозначной разрешимости задачи Коши (3.1.7), являющиеся аналогами теоремы Пикара, Коши.

Ÿ3.2. Непродолжаемые решения

Âпредыдущем параграфе мы сформулировали теоремы Пеано

èПикара, в которых изучается вопрос о локальном существовании решения задачи Коши для системы дифференциальных уравне-

íèé

 

dy

 

 

 

<

 

 

 

 

 

 

 

8

dt

= f(t; y);

(3:2:1)

:yjt=t0 = y0:

Âэтом параграфе мы рассмотрим вопрос о существовании решения задачи (3.2.1), определенного на максимально широком ин-

тервале (t ; t+). Для этой цели мы дадим некоторые новые определения.

58

Определение. Пусть y(t) è ye(t) решения системы

dy

= f(t; y);

f(t; y) 2 C(G);

(3:2:2)

dt

определенные на промежутках a; b è a; b соответственно, при

ýòîì a; b a; b è y(t) y(t) ïðè t 2 a; b . Тогда решение y(t)

называется продолжением решения

 

 

ñ

e

íà

e e

, а решение

 

 

 

 

 

 

 

 

 

 

 

 

 

ye(t)

 

 

 

y(t) e a; b

 

 

a; b

 

e

y(t)

сужением решения

íà

 

a; b

 

.

 

 

e e

 

 

 

 

 

 

 

 

 

Определение. Решениеey(t) системы дифференциальных урав-

нений (3.2.2), определенное на промежутке

t ; t+ , называется

непродолжаемым, если не существует решения системы (3.2.2) с более широкой областью определения, сужение которого на t ; t+ совпадало бы с y(t).

Упражнение. Нарисовать интегральные кривые в плоскости переменных (t; y) для задачи Коши

<

dy

 

8 dt

= "y + y2; " > 0;

: yjt=0 = y0:

Определить зависимость интервала (t ; t+) от параметров " è y0. В дальнейшем будем предполагать, что вектор-функция f(t; y)

удовлетворяет условиям теоремы Пикара.

Теорема 3.2.1. Для любой начальной точки (t0; y0) 2 G задача Коши (3.2.1) имеет единственное непродолжаемое решение y(t), определенное на некотором интервале (t ; t+).

Доказательство этой теоремы можно прочитать в [10, 12]. Следующие две теоремы позволяют получить наглядное объ-

яснение того, что интервал (t ; t+) определения непродолжаемого решения y(t) задачи Коши (3.2.1) является наиболее широким интервалом, на котором может существовать решение этой задачи.

59

Теорема 3.2.2. Пусть y(t) непродолжаемое решение задачи Коши (3.2.1) с интервалом определения (t ; t+). Тогда для любого

компакта K G существует отрезок [t1; t2] (t

; t+) такой,

÷òî (t; y(t)) не принадлежит K ïðè t ̸2[t1; t2].

 

Теорема 3.2.3. Пусть область G является полосой

G = f(t; y) 2 Rn+1 : t 2 (a; b); y 2 Rng

(3:2:3)

è y(t) непродолжаемое решение задачи Коши (3.2.1) с интервалом определения (t ; t+). Тогда

à) ëèáî t+ = b, ëèáî t+ < b è y(t) ! 1 ïðè t ! t+; á) ëèáî t = a, ëèáî t > a è y(t) ! 1 ïðè t ! t .

Доказательство этой теоремы можно прочитать в [12]. Приведем теперь две теоремы о глобальном существовании ре-

шения задачи Коши.

Теорема 3.2.4 (A. Уинтнер). Пусть область G является полосой (3.2.3). Предположим, что вектор-функция f(t; y) удовлетворяет условию

f(t; y) φ(t)g( y );

ãäå φ(t) > 0 непрерывная функция на [a; b], g(u) > 0 непрерывная функция на [0; 1), при этом для любого > 0 > 0

1

1 du = 1: g(u)

Тогда для любой начальной точки (t0; y0) 2 G непродолжаемое решение задачи Коши (3.2.1) определено на всем интервале (a; b).

Доказательство этой теоремы можно прочитать в [12]. Рассмотрим задачу Коши для систем дифференциальных урав-

нений в случае, когда вектор-функция, стоящая в правой части,

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]