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

main_du_ast

.pdf
Скачиваний:
54
Добавлен:
11.02.2015
Размер:
613.31 Кб
Скачать

Определим, при каком T имеем x(T ) = x20 .

x20 = x0e k(T t0);

k(T t0) = ln 2;

T = t0 + k1 ln 2 период полураспада не зависит от начального количества вещества:

Пример 1.17. Уравнение движения материальной точки под действием силы, приложенной вдоль прямой, получается из второго закона Ньютона:

mx• = f(x; t):

Пример 1.18. Скорость распространения бактерий.

Р е ш е н и е. Известно, что скорость распространения бактерий пропорциональна их количеству:

 

_

(0)

=

N0:

N

 

N(t)

=

kN(t);

dNN = kdt;

ln jNj = kt + ln c;

N = cekt;

N0 = N(0) = c;

N = N0ekt численность популяции растет экспоненциально.

Пример 1.19. Движение материальной точки под действием силы тяжести по вертикальной прямой.

Р е ш е н и е. Пусть известны скорость и положение точки в начальный момент времени, то есть выполняются условия x(t0) = x0 и x(t0) = v0. В силу второго закона Ньютона имеем

x•(t) = g; x(t) = gt + c1;

x(t) = 12gt2 + c1t + c2;

v0 = x(t0) = gt0 + c1 ) c1 = v0 gt0;

x0 = x(t0) = 12gt20 + (v0 gt0)t0 + c2 ) c2 = x0 + 12gt20 v0t0; x(t) = 12gt2 + (v0 gt0)t + x0 + 12gt20 v0t0:

Пример 1.20. Малые колебания математического маятника.

Р е ш е н и е. Шарик массы m закреплен на конце невесомой нерастяжимой нити. Отклонение нити от положения равновесия зада¸тся переменной x. На шарик действует сила тяжести F = mg и сила натяжения нити. Закон движения в проекции на касательную:

m

d2(lx)

= mg sin x;

dt2

x• = gl sin x:

Если колебания малые, то sin x x, и уравнение принимает вид

x• = gl x:

21

Пример 1.21. Уравнение семейства окружностей.

Р е ш е н и е. Окружность с центром в точке (x0; y0) радиуса R задается уравнением

(x x0)2 + (y y0)2 = R2:

Трижды продифференцируем данное уравнение по x, считая y функцией от x:

2(x x0) + 2(y y0)y0 = 0;

1 + (y0)2 + (y y0)y00 = 0 ) y y0 = 1 +y(00y0)2

2y0y00 + y0y00 + (y y0)y000 = 0 ) 3y0(y00)2 1 + (y0)2 y000 = 0:

2Существование и единственность решения задачи Коши. Продолжение решений. Непрерывная зависимость решения от начальных условий, правой части и параметра.

2.1Теорема существования и единственности решения задачи Коши

Рассмотрим задачу Коши:

y(x0)0

=

y0

(2.1)

y

=

f(x; y);

 

Теорема 2.1 ( Пикара (существования и единственности решения)). Пусть функция f(x; y)

определена и непрерывна по совокупности переменных и удовлетворяет условию Липшица по y

в прямоугольнике

 

=

f

(x y) :

j

x

 

x

0j

a;

j

y

 

y

0j

b

g. Тогда существует

единственное

 

 

 

 

 

 

 

 

 

 

 

b

решение задачи Коши (2.1), определенное на Uh(x0)

=

fx0

h; x0 + hg, где h

= minfa;

 

g;

M

M = max f(x; y);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство. I. Существование.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Запишем интегральное уравнение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x) = y0 + Zx0

f(s; y(s) ds:

(2.2)

Определение 2.1. Решением уравнения (2.2) называется непрерывная функция, обращающая это уравнение в тождество.

Лемма 2.1 (Об интегральном уравнении). Функция y(x) является решением уравнения (2.2) тогда и только тогда, когда она является решением задачи (2.1)

Доказательство утверждения. Пусть y(x) решение задачи (2.1), тогда оно удовлетворяет тождеству

y0 = f(x; y(x)):

Проинтегрируем:

 

x

 

y(x) y(x0) = Zx0

f(s; y(s) ds:

Функция y(x) непрерывна. Учитывая начальное условие задачи (2.1), получаем, что y(x) решение уравнения (2.2).

Пусть y(x) решение уравнения (2.2). Функция f(x; y) непрерывна по условию, функция y(x)по определению решения. Продифференцируем (2.2):

y0 = f(x; y(x)):

Кроме того, y(x0) = y0. Утверждение доказано.

22

Далее будем рассматривать решение интегрального уравнения (2.2). Метод последовательных приближений Пусть

x

 

x

 

y1(x) = y0 + Zx0

f(s; y0) ds;

y2(x) = y0 + Zx0

f(s; y1(s)) ds:

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

x

 

ym(x) = y0 + Zx0

f(s; ym 1(s)) ds:

Предположим, что эта последовательность сходится к y(x), и разрешен предельный переход под знаком интеграла и функции f. Тогда при m ! 1 получим

x

 

 

y(x) = y0 + Zx0

f(s; y(s)) ds;

(2.3)

то есть y решение уравнения (2.2), а, следовательно, и задачи Коши (2.1). Непрерывность функции y(x) будет доказана позже.

1) Докажем, что все функции последовательности ym(x) определены и непрерывны в Uh(x0)

и не выходят за , то есть jyk(x) y0j b; k = 1; 2; ::: при jx x0j < h; где h = minfa; Mb g; M = sup(x;y)2 jf(x; y)j. (Для всех точек (x; y) 2 выполняется f(x; y) 2 C( ), и, так как

замкнуто и ограничено, то M конечно.)

Докажем это по индукции. Основание индукции k = 1: Рассмотрим функцию y1(x). Она непрерывна при jx x0j a, так как f(x; y) непрерывна, а интеграл с переменным верхним пределом - непрерывная функция на том же отрезке. Если точка (x; y) 2 , то для нее выполняются условия

x

x

 

jy1(x) y0j = j Zx0 f(s; y0) dsj j

Zx0 jf(s; y0)j dsj Mjx x0j b

 

при условии, что jx x0j < h.

 

 

Далее, пусть это утверждение справедливо для yn 1(x), то есть выполняется

 

jyn 1(x) y0j Mjx x0j b:

 

 

Для функции yn(x) получим

 

 

x

 

 

jyn(x) y0j j Zx0 jf(s; yn 1(s))j dsj:

 

Так как (s; yn 1(s)) 2 , то jf(s; yn 1(s))j M, поэтому

 

jyn(x) y0j Mjx x0j b:

(2.4)

2) Докажем, что последовательность функций fyn(x)g1n=1 сходится равномерно к y(x) (yn(x) y(x)) на отрезке jx x0j h.

Построим такой ряд, чтобы последовательность fyn(x)g1n=1 являлась последовательностью его

частичных сумм:

(2.5)

y0(x) + (y1(x) y0(x)) + + (yn(x) yn 1(x)) + : : : :

Докажем равномерную сходимость этого ряда, используя признак Вейерштрасса. Для этого докажем, что для любого n 2 N

jyn(x) yn 1(x)j MLn 1jx x0jn=n! MLn 1hn=n!:

Доказательство будем проводить методом индукции. Пусть n = 1. Тогда, согласно (2.4),

jy1(x) y0(x)j Mjx x0j:

Пусть для n = k неравенство доказано. Докажем его для n = k + 1:

 

x

 

(s))j ds

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

jyk+1(x) yk(x)j

Zx0

jf(s; yk(s)) f(s; yk 1

 

Zx0 Ljyk(s) yk 1(s)j ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

k

1

 

x

 

x0

 

k

 

 

 

L

k

 

x

 

x0

 

k+1

 

 

 

 

 

 

 

 

ML

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

L

Zx0

 

 

j

 

 

 

j

 

=

 

 

 

 

j

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k!

 

 

 

 

 

 

 

 

(k + 1)!

 

 

23

Из того, что jx x0j h, получаем требуемое неравенство: jyn(x) yn 1(x)j MLn 1hn=n!. Ряд PLn 1hn=n! сходится (признак Даламбера). Таким образом, по признаку Вейерштрасса ряд (2.5) сходится равномерно, и его сумма y(x) является непрерывной функцией.

3) Докажем, что функция y(x) решение задачи Коши (2.1). Для этого сначала докажем, что эта функция является решением интегрального уравнения (2.3).

Используя условие Липшица и равномерную сходимость последовательности fyn(x)g к y(x), получим, что для любого " > 0 найдется такое N 2 N, что при всех n N будет выполняться неравенство

x

x

 

 

 

x

 

 

 

 

 

 

Zx0

f(s; yn(s)) ds Zx0

f(s;

 

(s)) ds

Zx0 jf(s; yn(s)) f(s;

 

(s))j ds

y

y

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L Z

 

jyn(s)

 

(s)j ds L" jx x0j L"h:

 

 

 

 

 

 

y

x0

Таким образом, мы обосновали законность предельного перехода под знаком интеграла, откуда следует, что y(x) решение интегрального уравнения (2.3), а следовательно, и задачи Коши 2.1.

II. Единственность.

Пусть y1(x) и y2(x) решения задачи Коши (2.1), тогда они могут быть записаны в следующем виде:

 

 

 

x

 

 

 

 

x

 

 

 

 

 

 

 

 

 

y1(x) = y0 + Zx0

f(s; y1(s)) ds;

y2(x) = y0 + Zx0

f(s; y2(s)) ds:

 

 

 

 

 

Получим оценку для их разности:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 jy1(x) y2(x)j

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

Zx0 jf(s; y1(s)) f(s; y2(s))jds

L Zx0 jy1(s) y2(s)jds

 

 

 

 

 

 

 

 

 

 

 

 

 

sup

 

y1(x)

 

(x)

 

2

M

 

0;

 

 

 

 

 

 

 

Lh

 

y2

 

Lh

 

 

 

 

 

 

 

 

 

 

j

 

 

 

j

 

 

!

 

откуда следует, что

 

 

 

 

 

 

 

jx x0jh

 

 

 

 

 

 

 

 

 

 

(1 Lh)

 

jy1(x) y2(x)j 0

 

 

 

 

 

 

 

 

 

 

 

 

sup

 

 

 

 

 

 

 

 

 

 

 

 

 

jx x0jh

 

 

 

 

 

 

 

 

 

 

 

 

и, выбирая h так, чтобы h <

1

, получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sup

jy1(x) y2(x)j 0:

 

 

 

 

 

 

 

 

 

 

 

 

 

jx x0jh

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, y1(x) y2(x).

2.1.1Лемма Гронуолла.

Лемма 2.2 (Гронуолла). Пусть функция u(x) 2 C[x0; ); u(x) 0. Пусть a 0;

неравенство

u(x) определена на

x 2 [x0; ); 1,

b(x) 2 C[x0; ); b(x)

0, такие, что выполнено

x

 

 

u(x) a +xZ0

b(t)u(t)dt:

Тогда

x

 

 

:

u(x) aeR0

 

 

b(t)dt

 

x

 

x

R

Доказательство. 1. Пусть a > 0. Обозначим через v(x) = a+ b(t)u(t)dt. Тогда v(x0) = a; v(x) > 0.

x0

Очевидно, что

v0(x) = b(x)u(x) b(x)v(x):

Разделим обе части неравенства на v(x), имеем

v0(x) b(x): v(x)

24

Проинтегрируем неравенство на [x0; x]. Получим

x

 

 

 

x

b(t)dt;

xZ0

v0(t) dt xZ0

v (t)

 

 

 

 

 

 

 

 

 

 

x

 

ln jv(x)j ln jv(x0)j xZ0

b(t)dt:

Тогда

 

 

 

 

x

 

u(x)

v(x)

 

 

:

 

aeR0

 

 

 

 

 

b(t)dt

 

 

 

x

 

 

 

 

 

 

 

 

2. Пусть a = 0. Для любого a1 > a лемма доказана, то есть выполнено

x

R

b1(t)dt

u(x) a1ex0

для любого a1 > 0 и для любого x 2 [x0; ). Таким образом, при фиксированном x 2 [x0; ), выбирая a1 сколь угодно малыми, получим u(x) 0, но u(x) 0 по условию леммы, следовательно, u(x) = 0,

следовательно,

x

Z

u(x) b(t)u(t)dt:

x0

Поэтому u(x) 0. Лемма доказана.

Предъявим второе доказательство единственности решения задачи Коши.

Доказательство. Если предположим, что существуют два решения y(x) и z(x), рассмотрим модуль их разности

 

y(x)

 

z(x)

 

x

 

f(t; y(t))

 

f(t; z(t))

dt

 

L x

 

y(t)

 

z(t)

dt

:

j

 

 

 

j

Z

j

 

 

j

 

 

 

Z

j

 

 

j

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, для функции u(x) = jy(x) z(x)j выполнено условие леммы Гронуолла

(в нашем

случае a = 0; b(x) = L). Следовательно, u(x) 0, то есть y(x) = z(x).

 

Единственность доказана.

 

Замечание 2.1. В этом доказательстве единственности достаточно, чтобы существовала положительная непрерывная функция b(x), такая, что

jf(x; y(x)) f(x; z(x))j b(x)jy(x) z(x)j:

Замечание 2.2. Для y0(x) = b(x)y + a(x) при условии непрерывности a(x); b(x) выполняются условия теоремы существования и единственности в области непрерывности коэффициентов уравнения. Таким образом, решение задачи Коши для линейного уравнения

y0(x) + p(x)y = q(x); y(x0) = y0;

где функции p(x); q(x) непрерывны, единственно.

Замечание 2.3. В лемме Гронуолла можно заменить интервал [x0; ) интервалом ( ; x0], где 1. Объединяя утверждения для этих случаев, получим лемму для x 2 (x0 h; x0 + h) для любого h 0.

Замечание 2.4. Отметим, что выполнение условий этой теоремы на всей числовой прямой не гарантирует продолжаемости решения на всю прямую (y0 = y2 2y + 1).

25

2.2Теорема о продолжении решения задачи Коши.

Теорема 2.2. Пусть

y0 = f(x; y) y(x0) = y0;

причем f(x; y); @f@y (x; y) 2 C(G); G R2. Тогда решение задачи Коши можно продолжить так, что оно подойдет сколь угодно близко к границе @G или уйдет сколько угодно далеко от начала координат (x0; y0).

Определение 2.2. Функция y~(x) называется продолжением решения y(x), определенного на интервале ( ; ), на интервал ( 1; 1) ( ; ), если y~(x) – решение задачи Коши, и y~(x) = y(x) на

( ; ).

Докажем вспомогательную лемму.

Лемма 2.3. Пусть G1 – область с границей @G1, G2 – область с границей @G2. Тогда @(G1 \ G2)

@G1 [ @G2.

Доказательство леммы. Если x 2 @(G1 \ G2), то существуют последовательности fxng; fyng такие, что

xn ! x; xn 2 G1 \ G2 yn ! y; yn 62G1 \ G2:

Таким образом,

xn 2 G1 xn 2 G2;

и yn 62G1 или yn 62G2. В первом случае существует бесконечное число yn 62G1, во втором – бесконечное число yn 62G2. Переходим к подпоследовательности и получаем, что в одном случае существуют fxng 2 G1; fyng 62G1 такие, что xn ! x; yn ! x; поэтому x 2 @G1, в другом случае существуют последовательности fxng 2 G2; fyng 62G2, что xn ! x; yn ! x. Следовательно, x 2 @G2. Итак, x 2 @G1 [ @G2.

Лемма доказана.

Доказательство теоремы о продолжении. Рассмотрим три случая: 1. Область G ограничена и замкнута.

Возьмем произвольную точку (x0; y0) 2 G. По локальной теореме существования и единственности

существует такое число

h

0,

что решение задачи Коши существует и единственно на интервале (x

 

 

 

r0

 

 

 

0

h0; x0 + h0). Выберем h0 =

p

 

 

, где r0 = ((x0; y0); @G) ; M = max

f(x; y)

.

 

2

+1

 

 

 

 

 

M

G j

j

 

 

Далее возьмем точку x1 так, чтобы x0 < x1 < x0 + h0 (например, x1 = x0 + 10099 h0). Тогда в силу локальной теоремы существования и единственности решения, решение задачи Коши, проходящее

через точку x1, существует и единственно на интервале (x1 h1; x1 + h1), где h1 =

p

r1

; r1

=

 

 

M2+1

 

((x1; y(x1)); @G), причем эти решения совпадают на общей области определения. Далее берем точку x2 = x1 + h1; 0 < < 1, тогда существует решение, проходящее через точку x2, определенное на

[x2 h2; x2 + h2], где h2 =

p

r2

; r2

= ((x2; y(x2)); @G).

 

 

 

 

 

 

 

 

 

 

M2+1

 

 

 

 

 

Таким образом, построили последовательность точек xk+1 = xk + dk; dk =

p

rk

; rk

=

M2+1

((xk; yk); @G) : Так как G ограничена, то xk ограничена и монотонно возрастает, значит, существует lim!1 xk = b. При этом решение продолжалось на каждый отрезок [xk; xk+1]; а следовательно и

Так как

S

j

 

j

 

j

j

на их объединение

[xk; xk+1] = [x0

; b):

 

 

 

f(x; y) ограничена и

y0

 

=

 

f(x; y) , то y(x) удовлетворяет условию Липшица, где L =

max jy0j = max f(x; y) = M. Таким образом, к последовательности y(xn + hn) применим критерий

GG

Коши сходимости последовательности limx!b 0 y(x) = y . Тогда так как y(x) непрерывна на [x0; b]

и pk(xk; yk) ! p (b; y(b)): Так как xk+1 = xk + d1 + ::: + dk ! b; k ! 1; то ряд

dk сходится, а

d

0; k

! 1

;

откуда

r

k !

0:

Если предположим, что

p

@G

, то

существует " > 0, что

значит, k !

 

 

 

 

 

62

 

P

(p ; @G) > 2"

, но т.к.

p = lim p

k, то 9

N > 0 : n > Np

k 2

U

(p )

– противоречие. Таким образом,

 

 

 

k!1

 

 

8

"

 

 

p 2 @G:

2. Пусть G произвольная область, f(x; y) ограничена в G. Тогда сначала возьмем h0 = p r0 ,

M2+1

и в силу локальной теоремы существования и единственности решение существует и единственно

26

на [x

 

h

 

; x

+ h

 

]

, далее берем точку

x

= x

 

+ h

 

,

тогда существует

единственное решение,

0

 

0

0

 

0

 

1

 

0

 

0

 

r1

 

проходящее через точку x1 и определенное на [x1 h1; x1 + h1], где h1 =

p

 

 

и т.д.

M2+1

 

Тогда либо решение пересечет границу @G, либо уйдет сколь угодно далеко от начала координат (пройдет сколь угодно близко к границе).

3. Пусть G произвольная область, условие ограниченности на f(x; y); @f@y (x; y) не требуем. Тогда рассмотрим круг K с центром в точке (x0; y0) радиуса 2R и отступим от границы области G внутрь на ", получим @G" [ @K. Рассмотрим область G \ K. Тогда по лемме @(G \ K) @G [ @K. G \ K

– замкнутая ограниченная область. Тогда по доказанному решение подойдет сколь угодно близко к границе области G \ K, т.е. оно либо подойдет сколь угодно близко к @K и тогда удовлетворяет условию jy(x)j R (т.к. R можно выбрать произвольно, то это означает, что в этом случае решение уйдет сколь угодно далеко от начала координат (x0; y0)), либо подойдет сколь угодно близко к @G (внутри K), т. е. сколь угодно близко к границе области G (можно доказать, что hn ! 0; n ! 1, следовательно, rn ! 0; n ! 1, т.е. (xn + hn; y(xn + hn)) подойдет сколь угодно близко к границе области G).

Теорема доказана.

Следствие 2.1. Если f(x; y) определена и непрерывна вместе с @f@y в полосе fjx x0j a; 1 < y < +1g, то решение задачи Коши с начальными данными y(x0) = y0 либо пересечет прямые, x = x0 + a; x = x0 a, либо уйдет на бесконечность на (x0; x0 + a) (x0 a; x0).

Замечание 2.5. При доказательстве теоремы о продолжении решения мы изучали продолжение на отрезок [x0; x0 + h], т.е. вправо. Продолжаемость влево доказывается аналогично.

Пример 2.1 (уравнения, все решения которого непродолжаемые). y0 = y2 + 1

dy

y2 + 1 = dx arctg y = x + C y = tg(x + C):

Любое решение уравнения продолжаемо лишь на интервал длины .

Пример 2.2. Рассмотрим уравнение

y0 = x3 y3

Р е ш е н и е. y0 = 0 , x = y, y0 > 0 при y < x, решение возрастает, y0 < 0 при y > x, решение убывает.

Если решение начинается в точке x0 выше прямой y = x, то оно обязательно пересечет эту прямую и войдет в область y < x. Рассмотрим любую точку этой области (x0; y(x0)) и область (x0; y(x0)) < y < x. Тогда в силу поля направлений решение не пересекает границу области, следовательно, оно определено при сколь угодно больших x.

2.2.1Достаточные условия продолжаемости решения на весь интервал

Рассмотрим задачу Коши

y0 = f(x; y)

y(x0) = y0; x0 2 ( ; ); 1 < +1

Теорема 2.3. Пусть

; x0

2 ( 2; ); ;

 

y(0 x0) = y0

;

y = f(x; y); (x; y)

где f(x; y) 2 C( ; ) \ Lipy( ; ); ; = f 1 l < x < +1; y 2 Rg (можно считать, что

fy0 2 C( ; )),

(2.6)

jf(x; y)j a(x)jyj + b(x);

где a(x); b(x) – непрерывные функции на ( ; ).

Тогда решение y(x) задачи Коши продолжается на весь интервал ( ; ).

27

Докажем, что для любого отрезка [ 1; 1] ( ; ); x0 2 [ 1; 1]; решение может быть продолжено на [ 1; 1]. На этом отрезке a(x); b(x) – непрерывны, следовательно существуют такие k и m, что ja(x)j k; jb(x)j m:

Для доказательства используется следующая

Лемма 2.4 (Лемма о дифференциальном неравенстве.). Пусть z(x) удовлетворяет условию на

( ; ); 1 < +1:

jz0(x)j kjz(x)j + m;

где k 0; m > 0, и пусть jz(x)j r0.

Тогда, если k = 0, то jz(x)j r0 + mjx x0j. Если k > 0, то jz(x)j r0ekjx x0j + mk (ekjx x0j 1).

x

Доказательство леммы. Если k = 0, то jz0(x)j m. Справедливо z(x) = z(x0) + R z0(t)dt. Поэтому

x0

имеем

x

Z

jz(x)j jz(x0)j + j jz0(t)jdtj;

x0

jz(x)j r0 + mjx x0j:

Пусть k 0.

jzj2 = z z:

Заметим, что для дифференцируемости jzj надо требовать jzj =6 0. Рассмотрим точку x1 > x0 такую, что z(x1) 6= 0. Положим x = x0, если z(x) 6= 0 на (x0; x1).

Если существует точка x 2 (x0; x1) такая, что z(x) = 0. то возьмем в качестве точки x = supfx : x 2 (x0; x1); z(x) = 0g. Тогда z(x ) = 0, более того, jz(x )j r0. Тогда на (x ; x1) существует jzj0. Продифференцируем

тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

jzj2 = z z;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2jzjjzj0 = 2zz0 2jzjjz0j;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

= 0

на

(x ; x )

 

 

 

z

0

j

z0

j.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда, так как j

j 6

 

 

 

1

 

, получим j j

 

 

 

 

0

 

 

 

 

z0

 

 

 

k z

+ m; то есть r0

(x)

 

kr(x) + m:

Обозначим j

z(x)

j

= r(x)

 

 

 

 

 

условия леммы

z

j

j

j

 

 

 

 

 

 

 

. Тогда изm

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

j

j

 

 

 

 

 

Рассмотрим '(x) = e kx

r(x) +

k

. Докажем, что '(x) монотонно не возрастает на (x ; x ).

Действительно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

'0(x) = ke kx r(x) +

 

+ e kxr0(x)

ke kx r(x) +

 

 

 

 

 

 

 

 

+ e kx (kr(x) + m) = 0:

k

 

 

k

Следовательно, '0(x) 0 и '(x1) '(x ). Таким образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e kx1 r(x1) +

 

 

e kx r(x ) +

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

k

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r(x1) ek(x1 x ) r(x ) +

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

k

 

 

 

 

 

 

 

Заметим, что x1 x x1 x0 и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

m

 

 

 

 

 

 

 

 

 

r(x ) = jz(x )j r0 ek(x1 x0) r0 +

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

k

k

 

 

 

Тогда это неравенство выполнено в произвольной точке x1 > x0. Можно записать

 

 

 

 

 

 

 

 

 

 

 

r(x) ek(x x0) r0 +

m

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

; x > x0:

 

 

 

 

 

 

 

 

 

 

 

 

k

k

 

 

 

Заменой x ! x;

 

x0 ! x0 можно получить аналогичную оценку для x < x0. Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

r(x) = jz(x)j ekjx x0j r0 +

 

= r0ekjx x0j +

 

 

 

 

ekjx x0j 1 :

 

 

 

 

k

 

 

k

 

 

 

Лемма доказана.

28

Доказательство теоремы. Возьмем [ 1; 1] ( ; ). Тогда существуют m 0; k 0 такие, что

j (x)j k; jb(x)j m:

Пусть jy(x0)j = d, обозначим R = max(r0ekjx x0j + mk (ekjx x0j 1)). Тогда из (2.6) следует, что

jy0(x)j = jf(x; y)j kjyj + m на [ 1; 1]. Поэтому jy(x)j r0ekjx x0j + mk (ekjx x0j 1) R. Рассмотрим полосу jyj R + 1, т. е. (R + 1) y R + 1. Получим замкнутый прямоугольник

1; 1;R = f(x; y) : x 2 [ 1; 1]; (R + 1) y R + 1g:

Таким образом, по теореме о продолжении решения в замкнутой ограниченной области решение выйдет на границу 1; 1;R. При этом решение не может выйти на верхнее и нижнее основания y = (R + 1), так как jy(x)j R, следовательно, решение пересекает прямые x = 1; x = 1. Берем последовательности f ng; f ng : n ! ; n ! ; n ! 1. Тогда решение может быть продолжено на любой отрезок [ n; n], а, следовательно, на их объединение, то есть на ( ; ).

Теорема доказана.

2.3Теорема о непрерывной зависимости решения от начальных условий и правой части уравнения.

Теорема 2.4. Пусть задача Коши

(

y0 = f(x; y)

y0 = y(x0)

имеет решение (x), определенное на интервале x0 2 [ ; ], и 1 < +1, причем функции f(x; y) и fy0(x; y) непрерывны в некоторой замкнутой окресности U графика y = (x), x 2 [ ; ].

Тогда для всех " > 0 существует такое > 0 что, для всех решений z(x) другой задачи

 

 

 

 

 

 

 

(z00

= z(x0)

 

 

 

 

 

 

 

 

 

 

z

= g(x; z)

 

 

 

такой что

g(x; z); g0

(x; z)

2

C(U)

(

g(x; z)

- непрерывна и удовлетворяет условию Липшица по

z

)

z

 

 

 

 

и

 

 

 

 

 

jg(x; z) f(x; z)j

(2.7)

 

 

 

 

 

 

и jz0 y0j ; имеем jz(x) (x)j < ", причем z(x) продолжается на [ ; ]: (Малое изменение правой части и начальных условий приводит к малым изменениям решения)

Доказательство. Возьмем U = f(x; y) : x0; x 2 [ ; ]; jy (x)j g. Из условия Липшица и условия (2.7) имеем,

jz0 y0j = jf(x; y) g(x; z)j = jf(x; y) f(x; z)+f(x; z) g(x; z)j jf(x; y) f(x; z)j+jf(x; z) g(x; z)j Ljy zj+ :

Обозначим через u(x) = z(x) (x), тогда ju0(x)j Lju(x)j + на [ 0; 0] [ ; ]. С исполь-

зованием дифференциального неравенства получим (так как jz0 y0j , то есть ju(x0)j ) ju0(x)j Lju(x)j + ; L 0,

Тогда

(

ju(x)j

+ jx x0j; L = 0

ju(x)j

eLjx x0j +

 

eLjx x0j 1 j; L > 0;

L

пока z(x) лежит в U.

Зафиксируем " > 0. Выберем "1 = min("; ) тогда выберем таким образом, что

8

L xj x0

x

0j

 

L x x0

1

;

L = 0

 

+ x

 

 

 

 

 

< "

< e j j

+ L

 

e j j 1 < "1;

L > 0:

:

 

 

 

 

 

 

 

 

 

 

29

Можно записать, что

 

 

 

 

 

 

 

8

Ls

Ls

1

; L = 0

 

 

 

 

 

 

 

 

+ s

 

 

< "

 

 

 

 

 

 

 

< e +

 

(e 1) < "1; L > 0;

 

 

 

 

 

 

 

L

где s = min(

x1

 

x0

; x0

 

x2

). :

 

 

 

 

 

j

 

 

j j

 

j

 

 

 

 

 

Тогда u(x) может быть продолжена до границы области U1 = f(x; y) : x 2 [ 0; 0]; jy (x)j g: При этом решение не может выйти через верхнюю и нижнюю границы U1, (так как на них ju(x)j = "1, a ju(x)j < "1), а значит пересекает прямые x = 0; x = 0, и решение u(x) (а значит и z(x)) может быть продолжено за [ 0; 0].

Эти рассуждения можно повторить для всех [ 0; 0] [ ; ]. Таким образом, получим что u(x)(и z(x)) может быть продолжено на весь отрезок [ ; ].

При этом ju(x)j = jz(x) (x)j < "

2.3.1Теорема о непрерывной зависимости решения от параметра.

Рассмотрим задачу Коши, зависящую от параметра .

(

y0

= f(x; y; )

 

(2.8)

y(x0) = a( ):

Теорема 2.5. Пусть при = 0 решение (x; 0) задачи (2.8) существует на отрезке [x1; x2]. Пусть функции f(x; y; ); fy(x; y; ) непрерывны по совокупности переменных на множестве V = f(x; y; ) : x 2 [x1; x2]; jy (x; 0)j ; : j 0j 1g, a( ) непрерывна в U:

Тогда существует такое > 0, что решение y(x; ) задачи (2.8) непрерывно по ; j 0j < .

3Уравнения первого порядка, не разрешенные относительно первой производной.

Уравнением, не разрешенным относительно первой производной, называется уравнение вида

f(x; y; y0) = 0:

Пример 3.1. Уравнение (y0)2 9x2 = 0 можно свести к двум уравнениям: y0 = 3x и y0 = 3x.

Р е ш е н и е. Их решения

y = 32x2 + c и y = c 32x2:

Через каждую точку плоскости x; y проходит не менее двух решений: по одному из каждого семейства решений.

Теорема 3.1 (существования и единственности для уравнений первого порядка, не разрешенных относительно первой производной). Рассмотрим уравнение

 

 

 

 

 

 

 

 

f(x; y; y0) = 0:

 

 

(3.1)

Пусть

f

2

C1

в области

D

и в точке

(x

; y

; y0

)

2

D

имеем

f = 0

@f=@y0 = 0

.

 

 

 

0

0

0

 

 

,

6

Тогда на любом достаточно малом отрезке [x0 d; x0 + d] существует единственное решение

уравнения (3.1), удовлетворяющее условиям

 

y(x0) = y0; y0(x0) = y00 :

(3.2)

Доказательство. По теореме о неявной функции в некоторой окрестности U точки (x0; y0) существует единственная непрерывная функция g(x; y), удовлетворяющая условиям

f(x; y; g(x; y))

 

0; g(x ; y ) = y0

;

(3.3)

 

0 0

0

 

30