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

Дифференциальные уравнения 2

.pdf
Скачиваний:
79
Добавлен:
18.02.2016
Размер:
515.6 Кб
Скачать

§ 25. Существование и единственность решения задачи Коши.

Нормальная система обыкновенных дифференциальных урав-

нений первого порядка имеет вид

 

 

 

 

 

 

dxi

 

 

 

 

 

 

 

 

 

= f (t, x

, . . . , x

n

), i = 1, n,

(25.1)

 

 

 

 

 

dt

i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где функции fi : G → R непрерывны на области G Rn+1, i = 1, n. В векторном виде система (25.1) записывается следующим образом:

dx

= f(t, x),

(25.2)

dt

 

 

где x = (x1, . . . , xn)T , (t, x) = (t, x1, . . . , xn)T , f = (f1, . . . , fn)T . Напомним (§2), что решением дифференциальной системы (25.2) мы называем такую гладкую векторную функцию

x : I → Rn, I R, что

 

 

dx(t)

= f(t, x(t)), t I,

(25.3)

dt

и (t, x(t)) G, t I. Если, кроме того, x(t0) = x0, то функция x разрешает задачу Коши

dx(t)

= f(t, x(t)), (t, x) G, x(t0) = x0, x0 I. (25.4)

dt

Теорема 25.1 (интегральный критерий). Для того, чтобы функция x : I → Rn была решением задачи Коши (25.4), необходимо и достаточно выполнение векторного ин-

тегрального тождества

 

 

x(t) = x0 + 0t

f(τ, x(τ))dτ, t I.

(25.5)

t

 

 

111

Доказательство. Необходимость. Пусть функция x является решением задачи Коши (25.4). Тогда выполняется тождество (25.3). Интегрируя его в пределах от t0 до t, получаем соотношение

t

x(t) − x(t0) = f(τ, x(τ))dτ, t I.

t0

А оно в силу x(t0) = x0 равносильно соотношению (25.5). Достаточность. При гладкой функции x суперпозиция

(сложная функция) f(τ, x(τ)) непрерывна по τ. Поэтому векторный интеграл, а с ним правая и левая части соотношения (25.5) дифференцируемы по t. Дифференцируя (25.5) по t, приходим к (25.3). Это означает, что функция x есть решение системы (25.2). Кроме того, из (25.5) вытекает, что x(t0) = x0 + 0 = x0. Таким образом, x есть решение задачи Коши (25.5). Теорема 25.1 доказана.

Будем говорить, что правая часть f системы (25.2) удовлетворяет условию Липшица по переменной x на области

G, если

||f(t, x)−f(t, y)|| 6 L||x−y||, (t, x) G, (t, y) G, (25.6)

где L есть некоторая постояная (постоянная Липшица).

Лемма 25.1 (об условии Липшица). Если на выпуклом по x множестве V Rn+1 существует ограниченная по норме матрица

Df(t, x),

Dx

то есть

DfDx

 

6 K, (t, y) V

 

 

(t, x)

 

 

 

 

 

 

 

 

 

 

(в частности, если матрица Df(t, x) непрерывна на выпук-

112 Dx

лом компактном множестве V ), то функция f удовлетворяет условию Липшица по x на V .

Теорема 25.2 (Пикара–Линделефа). Пусть функция f непрерывна на области G Rn+1 и удовлетворяет условию Липшица по x в некоторой окрестности U G точки (t0, x0). Тогда задача Коши (25.4) локально однозначно разрешима. При этом ее решение x определено, по крайней мере, на отрезке [t0 − h, t0 + h], h > 0 (где h достаточно мало), и может быть построено по методу последовательных приближений

x0(t) = x0, xk(t) = x0 + 0t

f(τ, xk−1(τ))dτ,

 

 

 

t

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

k−1

 

x(t) = lim xk

(t) = x

0

+

k

(t) − x

(t)).

k +

 

(x

 

→ ∞

 

 

 

 

 

 

 

k=1

113

§ 26. Сравнение решений и продолжимость.

На множестве G Rn+1 рассмотрим дифференциальную

систему

 

 

 

dx

= f(t, x).

(26.1)

 

dt

 

 

 

Возьмем решение x : I → R дифференциальной системы (26.1). Каждое сужение x1 данного решения x на промежуток I1 I также является решением этой дифференциальной системы. При этом x будем называть продолжением решения x1. Если же решение x не служит сужением ни одного решения (т.е. если для любого решения x данной дифференциальной системы из x x следует x = x ), то решение x

будем называть непродолжаемым или продолженным.

Теорема 26.1. Каждое решение дифференциальной системы (26.1) является сужением хотя бы одного продолженного решения этой системы.

Далее будем рассматривать дифференциальную систему (26.1) в предположении, что G Rn+1 есть область.

Теорема 26.2 (критерий продолжимости решения).

Если вектор–функция f непрерывна на G, то для продолжимости решения x : (α, β) Rn дифференциальной системы (26.1) вправо от β (влево от α) необходимо и достаточно,

чтобы существовал левый предел x(β − 0) =

t β 0

и

 

lim x(t)

 

 

→ −

 

(β, x(β − 0)) G (существовал правый предел x(α + 0) =

lim x(t) и (α, x(α + 0)) G)).

t→α+0

Будем говорить, что решение x : (α, β) Rn дифференциальной системы (26.1) стремится к границе ∂G области G при t → β − 0, если для любого наперед заданного компактного подмножества G0 G существует β0, α < β0 < β, такое, что при всех t (β0, β) точка (t, x(t)) не принадлежит

G0.

114

Аналогичным образом определяется стремление решения x : (α, β) Rn дифференциальной системы (26.1) стремится к границе ∂G области G при t → α + 0.

Теорема 26.3 (критерий непродолжимости решения).

Если вектор–функция f непрерывна на области G, то решение x : (α, β) Rn дифференциальной системы (26.1) непродолжимо вправо (влево) тогда и только тогда, когда оно стремится к границе ∂G при t → β − 0 (при t → α + 0).

В частности, если область G ограниченная, то график непродолжимого решения целиком не лежит ни в одном замкнутом подмножестве из G.

Теорема 26.4 (теорема сравнения Чаплыгина). Пусть вещественные функции f1 и f2 определены в области G R2 и удовлетворяют неравенству

f1(t, x) < f2(t, x), (t, x) G.

Тогда, если функции x1[s, β] R, x2[s, β] R, являются, соответственно, решениями дифференциальных уравнений

dx

= f1(t, x),

dx

= f2(t, x),

dt

dt

 

 

и удовлетворяют условию x1(s) = x2(s), то

x1(t) < x2(t), t (s, β].

Теорема 26.5 (о промежутке существования непродолжимого решения). Пусть вектор–функция f непрерывна в полосе G =< α, β > ×Rn и удовлетворяет неравенству

|f(t, x)| < Φ(|x|), (t, x) G,

где вещественная функция Φ непрерывна, Φ > 0, если u > 0

115

и

 

 

+

 

 

 

 

 

 

u

du

= +∞.

 

 

 

 

 

 

 

 

Φ(u)

 

 

 

0

 

 

 

 

Тогда всякое непродолжимое решение задачи Коши

 

 

dx

= f(t, x),

x(t0) = x0,

(26.2)

 

 

 

dt

определено на промежутке < α, β >.

В частности, если вектор–функция f непрерывна в полупространстве t > t0 пространства Rn+1 и удовлетворяет нера-

венству

|f(t, x)| 6 M|x|, t > t0, x Rn,

то всякое решение дифференциальной системы (26.1) продолжимо на бесконечный промежуток [t0, +).

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

dx

= A(t)x + f(t),

(26.3)

dt

 

 

с непрерывными по t на I =< α, β > матричной функцией A и векторной функцией f.

Нетрудно видеть, что для произвольного отрезка [a, b] I существует постоянная M > 0 такая, что

|A(t)x + f(t)| 6 M(|x| + 1), t [a, b], x Rn.

Из теоремы о промежутке существования непродолжимого решения следует, что всякое решение линейной дифференциальной системы (26.3), определенное при t = t0, t0 [a, b], продолжимо на отрезок [a, b].

И на основании произвола в выборе отрезка [a, b] I заключаем, что всякое решение линейной дифференциальной системы продолжимо на промежуток I задания данной дифференциальной системы.

116

В заключение отметим, что продолжимость всех решений линейной дифференциальной системы с непрерывными коэффициентами на промежуток задания дифференциальной системы является характерным свойством линейных дифференциальных систем. Нетрудно привести пример нелинейной дифференциально системы, у которой непродолжимые решения имеют разные множества задания (У–13).

117

Dx
118

§ 27. Общий интеграл.

Определение 27.1. Гладкую функцию H(t, x) : G → R будем называть первым интегралом дифференциальной системы (26.1), если для любого ее решения φ(t) :< a, b >→ Rn имеет место тождество H(t, φ(t)) ≡ const, t < a, b >.

Уравнение H(t, x) = C определяет множество Γ в пространстве Otx1 . . . xn. Из определения 27.1 следует, что интегральные кривые, проходящие через точки множества Γ, принадлежат Γ. Множество Γ является n–мерной поверхностью (гиперповерхностью), образованной интегральными кривыми. Такие поверхности будем называть интегральными.

Пусть H(x, y) есть гладкая на области G функция. На множестве таких функций определим линейный оператор D:

DH =

∂H

+

DH

f

(27.1)

∂t

 

 

 

Dx

 

со значениями в множестве непрерывных на G функций. Функцию DH будем называть производной функции H в силу системы (26.1). Данное название объясняется формулой

DH(t, φ(t)) =

d

 

H(t, φ(t)),

(27.2)

dt

 

 

 

где φ(t) есть любое решение системы (26.1). Для доказательства формулы (27.2) достаточно продифференцировать H(t, φ(t)) по t. Из (27.2) вытекает утверждение.

Теорема 27.1. Для того, чтобы гладкая функция H : G → R была первым интегралом системы (26.1), необходимо и достаточно, чтобы

DH = 0, (t, x) G.

Определение 27.2. Первые интегралы H1, . . . , Hm : G0

R будем называть независимыми, если ранг rankD(H1, . . . , Hm) = m

в каждой точке области G0.

Имеет место следующее утверждение.

Теорема 27.2. Пусть матрица Df(t, x) непрерывна на

Dx

области G. Тогда:

1)в некоторой окрестности точки (t0, x0) существуют n независимых первых интегралов системы (26.1);

2)если H1, . . . , Hn есть независимые первые интегралы системы (26.1), определенные в некоторой окрестности точки (t0, x0) G, то система

(H1(t, x), . . . , Hn(t, x)) = (H1(t0, x0), . . . , Hn(t0, x0)) (27.3)

имеет единственное решение y = φ(x), являющееся решением задачи Коши (25.4).

Определение 27.3. Если H1, . . . , Hn : G0 R есть независимые первые интегралы, то соотношение

H1(t, x) = C1, . . . , Hn(t, x) = Cn

будем называть общим интегралом системы (26.1) на области G0.

Согласно утверждению 2 теоремы 27.2 общий интеграл определяет как неявную функцию любое решение системы (26.1) на области G0. Для получения же общего интеграла достаточно разрешить формулу общего решения относительно произвольных постоянных.

Теорема 27.3 (об интегрируемой комбинации). Пусть выражение

φ1dx1 + . . . + φndxn

является полным дифференциалом некоторой функции Φ и пусть

φ1f1 + . . . + φnfn = 0, (t, x) G.

119

Тогда Φ : G → R есть первый интеграл системы (26.1).

Доказательство. Так как

dΦ = φ1dx1 + . . . + φndxn,

то

Φ

∂xk = φk, k = 1, n.

Поэтому

Φ

 

Φ

+ + φnfn = 0, (t, x) G.

 

f1

+ . . . +

 

fn = φ1f1

∂x1

∂xn

Теперь на основании теоремы 27.1 приходим к нашему утверждению.

Рассмотрим теперь автономную систему

 

 

 

 

dx

= f(x).

 

 

 

 

(27.4)

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Представим ее в виде

 

 

 

 

 

 

 

 

 

 

 

 

dx1

= . . . =

dxn

 

= dt.

 

 

 

 

 

 

 

 

 

 

f1(x)

 

 

 

 

 

 

 

 

 

 

fn(x)

 

Тогда систему

 

 

 

 

 

 

 

 

 

 

 

 

 

dx1

 

= . . . =

dxn

 

(27.5)

 

 

 

 

 

f

 

(x)

 

 

f (x)

 

 

n

 

1

 

 

 

 

 

 

 

 

 

 

будем называть системой в симметричной форме. Отметим, что система (27.5) описывает геометрические свойства решений системы (27.4). Каждое дифференциальное уравне-

ние

dxk = dt, k = 1, n, fk(x)

выявляет характер параметризации интегральных кривых. Для представления неавтономной системы (26.1) в симметрической форме необходимо сначала перейти к соответствующей ей автономной системе (как мы указывали ранее).

120