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

книги / Обыкновенные дифференциальные уравнения

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
2.7 Mб
Скачать

откуда,

допуская, что значения yn 1 D при

 

x x0

 

h , получим

 

 

 

yn y0

 

Mh b , т.е. n-е приближение также не выходит из той же области.

 

 

Таким образом, методом полной индукции доказано, что ни одно из последовательных приближений не выйдет из области D, если

x0 h x x0 + h .

Теперь нужно показать, что существует предел последовательности yn

,

Y ( x) = lim yn ( x)

,

(1.47)

n → ∞

и что предельная функция удовлетворяет уравнению (1.42) и начальным условиям. Для того чтобы показать существование предела (1.47) последовательности {yn}, достаточно доказать сходимость ряда:

y0 + (y1 y0 ) + (y2 y1 ) + ... + (yn yn1 ) + ... ;

(1.48)

действительно, n-я частная сумма этого ряда Sn = yn , и если мы докажем, что ряд (1.48) сходится к функции Y (x) , то этим самым и будет доказано соотношение (1.47).

Оценим абсолютные величины членов ряда (1.18):

 

 

 

 

 

 

 

 

 

 

 

 

y1 y0

 

=

 

x f (x, y0 )dx

 

M

 

x1 x0

 

;

 

 

 

 

 

 

 

 

(1.49, a)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

далее

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 2 y1

 

=

 

 

{f (x , y1 ) f (x , y 0 )}dx

 

 

M

 

 

 

 

f

(x , y1 ) f (x , y 0 )

 

dx

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На основании условия Липшица (1.44) подынтегральная функция

удовлетворяет неравенству

 

 

 

 

 

f (x, y1 ) f (x, y0 )

 

N

 

y1 y0

 

 

; подставляя сюда

 

 

 

 

 

 

 

 

только что найденную оценку для

 

 

 

 

y1 y 0

 

 

 

 

, находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

MN

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2 y1

 

 

N

M

 

x x0

 

dx

=

 

 

x x0

 

 

 

 

 

 

 

(1.49, б)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

Аналогично получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y3 y2

 

 

f (x, y2 ) f

(x, y1 )

 

dx

N

 

 

y2 y1

 

dx

MN

 

 

 

 

(x x0 )2 dx

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 2

 

x0

 

и, наконец,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y3 y2

 

MN 3

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1,49, в)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x0

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Докажем, что для любого целого положительного n справедливо неравенство

31

yn yn 1

 

MN n 1

 

 

x x0

 

n

 

 

 

 

 

 

n!

 

 

 

.

(1,49, n)

 

 

 

Для этого применим метод полной индукции: допустив, что неравенство (1.49, n) справедливо для n, докажем его справедливость для n + 1. Мы имеем:

 

 

 

x

 

x

yn +1 yn

 

 

f (x, yn ) f (x, yn 1 )

 

dx

N

 

yn yn 1

 

dx

 

 

 

 

 

 

 

 

x0

 

x0

(последнее неравенство получаем на основании условия Липшица).

Подставив в последний интеграл вместо

 

y n y n 1

 

выражение (1.49n) (отчего

 

 

неравенство усилится), получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

MN

n 1

 

x

 

 

 

n dx

 

MN

n

 

n +1 ,

 

 

 

 

 

 

 

 

 

 

yn +1 yn

 

 

 

 

x x0

 

=

 

 

 

x x0

 

 

 

 

 

 

 

 

 

 

n!

(n +

1)!

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. неравенство (1.49n) справедливо и для n + 1. А поскольку данное неравенство оправдывается для n = 2 и n = 3, то оно верно для каждого натурального числа n. Заменяя еще x x 0 его наибольшим допустимым

значением h, мы приходим к заключению, что каждый член ряда (1.48) не превосходит соответствующего члена числового ряда с положительными членами:

 

 

 

 

 

 

 

 

 

Nh 2

 

 

N

2 h3

 

 

N n 1hn

 

 

 

y0

 

 

+ Mh +

M

 

 

 

+ M

 

 

+ ... + M

 

 

 

 

 

+ ... .

(1.50)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2!

 

 

 

3!

 

 

 

 

n!

Применяя к последнему ряду признак Даламбера, находим:

 

 

lim

 

u

n+1

= lim

 

MN nhn+1 n!

= lim

Nh

 

= 0

< 1.

 

 

 

un

 

 

 

 

 

 

n + 1

 

n→ ∞

 

 

n→ ∞ (n + 1)!MN n1hn

n

 

 

 

 

Итак, ряд (1.50) сходится. Поскольку все члены ряда (1.48) по абсолютной величине меньше членов числового ряда (1.50), то ряд (1.48) не только сходится, но и, в силу критерия Вейерштрасса, сходится равномерно для всех х, удовлетворяющих условию x x0 h. Каждый член ряда (1.48) есть непрерывная функция верхнего предела. Следовательно,

Y = lim y n

n → ∞

существует и является непрерывной функцией от x.

Докажем, что полученная таким образом функция Y (x) удовлетворяет данному уравнению (очевидно, что начальное условие удовлетворено, так

как Y ( x0 ) = lim yn ( x) = y0 ).

n→ ∞

Возьмем равенство (1.46n):

x

yn (x) = y0 + f (x, yn1 )dx

x0

и перейдем к пределу при n → ∞ .

32

Благодаря равномерной непрерывности функции f (x, y )

по y в области

D мы для любого заренее заданного положительного числа ε

можем найти

такое δ > 0, что неравенство

 

| f (x, y) f (x, y′′) |< ε

 

будет выполнено для тех пар точек (x, y) и (x, y) области D, для которых выполняется неравенство y′ − y′′ < δ (в силу условия Липшица достаточно

взять δ = ε ). Далее из равномерности стремления последовательности yn

N

к пределу возникает возможность для выбранного δ так подобрать натуральное число n0, чтобы при n – 1 > n0 для всех значений x в интервале (x0 h, x0 + h) имело место неравенство

 

 

 

 

 

yn 1 ( x) Y ( x)

 

< δ .

 

 

 

 

 

 

 

 

Сопоставляя оба эти неравенства, мы получаем при n 1 > n0 :

 

 

 

 

 

f ( x, y n 1 ( x)) f ( x,Y ( x))

 

< ε .

 

 

 

 

 

 

Отсюда следует:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

x

 

 

x

 

 

 

 

 

 

 

 

 

 

f ( x, yn1)dx f (x,Y )dx

 

f ( x, yn1) f ( x,Y )

 

dx

< ε h .

 

 

 

 

 

 

x0

 

 

x0

 

 

x0

 

 

Пользуясь произвольностью числа ε, находим

 

 

 

 

 

lim

x f (x, yn 1 )dx = x f (x,Y )dx .

 

 

 

n → ∞

x0

 

 

 

x0

 

Таким образом, переходя к пределу в равенстве (1.46n) при

n → ∞ ,

получаем тождество

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y = y0 + x f (x,Y )dx,

(1.51)

 

 

 

 

 

 

x0

 

т.е. функция Y (x)

удовлетворяет интегральному уравнению (1.45). Далее

функция Y допускает производную по x, так как в правой части тождества (1.51) стоит интеграл от непрерывной функции, допускающий производную по верхнему пределу. Дифференцируя (1.51), получаем

dY = f (x,Y ), dx

т.е. Y(x) удовлетворяет данному уравнению (1.42). Первая часть теоремы доказана.

2. Единственность решения задачи Коши (1.42), (1.43).

Допустим, что кроме решения Y(x) существует еще другое решение Z(x), удовлетворяющее тому же начальному условию: Z (x0 ) = y0 . Без ограничения общности можно предположить, что значения x, для которых

33

Y (x) ≠ Z(x) , находятся вправо от x0 в любой близости от точки x0 (иначе за точку x0 мы взяли бы ту точку, в любой близости которой значения Y (x) и Z(x) перестают быть равными, или заменили бы x на – x ).

Покажем, что наше допущение приводит к противоречию. Возьмем

некоторое малое постоянное число

ε >0 ; по предположению, в замкнутом

интервале x0 x

x0 + ε не всюду

Y = Z ; следовательно, положительная

функция

 

Y ( x) Z ( x)

 

достигает в этом интервале в некоторой точке ξ своего

 

 

наибольшего значения θ >0, причем не может быть ξ = x0 , так как при x = x0

обе функции Y

и Z

равны между собой. Мы имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y ( x) = y0 + x f (x,Y )dx,

Z ( x) = y0 + x f (x, Z )dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

Вычтем

оба

тождества

одно

из

другого, придавая

x значение ξ

и используя при оценке разности условие Липшица:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y (ξ ) Z (ξ )

 

= θ ≤

 

f ( x,Y ) f ( x, Z )

 

dx N

 

Y ( x) Z ( x)

 

dx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

Мы только усилим последнее неравенство, если под интегралом заменим

разность

 

 

Y ( x) Z ( x)

 

 

 

ее

наибольшим

значением

θ

и

если

интеграл

 

 

 

 

распространим

 

на

 

промежуток (x0 , x0 + ε )

вместо промежутка (x0 ,ξ ). Получим

 

 

θ ≤ εθ

откуда

по

 

предположению

θ > 0

,

находим

что

1

Nε ,

что

 

N ,

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

противоречиво,

поскольку ε можно выбрать сколь угодно малым,

например

ε <

1

 

Следовательно, допущение существования второго решения z,

N .

 

удовлетворяющего

тому

же

начальному

условию,

что

 

 

 

Y,

приводит

кпротиворечию.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, при выполнении условий теоремы существует единственное решение задачи Коши (1.42), (1.43).

Продолжение решения задачи Коши (1.42), (1.43)

Существование решения задачи Коши (1.42), (1.43) доказано только

для отрезка I = [x0 h, x0 + h].

Покажем,

что найденное решение может быть

продолжено на всю область Q,

где функция

f (x, y )

определена, непрерывна

и удовлетворяет условию Липшица по y.

 

 

 

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

 

x

(1)

= x

0

 

+ h

 

 

 

 

 

0

 

 

,

 

 

 

y(x0(1) )= y0(1) , причем (x0(1) , y0(1) )

Q .

Найдем прямоугольник D :

 

x x(1)

 

 

a ,

 

 

 

 

 

 

 

 

 

1

 

 

 

0

 

 

1

 

 

 

 

 

 

y y(1)

 

b : D

Q.

 

 

 

 

 

 

 

 

 

 

0

 

 

1

1

 

 

34

Пусть M1 = max f (x, y) , тогда по теореме Коши – Пикара существует

 

 

D1

y(x

(1) )= y

 

 

 

 

 

 

 

 

 

(1)

 

 

 

(1.52)

 

 

 

 

 

0

0 ,

I1

= x0

h1; x0

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

+ h1 , где

 

 

 

 

 

 

 

 

(1)

(1)

 

 

 

b1

 

 

 

 

 

 

 

 

h1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= min a1,

M1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Середина отрезка I1 совпадает с концом I в этой точке? оба решения

совпадают и равны y0(1) . В силу единственности оба решения будут также

совпадать в общей части отрезков I и I1. Таким образом, найденное решение является продолжением ранее полученного на отрезке I решения на вторую половину отрезка I1.

Пусть теперь x0(2) = x0(1) + h1, y(x0(2) )= y0(2) и точка (x0(2) , y0(2) ) Q. Мы можем по начальному условию определить решение на отрезке I2, имеющего общую половину с I1, т.е. продолжить решение на вторую половину отрезка I2 и т.д.

Аналогичное рассуждение может быть проведено и для убывающих значений x.

Таким образом, можно продолжить решение до границы области Q. Пример. С помощью метода последовательного приближения найти

решение задачи Коши

y′ = y, y(0) =1.

x

Эта задача Коши эквивалентна уравнению y(x) =1+ y(τ )dτ .

0

На его основе строим последовательные приближения: y0 =1,

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y1 =1 + dτ = 1+ x,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

x

2

 

 

 

 

 

 

 

y2 =1 + (1 + τ )dτ = 1+ x+

 

,

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

y

=1

+

1 + τ +

τ

 

dτ =

1+

 

x+

x

+

x

 

,

x

2

 

3

 

 

 

 

 

 

 

 

 

 

2

 

 

 

3

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

3!

 

 

 

 

 

2

 

 

 

 

 

 

 

...

yn =1 + x + x2 + ... + xn , 2! n!

....

Очевидно, что lim yn (x) = ex = y(x) .

n→ ∞

35

Теорема существования решения задачи Коши (1.42), (1.43) (Теорема Пеано)

 

 

 

Пусть

 

функция

f (x , y )

непрерывна

 

в

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

D :

 

 

x x0

 

a,

 

y y0

 

b . Пусть, далее,

M = max

 

f (x, y)

 

. Тогда задача Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

(1.42), (1.43)

 

имеет по

крайней мере

одно

решение,

определенное для

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

x

[x0 h, x0 + h], где h = min a,

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

Для доказательства

этой теоремы

мы

применим

теорему Арцела и

принцип неподвижной точки Шаудера. Введем предварительно следующие определения.

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

 

 

 

 

 

f1(x), f2 (x), ..., fn (x), ...

 

называется равномерно ограниченной на [a,b], если существует число K > 0

такое, что

для всех

натуральных

n

и любого

x [a,b] справедливо

неравенство

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f n ( x)

 

K.

 

 

 

 

 

 

 

 

 

Определение. Представленная выше последовательность непрерывных

на [a,b] функций называется равностепенно непрерывной, если для любого

ε >0 существует δ (ε ) > 0

такое, что для любых x1, x2

[a,b], удовлетворяющих

неравенству

 

x1 x2

 

< δ , и любого натурального n выполняется неравенство

 

 

f n ( x1 ) f n ( x2 ) < ε .

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

Эта теорема является критерием компактности в пространстве C [a, b] непрерывных на [a, b] функций.

Ее доказательство приведено, например, в [11].

Теорема Шаудера. Если вполне непрерывный оператор F отображает ограниченное замкнутое выпуклое множество B банахова пространства Y на свою часть, то существует неподвижная точка этого отображения, т.е. такая точка y B , что Fy = y.

Доказательство этой теоремы, а также соответствующие определения также приведены в [11].

Перейдем теперь к доказательству теоремы Пеано, воспользовавшись теоремами Арцела и Шаудера.

Уравнение (1.52) вместе с начальным условием эквивалентно интегральному уравнению

36

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y ( x) = y0 + x f (τ , y (τ ))dτ

.

 

(1.53)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим оператор F, определенный равенством

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fy = y0 + x f (τ , y(τ ))dτ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

на шаре B :

 

 

 

y y0

 

 

 

b пространства C[x0 h, x0 + h]. Покажем, что оператор F

 

 

 

 

вполне непрерывен на этом шаре.

{yn (x)},

 

 

 

 

Прежде

всего, если последовательность

принадлежащая шару

 

 

y y 0

 

 

 

b, равномерно сходится к функции

y(x) , очевидно, принадлежащей

 

 

 

 

 

 

тому же

 

шару, то в силу непрерывности функции f (x, y )

имеем

f (x, y (x))

f (x, y(x))

равномерно на [x0 h, x0 + h] . Отсюда вследствие

n

возможности

предельного перехода под знаком интеграла при равномерной сходимости

получаем, что Fyn

Fy , т.е. оператор F непрерывеннашаре

 

 

 

y y 0

 

 

 

b .

 

 

 

 

Далее для любого элемента y(x) шара

 

 

 

 

y y 0

 

 

 

b имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fy( x)

 

 

y0

 

+

 

 

 

x f (τ , y(τ ))dτ

 

 

 

 

y0

 

+ Mh.

(1.54)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если x1 и x2

– две любые точки отрезка [x0 h, x0 + h] , то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fy( x2 ) Fy( x1 )

 

 

 

 

 

2 f (τ , y(τ ))dτ

M

 

x2 x1

 

.

(1.55)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Неравенства (1.54) и (1.55) в силу теоремы Арцела показывают, что

оператор F преобразует

шар

 

 

 

y y0

 

 

 

b

 

 

 

в

относительно

 

компактное

 

 

 

 

 

 

 

 

множество.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Покажем, наконец, что оператор F преобразует этот шар в себя. Действительно,

Fy ( x) y

 

= max

x

f (τ , y(τ ))dτ

Mh M

b

= b

0

 

 

x

 

 

 

M

.

 

 

x0

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, оператор F удовлетворяет всем условиям теоремы Шаудера. Поэтому существует неподвижная точка этого оператора, т.е. такая функция y(x) , что

y( x) y0 + xf (τ , y(τ ))dτ .

x0

Это равенство равносильно двум равенствам:

dy(x) = f (x, y(x)), y(x0 ) = y0, dx

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

37

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

2.1. Уравнения первого порядка n-й степени

 

Рассмотрим уравнение первого порядка

 

F (x, y, y) = 0 .

(2.1)

Будем предполагать, что левая часть уравнения (2.1) является

многочленом относительно y степени n, причем коэффициенты

этого

многочлена в некоторой области D есть непрерывные функции от x и y, допускающие по y непрерывные частные производные.

F (x, y, y)

A (x, y )yn + A

(x, y )yn 1 + ... + A (x, y) = 0 .

(2.2)

 

n

n 1

0

 

Уравнение вида (2.2) называется уравнением первого порядка n

степени относительно

y . Затем допустим, что в области D коэффициент

при старшей степени

y , т.е.

An (x, y ) , нигде не обращается в нуль. Тогда по

основной теореме алгебры уравнение (2.2) для всякой пары значений x, y в рассматриваемой области имеет n решений y (действительных или мнимых). По теореме о неявных функциях каждое из действительных решений является непрерывной функцией от x и y и имеет конечную частную

производную

 

y

, если для рассматриваемых значений x, y, y

имеем

F

0

 

 

 

 

 

 

 

 

 

 

y

 

y

 

(т.е. если для данных значений x и y уравнение (2.2) относительно

y

не

имеет

кратного корня). Таким образом, во всякой области, в которой

 

 

F

 

a > 0 ,

разрешенное относительно y уравнение (2.2)

удовлетворяет

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

условию Липшица.

 

 

 

 

 

 

 

 

 

Итак, пусть в области D уравнение (2.1) имеет k действительных

решений ( k

 

n ):

 

 

 

 

 

 

 

 

 

 

 

 

y′ = f1 (x, y ), y′ = f 2 (x, y ), ..., y′ = f k (x, y ).

 

(2.3)

Каждое из уравнений (2.3) есть уравнение такого вида, для которого доказано существование единственного решения (теорема Коши – Пикара), проходящего через точку (x0, y0); следовательно, уравнение (2.1) допускает в точности k интегральных кривых, проходящих через данную точку (x0, y0) области D.

Пример 1.

y2 + yy′ − x2 xy = 0.

38

Разрешая относительно y (или разлагая левую часть на множители), получаем два уравнения: 1) y′ = x , 2) y′ = − y x ; обе правые части однозначны и непрерывны во всей плоскости xy; соответственно, имеем два решения

1) y = x2 + C ; 2) y = Cex x +1.

2

Пример 2.

y2 + y2 1 = 0 .

 

 

Разрешая относительно

y,

 

 

:

 

y

= ±

1

y2

;

в полосе

1

<

y

<

1

 

 

 

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

плоскости xOy оба рассматриваемых значения

y действительны и различны,

при

 

y

 

> 1 значения

y

мнимы,

и при

y = ±1

 

два корня уравнения для

y

 

 

 

совпадают. Рассмотрим

полосу

1 < y < +1;

через

 

каждую точку (x0 ,

y0 )

проходят два решения

y

:

для одного из них значение производной в этой точке

0

 

y′ = + 1 y2 , для другого

y′ = −

1 y2

. В

 

данном

 

случае

в уравнении,

0

0

 

 

 

 

0

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

разрешенном относительно

y,

переменные разделяются и мы имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

:

 

 

 

 

 

 

 

 

 

 

 

1)

 

 

dy

 

= dx

; 2)

 

 

dy

 

= −dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 y2

 

 

1y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получаем решения: 1) arcsin y = x + C , 2)

arcsin y = −x + C . Откуда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = sin(x + C),

 

 

< x + C <

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = sin(x + C),

 

 

 

π

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

< −x + C <

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возьмем

начальные условия

 

 

 

x = x0 , y = y0 ; для

 

определения

произвольного постоянного C получаем уравнения:

 

 

 

 

y0 = sin (± x0 + C ),

± x0 + C = arcsin y0 ,

C = ± x0 + arcsin y0 .

Решения, удовлетворяющие нашему начальному условию:

y = sin[(x x0 ) + arcsin y0 ],

π

 

< x x0

+ arcsin y0 <

π

;

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

y = sin[(x x0 ) + arcsin y0

],

π

< −x + x0 + arcsin y0 <

π

.

 

 

 

 

 

 

2

 

 

2

 

Первое

семейство

решений

представляет

семейство

дуг синусоид,

взятых от точки минимума до точки максимума (исключая концы); все кривые получаются из одной, например

y = sin x,

π

< x <

π

,

 

 

2

2

 

39

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

y = sin(x),

π

< x <

π

.

 

 

2

2

 

Через каждуюточку полосы 1 < y < 1 проходитпооднойкривойкаждого семейства – возрастающая криваяпервогосемейства иубывающая– второго.

Расширим теперь область, в которой рассматриваются кривые, добавив обе граничные прямые y = ±1 , на которых условия Липшица не выполняются

(у нас

f (x, y) = ± 1y2 ; производная f =

 

y

2

становится бесконечной

 

y

1

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при y

±1 ), но правые части уравнений остаются непрерывными. Теперь нам

 

 

 

 

 

 

 

π

π

 

нет необходимости ограничивать изменение аргумента пределами

 

,

 

.

2

2

 

 

 

 

 

 

 

 

 

Замечая, далее, что sin (x + C ) = sin (x + π + C ), мы можем представить общее

решение одной формулой:

y = sin (x + C ), − ∞ < x < +∞ .

Это семейство полных синусоид, получаемых из одной параллельным переносом вдоль оси x; чтобы получить каждую один раз, надо ограничить область изменения произвольного постоянного С, например 0 C≤ π2 . В точках

максимума

одной полной синусоиды кривая первого семейства

переходит

вкривую

второго, в точке минимума имеет место обратный

переход

(см. пунктирные кривые на рис. 2.1). Но эти переходы совершаются на прямых y = ±1 , где условия Липшица не выполнены; поэтому, применяя теорему Коши– Пикара, мы не можем связать убывающую и возрастающую дугу синусоиды в одну интегральную кривую; теорема существования дает две отдельные ветви, проходящиечерезкаждуюточкуполосы(сплошныелиниинарис. 2.1).

y

I

0

x

 

II

Рис. 2.1

40

Соседние файлы в папке книги